Calculations of depth-dose distributions, cross sections and momentum loss

The impact of DNA double-strand break repair pathways throughout the carbon ion spread-out Bragg peak beam

Following carbon ion beam irradiation in mammalian cells, such as used in carbon ion radiotherapy (CIRT), it has been suggested that the balance between whether nonhomologous end joining (NHEJ) or homologous recombination (HR) is utilized depends on the DNA double-strand break (DSB) complexity. Here, we quantified DSB distribution and identified the importance of each DSB repair pathway at increasing depths within the carbon ion spread-out Bragg peak (SOBP) beam range. Chinese hamster ovary (CHO) cell lines were irradiated in a single biological system capable of incorporating the full carbon ion SOBP beam range. Cytotoxicity and DSB distribution/repair kinetics were examined at increasing beam depths using cell survival as an endpoint and γ-H2AX as a surrogate marker for DSBs. We observed that proximal SOBP had the highest number of total foci/cell and lowest survival, while distal SOBP had the most dense tracks. Both NHEJ- and HR-deficient CHO cells portrayed an increase in radiosensitivity throughout the full carbon beam range, although NHEJ-deficient cells were the most radiosensitive cell line from beam entrance up to proximal SOBP and demonstrated a dose-dependent decrease in ability to repair DSBs. In contrast, HR-deficient cells had the greatest ratio of survival fraction at entrance depth to the lowest survival fraction within the SOBP and demonstrated a linear energy transfer (LET)-dependent decrease in ability to repair DSBs. Collectively, our results provide insight into treatment planning and potential targets to inhibit, as HR was a more beneficial pathway to inhibit than NHEJ to enhance the cell killing effect of CIRT in targeted tumor cells within the SOBP while maintaining limited unwanted damage to surrounding healthy cells.

Calculating dose-averaged linear energy transfer in an analytical treatment planning system for carbon-ion radiotherapy

BACKGROUND

Compelling evidence shows the association between the relative biological effectiveness (RBE) of carbon-ion radiotherapy (CIRT) and the dose averaged linear energy transfer (LETd). However, the ability to calculate the LETd in commercially available treatment planning systems (TPS) is lacking.

PURPOSE

This study aims to develop a method of calculating the LETd of CIRT plans that could be robustly carried out in RayStation (V10B, Raysearch, Sweden).

METHODS

The calculation used the fragment spectra in RayStation for the CIRT treatment planning. The dose-weighted averaging procedure was supported by the microdosimetric kinetic model (MKM). The MKM-based pencil beam dose engine (PBA, v4.2) for calculating RBE-weighted doses was reformulated to become a LET-weighted calculating engine. A separate module was then configured to inversely calculate the LETd from the absorbed dose of a plan and the associated fragment spectra. In this study, the ion and energy-specific LET table in the LETd module was further matched with the values decoded from the baseline data of the Syngo TPS (V13C, Siemens, Germany). The LETd distributions of several monoenergetic and modulated beams were calculated and validated against the values derived from the Syngo TPS and the published data.

RESULTS

The differences in LETds of the monoenergetic beams between the new method and the traditional method were within 3% in the entrance and Bragg-peak regions. However, a larger difference was observed in the distal region. The results of the modulated beams were in good agreement with the works from the published literature.

CONCLUSIONS

The method presented herein reformulates the MKM dose engine in the RayStation TPS to inversely calculate LETds. The robustness and accuracy were demonstrated.

Effective radiotherapeutic treatment intensification in patients with pancreatic cancer: higher doses alone, higher RBE or both?

Pancreatic cancer, especially in case of locally advanced stage has a poor prognosis. Radiotherapy in general can lead to tumor volume reduction, but further improvements, such as ion beam therapy have to be promoted in order to enable dose escalation, which in turn results in better local control rates and downsizing of the tumor itself. Ion beam therapy with its highly promising physical properties is also accompanied by distinct inter- and intrafractional challenges in case of robustness. First clinical results are promising, but further research in motion mitigation and biological treatment planning is necessary, in order to determine the best clinical rationales and conditions of ion beam therapy of pancreatic cancer. This review summarizes the current knowledge and studies on ion beam therapy of pancreatic cancer.

