Kenta Takada

Erratum to ‘Microdosimetric calculation of relative biological effectiveness for design of therapeutic proton beams’

The authors attempt to establish the relative biological effectiveness (RBE) calculation for designing therapeutic proton beams on the basis of microdosimetry. The tissue-equivalent proportional counter (TEPC) was used to measure microdosimetric lineal energy spectra for proton beams at various depths in a water phantom. An RBE-weighted absorbed dose is defined as an absorbed dose multiplied by an RBE for cell death of human salivary gland (HSG) tumor cells in this study. The RBE values were calculated by a modified microdosimetric kinetic model using the biological parameters for HSG tumor cells. The calculated RBE distributions showed a gradual increase to about 1cm short of a beam range and a steep increase around the beam range for both the mono-energetic and spread-out Bragg peak (SOBP) proton beams. The calculated RBE values were partially compared with a biological experiment in which the HSG tumor cells were irradiated by the SOBP beam except around the distal end. The RBE-weighted absorbed dose distribution for the SOBP beam was derived from the measured spectra for the mono-energetic beam by a mixing calculation, and it was confirmed that it agreed well with that directly derived from the microdosimetric spectra measured in the SOBP beam. The absorbed dose distributions to planarize the RBE-weighted absorbed dose were calculated in consideration of the RBE dependence on the prescribed absorbed dose and cellular radio-sensitivity. The results show that the microdosimetric measurement for the mono-energetic proton beam is also useful for designing RBE-weighted absorbed dose distributions for range-modulated proton beams.

Proton beam irradiation stimulates migration and invasion of human U87 malignant glioma cells

Migration and invasion of malignant glioma play a major role in tumor progression and can be increased by low doses of gamma or X-ray irradiation, especially when the migrated tumor cells are located at a distance from the main tumor mass or postoperative cavity and are irradiated in fractions. We studied the influence of proton beam irradiation on migration and invasion of human U87 malignant glioma (U87MG) cells. Irradiation at 4 and 8 Gy increased cell migration by 9.8% (±4, P = 0.032) and 11.6% (±6.6, P = 0.031) and invasion by 45.1% (±16.5, P = 0.04) and 40.5% (±12.7, P = 0.041), respectively. After irradiation at 2 and 16 Gy, cell motility did not differ from that at 0 Gy. We determined that an increase in proton beam irradiation dose to over 16 Gy might provide tumor growth control, although additional specific treatment might be necessary to prevent the potentially increased motility of glioma cells during proton beam therapy.

Verification of nuclear data for the Tsukuba plan, a newly developed treatment planning system for boron neutron capture therapy.

Various verifications were performed to apply JENDL-4.0 as nuclear data for a newly developed treatment planning system with a homogeneous or precise human-like phantom. The nitrogen dose calculated by JENDL-4.0 differed slightly from that calculated by ENDF/B-VII.0. However, the total weighted dose-based dose volume histogram in the boron neutron capture therapy (BNCT) treatment for brain tumors calculated by JENDL-4.0 was in good agreement with the results of the ENDF/B-VII.0 calculation. Therefore, calculation with JENDL-4.0 can be applied to the BNCT dose calculation.

Neutron spectral fluence measurements using a Bonner sphere spectrometer in the development of the iBNCT accelerator-based neutron source

The neutron spectral fluence of an accelerator-based neutron source facility for boron neutron capture therapy (BNCT) based on a proton linac and a beryllium target was evaluated by the unfolding method using a Bonner sphere spectrometer (BSS). A 3He-proportional-counter-based BSS was used with weak beam during the development of the facility. The measured epithermal neutron spectra were consistent with calculations. The epithermal neutron intensity at the beam port was estimated and the results gave a numerical target for the enhancement of the proton beam intensity and will be used as reference data for measurements performed after the completion of the facility.

Influence of echo time in quantitative proton MR spectroscopy using LCModel

OBJECTIVE The objective of this study was to elucidate the influence on quantitative analysis using LCModel with the condition of echo time (TE) longer than the recommended values in the spectrum acquisition specifications. METHODS A 3T magnetic resonance system was used to perform proton magnetic resonance spectroscopy. The participants were 5 healthy volunteers and 11 patients with glioma. Data were collected at TE of 72, 144 and 288ms. LCModel was used to quantify several metabolites (N-acetylaspartate, creatine and phosphocreatine, and choline-containing compounds). The results were compared with quantitative values obtained by using the T2-corrected internal reference method. RESULTS In healthy volunteers, when TE was long, the quantitative values obtained using LCModel were up to 6.8-fold larger (p<0.05) than those obtained using the T2-corrected internal reference method. The ratios of the quantitative values obtained by the two methods differed between metabolites (p<0.05). In patients with glioma, the ratios of quantitative values obtained by the two methods tended to be larger at longer TE, similarly to the case of healthy volunteers, and large between-individual variation in the ratios was observed. CONCLUSIONS In clinical practice, TE is sometimes set longer than the value recommended for LCModel. If TE is long, LCModel overestimates the quantitative value since it cannot compensate for signal attenuation, and this effect is different for each metabolite and condition. Therefore, if TE is longer than recommended, it is necessary to account for the possibly reduced reliability of quantitative values calculated using LCModel.

Validation of the physical and RBE-weighted dose estimator based on PHITS coupled with a microdosimetric kinetic model for proton therapy

