Toshiyuki Toshito

Treatment planning for a scanned carbon beam with a modified microdosimetric kinetic model

We describe a method to calculate the relative biological effectiveness in mixed radiation fields of therapeutic ion beams based on the modified microdosimetric kinetic model (modified MKM). In addition, we show the procedure for integrating the modified MKM into a treatment planning system for a scanned carbon beam. With this procedure, the model is fully integrated into our research version of the treatment planning system. To account for the change in radiosensitivity of a cell line, we measured one of the three MKM parameters from a single survival curve of the current cells and used the parameter in biological optimization. Irradiation of human salivary gland tumor cells was performed with a scanned carbon beam in the Heavy Ion Medical Accelerator in Chiba (HIMAC), and we then compared the measured depth-survival curve with the modified MKM predicted survival curve. Good agreement between the two curves proves that the proposed method is a candidate for calculating the biological effects in treatment planning for ion irradiation.

Evaluation of hybrid depth scanning for carbon-ion radiotherapy

PURPOSE In radiotherapy with a scanned carbon-ion beam, its Bragg peak is shifted along the depth direction either by inserting the range shifter plates or by changing the beam-extraction energy of a synchrotron. In the former technique (range shifter scanning: RS), the range shifter plates broaden the beam size and produce secondary fragments through nuclear reactions. In the latter technique (active-energy scanning: ES), it may take several seconds to change the beam energy depending on the synchrotron operation cycle, leading to a long treatment time. The authors propose a hybrid depth scan technique (hybrid scanning: HS), where several beam energies are used in conjunction with the range shifter plates for finer range shift. In this study, HS is evaluated from the viewpoints of dose distribution and treatment time. METHODS Assuming realistic accelerator and beam-delivery systems, the authors performed computer simulations using GEANT4 Monte Carlo code for beam modeling and a treatment planning system to evaluate HS. Three target volumes with the same dimensions of 60 × 60 × 60 mm(3) were generated at depths of 45, 85, and 125 mm in water phantom, and uniform clinical dose was planned for these targets. The sizes of lateral dose falloff and the peak to plateau ratio defined as the ratio of the clinical dose averaged over the target to the clinical dose at the entrance as well as the treatment time were compared among the three depth scan techniques. RESULTS The sizes of lateral dose falloffs at the center of SOBP are 11.4, 8.5, and 5.9 mm for the three targets in RS, while they are 5.7, 4.8, and 4.6 mm in ES and 6.6, 5.7, and 5.0 mm in HS, respectively. The peak to plateau ratios are 1.39, 1.96, and 2.15 in RS, while they are 1.48, 2.04, and 2.19 in ES and 1.47, 2.03, and 2.18 in HS, respectively. The treatment times are 128.7, 128.6, and 128.6 s in ES, while they are 61.2, 54.6, and 47.8 s in RS and 43.2, 44.1, and 44.7 s in HS, respectively. The multiple scattering and the nuclear reaction by range shifter degraded the beam qualities such as lateral dose falloff and peak to plateau ratio, which was especially pronounced for the shallow target in RS. The depth scan timing was limited by accelerator cycle in ES. That increased the treatment time by a few times. CONCLUSIONS This study revealed that HS can provide dose distributions with steeper lateral dose falloffs and higher peak to plateau ratio comparing to RS and comparable to ES. In addition, the treatment time can be considerably reduced in HS compared to ES.

Compatibility of the repairable-conditionally repairable, multi-target and linear-quadratic models in converting hypofractionated radiation doses to single doses

