Assessment of Cystamine's Radioprotective/Antioxidant Ability under High-Dose-Rate Irradiation: A Monte Carlo Multi-Track Chemistry Simulation Study.

(1) Background: cystamine and its reduced form cysteamine have radioprotective/antioxidant effects in vivo. In this study, we use an in vitro model system to examine the behavior of cystamine towards the reactive primary species produced during the radiolysis of the Fricke dosimeter under high dose-rate irradiation conditions. (2) Methods: our approach was to use the familiar radiolytic oxidation of ferrous to ferric ions as an indicator of the radioprotective/antioxidant capacity of cystamine. A Monte Carlo computer code was used to simulate the multi-track radiation-induced chemistry of aerated and deaerated Fricke-cystamine solutions as a function of dose rate while covering a large range of cystamine concentrations. (3) Results: our simulations revealed that cystamine provides better protection at pulsed dose rates compared to conventional, low-dose-rate irradiations. Furthermore, our simulations confirmed the radical-capturing ability of cystamine, clearly indicating the strong antioxidant profile of this compound. (4) Conclusion: assuming that these findings can be transposable to cells and tissues at physiological pH, it is suggested that combining cystamine with FLASH-RT could be a promising approach to further enhance the therapeutic ratio of cancer cure.

A Computer Modeling Study of Water Radiolysis at High Dose Rates. Relevance to FLASH Radiotherapy.

"FLASH radiotherapy" is a new method of radiation treatment by which large doses of radiation are delivered at high dose rates to tumors almost instantaneously (a few milliseconds), paradoxically sparing healthy tissue while preserving anti-tumor activity. To date, no definitive mechanism has been proposed to explain the different responses of the tumor and normal tissue to radiation. As a first step, and given that living cells and tissues consist mainly of water, we studied the effects of high dose rates on the transient yields (G values) of the radical and molecular species formed in the radiolysis of deaerated/aerated water by irradiating protons, using Monte Carlo simulations. Our simulation model consisted of two steps: 1. The random irradiation of a right circular cylindrical volume of water, embedded in nonirradiated bulk water, with single and instantaneous pulses of N 300-MeV incident protons ("linear energy transfer" or LET ∼ 0.3 keV/µm) traveling along the axis of the cylinder; and 2. The development of these N proton tracks, which were initially contained in the irradiated cylinder, throughout the solution over time. The effect of dose rate was studied by varying N, which was calibrated in terms of dose rate. For this, experimental data on the yield G(Fe3+) of the super-Fricke dosimeter as a function of dose rate up to ∼1010 Gy/s were used. Confirming previous experimental and theoretical studies, significant changes in product yields were found to occur with increasing dose rate, with lower radical and higher molecular yields, which result from an increase in the radical density in the bulk of the solution. Using the kinetics of the decay of hydrated electrons, a critical time (τc), which corresponds to the "onset" of dose-rate effects, was determined for each value of N. For the cylindrical irradiation model, τc was inversely proportional to the dose rate. Moreover, the comparison with experiments with pulsed electrons underlined the importance of the geometry of the irradiation volume for the estimation of τc. Finally, in the case of aerated water radiolysis, we calculated the yield of oxygen consumption and estimated the corresponding concentration of consumed (depleted) oxygen as a function of time and dose rate. It was shown that this concentration increases substantially with increasing dose rate in the time window ∼1 ns-10 µs, with a very pronounced maximum around 0.2 µs. For high-dose-rate irradiations (>109 Gy/s), a large part of the available oxygen (∼0.25 mM for an air-saturated solution) was found to be consumed. This result, which was obtained on a purely water radiation chemistry basis, strongly supports the hypothesis that the normal tissue-sparing effect of FLASH stems from temporary hypoxia due to oxygen depletion induced by high-dose-rate irradiation.

Ultra-High Dose-Rate, Pulsed (FLASH) Radiotherapy with Carbon Ions: Generation of Early, Transient, Highly Oxygenated Conditions in the Tumor Environment.