Abstract The microdosimetric kinetic model (MKM) is widely used for estimating relative biological effectiveness (RBE)-weighted doses for various radiotherapies because it can determine the surviving fraction of irradiated cells based on only the lineal energy distribution, and it is independent of the radiation type and ion species. However, the applicability of the method to proton therapy has not yet been investigated thoroughly. In this study, we validated the RBE-weighted dose calculated by the MKM in tandem with the Monte Carlo code PHITS for proton therapy by considering the complete simulation geometry of the clinical proton beam line. The physical dose, lineal energy distribution, and RBE-weighted dose for a 155 MeV mono-energetic and spread-out Bragg peak (SOBP) beam of 60 mm width were evaluated. In estimating the physical dose, the calculated depth dose distribution by irradiating the mono-energetic beam using PHITS was consistent with the data measured by a diode detector. A maximum difference of 3.1% in the depth distribution was observed for the SOBP beam. In the RBE-weighted dose validation, the calculated lineal energy distributions generally agreed well with the published measurement data. The calculated and measured RBE-weighted doses were in excellent agreement, except at the Bragg peak region of the mono-energetic beam, where the calculation overestimated the measured data by ~15%. This research has provided a computational microdosimetric approach based on a combination of PHITS and MKM for typical clinical proton beams. The developed RBE-estimator function has potential application in the treatment planning system for various radiotherapies.

Development of Monte Carlo based real-time treatment planning system with fast calculation algorithm for boron neutron capture therapy.

PURPOSE We simulated the effect of patient displacement on organ doses in boron neutron capture therapy (BNCT). In addition, we developed a faster calculation algorithm (NCT high-speed) to simulate irradiation more efficiently. METHODS We simulated dose evaluation for the standard irradiation position (reference position) using a head phantom. Cases were assumed where the patient body is shifted in lateral directions compared to the reference position, as well as in the direction away from the irradiation aperture. For three groups of neutron (thermal, epithermal, and fast), flux distribution using NCT high-speed with a voxelized homogeneous phantom was calculated. The three groups of neutron fluxes were calculated for the same conditions with Monte Carlo code. These calculated results were compared. RESULTS In the evaluations of body movements, there were no significant differences even with shifting up to 9mm in the lateral directions. However, the dose decreased by about 10% with shifts of 9mm in a direction away from the irradiation aperture. When comparing both calculations in the phantom surface up to 3cm, the maximum differences between the fluxes calculated by NCT high-speed with those calculated by Monte Carlo code for thermal neutrons and epithermal neutrons were 10% and 18%, respectively. The time required for NCT high-speed code was about 1/10th compared to Monte Carlo calculation. CONCLUSIONS In the evaluation, the longitudinal displacement has a considerable effect on the organ doses. We also achieved faster calculation of depth distribution of thermal neutron flux using NCT high-speed calculation code.

DEVELOPMENT OF A MULTIMODAL MONTE CARLO BASED TREATMENT PLANNING SYSTEM

To establish boron neutron capture therapy (BNCT), the University of Tsukuba is developing a treatment device and peripheral devices required in BNCT, such as a treatment planning system. We are developing a new multimodal Monte Carlo based treatment planning system (developing code: Tsukuba Plan). Tsukuba Plan allows for dose estimation in proton therapy, X-ray therapy and heavy ion therapy in addition to BNCT because the system employs PHITS as the Monte Carlo dose calculation engine. Regarding BNCT, several verifications of the system are being carried out for its practical usage. The verification results demonstrate that Tsukuba Plan allows for accurate estimation of thermal neutron flux and gamma-ray dose as fundamental radiations of dosimetry in BNCT. In addition to the practical use of Tsukuba Plan in BNCT, we are investigating its application to other radiation therapies.

Investigation of the neutron spectrum measurement method for dose evaluation in boron neutron capture therapy

In boron neutron capture therapy, it is important to evaluate the dose administered to a patients body outside the tumour area. The exposure dose is evaluated by calculation; however, the calculated value must be validated using a measured value. The dose evaluations based on the measured neutron spectrum are investigated. Multi-foil activation, combined with a LiCaAlF6 scintillation detector and an imaging plate, is proposed as a measurement method. The proposed method can measure the neutron spectrum at various points quickly.

Effect of biological factors on successful measurements with skeletal-muscle 1 H-MRS

Background Our purpose in this study was to clarify whether differences in subject group attributes could affect data acquisition in proton magnetic resonance spectroscopy (1H-MRS). Methods Subjects without diabetes mellitus (DM) were divided into two groups (group A, in their 20s; group B, 30–60 years old). Subjects with DM formed group C (30–60 years old). The numbers of subjects were 19, 27, and 22 for group A, B, and C respectively. For all subjects, 1H-MRS measurements were taken of the soleus muscle (SOL) and the anterior tibial muscle (AT). We defined the success of the measurements by the detection of intramyocellular lipids. Moreover, we also measured the full width at half maximum of the water peaks for all subjects. Results The success rate was significantly higher for the AT (100%) than for the SOL (81.6%) (P<0.01). For the SOL, the success rate was 100% in group A, 85.2% in group B, and 77.3% in group C. There was a significant difference (P<0.05) between groups A and B, as well as between groups A and C. In all subjects, there was a significant difference (P<0.01) in the full width at half maximum (Hz) of the water peak between the AT and SOL measurements. Conclusion We conclude that differences in the age and DM history of subjects could affect the probability of successful 1H-MRS data acquisition.