Carbon ion fragmentation effects on the nanometric level behind the Bragg peak depth

In this study, fragmentation yields of carbon therapy beams are estimated using the Geant4 simulation toolkit version 9.5. Simulations are carried out in a step-by-step mode using the Geant4-DNA processes for each of the major contributing fragments. The energy of the initial beam is taken 400 MeV amu(-1) as this is the highest energy, which is used for medical accelerators and this would show the integral role of secondary contributions in radiotherapy irradiations. The obtained results showed that 64% of the global dose deposition is initiated by carbon ions, while up to 36% is initiated by the produced fragments including all their isotopes. The energy deposition clustering yields of each of the simulated fragments are then estimated using the DBSCAN clustering algorithm and they are compared to the yields of the incident primary beam.

Dependence of simulated positron emitter yields in ion beam cancer therapy on modeling nuclear fragmentation

In ion beam cancer therapy, range verification in patients using positron emission tomography (PET) requires the comparison of measured with simulated positron emitter yields. We found that (1) changes in modeling nuclear interactions strongly affected the positron emitter yields and that (2) Monte Carlo simulations with SHIELD-HIT10Areasonably matched the most abundant PET isotopes (11)C and (15)O. We observed an ion-energy (i.e., depth) dependence of the agreement between SHIELD-HIT10Aand measurement. Improved modeling requires more accurate measurements of cross-section values.

Measurement of electron density in dual-energy x-ray CT with monochromatic x rays and evaluation of its accuracy

Information on electron density is important for radiotherapy treatment planning in order to optimize the dose distribution in the target volume of a patient. At present, the electron density is derived from a computed tomography (CT) number measured in x-ray CT scanning; however, there are uncertainties due to the beam hardening effect and the method by which the electron density is converted from the CT number. In order to measure the electron density with an accuracy of +/-1%, the authors have developed dual-energy x ray CT using monochromatic x rays. They experimentally proved that the measured linear attenuation coefficients were only a few percent lower than the theoretical ones, which led to an accuracy within 2% for the electron density. There were three factors causing inaccuracy in the linear attenuation coefficient and the electron density: the influence of scattered radiation, the nonlinearity in the detector response function, and a theoretical process to derive the electron density from the linear attenuation coefficients. The linear attenuation coefficients of water were experimentally proved to differ by 1%-2% from the theoretical one even when the scattering effect was negligible. The nonlinearity of the response function played an important role in correcting the difference in the linear attenuation coefficient. Furthermore, the theoretical process used for deriving the electron density from the linear attenuation coefficients introduces about 0.6% deviation from the theoretical value into the resultant electron density. This deviation occurs systematically so that it can be corrected. The authors measured the electron densities for seven samples equivalent to soft tissue in dual-energy x-ray CT, and finally obtained them with an accuracy of around +/-1%.

Effectiveness of monoenergetic and spread-out bragg peak carbon-ions for inactivation of various normal and tumour human cell lines

This work aimed at measuring cell-killing effectiveness of monoenergetic and Spread-Out Bragg Peak (SOBP) carbon-ion beams in normal and tumour cells with different radiation sensitivity. Clonogenic survival was assayed in normal and tumour human cell lines exhibiting different radiosensitivity to X- or gamma-rays following exposure to monoenergetic carbon-ion beams (incident LET 13-303 keV/microm) and at various positions along the ionization curve of a therapeutic carbon-ion beam, corresponding to three dose-averaged LET (LET(d)) values (40, 50 and 75 keV/microm). Chinese hamster V79 cells were also used. Carbon-ion effectiveness for cell inactivation generally increased with LET for monoenergetic beams, with the largest gain in cell-killing obtained in the cells most radioresistant to X- or gamma-rays. Such an increased effectiveness in cells less responsive to low LET radiation was found also for SOBP irradiation, but the latter was less effective compared with monoenergetic ion beams of the same LET. Our data show the superior effectiveness for cell-killing exhibited by carbon-ion beams compared to lower LET radiation, particularly in tumour cells radioresistant to X- or gamma-rays, hence the advantage of using such beams in radiotherapy. The observed lower effectiveness of SOBP irradiation compared to monoenergetic carbon beam irradiation argues against the radiobiological equivalence between dose-averaged LET in a point in the SOBP and the corresponding monoenergetic beams.