We investigated the applicability of the repairable-conditionally repairable (RCR) model and the multi-target (MT) model to dose conversion in high-dose-per-fraction radiotherapy in comparison with the linear-quadratic (LQ) model. Cell survival data of V79 and EMT6 single cells receiving single doses of 2–12 Gy or 2 or 3 fractions of 4 or 5 Gy each, and that of V79 spheroids receiving single doses of 5–26 Gy or 2–5 fractions of 5–12 Gy, were analyzed. Single and fractionated doses to actually reduce cell survival to the same level were determined by a colony assay. Single doses used in the experiments and surviving fractions at the doses were substituted into equations of the RCR, MT and LQ models in the calculation software Mathematica, and each parameter coefficient was computed. Thereafter, using the coefficients and the three models, equivalent single doses for the hypofractionated doses were calculated. They were then compared with actually-determined equivalent single doses for the hypofractionated doses. The equivalent single doses calculated using the RCR, MT and LQ models tended to be lower than the actually determined equivalent single doses. The LQ model seemed to fit relatively well at doses of 5 Gy or less. At 6 Gy or higher doses, the RCR and MT models seemed to be more reliable than the LQ model. In hypofractionated stereotactic radiotherapy, the LQ model should not be used, and conversion models incorporating the concept of the RCR or MT models, such as the generalized linear-quadratic models, appear to be more suitable.

Luminescence imaging of water during proton-beam irradiation for range estimation

PURPOSE Proton therapy has the ability to selectively deliver a dose to the target tumor, so the dose distribution should be accurately measured by a precise and efficient method. The authors found that luminescence was emitted from water during proton irradiation and conjectured that this phenomenon could be used for estimating the dose distribution. METHODS To achieve more accurate dose distribution, the authors set water phantoms on a table with a spot scanning proton therapy system and measured the luminescence images of these phantoms with a high-sensitivity, cooled charge coupled device camera during proton-beam irradiation. The authors imaged the phantoms of pure water, fluorescein solution, and an acrylic block. RESULTS The luminescence images of water phantoms taken during proton-beam irradiation showed clear Bragg peaks, and the measured proton ranges from the images were almost the same as those obtained with an ionization chamber. Furthermore, the image of the pure-water phantom showed almost the same distribution as the tap-water phantom, indicating that the luminescence image was not related to impurities in the water. The luminescence image of the fluorescein solution had ∼3 times higher intensity than water, with the same proton range as that of water. The luminescence image of the acrylic phantom had a 14.5% shorter proton range than that of water; the proton range in the acrylic phantom generally matched the calculated value. The luminescence images of the tap-water phantom during proton irradiation could be obtained in less than 2 s. CONCLUSIONS Luminescence imaging during proton-beam irradiation is promising as an effective method for range estimation in proton therapy.

A proton therapy system in Nagoya Proton Therapy Center

The purpose of this paper is to describe an outline of a proton therapy system in Nagoya Proton Therapy Center (NPTC). The NPTC has a synchrotron with a linac injector and three treatment rooms: two rooms are equipped with a gantry and the other one is equipped with a fixed horizontal beamline. One gantry treatment room has a pencil beam scanning treatment delivery nozzle. The other two treatment rooms have a passive scattering treatment delivery nozzle. In the scanning treatment delivery nozzle, an energy absorber and an aperture system to treat head and neck cancer have been equipped. In the passive treatment delivery nozzle, a multi-leaf collimator is equipped. We employ respiratory gating to treat lung and liver cancers for passive irradiation. The proton therapy system passed all acceptance tests. The first patient was treated on February 25, 2013, using passive scattering fixed beams. Respiratory gating is commonly used to treat lung and liver cancers in the passive scattering system. The MLCs are our first choice to limit the irradiation field. The use of the aperture for scanning irradiation reduced the lateral fall off by half or less. The energy absorber and aperture system in scanning delivery is beneficial to treat head and neck cancer.

Luminescence imaging of water during carbon-ion irradiation for range estimation: Luminescence imaging of water during carbon-ion irradiation

PURPOSE The authors previously reported successful luminescence imaging of water during proton irradiation and its application to range estimation. However, since the feasibility of this approach for carbon-ion irradiation remained unclear, the authors conducted luminescence imaging during carbon-ion irradiation and estimated the ranges. METHODS The authors placed a pure-water phantom on the patient couch of a carbon-ion therapy system and measured the luminescence images with a high-sensitivity, cooled charge-coupled device camera during carbon-ion irradiation. The authors also carried out imaging of three types of phantoms (tap-water, an acrylic block, and a plastic scintillator) and compared their intensities and distributions with those of a phantom containing pure-water. RESULTS The luminescence images of pure-water phantoms during carbon-ion irradiation showed clear Bragg peaks, and the measured carbon-ion ranges from the images were almost the same as those obtained by simulation. The image of the tap-water phantom showed almost the same distribution as that of the pure-water phantom. The acrylic block phantoms luminescence image produced seven times higher luminescence and had a 13% shorter range than that of the water phantoms; the range with the acrylic phantom generally matched the calculated value. The plastic scintillator showed ∼15 000 times higher light than that of water. CONCLUSIONS Luminescence imaging during carbon-ion irradiation of water is not only possible but also a promising method for range estimation in carbon-ion therapy.

