Design of spread-out Bragg peaks in hadron therapy with oxygen ions

Aim

Design of a numerical method for creating spread-out Bragg peak (SOBP) and evaluation of the best parameter in Bortfeld Model to this aim in oxygen ion therapy.

Accurate proton treatment planning for pencil beam crossing titanium fixation implants

PurposeTo present a planning strategy for proton pencil-beam scanning when titanium implants need to be crossed by the beam.MethodsWe addressed three issues: the implementation of a CT calibration curve to assign to titanium the correct stopping power; the effect of artefacts on CT images and their reduction by a dedicated algorithm; the differences in dose computation depending on the dose engine, pencil-beam vs Monte-Carlo algorithms. We performed measurement tests on a simple cylinder phantom and on a real implant. These phantoms were irradiated with three geometries (single spots, uniform mono-energetic layer and uniform box), measuring the exit dose either by radio-chromic film or multi-layer ionization chamber. The procedure was then applied on two patients treated for chordoma.ResultsWe had to set in the calibration curve a mass density equal to 4.37 g/cm³ to saturated Hounsfield Units, in order to have the correct stopping power assigned to titanium in TPS. CT artefact reduction algorithm allowed a better reconstruction of the shape and size of the implant. Monte-Carlo resulted accurate in computing the dose distribution whereas the pencil-beam algorithm failed due to sharp density interfaces between titanium and the surrounding material. Finally, the treatment plans obtained on two patients showed the impact of the dose engine algorithm, with 10–20% differences between pencil-beam and Monte-Carlo in small regions distally to the titanium screws.ConclusionThe described combination of CT calibration, artefacts reduction and Monte-Carlo computation provides a reliable methodology to compute dose in patients with titanium implants.     

Out-of-field doses from secondary radiation produced in proton therapy and the associated risk of radiation-induced cancer from a brain tumor treatment

PurposeTo determine out-of-field doses produced in proton pencil beam scanning (PBS) therapy using Monte Carlo simulations and to estimate the associated risk of radiation-induced second cancer from a brain tumor treatment.MethodsSimulations of out-of-field absorbed doses were performed with MCNP6 and benchmarked against measurements with tissue-equivalent proportional counters (TEPC) for three irradiation setups: two irradiations of a water phantom using proton energies of 78–147 MeV and 177–223 MeV, and one brain tumor irradiation of a whole-body phantom. Out-of-field absorbed and equivalent doses to organs in a whole-body phantom following a brain tumor treatment were subsequently simulated and used to estimate the risk of radiation-induced cancer. Additionally, the contribution of absorbed dose originating from radiation produced in the nozzle was calculated from simulations.ResultsOut-of-field absorbed doses to the TEPC ranged from 0.4 to 135 µGy/Gy. The average deviation between simulations and measurements of the water phantom irradiations was about 17%. The absorbed dose contribution from radiation produced in the nozzle ranged between 0 and 70% of the total dose; the contribution was however small in absolute terms. The absorbed and equivalent doses to the organs ranged between 0.2 and 60 µGy/Gy and 0.5–151 µSv/Gy. The estimated lifetime risk of radiation-induced second cancer was approximately 0.01%.ConclusionsThe agreement of out-of-field absorbed doses between measurements and simulations was good given the sources of uncertainties. Calculations of out-of-field organ doses following a brain tumor treatment indicated that proton PBS therapy of brain tumors is associated with a low risk of radiation-induced cancer.     

Medical physics aspects of the synchrotron radiation therapies: Microbeam radiation therapy (MRT) and synchrotron stereotactic radiotherapy (SSRT)