Specification of Carbon Ion Dose at the National Institute of Radiological Sciences (NIRS)

The clinical dose distributions of therapeutic carbon beams, currently used at NIRS HIMAC, are based on in-vitro Human Salivary Gland (HSG) cell survival response and clinical experience from fast neutron radiotherapy. Moderate radiosensitivity of HSG cells is expected to be a typical response of tumours to carbon beams. At first, the biological dose distribution is designed so as to cause a flat biological effect on HSG cells in the spread-out Bragg peak (SOBP) region. Then, the entire biological dose distribution is evenly raised in order to attain a RBE (relative biological effectiveness) = 3.0 at a depth where dose-averaged LET (linear energy transfer) is 80 keV/mum. At that point, biological experiments have shown that carbon ions can be expected to have a biological effect identical to fast neutrons, which showed a clinical RBE of 3.0 for fast neutron radiotherapy at NIRS. The resulting clinical dose distribution in this approximation is not dependent on dose level, tumour type or fractionation scheme and thus reduces the unknown parameters in the analysis of the clinical results. The width SOBP and the clinical / physical dose at the center of SOBP specify the dose distribution. The clinical results analysed in terms of TCP were found to show good agreement with the expected RBE value at higher TCP levels. The TCP analysis method was applied for the prospective dose estimation of hypofractionation.

Irradiation System for HIMAC

Clinical trials of carbon radiotherapy started at HIMAC in 1994 using three treatment rooms and four beam ports, two horizontal and two vertical. The broad beam method was adopted to make a three-dimensionally uniform field at an isocenter. A spot beam extracted from an accelerator was laterally spread out by using a pair of wobbler magnets and a scatterer. A bar ridge filter modulated the beam energy to obtain the spread out Bragg peak (SOBP). The SOBP was designed to be flat in terms of the biological dose based on the consideration that the field consisted of various beams with different LET. Finally, the field of 20 cm in diameter with +/- 2.5% uniformity was formed at the isocenter. The width of the maximum SOBP was 15 cm. When treating the lung or liver, organs that move due to breathing, the beam was irradiated only during the expiration period in a respiration-gated irradiation method. This reduced the treatment margin of the moving target. In order to prevent normal tissues adjacent to the target volume from irradiation by an unwanted dose, a layer-stacking method was developed. In this method, thin SOBP layers which have different ranges were piled up step by step from the distal end to the entrance of the target volume. At the same time, a multi-leaf collimator was used to change the aperture shape to match the shape of each layer to the cross-sectional shape of the target. This method has been applied to rather large volume cancers including bone and soft-tissue cancers. Only a few serious problems in the irradiation systems have been encountered since the beginning of the clinical trials. Overall the systems have been working stably and reliably.

The potential application of β-delayed particle decay beam9C in cancer therapy

A radioactive ion beam like 9C serves as a double radiation source and may be useful in cancer treatment, where the essential irradiation comes from the external beam itself and the extra one is due to the low-energy particles emitted internally during the decay of 9C. Based on the microdosimetric specific energy spectrum in cell nuclei, a model to evaluate the biological effect induced by the internally emitted particles from a beta-delayed particle decay beam has been developed. In this paper, using this model the additional contributions to the cell-killing effect due to the emitted particles from stopping 9C ions were incorporated in the design of spread-out Bragg peaks (SOBP) for radioactive 9C beams. For this purpose, a simulated annealing algorithm was employed to optimize the superposing weighting fractions of all monoenergetic beams so that a uniform cell survival level could be realized across the SOBP within an acceptable deviation of 5%. SOBPs with different widths and at different cell survival levels were designed for both therapeutic 9C and 12C beams for comparison. The potential use of the 9C beam in radiotherapy compared to the 12C beam, which is commonly adopted in the practices of current heavy-ion therapy, is shown systematically in terms of the distributions of biological effective dose and cell survival along the beam penetration.

