Double scattering and pencil beam scanning Monte Carlo workflows for proton therapy retrospective studies on radiation-induced toxicities

PURPOSE

Monte Carlo (MC) simulations can be used to accurately simulate dose and linear energy transfers (LET) distributions, thereby allowing for the calculation of the relative biological effectiveness (RBE) of protons. We present hereby the validation and implementation of a workflow for the Monte Carlo modelling of the double scattered and pencil beam scanning proton beamlines at our institution.

METHODS

The TOPAS/Geant4 MC model of the clinical nozzle has been comprehensively validated against measurements. The validation also included a comparison between simulated clinical treatment plans for four representative patients and the clinical treatment planning system (TPS). Moreover, an in-house tool implemented in Python was tested to assess the variable RBE-weighted dose in proton plans, which was illustrated for a patient case with a developing radiation-induced toxicity.

RESULTS

The simulated range and modulation width closely matches the measurements. Gamma-indexes (3%/3mm 3D), which compare the TPS and MC computations, showed a passing rate superior to 98%. The calculated RBE-weighted dose presented a slight increase at the necrosis location, within the PTV margins. This indicates the need for reporting on the physical and biological effects of irradiation in high dose regions, especially at the healthy tissues and increased LET distributions location.

CONCLUSION

The results demonstrate that the Monte Carlo method can be used to independently validate a TPS calculation, and to estimate LET distributions. The features of the in-house tool can be used to correlate LET and RBE-weighted dose distributions with the incidence of radiation-induced toxicities following proton therapy treatments.

Secondary neutron dose contribution from pencil beam scanning, scattered and spatially fractionated proton therapy

The Orsay Proton therapy Center (ICPO) has a long history of intracranial radiotherapy using both double scattering (DS) and pencil beam scanning (PBS) techniques, and is actively investigating a promising modality of spatially fractionated radiotherapy using proton minibeams (pMBRT). This work provides a comprehensive comparison of the organ-specific secondary neutron dose due to each of these treatment modalities, assessed using Monte Carlo (MC) algorithms and measurements. A MC model of a universal nozzle was benchmarked by comparing the neutron ambient dose equivalent,H*(10), in the gantry room with measurements obtained using a WENDI-II counter. The secondary neutron dose was evaluated for clinically relevant intracranial treatments of patients of different ages, in which secondary neutron doses were scored in anthropomorphic phantoms merged with the patients' images. The MC calculatedH*(10) values showed a reasonable agreement with the measurements and followed the expected tendency, in which PBS yields the lowest dose, followed by pMBRT and DS. Our results for intracranial treatments show that pMBRT yielded a higher secondary neutron dose for organs closer to the target volume, while organs situated furthest from the target volume received a greater quantity of neutrons from the passive scattering beam line. To the best of our knowledge, this is the first study to compare MC secondary neutron dose estimates in clinical treatments between these various proton therapy modalities and to realistically quantify the secondary neutron dose contribution of clinical pMBRT treatments. The method established in this study will enable epidemiological studies of the long-term effects of intracranial treatments at ICPO, notably radiation-induced second malignancies.

Future technological developments in proton therapy - A predicted technological breakthrough

In the current spectrum of cancer treatments, despite high costs, a lack of robust evidence based on clinical outcomes or technical and radiobiological uncertainties, particle therapy and in particular proton therapy (PT) is rapidly growing. Despite proton therapy being more than fifty years old (first proposed by Wilson in 1946) and more than 220,000 patients having been treated with in 2020, many technological challenges remain and numerous new technical developments that must be integrated into existing systems. This article presents an overview of on-going technical developments and innovations that we felt were most important today, as well as those that have the potential to significantly shape the future of proton therapy. Indeed, efforts have been done continuously to improve the efficiency of a PT system, in terms of cost, technology and delivery technics, and a number of different developments pursued in the accelerator field will first be presented. Significant developments are also underway in terms of transport and spatial resolution achievable with pencil beam scanning, or conformation of the dose to the target: we will therefore discuss beam focusing and collimation issues which are important parameters for the development of these techniques, as well as proton arc therapy. State of the art and alternative approaches to adaptive PT and the future of adaptive PT will finally be reviewed. Through these overviews, we will finally see how advances in these different areas will allow the potential for robust dose shaping in proton therapy to be maximised, probably foreshadowing a future era of maturity for the PT technique.