Stereotactic Synchrotron Radiotherapy (SSRT) and Microbeam Radiation Therapy (MRT) are both novel approaches to treat brain tumor and potentially other tumors using synchrotron radiation. Although the techniques differ by their principles, SSRT and MRT share certain common aspects with the possibility of combining their advantages in the future. For MRT, the technique uses highly collimated, quasi-parallel arrays of X-ray microbeams between 50 and 600 keV. Important features of highly brilliant Synchrotron sources are a very small beam divergence and an extremely high dose rate. The minimal beam divergence allows the insertion of so called Multi Slit Collimators (MSC) to produce spatially fractionated beams of typically ∼25–75 micron-wide microplanar beams separated by wider (100–400 microns center-to-center(ctc)) spaces with a very sharp penumbra. Peak entrance doses of several hundreds of Gy are extremely well tolerated by normal tissues and at the same time provide a higher therapeutic index for various tumor models in rodents. The hypothesis of a selective radio-vulnerability of the tumor vasculature versus normal blood vessels by MRT was recently more solidified.SSRT (Synchrotron Stereotactic Radiotherapy) is based on a local drug uptake of high-Z elements in tumors followed by stereotactic irradiation with 80 keV photons to enhance the dose deposition only within the tumor. With SSRT already in its clinical trial stage at the ESRF, most medical physics problems are already solved and the implemented solutions are briefly described, while the medical physics aspects in MRT will be discussed in more detail in this paper.     

Beam and tissue factors affecting Cherenkov image intensity for quantitative entrance and exit dosimetry on human tissue

This study's goal was to determine how Cherenkov radiation emission observed in radiotherapy is affected by predictable factors expected in patient imaging. Factors such as tissue optical properties, radiation beam properties, thickness of tissues, entrance/exit geometry, curved surface effects, curvature and imaging angles were investigated through Monte Carlo simulations. The largest physical cause of variation of the correlation ratio between of Cherenkov emission and dose was the entrance/exit geometry (?50%). The largest human tissue effect was from different optical properties (?45%). Beyond these, clinical beam energy varies the correlation ratio significantly (?20% for X‐ray beams), followed by curved surfaces (?15% for X‐ray beams and ?8% for electron beams), and finally, the effect of field size (?5% for X‐ray beams). Other investigated factors which caused variations less than 5% were tissue thicknesses and source to surface distance. The effect of non‐Lambertian emission was negligible for imaging angles smaller than 60 degrees. The spectrum of Cherenkov emission tends to blue‐shift along the curved surface. A simple normalization approach based on the reflectance image was experimentally validated by imaging a range of tissue phantoms, as a first order correction for different tissue optical properties.

     

Minibeam radiation therapy at a conventional irradiator: Dose-calculation engine and first tumor-bearing animals irradiation

Purpose: Minibeam radiation therapy (MBRT) is a novel therapeutic strategy, whose exploration was hindered due to its restriction to large synchrotrons. Our recent implementation of MBRT in a wide-spread small animal irradiator offers the possibility of performing systematic radiobiological studies. The aim of this research was to develop a set of dosimetric tools to reliably guide biological experiments in the irradiator.Methods: A Monte Carlo (Geant4)-based dose calculation engine was developed. It was then benchmarked against a series of dosimetric measurements performed with gafchromic films. Two voxelized rat phantoms (ROBY, computer tomography) were used to evaluate the treatment plan of F98 tumor-bearing rats. The response of a group of 7 animals receiving a unilateral irradiation of 58 Gy was compared to a group of non-irradiated controls.Results: The good agreement between calculations and the experimental data allowed the validation of the dose-calculation engine. The latter was first used to compare the dose distributions in computer tomography images of a rat’s head and in a digital model of a rat’s head (ROBY), obtaining a good general agreement. Finally, with respect to the in vivo experiment, the increase of mean survival time of the treated group with respect to the controls was modest but statistically significant.Conclusions: The developed dosimetric tools were used to reliably guide the first MBRT treatments of intracranial glioma-bearing rats outside synchrotrons. The significant tumor response obtained with respect to the non-irradiated controls, despite the heterogenous dose coverage of the target, might indicate the participation of non-targeted effects.     

A review on gold nanoparticles radiosensitization effect in radiation therapy of cancer