Evaluation of the influence of double and triple Gaussian proton kernel models on accuracy of dose calculations for spot scanning technique.

PURPOSE The main purpose in this study was to present the results of beam modeling and how the authors systematically investigated the influence of double and triple Gaussian proton kernel models on the accuracy of dose calculations for spot scanning technique. METHODS The accuracy of calculations was important for treatment planning software (TPS) because the energy, spot position, and absolute dose had to be determined by TPS for the spot scanning technique. The dose distribution was calculated by convolving in-air fluence with the dose kernel. The dose kernel was the in-water 3D dose distribution of an infinitesimal pencil beam and consisted of an integral depth dose (IDD) and a lateral distribution. Accurate modeling of the low-dose region was important for spot scanning technique because the dose distribution was formed by cumulating hundreds or thousands of delivered beams. The authors employed a double Gaussian function as the in-air fluence model of an individual beam. Double and triple Gaussian kernel models were also prepared for comparison. The parameters of the kernel lateral model were derived by fitting a simulated in-water lateral dose profile induced by an infinitesimal proton beam, whose emittance was zero, at various depths using Monte Carlo (MC) simulation. The fitted parameters were interpolated as a function of depth in water and stored as a separate look-up table. These stored parameters for each energy and depth in water were acquired from the look-up table when incorporating them into the TPS. The modeling process for the in-air fluence and IDD was based on the method proposed in the literature. These were derived using MC simulation and measured data. The authors compared the measured and calculated absolute doses at the center of the spread-out Bragg peak (SOBP) under various volumetric irradiation conditions to systematically investigate the influence of the two types of kernel models on the dose calculations. RESULTS The authors investigated the difference between double and triple Gaussian kernel models. The authors found that the difference between the two studied kernel models appeared at mid-depths and the accuracy of predicting the double Gaussian model deteriorated at the low-dose bump that appeared at mid-depths. When the authors employed the double Gaussian kernel model, the accuracy of calculations for the absolute dose at the center of the SOBP varied with irradiation conditions and the maximum difference was 3.4%. In contrast, the results obtained from calculations with the triple Gaussian kernel model indicated good agreement with the measurements within ±1.1%, regardless of the irradiation conditions. CONCLUSIONS The difference between the results obtained with the two types of studied kernel models was distinct in the high energy region. The accuracy of calculations with the double Gaussian kernel model varied with the field size and SOBP width because the accuracy of prediction with the double Gaussian model was insufficient at the low-dose bump. The evaluation was only qualitative under limited volumetric irradiation conditions. Further accumulation of measured data would be needed to quantitatively comprehend what influence the double and triple Gaussian kernel models had on the accuracy of dose calculations.

Spot Scanning and Passive Scattering Proton Therapy: Relative Biological Effectiveness and Oxygen Enhancement Ratio in Cultured Cells.