Metallic implants and CT artefacts in the CTV area: Where are we in 2020?

Radiation therapy (RT) is one of the main modalities of cancer treatment worldwide with computed tomography (CT), as the most commonly used imaging method for treatment planning system (TPS). Image reconstruction errors may greatly affect all the radiation therapy planning process, such as target delineation, dose calculation and delivery, particularly with particle therapy. Metallic implants, such as hip and spinal implants, and dental filling significantly deteriorate image quality. These hardware structures are often very complex in geometry leading to geometric complex artefacts in the clinical target volume (CTV) area, rendering the delineation of CTV challenging. In our review, we focus on the methods to overcome artefact consequences on CTV delineation: 1- medical approaches anticipating issues associated with imaging artefacts during preoperative multidisciplinary discussions while following standard recommendations; 2- common metal artefact reduction (MAR) methods such as manually override artefact regions, ballistics avoiding beam paths through implanted materials, megavoltage-CT (MVCT); 3- prospects with radiolucent implants, MAR algorithms and various methods of dual energy computed tomography (DECT). Despite substantial and broad evidence for their benefits, there is still no universal solution for cases involving implanted metallic devices. There is still a high need for research efforts to adapt technologies to our issue: "how do I accurately delineate the ideal CTV in a metal artefact area?"

Exploiting the full potential of proton therapy: An update on the specifics and innovations towards spatial or temporal optimisation of dose delivery

Prescription and delivery of protons are somewhat different compared to photons and may influence outcomes (tumour control and toxicity). These differences should be taken into account to fully exploit the clinical potential of proton therapy. Innovations in proton therapy treatment are also required to widen the therapeutic window and determine appropriate populations of patients that would benefit from new treatments. Therefore, strategies are now being developed to reduce side effects to critical normal tissues using alternative treatment configurations and new spatial or temporal-driven optimisation approaches. Indeed, spatiotemporal optimisation (based on flash, proton minibeam radiation therapy or hypofractionated delivery methods) has been gaining some attention in proton therapy as a mean of improving (biological and physical) dose distribution. In this short review, the main differences in planning and delivery between protons and photons, as well as some of the latest developments and methodological issues (in silico modelling) related to providing scientific evidence for these new techniques will be discussed.

A diamond guard ring microdosimeter for ion beam therapy

A single crystal chemical vapor deposition diamond-based microdosimeter prototype featuring an array of micro-sensitive volumes (μSVs) and surrounded by a so-called guard ring (GR) electrode has been fabricated using various microfabrication techniques available at Diamond Sensors Laboratory of CEA, Saclay. The GR microdosimeter was irradiated by a raster scanning method with 2 MeV proton microbeams. The charge transport properties of the GR sensor were determined with sub-micron spatial resolution by measuring the charge collection efficiency (CCE), the μSV geometry, and the pulse-height spectra. The response of the microdosimeter showed a well-defined and homogeneously active μSV. Appropriate biasing of the μSV structures led toward a full CCE for protons with lineal energies of ∼46 keV/μm. This shows the GR microdosimeter's great potential for applications in microdosimetry in clinical beam conditions.

First proton minibeam radiation therapy treatment plan evaluation

Proton minibeam radiation therapy (pMBRT) is a novel dose delivery method based on spatial dose fractionation. pMBRT has been shown to be promising in terms of reduced side effects and superior tumour control in high-grade glioma-bearing rats compared to standard irradiation. These findings, together with the recent optimized implementation of pMBRT in a clinical pencil beam scanning system, have triggered reflection on the possible application to patient treatments. In this context, the present study was designed to conduct a first theoretical investigation of the clinical potential of this technique. For this purpose, a dedicated dose engine was developed and used to evaluate two clinically relevant patient treatment plans (high-grade glioma and meningioma). Treatment plans were compared with standard proton therapy plans assessed by means of a commercial treatment planning system (ECLIPSE-Varian Medical systems) and Monte Carlo simulations. A multislit brass collimator consisting of 0.4 mm wide slits separated by a centre-to-centre distance of 4 or 6 mm was placed between the nozzle and the patient to shape the planar minibeams. For each plan, spread-out Bragg peaks and homogeneous dose distributions (±7% dose variations) can be obtained in target volumes. The Peak-to-Valley Dose Ratios (PVDR) were evaluated between 9.2 and 12.8 at a depth of 20 mm for meningioma and glioma, respectively. Dose volume histograms (DVHs) for target volumes and organs at risk were quantitatively compared, resulting in a slightly better target homogeneity with standard PT than with pMBRT plans, but similar DVHs for deep-seated organs-at-risk and lower average dose for shallow organs. The proposed delivery method evaluated in this work opens the way to an effective treatment for radioresistant tumours and will support the design of future clinical research.