It is well known that molecular oxygen is a product of the radiolysis of water with high-linear energy transfer (LET) radiation, which is distinct from low-LET radiation wherein O2 radiolytic yield is negligible. Since O2 is a powerful radiosensitizer, this fact is of practical relevance in cancer therapy with energetic heavy ions, such as carbon ions. It has recently been discovered that large doses of ionizing radiation delivered to tumors at very high dose rates (i.e., in a few milliseconds) have remarkable benefits in sparing healthy tissue while preserving anti-tumor activity compared to radiotherapy delivered at conventional, lower dose rates. This new method is called "FLASH radiotherapy" and has been tested using low-LET radiation (i.e., electrons and photons) in various pre-clinical studies and recently in a human patient. Although the exact mechanism(s) underlying FLASH are still unclear, it has been suggested that radiation delivered at high dose rates spares normal tissue via oxygen depletion. In addition, heavy-ion radiation achieves tumor control with reduced normal tissue toxicity due to its favorable physical depth-dose profile and increased radiobiological effectiveness in the Bragg peak region. To date, however, biological research with energetic heavy ions delivered at ultra-high dose rates has not been performed and it is not known whether heavy ions are suitable for FLASH radiotherapy. Here we present the additive or even synergistic advantages of integrating the FLASH dose rates into carbon-ion therapy. These benefits result from the ability of heavy ions at high LET to generate an oxygenated microenvironment around their track due to the occurrence of multiple (mainly double) ionization of water. This oxygen is abundant immediately in the tumor region where the LET of the carbon ions is very high, near the end of the carbon-ion path (i.e., in the Bragg peak region). In contrast, in the "plateau" region of the depth-dose distribution of ions (i.e., in the normal tissue region), in which the LET is significantly lower, this generation of molecular oxygen is insignificant. Under FLASH irradiation, it is shown that this early generation of O2 extends evenly over the entire irradiated tumor volume, with concentrations estimated to be several orders of magnitude higher than the oxygen levels present in hypoxic tumor cells. Theoretically, these results indicate that FLASH radiotherapy using carbon ions would have a markedly improved therapeutic ratio with greater toxicity in the tumor due to the generation of oxygen at the spread-out Bragg peak.

Modeling the radiolysis of supercritical water by fast neutrons: density dependence of the yields of primary species at 400°c.

A reliable understanding of radiolysis processes in supercritical water (SCW)-cooled reactors is crucial to developing chemistry control strategies that minimize the corrosion and degradation of materials. However, directly measuring the chemistry in reactor cores is difficult due to the extreme conditions of high temperature and pressure and mixed neutron and gamma-radiation fields, which are incompatible with normal chemical instrumentation. Thus, chemical models and computer simulations are an important route of investigation for predicting the detailed radiation chemistry of the coolant in a SCW reactor and the consequences for materials. Surprisingly, information on the fast neutron radiolysis of water at high temperatures is limited, and even more so for fast neutron irradiation of SCW. In this work, Monte Carlo simulations were used to predict the G values for the primary species e(-)aq, H(•), H2, (•)OH and H2O2 formed from the radiolysis of pure, deaerated SCW (H2O) by 2 MeV monoenergetic neutrons at 400°C as a function of water density in the range of ∼0.15-0.6 g/cm(3). The 2 MeV neutron was taken as representative of a fast neutron flux in a reactor. For light water, the moderation of these neutrons after knock-on collisions with water molecules generated mostly recoil protons of 1.264, 0.465, 0.171 and 0.063 MeV. Neglecting oxygen ion recoils and assuming that the most significant contribution to the radiolysis came from these first four recoil protons, the fast neutron yields were estimated as the sum of the G values for these protons after appropriate weightings were applied according to their energy. Calculated yields were compared with available experimental data and with data obtained for low-LET radiation. Most interestingly, the reaction of H(•) atoms with water was found to play a critical role in the formation yields of H2 and (•)OH at 400°C. Recent work has underscored the potential importance of this reaction above 200°C, but its rate constant is still controversial.

Calculation of the yields for the primary species formed from the radiolysis of liquid water by fast neutrons at temperatures between 25-350°C.