In the recent years, application of nanoparticles in diagnosis and treatment of cancer has been the issue of extensive research. Among these studies some have focused on the dose enhancement effect of gold nanoparticles (GNPs) in radiation therapy of cancer. On the other hand, some studies indicated energy dependency of dose enhancement effect, and the others have studied the GNP size effect in association with photon energy. However, in some aspects of GNP-based radiotherapy the results of recent studies do not seem very conclusive in spite of relative agreement on the basic physical interaction of photoelectric between GNPs and low energy photons. The main idea behind the GNP dose enhancement in some studies is not able to explain the results especially in recent investigation on cell lines and animal models radiation therapy using GNPs. In the present article the results of the available reports and articles were analyzed and compared and the final status of the GNP-RT was discussed.     

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

PurposeThis study focuses on the configuration and validation of an analytical model predicting leakage neutron doses in proton therapy.MethodsUsing 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.ResultsFollowing 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.ConclusionAnalytical models represent a promising solution that substitutes for time-consuming MC calculations when assessing doses to healthy organs.     

Modifying the tumour microenvironment and reverting tumour cells: New strategies for treating malignant tumours

The tumour microenvironment (TME) plays a pivotal role in tumour fate determination. The TME acts together with the genetic material of tumour cells to determine their initiation, metastasis and drug resistance. Stromal cells in the TME promote the growth and metastasis of tumour cells by secreting soluble molecules or exosomes. The abnormal microenvironment reduces immune surveillance and tumour killing. The TME causes low anti‐tumour drug penetration and reactivity and high drug resistance. Tumour angiogenesis and microenvironmental hypoxia limit the drug concentration within the TME and enhance the stemness of tumour cells. Therefore, modifying the TME to effectively attack tumour cells could represent a comprehensive and effective anti‐tumour strategy. Normal cells, such as stem cells and immune cells, can penetrate and disrupt the abnormal TME. Reconstruction of the TME with healthy cells is an exciting new direction for tumour treatment. We will elaborate on the mechanism of the TME to support tumours and the current cell therapies for targeting tumours and the TME—such as immune cell therapies, haematopoietic stem cell (HSC) transplantation therapies, mesenchymal stem cell (MSC) transfer and embryonic stem cell‐based microenvironment therapies—to provide novel ideas for producing breakthroughs in tumour therapy strategies.     

Optimization of an ultra-fast silicon detector for proton and carbon beams using GEANT4 Monte Carlo toolkit

AimThe purpose of this study was to investigate the crosstalk effects between adjacent pixels in a thin silicon detector with 50 um thickness.BackgroundThere are some limitations in the applications of detectors in hadron therapy. So it is necessary to have a detector with concurrent excellent time and resolution. In this work, the GEANT4 toolkit was applied to estimate the best value for energy cutoff in the thin silicon detector in order to optimize the detector.Materials and MethodsGEANT4 toolkit was applied to simulate the transport and interactions of particles. Calculations were performed for a thin silicon detector (2 cm × 2 cm×0.005 cm) irradiated by proton and carbon ion beams. A two-dimensional array of silicon pixels in the x-y plane with 100 um × 100 um × 50 um dimensions build the whole detector. In the end, the ROOT package is used to interpret and analyze the resultsResultsIt is seen that by the presence of energy cutoff, pixels with small deposited energy are ignored. The best values for energy cutoff are 0.01 MeV and 0.7 MeV for proton and carbon ion beams, respectively. By applying these energy cutoff values, efficiency and purity values are maximized and also minimum output errors are achieved.ConclusionsThe results are reasonable, good and useful to optimize the geometry of future silicon detectors in order to be used as beam monitoring in hadron therapy applications.     

Use of Monte Carlo simulations for improved facility fit planning in downstream biomanufacturing and technology transfer