Proton radiography with a commercial range telescope detector using dedicated post processing methods

Proton transmission imaging uses protons with high enough energy to fully traverse the phantom/patient and to be captured in a suitable detector placed behind it. The measured residual energy or residual range provide a direct estimate of the water equivalent thickness (WET) of the image volume. Requirements for proton imaging to be exploitable in clinical practice include: sufficient WET accuracy and integrability into the treatment room and the clinical workflow, as well as an acceptably low dose to the patient and a sufficient spatial resolution. In this work, we report on experiments performed at the Institut Curie-Proton therapy center in Orsay (IC-CPO), France, using a commercial range telescope commonly employed for quality assurance measurements. The purpose was to keep the experimental set-up as simple as possible and to achieve nonetheless high WET accuracy radiographies by developing and applying dedicated post processing methods. We explain these methods in detail and discuss their performance. We assess the WET accuracy based on two different reference phantoms: a CIRS electron density phantom with tissue equivalent inserts and a homogeneous step phantom. We find an agreement between the measured and the reference WET values of 0.2-0.5 mm. The lowest investigated dose was 10 mGy per acquisition. It could be lowered by modifying the irradiation plan and lowering the beam current, though the latter would impose slight optimisations of the detector hardware. Our work suggests that proton radiographies with good WET accuracy can be obtained with a reasonable experimental effort that would facilitate integration into clinical routine.

Mechanisms of phosphene generation in ocular proton therapy as related to space radiation exposure

Particle therapy provides an opportunity to study the human response to space radiation in ground-based facilities. On this basis, a study of light flashes analogous to astronauts' phosphenes reported by patients undergoing ocular proton therapy has been undertaken. The influence of treatment parameters on phosphene generation was investigated for 430 patients treated for a choroidal melanoma at the proton therapy centre of the Institut Curie (ICPO) in Orsay, France, between 2008 and 2011. 60% of them report light flashes, which are predominantly (74%) blue. An analysis of variables describing the patient's physiology, properties of the tumour and dose distribution shows that two groups of tumour and beam variables are correlated with phosphene occurrence. Physiology is found to have no influence on flash triggering. Detailed correlation study eventually suggests a possible twofold mechanism of phosphene generation based on (i) indirect Cerenkov light in the bulk of the eye due to nuclear interactions and radioactive decay and (ii) direct excitation of the nerve fibres in the back of the eye and/or radical excess near the retina.

Influence of beam incidence and irradiation parameters on stray neutron doses to healthy organs of pediatric patients treated for an intracranial tumor with passive scattering proton therapy

PURPOSE

In scattering proton therapy, the beam incidence, i.e. the patient's orientation with respect to the beam axis, can significantly influence stray neutron doses although it is almost not documented in the literature.

METHODS

MCNPX calculations were carried out to estimate stray neutron doses to 25 healthy organs of a 10-year-old female phantom treated for an intracranial tumor. Two beam incidences were considered in this article, namely a superior (SUP) field and a right lateral (RLAT) field. For both fields, a parametric study was performed varying proton beam energy, modulation width, collimator aperture and thickness, compensator thickness and air gap size.

RESULTS

Using a standard beam line configuration for a craniopharyngioma treatment, neutron absorbed doses per therapeutic dose of 63μGyGy(-1) and 149μGyGy(-1) were found at the heart for the SUP and the RLAT fields, respectively. This dose discrepancy was explained by the different patient's orientations leading to changes in the distance between organs and the final collimator where external neutrons are mainly produced. Moreover, investigations on neutron spectral fluence at the heart showed that the number of neutrons was 2.5times higher for the RLAT field compared against the SUP field. Finally, the influence of some irradiation parameters on neutron doses was found to be different according to the beam incidence.