Monte Carlo simulations were used to calculate the yields for the primary species (e(-)aq, H(•), H2, (•)OH and H2O2) formed from the radiolysis of neutral liquid water by mono-energetic 2 MeV neutrons at temperatures between 25-350°C. The 2 MeV neutron was taken as representative of a fast neutron flux in a reactor. For light water, the moderation of these neutrons generated elastically scattered recoil protons of ∼1.264, 0.465, 0.171 and 0.063 MeV, which at 25°C, had linear energy transfers (LETs) of ∼22, 43, 69 and 76 keV/µm, respectively. Neglecting the radiation effects due to oxygen ion recoils and assuming that the most significant contribution to the radiolysis came from these first four recoil protons, the fast neutron yields could be estimated as the sum of the yields for these protons after allowance was made for the appropriate weightings according to their energy. Yields were calculated at 10(-7), 10(-6) and 10(-5) s after the ionization event at all temperatures, in accordance with the time range associated with the scavenging capacities generally used for fast neutron radiolysis experiments. The results of the simulations agreed reasonably well with the experimental data, taking into account the relatively large uncertainties in the experimental measurements, the relatively small number of reported radiolysis yields, and the simplifications included in the model. Compared with data obtained for low-LET radiation ((60)Co γ rays or fast electrons), our computed yields for fast neutron radiation showed essentially similar temperature dependences over the range of temperatures studied, but with lower values for yields of free radicals and higher values for molecular yields. This general trend is a reflection of the high-LET character of fast neutrons. Although the results of the simulations were consistent with the experiment, more experimental data are required to better describe the dependence of radiolytic yields on temperature and to test more thoroughly our modeling calculations.

Effect of temperature on the low-linear energy transfer radiolysis of the ceric-cerous sulfate dosimeter: a Monte Carlo simulation study.

The stochastic modeling of the (60)Co γ/fast-electron radiolysis of the ceric-cerous chemical dosimeter has been performed as a function of temperature from 25-350°C. The system used is a dilute solution of ceric sulfate and cerous sulfate in aqueous 0.4 M sulfuric acid. In this system, H(•) (or HO2(•) in the presence of dissolved oxygen) and H2O2 produced by the radiolytic decomposition of water both reduce Ce(4+) ions to Ce(3+) ions, while (•)OH radicals oxidize the Ce(3+) present in the solution back to Ce(4+). The net Ce(3+) yield is given by G(Ce(3+)) = g(H(•)) + 2 g(H2O2) - g((•)OH), where the primary (or "escape") yields of H(•), H2O2 and (•)OH are represented by lower case g's. At room temperature, G(Ce(3+)) has been established to be 2.44 ± 0.8 molecules/100 eV. In this work, we investigated the effect of temperature on the yield of Ce(3+) and on the underlying chemical reaction kinetics using Monte Carlo track chemistry simulations. The simulations showed that G(Ce(3+)) is time dependent, a result of the differences in the lifetimes of the reactions that make up the radiolysis mechanism. Calculated G(Ce(3+)) values were found to decrease almost linearly with increasing temperature up to about 250°C, and are in excellent agreement with available experimental data. In particular, our calculations confirmed previous estimated values by Katsumura et al. (Radiat Phys Chem 1988; 32:259-63) showing that G(Ce(3+)) at ∼250°C is about one third of its value at room temperature. Above ∼250°C, our model predicted that G(Ce(3+)) would drop markedly with temperature until, instead of Ce(4+) reduction, Ce(3+) oxidation is observed. This drop is shown to occur as a result of the reaction of hydrogen atoms with water in the homogeneous chemical stage.

Density dependence of the yield of hydrated electrons in the low-LET radiolysis of supercritical water at 400 °C: influence of the geminate recombination of subexcitation-energy electrons prior to thermalization.