Biologics manufacturing is capital and consumable intensive with need for advanced inventory planning to account for supply chain constraints. Early-stage process design and technology transfer are often challenging due to limited information on process variability regarding bioreactor titer, process yield, and product quality. Monte Carlo (MC) methods offer a stochastic modeling approach for process optimization where probabilities of occurrence for process inputs are incorporated into a deterministic model to simulate more likely scenarios for process outputs. In this study, we explore MC simulation-based design of a monoclonal antibody downstream manufacturing process. We demonstrate that this probabilistic approach offers more representative outcomes over the conventional worst-case approach where the theoretical minimum and maximum values of each process parameter are used without consideration for their probability of occurrence. Our work demonstrates case studies on more practically sizing unit operations to improve consumable utilization, thereby reducing manufacturing costs. We also used MC simulations to minimize process cadence by constraining the number of cycles per unit operation to fit facility preferences. By factoring in process uncertainty, we have implemented MC simulation-based facility fit analyses to efficiently plan for inventory when accounting for process constraints during technology transfer from lab-scale to clinical or commercial manufacturing.     

Intelligence-guided beam angle optimization in treatment planning of intensity-modulated radiation therapy

An intelligence guided approach based on fuzzy inference system (FIS) was proposed to automate beam angle optimization in treatment planning of intensity-modulated radiation therapy (IMRT). The model of FIS is built on inference rules in describing the relationship between dose quality of IMRT plan and irradiated region of anatomical structure. Dose quality of IMRT plan is quantified by the difference between calculated and constraint doses of the anatomical structures in an IMRT plan. Irradiated region of anatomical structure is characterized by the metric, covered region of interest, which is the region of an anatomical structure under radiation field while beam’s eye-view is conform to target volume. Initially, an IMRT plan is created with a single beam. The dose difference is calculated for the input of FIS and the output of FIS is obtained with processing of fuzzy inference. Later, a set of candidate beams is generated for replacing the current beam. This process continues until no candidate beams is found. Then the next beam is added to the IMRT plan and optimized in the same way as the previous beam. The new beam keeps adding to the IMRT plan until the allowed beam number is reached. Two spinal cases were investigated in this study. The preliminary results show that dose quality of IMRT plans achieved by this approach is better than those achieved by the default approach with equally spaced beam setting. It is effective to find the optimal beam combination of IMRT plan with the intelligence-guided approach.     

Dosimetric impact of different multileaf collimators on prostate intensity modulated treatment planning

AimThe main purpose of this study is to perform a dosimetric comparison on target volumes and organs at risks (OARs) between prostate intensity modulated treatment plans (IMRT) optimized with different multileaf collimators (MLCs).BackgroundThe use of MLCs with a small leaf width in the IMRT optimization may improve conformity around the tumor target whilst reducing the dose to normal tissues.Materials and methodsTwo linacs mounting MLCs with 5 and 10 mm leaf-width, respectively, implemented in Pinnacle³ treatment planning system were used for this work. Nineteen patients with prostate carcinoma undergoing a radiotherapy treatment were enrolled. Treatment planning with different setup arrangements (7 and 5 beams) were performed for each patient and each machine. Dose volume histograms (DVHs) cut-off points were used in the treatment planning comparison.ResultsComparable planning target volume (PTV) coverage was obtained with 7- and 5-beam configuration (both with 5 and 10 mm MLC leaf-width). The comparison of bladder and rectum DVH cut-off points for the 5-beam arrangement shows that 52.6% of the plans optimized with a larger leaf-width did not satisfy at least one of the OARs’ constraints. This percentage is reduced to 10.5% for the smaller leaf-width. If a 7-beam arrangement is used the value of 52.6% decreases to 21.1% while the value of 10.5% remains unchanged.ConclusionMLCs collimators with different widths and number of leaves lead to a comparable prostate treatment planning if a proper adjustment is made of the number of gantry angles.     

Localized interleukin‐12 delivery for immunotherapy of solid tumours