A model to evaluate the biological effect induced by the emitted particles from a -delayed particle decay beam

Due to their favourable properties such as high dose localization and high RBE heavy-ion beams have attracted increasing interest in cancer treatment. Efforts to exploit these advantages to the maximum extent in cancer therapy have never been given up. A new idea of applying a radioactive ion beam with beta-delayed particle decay such as 9C or 8B to cancer therapy is put forward in this paper. A model to evaluate the biological effect in terms of cell survival induced by the emitted particles from the decays of the stopped ions has been established. Because of the difference of the internally emitted particle irradiation from the external ion beam, the microdosimetric quantity such as specific energy is applied to evaluate the cell surviving effect induced by the emitted particles from the decays of the radioactive ions. Within the framework of this model, the cell-killing effects resulting from the emitted particles were calculated under different conditions. Finally, the potential application of the radioactive ion beam 9C in cancer therapy is demonstrated.

Treatment planning for the layer-stacking irradiation system for three-dimensional conformal heavy-ion radiotherapy

We have upgraded a heavy-ion radiotherapy treatment-planning system to adapt for the layer-stacking irradiation method, which is to conform a variable spread-out Bragg peak to a target volume by means of dynamic control of the conventional beam-modifying devices. The biophysical model, the beam-setup logic, and the dose-calculation algorithm implemented for the layer-stacking method are described and the expected clinical usability is discussed. The layer-stacking method was integrated in perfect accordance with the ongoing conventional treatments so that the established protocols, which are the clinically optimized dose fractionation schemes, will still be valid. On the other hand, a simulation study indicated a substantial improvement of dose distribution with the layer-stacking method though the significance may depend on the size, shape, and location of the tumor. The completed treatment system will provide an option for improved conformal radiotherapy without interfering with the conventional method and we expect a gradual expansion of the clinical cases applicable to the layer-stacking method.

Biophysical characteristics of HIMAC clinical irradiation system for heavy-ion radiation therapy

PURPOSE

The irradiation system and biophysical characteristics of carbon beams are examined regarding radiation therapy.

METHODS AND MATERIALS

An irradiation system was developed for heavy-ion radiotherapy. Wobbler magnets and a scatterer were used for flattening the radiation field. A patient-positioning system using X ray and image intensifiers was also installed in the irradiation system. The depth-dose distributions of the carbon beams were modified to make a spread-out Bragg peak, which was designed based on the biophysical characteristics of monoenergetic beams. A dosimetry system for heavy-ion radiotherapy was established to deliver heavy-ion doses safely to the patients according to the treatment planning. A carbon beam of 80 keV/microm in the spread-out Bragg peak was found to be equivalent in biological responses to the neutron beam that is produced at cyclotron facility in National Institute Radiological Sciences (NIRS) by bombarding 30-MeV deuteron beam on beryllium target. The fractionation schedule of the NIRS neutron therapy was adapted for the first clinical trials using carbon beams.

RESULTS

Carbon beams, 290, 350, and 400 MeV/u, were used for a clinical trial from June of 1994. Over 300 patients have already been treated by this irradiation system by the end of 1997.

Relationship between CT number and electron density, scatter angle and nuclear reaction for hadron-therapy treatment planning

The precise conversion of CT numbers to their electron densities is essential in treatment planning for hadron therapy. Although some conversion methods have already been proposed, it is hard to check the conversion accuracy during practical therapy. We have estimated the CT numbers of real tissues by a calculational method established by Mustafa and Jackson. The relationship between the CT numbers and the electron densities was investigated for various body tissues as well as some tissue-equivalent materials used for a conversion to check the accuracy of the current conversion methods. The result indicates a slight disagreement at the high-CT-number region. A precise estimation of the multiple scattering, nuclear reaction and range straggling of incident particles has been considered as being important to realize higher-level conformal therapy in the future. The relationship between these parameters and the CT numbers was also investigated for tissues and water. The result shows that it is sufficiently practical to replace these parameters for real tissues with those for water by adjusting the density.