Monte Carlo simulations were used to calculate the yield of hydrated electrons (eaq(-)) in the low-linear energy transfer radiolysis of supercritical water at 400 °C as a function of water density over the range of ~0.15 to 0.6 g cm(-3). Very good agreement was found between our calculations and picosecond pulse radiolysis experimental data at ~60 ps and 1 ns at high density (>0.35 g cm(-3)). At densities lower than ~0.35 g cm(-3), our eaq(-) yields were lower than the experimental data, especially at ~60 ps. However, if we incorporated into the simulations a prompt geminate electron-cation (H2O˙(+)) recombination (prior thermalization of the electron) that decreased as the density decreased, our computed eaq(-) yields at ~60 ps and 1 ns compared fairly well with the experimental data for the entire density range studied.

Time-dependent yield of the hydrated electron in subcritical and supercritical water studied by ultrafast pulse radiolysis and Monte-Carlo simulation.

Fast kinetics and time-dependent yields of the hydrated electron (e(-)(aq)) in pure water under conditions of high temperature and pressure up to the supercritical region were investigated by picosecond and nanosecond pulse radiolysis experiments. More significant decays at short times followed by plateau components at longer times were observed with increasing temperature, suggesting faster spur reaction processes. In supercritical water, it was also found that the e(-)(aq) yields strongly depend on the pressure (density). Comparison of these measurements with Monte-Carlo computer simulations allowed us to identify spur reactions of e(-)(aq) that occur predominantly at high temperatures and also to provide new key information on certain spur model parameters. In particular, the experimental time-dependent e(-)(aq) yields were best reproduced if the electron thermalization distance decreases with increasing temperature. This "shrinkage" of spur sizes at high temperatures was attributed to an increase in the scattering cross sections of subexcitation electrons, likely originating from a decrease in the degree of structural order of water molecules as the temperature is increased.

Utilization of the ferrous sulfate (Fricke) dosimeter for evaluating the radioprotective potential of cystamine: experiment and Monte Carlo simulation.

Cystamine, an organic disulfide (RSSR), is among the best of the known radiation-protective compounds and has been used to protect normal tissues in clinical radiation therapy. Recently, it has also proved to be beneficial in the treatment of disorders of the central nervous system in animal models. However, the underlying mechanism of its action at the chemical level is not yet well understood. The present study aims at using the ferrous sulfate (Fricke) dosimeter to quantitatively evaluate, both experimentally and theoretically, the radioprotective potential of this compound. The well-known radiolysis of the Fricke dosimeter by (60)Co γ rays or fast electrons, based on the oxidation of ferrous ions to ferric ions by the oxidizing species (•)OH, HO(2)(•), and H(2)O(2) produced in the radiolytic decomposition of water, forms the basis for our method. The presence of cystamine in Fricke dosimeter solutions during irradiation prevents the radiolytic oxidation of Fe(2+) and leads to decreased ferric yields (or G values). The observed decrease in G(Fe(3+)) increases upon increasing the concentration of the disulfide compound over the range 0-0.1 M under both aerated and deaerated conditions. To help assess the basic radiation-protective mechanism of this compound, a full Monte Carlo computer code is developed to simulate in complete detail the radiation-induced chemistry of the studied Fricke/cystamine solutions. Benefiting from the fact that cystamine is reasonably well characterized in terms of radiation chemistry, this computer model proposes reaction mechanisms and incorporates specific reactions describing the radiolysis of cystamine in aerated and deaerated Fricke solutions that lead to the observable quantitative chemical yields. Results clearly indicate that the protective effect of cystamine originates from its radical-capturing ability, which allows this compound to act by competing with the ferrous ions for the various free radicals--especially (•)OH radicals and H(•) atoms--formed during irradiation of the surrounding water. Most interestingly, our simulation modeling also shows that the predominant pathway in the oxidation of cystamine by (•)OH radicals involves an electron-transfer mechanism, yielding RSSR(•+) and OH(-). A very good agreement is found between calculated G(Fe(3+)) values and experiment. This study concludes that Monte Carlo simulations represent a very efficient method for understanding indirect radiation damage at the molecular level.

Temperature dependence of the Fricke dosimeter and spur expansion time in the low-LET high-temperature radiolysis of water up to 350 °C: a Monte-Carlo simulation study.