CONCLUSION

Beam incidence was thus found to induce large variations in stray neutron doses, proving that this parameter could be optimized to enhance the radiation protection of the patient.

An empirical model for calculation of the collimator contamination dose in therapeutic proton beams

Collimators are used as lateral beam shaping devices in proton therapy with passive scattering beam lines. The dose contamination due to collimator scattering can be as high as 10% of the maximum dose and influences calculation of the output factor or monitor units (MU). To date, commercial treatment planning systems generally use a zero-thickness collimator approximation ignoring edge scattering in the aperture collimator and few analytical models have been proposed to take scattering effects into account, mainly limited to the inner collimator face component. The aim of this study was to characterize and model aperture contamination by means of a fast and accurate analytical model. The entrance face collimator scatter distribution was modeled as a 3D secondary dose source. Predicted dose contaminations were compared to measurements and Monte Carlo simulations. Measurements were performed on two different proton beam lines (a fixed horizontal beam line and a gantry beam line) with divergent apertures and for several field sizes and energies. Discrepancies between analytical algorithm dose prediction and measurements were decreased from 10% to 2% using the proposed model. Gamma-index (2%/1 mm) was respected for more than 90% of pixels. The proposed analytical algorithm increases the accuracy of analytical dose calculations with reasonable computation times.

Analytical linear energy transfer model including secondary particles: calculations along the central axis of the proton pencil beam