Interleukin (IL)‐12 is the key cytokine in the initiation of a Th1 response and has shown promise as an anti‐cancer agent; however, clinical trials involving IL‐12 have been unsuccessful due to toxic side‐effects. To address this issue, lentiviral vectors were used to transduce tumour cell lines that were injected as an autologous tumour cell vaccine. The focus of the current study was to test the efficacy of this approach in a solid tumour model. SCCVII cells that were transduced to produce IL‐12 at different concentrations were then isolated. Subcutaneous injection of parental SCCVII cells results in tumour development, while a mixture of IL‐12‐producing and non‐producing cells results in tumour clearance. Interestingly, when comparing mice injected a mixture of SCCVII and either high IL‐12‐producing tumour cells or low IL‐12‐producing tumour cells, we observed that mixtures containing small amounts of high producing cells lead to tumour clearance, whereas mixtures containing large amounts of low producing cells fail to elicit protection, despite the production of equal amounts of total IL‐12 in both mixtures. Furthermore, immunizing mice with IL‐12‐producing cells leads to the establishment of both local and systemic immunity against challenge with SCCVII. Using depletion antibodies, it was shown that both CD4⁺ and CD8⁺ cells are crucial for therapy. Lastly, we have established cell clones of other solid tumour cell lines (RM‐1, LLC1 and moto1.1) that produce IL‐12. Our results show that the delivery of IL‐12 by cancer cells is an effective route for immune activation.     

FLUKA simulation of target fragmentation in proton therapy

In proton therapy, secondary fragments are created in nuclear interactions of the beam with the target nuclei. The secondary fragments have low kinetic energies and high atomic numbers as compared to primary protons. Fragments have a high LET and deposit all their energy close to the generation point. For their characteristics, secondary fragments can alter the dose distribution and lead to an increase of RBE for the same delivered physical dose. Moreover, the radiobiological impact of target fragmentation is significant mostly in the region before the Bragg peak, where generally healthy tissues are present, and immediately after Bragg peak. Considering the high biological impact of those particles, especially in the case of healthy tissues or organs at risk, the inclusion of target fragmentation processes in the dose calculation of a treatment planning system can be relevant to improve the treatment accuracy and for this reason it is one of the major tasks of the MoVe IT project.In this study, Monte Carlo simulations were employed to fully characterize the mixed radiation field generated by target fragmentation in proton therapy. The dose averaged LET has been evaluated in case of a Spread Out Bragg Peak (SOBP). Starting from LET distribution, RBE has been evaluated with two different phenomenological models. In order to characterize the mixed radiation field, the production cross section has been evaluated by means of the FLUKA code. The future development of present work is to generate a MC database of fragments fluence to be included in TPS.     

Static and rotational intensity modulated techniques for head-neck cancer radiotherapy: A planning comparison

PurposeTo compare helical Tomotherapy (HT), two volumetric-modulated arc techniques and conventional fixed-field intensity modulated techniques (S-IMRT) for head-neck (HN) cancers.Methods and materialsEighteen HN patients were considered. Four treatment plans were generated for each patient: HT, S-IMRT optimised with Eclipse treatment planning system and two volumetric techniques using Elekta–Oncentra approach (VMAT) and Varian-RapidArc (RA), using two full arcs. All techniques were optimised to simultaneously deliver 66Gy to PTV1 (GTV and enlarged nodes) and 54Gy to PTV2 (subclinical and electively treated nodes). Comparisons were assessed on several dosimetric parameters and, secondarily, on planned MUs and delivery time.ResultsConcerning PTV coverage, significantly better results were found for HT and RA. HT significantly improved the target coverage both compared to S-IMRT and VMAT. No significant differences were found between S-IMRT and volumetric techniques in terms of dose homogeneity. For OARs, all the techniques were able to satisfy all hard constraints; significantly better results were found for HT, especially in the intermediate dose range (15–30 Gy). S-IMRT reached a significantly better OARs sparing with respect to VMAT and RA. No significant differences were found for body mean dose, excepting higher values of V5–V10 for HT. A reduction of planned MUs and delivery treatment time was found with volumetric techniques.ConclusionsThe objectives of satisfying target coverage and sparing of critical structures were reached with all techniques. S-IMRT techniques were found more advantageous compared to RA and VMAT for OARs sparing. HT reached the best overall treatment plan quality.