Monte-Carlo simulations of the radiolysis of the ferrous sulfate (Fricke) dosimeter with low-linear energy transfer (LET) radiation (such as (60)Co γ-rays or fast electrons) have been performed as a function of temperature from 25 to 350 °C. The predicted yields of Fe(2+) oxidation are found to increase with increasing temperature up to ∼100-150 °C, and then tend to remain essentially constant at higher temperatures, in very good agreement with experiment. By using a simple method based on the direct application of the stoichiometric relationship that exists between the ferric ion yields so obtained G(Fe(3+)) and the sum {3 [g(e(-)(aq) + H˙) + g(HO(2)˙)] + g(˙OH) + 2 g(H(2)O(2))}, where g(e(-)(aq) + H˙), g(HO(2)˙), g(˙OH), and g(H(2)O(2)) are the primary radical and molecular yields of the radiolysis of deaerated 0.4 M H(2)SO(4) aqueous solutions, the lifetime (τ(s)) of the spur and its temperature dependence have been determined. In the spirit of the spur model, τ(s) is an important indicator for overlapping spurs, giving the time required for the changeover from nonhomogeneous spur kinetics to homogeneous kinetics in the bulk solution. The calculations show that τ(s) decreases by about an order of magnitude over the 25-350 °C temperature range, going from ∼4.2 × 10(-7) s at 25 °C to ∼5.7 × 10(-8) s at 350 °C. This decrease in τ(s) with increasing temperature mainly originates from the quicker diffusion of the individual species involved. Moreover, the observed dependence of G(Fe(3+)) on temperature largely reflects the influence of temperature upon the primary free-radical product yields of the radiolysis, especially the yield of H˙ atoms. Above ∼200-250 °C, the more and more pronounced intervention of the reaction of H˙ atoms with water also contributes to the variation of G(Fe(3+)), which may decrease or increase slightly, depending on the choice made for the rate constant of this reaction. All calculations reported herein use the radiolysis database of Elliot (Atomic Energy of Canada Limited) and Bartels (University of Notre Dame) that contains all the best currently available information on the rate constants, reaction mechanisms, and g-values in the range 20 to 350 °C.

High-LET ion radiolysis of water: oxygen production in tracks.

It is known that molecular oxygen is a product of the radiolysis of water with high-linear energy transfer (LET) radiation, a result that is of particular significance in radiobiology and of practical relevance in radiotherapy. In fact, it has been suggested that the radiolytic formation of an oxygenated microenvironment around the tracks of high-LET heavy ions is an important factor in their enhanced biological efficiency in the sense that this may be due to an "oxygen effect" by O(2) produced by these ions in situ. Using Monte Carlo track simulations of pure, deaerated water radiolysis by 4.8 MeV (4)He(2+) (LET approximately 94 keV/microm) and 24 MeV (12)C(6+) (LET approximately 490 keV/microm) ions, including the mechanism of multiple ionization of water, we have calculated the yields and concentrations of O(2) in the tracks of these irradiating ions as a function of time between approximately 10(-12) and 10(-5) s at 25 and 37 degrees C. The track oxygen concentrations obtained compare very well with O(2) concentrations estimated from the "effective" amounts of oxygen that are needed to produce the observed reduction in oxygen enhancement ratio (OER) with LET (assuming this decrease is attributable to the sole radiolytic formation of O(2) in the tracks). For example, for 24 MeV (12)C(6+) ions, the initial track concentration of O(2) is estimated to be more than three orders of magnitude higher than the oxygen levels present in normally oxygenated and hypoxic tumor regions as well as in normal human cells. Such results, which largely plead in favor of the "oxygen in the heavy-ion track" hypothesis, could explain at least in part the greater efficiency of high-LET radiation for cell inactivation (at equal radiation dose).

Water radiolysis with heavy ions of energies up to 28 GeV. 3. Measurement of G(MV*+) in deaerated methyl viologen solutions containing various concentrations of sodium formate and Monte Carlo simulation.