In proton therapy, the relative biological effectiveness (RBE) depends on various types of parameters such as linear energy transfer (LET). An analytical model for LET calculation exists (Wilkens' model), but secondary particles are not included in this model. In the present study, we propose a correction factor, L sec, for Wilkens' model in order to take into account the LET contributions of certain secondary particles. This study includes secondary protons and deuterons, since the effects of these two types of particles can be described by the same RBE-LET relationship. L sec was evaluated by Monte Carlo (MC) simulations using the GATE/GEANT4 platform and was defined by the ratio of the LET d distributions of all protons and deuterons and only primary protons. This method was applied to the innovative Pencil Beam Scanning (PBS) delivery systems and L sec was evaluated along the beam axis. This correction factor indicates the high contribution of secondary particles in the entrance region, with L sec values higher than 1.6 for a 220 MeV clinical pencil beam. MC simulations showed the impact of pencil beam parameters, such as mean initial energy, spot size, and depth in water, on L sec. The variation of L sec with these different parameters was integrated in a polynomial function of the L sec factor in order to obtain a model universally applicable to all PBS delivery systems. The validity of this correction factor applied to Wilkens' model was verified along the beam axis of various pencil beams in comparison with MC simulations. A good agreement was obtained between the corrected analytical model and the MC calculations, with mean-LET deviations along the beam axis less than 0.05 keV μm(-1). These results demonstrate the efficacy of our new correction of the existing LET model in order to take into account secondary protons and deuterons along the pencil beam axis.

Configuration and validation of an analytical model predicting secondary neutron radiation in proton therapy using Monte Carlo simulations and experimental measurements

PURPOSE

This study focuses on the configuration and validation of an analytical model predicting leakage neutron doses in proton therapy.

METHODS

Using Monte Carlo (MC) calculations, a facility-specific analytical model was built to reproduce out-of-field neutron doses while separately accounting for the contribution of intra-nuclear cascade, evaporation, epithermal and thermal neutrons. This model was first trained to reproduce in-water neutron absorbed doses and in-air neutron ambient dose equivalents, H*(10), calculated using MCNPX. Its capacity in predicting out-of-field doses at any position not involved in the training phase was also checked. The model was next expanded to enable a full 3D mapping of H*(10) inside the treatment room, tested in a clinically relevant configuration and finally consolidated with experimental measurements.

RESULTS

Following the literature approach, the work first proved that it is possible to build a facility-specific analytical model that efficiently reproduces in-water neutron doses and in-air H*(10) values with a maximum difference less than 25%. In addition, the analytical model succeeded in predicting out-of-field neutron doses in the lateral and vertical direction. Testing the analytical model in clinical configurations proved the need to separate the contribution of internal and external neutrons. The impact of modulation width on stray neutrons was found to be easily adjustable while beam collimation remains a challenging issue. Finally, the model performance agreed with experimental measurements with satisfactory results considering measurement and simulation uncertainties.

CONCLUSION

Analytical models represent a promising solution that substitutes for time-consuming MC calculations when assessing doses to healthy organs.

Use of gEUD for predicting ear and pituitary gland damage following proton and photon radiation therapy

OBJECTIVE

To determine the relationship between the dose to the inner ear or pituitary gland and radiation-induced late effects of skull base radiation therapy.

METHODS

140 patients treated between 2000 and 2008 were considered for this study. Hearing loss and endocrine dysfunction were retrospectively reviewed on pre- and post-radiation therapy audiometry or endocrine assessments. Two normal tissue complication probability (NTCP) models were considered (Lyman-Kutcher-Burman and log-logistic) whose parameters were fitted to patient data using receiver operating characteristics and maximum likelihood analysis. The method provided an estimation of the parameters of a generalized equivalent uniform dose (gEUD)-based NTCP after conversion of dose-volume histograms to equivalent doses.

RESULTS

All 140 patients had a minimum follow up of 26 months. 26% and 44% of patients experienced mild hearing loss and endocrine dysfunction, respectively. The fitted values for TD50 and γ50 ranged from 53.6 to 60.7 Gy and from 1.9 to 2.9 for the inner ear and were equal to 60.6 Gy and 4.9 for the pituitary gland, respectively. All models were ranked equal according to Akaike's information criterion.

CONCLUSION

Mean dose and gEUD may be used as predictive factors for late ear and pituitary gland late complications after skull base proton and photon radiation therapy.

ADVANCES IN KNOWLEDGE

In this study, we have reported mean dose effects and dose-response relationship of small organs at risk (partial volumes of the inner ear and pituitary gland), which could be useful to define optimal dose constraints resulting in an improved therapeutic ratio.

Secondary neutron doses received by paediatric patients during intracranial proton therapy treatments

This paper's goal is to assess secondary neutron doses received by paediatric patients treated for intracranial tumours using a 178 MeV proton beam. The MCNPX Monte Carlo model of the proton therapy facility, previously validated through experimental measurements for both proton and neutron dosimetry, was used. First, absorbed dose was calculated for organs located outside the clinical target volume using a series of hybrid computational phantoms for different ages and considering a realistic treatment plan. In general, secondary neutron dose was found to decrease as the distance to the treatment field increases and as the patient age increases. In addition, secondary neutron doses were studied as a function of the beam incidence. Next, neutron equivalent dose was assessed using organ-specific energy-dependent radiation weighting factors determined from Monte Carlo simulations of neutron spectra at each organ. The equivalent dose was found to reach a maximum value of ∼155 mSv at the level of the breasts for a delivery of 49 proton Gy to an intracranial tumour of a one-year-old female patient. Finally, a thorough comparison of the calculation results with published data demonstrated the dependence of neutron dose on the treatment configuration and proved the need for facility-specific and treatment-dependent neutron dose calculations.

Flake size-dependent cyto and genotoxic evaluation of graphene oxide on in vitro A549, CaCo2 and vero cell lines

This study was carried out by varying both graphene oxide (GO) concentration (10 μg/mL, 50 μg/mL, 100 μg/mL) and flakes sizes of 1320 nm and 130 nm. Characterization by scanning electron microscopy and Raman spectroscopy demonstrate that the area of GO flakes varies of one order of magnitude but their chemical structure remains unmodified. A 24-h cytotoxicity test showed, for A549, a loss in the viability, while the test exhibits overall a positive increase in the viability for CaCo2 and Vero. A 24-h comet assay shows a marked GO genotoxicity: for micrometer-sized GO flakes the genotoxicity is in positive correlation with the concentration, while for nanometer-sized GO flakes there was a high degree of genotoxicity at the lowest concentration tested.

Calibration of CT Hounsfield units for proton therapy treatment planning: use of kilovoltage and megavoltage images and comparison of parameterized methods

Proton beam range is of major concern, in particular, when images used for dose computations are artifacted (for example in patients with surgically treated bone tumors). We investigated several conditions and methods for determination of computed tomography Hounsfield unit (CT-HU) calibration curves, using two different conversion schemes. A stoichiometric methodology was used on either kilovoltage (kV) or megavoltage (MV) CT images and the accuracy of the calibration methods was evaluated. We then studied the effects of metal artifacts on proton dose distributions using metallic implants in rigid phantom mimicking clinical conditions. MV-CT images were used to evaluate relative proton stopping power in certain high density implants, and a methodology is proposed for accurate delineation and dose calculation, using a combined set of kV- and MV-CT images. Our results show good agreement between measurements and dose calculations or relative proton stopping power determination (<5%). The results also show that range uncertainty increases when only kV-CT images are used or when no correction is made on artifacted images. However, differences between treatment plans calculated on corrected kV-CT data and MV-CT data remained insignificant in the investigated patient case, even with streak artifacts and volume effects that reduce the accuracy of manual corrections.

Characterization of Bordetella pertussis strains of recent isolation

During the clinical trial conducted in Italy to evaluate the efficacy of new acellular pertussis vaccines, the most favorable conditions for the recovery and characterization of Bordetella pertussis strains, isolated from children with cough, were adopted. The nasopharyngeal aspirates were collected and sent to the laboratory in refrigerated conditions within 24 hours. Charcoal agar selective and non selective plates were used, and most of the isolates were recovered after 3-4 days of incubation. Confirmation of all suspected colonies included the use of biochemical tests and specific agglutination reaction with B. pertussis antiserum. Serotyping of fimbriae, susceptibility to erythromycin and DNA fingerprinting by Pulsed Field Gel Electrophoresis (PFGE), were performed to characterize B. pertussis isolates and to determine relatedness among different strains. Serotype 1,3 was the most represented in the bacterial population examined. A predominant pulsetype (PTA) characterized most of the isolates accounting for 71.4% of the strains examined. Eight subclones (23.5%) and three unrelated pulsetypes were also found. No resistant strains to erythromycin were detected.

Bordetella parapertussis infections

The rate of isolation of Bordetella parapertussis among children with cough during the follow-up of different clinical efficacy studies has been evaluated. In the Italian trial, a comparison of clinical characteristics between B. pertussis and B. parapertussis infections showed lower frequencies and shorter duration of typical symptoms of whooping cough such as paroxysmal coughing, whooping, and vomiting in the group of children affected with B. parapertussis infections. In about 70% of B. parapertussis infections, there was a two-fold increase of IgA or IgG anti-FHA from acute- and convalescent-phase serum specimens. The analysis of the distribution of B. parapertussis cases in children fully immunized with each pertussis vaccine suggested that vaccination is irrelevant in preventing B. parapertussis infection.

Immunity status to polioviruses among non-European union immigrants in Veneto Region (North-East Italy)

Widespread use of poliovirus vaccines has lead to a dramatic decline on the incidence of paralytic poliomyelitis in the Western Hemisphere. In Italy, the use of live attenuated vaccine (OPV) has lead to virtual eradication of poliomyelitis. However, wild poliovirus is still endemic in some less developed countries, and there is still a risk of importation of wild-virus poliomyelitis as result of immigration from these areas. Therefore, a serologic study of the immunologic status to poliovirus was carried out in 242 immigrants from less developed countries in the Veneto Region (North-East Italy). The overall prevalence of serum neutralizing antibodies was 98.3 percent for poliovirus type 1, 99.6 percent for type 2, and 95.9 percent for type 3. The modal titres were 1:128 (21.1 percent), 1:512 (26.4 percent), and 1:32 (22.3 percent) for poliovirus type 1, type 2, and type 3, respectively. No subject was simultaneously lacking neutralizing antibodies to all three polioviruses. The mean duration of stay in Italy before testing was 1.53 years (range 1 month to 7 years). These results show a good level of the immunity status, similar to those found in the Veneto general population. Although the risk of paralytic poliomyelitis in our population is very low, seronegative immigrants originating from areas of high endemicity could reintroduce wild poliovirus. It is therefore important to maintain a careful surveillance system on the importation of wild poliovirus and to plan vaccination programs for unprotected immigrants.