PURPOSE To determine the relative biological effectiveness (RBE), oxygen enhancement ratio (OER), and contribution of the indirect effect of spot scanning proton beams, passive scattering proton beams, or both in cultured cells in comparison with clinically used photons. METHODS AND MATERIALS The RBE of passive scattering proton beams at the center of the spread-out Bragg peak (SOBP) was determined from dose-survival curves in 4 cell lines using 6-MV X rays as controls. Survival of 2 cell lines after spot scanning and passive scattering proton irradiation was then compared. Biological effects at the distal end region of the SOBP were also investigated. The OER of passive scattering proton beams and 6 MX X rays were investigated in 2 cell lines. The RBE and OER values were estimated at a 10% cell survival level. The maximum degree of protection of radiation effects by dimethyl sulfoxide was determined to estimate the contribution of the indirect effect against DNA damage. All experiments comparing protons and X rays were made under the same biological conditions. RESULTS The RBE values of passive scattering proton beams in the 4 cell lines examined were 1.01 to 1.22 (average, 1.14) and were almost identical to those of spot scanning beams. Biological effects increased at the distal end of the SOBP. In the 2 cell lines examined, the OER was 2.74 (95% confidence interval, 2.56-2.80) and 3.08 (2.84-3.11), respectively, for X rays, and 2.39 (2.38-2.43) and 2.72 (2.69-2.75), respectively, for protons (P<.05 for both cells between X rays and protons). The maximum degree of protection was significantly higher for X rays than for proton beams (P<.05). CONCLUSIONS The RBE values of spot scanning and passive scattering proton beams were almost identical. The OER was lower for protons than for X rays. The lower contribution of the indirect effect may partly account for the lower OER of protons.

High resolution Cerenkov light imaging of induced positron distribution in proton therapy

PURPOSE In proton therapy, imaging of the positron distribution produced by fragmentation during or soon after proton irradiation is a useful method to monitor the proton range. Although positron emission tomography (PET) is typically used for this imaging, its spatial resolution is limited. Cerenkov light imaging is a new molecular imaging technology that detects the visible photons that are produced from high-speed electrons using a high sensitivity optical camera. Because its inherent spatial resolution is much higher than PET, the authors can measure more precise information of the proton-induced positron distribution with Cerenkov light imaging technology. For this purpose, they conducted Cerenkov light imaging of induced positron distribution in proton therapy. METHODS First, the authors evaluated the spatial resolution of our Cerenkov light imaging system with a (22)Na point source for the actual imaging setup. Then the transparent acrylic phantoms (100 × 100 × 100 mm(3)) were irradiated with two different proton energies using a spot scanning proton therapy system. Cerenkov light imaging of each phantom was conducted using a high sensitivity electron multiplied charge coupled device (EM-CCD) camera. RESULTS The Cerenkov lights spatial resolution for the setup was 0.76 ± 0.6 mm FWHM. They obtained high resolution Cerenkov light images of the positron distributions in the phantoms for two different proton energies and made fused images of the reference images and the Cerenkov light images. The depths of the positron distribution in the phantoms from the Cerenkov light images were almost identical to the simulation results. The decay curves derived from the region-of-interests (ROIs) set on the Cerenkov light images revealed that Cerenkov light images can be used for estimating the half-life of the radionuclide components of positrons. CONCLUSIONS High resolution Cerenkov light imaging of proton-induced positron distribution was possible. The authors conclude that Cerenkov light imaging of proton-induced positron is promising for proton therapy.

Optimization and verification of image reconstruction for a Compton camera towards application as an on-line monitor for particle therapy

Particle therapy is an advanced cancer therapy that uses a feature known as the Bragg peak, in which particle beams suddenly lose their energy near the end of their range. The Bragg peak enables particle beams to damage tumors effectively. To achieve precise therapy, the demand for accurate and quantitative imaging of the beam irradiation region or dosage during therapy has increased. The most common method of particle range verification is imaging of annihilation gamma rays by positron emission tomography. Not only 511-keV gamma rays but also prompt gamma rays are generated during therapy; therefore, the Compton camera is expected to be used as an on-line monitor for particle therapy, as it can image these gamma rays in real time. Proton therapy, one of the most common particle therapies, uses a proton beam of approximately 200 MeV, which has a range of ~ 25 cm in water. As gamma rays are emitted along the path of the proton beam, quantitative evaluation of the reconstructed images of diffuse sources becomes crucial, but it is far from being fully developed for Compton camera imaging at present. In this study, we first quantitatively evaluated reconstructed Compton camera images of uniformly distributed diffuse sources, and then confirmed that our Compton camera obtained 3 %(1 σ) and 5 %(1 σ) uniformity for line and plane sources, respectively. Based on this quantitative study, we demonstrated on-line gamma imaging during proton irradiation. Through these studies, we show that the Compton camera is suitable for future use as an on-line monitor for particle therapy.