Formation yields of methyl viologen cation radicals G(MV*+) (100 eV)(-1) have been measured in deaerated aqueous solutions of 0.25 mM methyl viologen (MV(2+)) containing various concentrations of formate anion (0.01-2 M) after irradiation with six different ion beams (4He(2+), 12C(6+), 20Ne(10+), 28Si(14+), 40Ar(18+) and 56Fe(26+) with incident energies varying from 0.6 to 28 GeV) provided by the Heavy Ion Medical Accelerator in Chiba (HIMAC) at the National Institute of Radiological Science (NIRS) in Japan. The sample solutions are irradiated at the incident energy of the ions using 1-cm irradiation cells. Corresponding LET values cover the range from 2.2 to 185 eV/nm. G(MV*+) increases with increasing formate concentration. In 4He(2+) radiolysis, it increases from 5.7 to 7.1 as the concentration of formate is increased from 0.01 to 2 M, while in 56Fe(26+) radiolysis, the MV*+ yield value changes from 2.2 to 4.1. The other values lie between the yields for 4He(2+) and 56Fe(26+). In addition, G(MV*+) decreases with increasing LET. In the case of 12C(6+) radiolysis, G(MV*+) increases with increasing energy of the carbon ions from 135 to 400 MeV/nucleon, i.e., with decreasing LET from 21 to 11 eV/nm. In parallel to the above measurements, Monte Carlo simulations of the radiolysis of the MV(2+)/formate solutions have been performed. Ionic strength effects on reactions between charged species are taken into account. To reproduce the experimental results, previously unreported reactions such as e(aq)(-) + MV*+, MV*(+) + *OH and *COO- + *OH have been introduced in the reaction scheme. After optimization, the rate constants of these latter two reactions are determined to be (3 +/- 0.5) x 10(10) and (5 +/- 0.5) x 10(10) M(-1) s(-1), respectively. By contrast, the reaction between e(aq)- and MV*+ is too slow to affect G(MV*+). On the basis of these calculations, characteristics of intratrack reactions induced by heavy-ion beams are discussed in reference not only to the scavenger method used for measurement of water decomposition product yields but also to the differences in the relative spatial distribution of the reactants as well as the places where their intratrack reaction occurs within the geometry of the ion track structure.

Effect of multiple ionization on the yield of H2O2 produced in the radiolysis of aqueous 0.4 M H2SO4 solutions by high-LET 12C6+ and 20Ne9+ ions.

Monte Carlo track structure simulations were performed to investigate the effect of multiple ionization of water on the primary (or "escape") (at approximately 10(-6) s) yield of hydrogen peroxide (G(H2O2)) produced in the radiolysis of deaerated 0.4 M H2SO4 solutions by 12C6+ and 20Ne9+ ions at high linear energy transfer (LET) up to approximately 900 keV/microm. It was found that, upon incorporating the mechanisms of double, triple and quadruple ionizations of water in the calculations, a quantitative agreement between theory and experiment can be obtained. The curve for G(H2O2) as a function of LET reaches a well-defined maximum of approximately 1.4 molecules/100 eV at approximately 180-200 keV/microm, in very good accord with the available experimental data. Our results also show that, for the highest LET values considered in this study, the H2O2 escape yields obtained in 0.4 M sulfuric acid solutions are about 45% greater in magnitude than those found in neutral water. Contrary to a recent assumption suggesting that the limiting value of G(H2O2) at infinite LET should be approximately 1 molecule/100 eV, somewhat similar for neutral and acidic water, our simulations show a clear decrease in the primary H2O2 yields with increasing LET at high LET, indicating that the question of the limiting value of G(H2O2) at very high LET for both neutral and acidic liquid water is still open.

Low-energy electron penetration range in liquid water.

Monte Carlo simulations of electron tracks in liquid water are performed to calculate the energy dependence of the electron penetration range at initial electron energies between 0.2 eV and 150 keV, including the subexcitation electron region (<7.3 eV). Our calculated electron penetration distances are compared with available experimental data and earlier calculations as well as with the results of simulations using newly reported amorphous ice electron scattering cross sections in the range approximately 1-100 eV.