Track Structure-Based Simulations on DNA Damage Induced by Diverse Isotopes.

Diverse isotopes such as 2H, 3He, 10Be, 11C and 14C occur in nuclear reactions in ion beam radiotherapy, in cosmic ray shielding, or are intentionally accelerated in dating techniques. However, only a few studies have specifically addressed the biological effects of diverse isotopes and were limited to energies of several MeV/u. A database of simulations with the PARTRAC biophysical tool is presented for H, He, Li, Be, B and C isotopes at energies from 0.5 GeV/u down to stopping. The doses deposited to a cell nucleus and also the yields per unit dose of single- and double-strand breaks and their clusters induced in cellular DNA are predicted to vary among diverse isotopes of the same element at energies < 1 MeV/u, especially for isotopes of H and He. The results may affect the risk estimates for astronauts in deep space missions or the models of biological effectiveness of ion beams and indicate that radiation protection in 14C or 10Be dating techniques may be based on knowledge gathered with 12C or 9Be.

Estimating long-term health risks after breast cancer radiotherapy: merging evidence from low and high doses.

In breast cancer radiotherapy, substantial radiation exposure of organs other than the treated breast cannot be avoided, potentially inducing second primary cancer or heart disease. While distant organs and large parts of nearby ones receive doses in the mGy-Gy range, small parts of the heart, lung and bone marrow often receive doses as high as 50 Gy. Contemporary treatment planning allows for considerable flexibility in the distribution of this exposure. To optimise treatment with regards to long-term health risks, evidence-based risk estimates are required for the entire broad range of exposures. Here, we thus propose an approach that combines data from medical and epidemiological studies with different exposure conditions. Approximating cancer induction as a local process, we estimate organ cancer risks by integrating organ-specific dose-response relationships over the organ dose distributions. For highly exposed organ parts, specific high-dose risk models based on studies with medical exposure are applied. For organs or their parts receiving relatively low doses, established dose-response models based on radiation-epidemiological data are used. Joining the models in the intermediate dose range leads to a combined, in general non-linear, dose response supported by data over the whole relevant dose range. For heart diseases, a linear model consistent with high- and low-dose studies is presented. The resulting estimates of long-term health risks are largely compatible with rate ratios observed in randomised breast cancer radiotherapy trials. The risk models have been implemented in a software tool PASSOS that estimates long-term risks for individual breast cancer patients.

Toxicity of internal mammary irradiation in breast cancer. Are concerns still justified in times of modern treatment techniques?

BACKGROUND: The purpose of this study was to estimate the additional risk of side effects attributed to internal mammary node irradiation (IMNI) as part of regional lymph node irradiation (RNI) in breast cancer patients and to compare it with estimated overall survival (OS) benefit from IMNI. MATERIAL AND METHODS: Treatment plans (n = 80) with volumetric modulated arc therapy (VMAT) were calculated for 20 patients (4 plans per patient) with left-sided breast cancer from the prospective GATTUM trial in free breathing (FB) and in deep inspiration breath hold (DIBH). We assessed doses to organs at risk ((OARs) lung, contralateral breast and heart) during RNI with and without additional IMNI. Based on the OAR doses, the additional absolute risks of 10-year cardiac mortality, pneumonitis, and secondary lung and breast cancer were estimated using normal tissue complication probability (NTCP) and risk models assuming different age and risk levels. RESULTS: IMNI notably increased the mean OAR doses. The mean heart dose increased upon IMNI by 0.2-3.4 Gy (median: 1.9 Gy) in FB and 0.0-1.5 Gy (median 0.4 Gy) in DIBH. However, the estimated absolute additional 10-year cardiac mortality caused by IMNI was <0.5% for all patients studied except 70-year-old high risk patients (0.2-2.4% in FB and 0.0-1.1% in DIBH). In comparison to this, the published oncological benefit of IMNI ranges between 3.3% and 4.7%. The estimated additional 10-year risk of secondary cancer of the lung or contralateral breast ranged from 0-1.5% and 0-2.8%, respectively, depending on age and risk levels. IMNI increased the pneumonitis risk in all groups (0-2.2%). CONCLUSION: According to our analyses, the published oncological benefit of IMNI outweighs the estimated risk of cardiac mortality even in case of (e.g., cardiac) risk factors during VMAT. The estimated risk of secondary cancer or pneumonitis attributed to IMNI is low. DIBH reduces the estimated additional risk of IMNI even further and should be strongly considered especially in patients with a high baseline risk.

First independent validation of the proton-boron capture therapy concept.

Boron has been suggested to enhance the biological effectiveness of proton beams in the Bragg peak region via the p + 11B → 3α nuclear capture reaction. However, a number of groups have observed no such enhancement in vitro or questioned its proposed mechanism recently. To help elucidate this phenomenon, we irradiated DU145 prostate cancer or U-87 MG glioblastoma cells by clinical 190 MeV proton beams in plateau or Bragg peak regions with or without 10B or 11B isotopes added as sodium mercaptododecaborate (BSH). The results demonstrate that 11B but not 10B or other components of the BSH molecule enhance cell killing by proton beams. The enhancement occurs selectively in the Bragg peak region, is present for boron concentrations as low as 40 ppm, and is not due to secondary neutrons. The enhancement is likely initiated by proton-boron capture reactions producing three alpha particles, which are rare events occurring in a few cells only, and their effects are amplified by intercellular communication to a population-level response. The observed up to 2-3-fold reductions in survival levels upon the presence of boron for the studied prostate cancer or glioblastoma cells suggest promising clinical applications for these tumour types.

Track structure based modelling of chromosome aberrations after photon and alpha-particle irradiation.

A computational model of radiation-induced chromosome aberrations in human cells within the PARTRAC Monte Carlo simulation framework is presented. The model starts from radiation-induced DNA damage assessed by overlapping radiation track structures with multi-scale DNA and chromatin models, ranging from DNA double-helix in atomic resolution to chromatin fibre loops, heterochromatic and euchromatic regions, and chromosome territories. The repair of DNA double-strand breaks via non-homologous end-joining is followed. Initial spatial distribution and complexity, diffusive motion, enzymatic processing, synapsis and ligation of individual DNA ends from the breaks are simulated. To enable scoring of different chromosome aberration types resulting from improper joining of DNA fragments, the repair module has been complemented by tracking the chromosome origin of the ligated fragments and the positions of centromeres. The modelled motion of DNA ends has sub-diffusive characteristics and corresponds to measured chromatin mobility within time-scales of a few hours. The calculated formation of dicentrics after photon and α-particle irradiation in human fibroblasts is compared to experimental data (Cornforth et al., 2002, Radiat Res 158, 43). The predicted yields of dicentrics overestimate the measurements by factors of five for γ-rays and two for α-particle irradiation. Nevertheless, the observed relative dependence on radiation dose is correctly reproduced. Calculated yields and size distributions of other aberration types are discussed. The present work represents a first mechanistic approach to chromosome aberrations and their kinetics, combining full track structure simulations with detailed models of chromatin and accounting for the kinetics of DNA repair.

Mechanistic modelling suggests that the size of preneoplastic lesions is limited by intercellular induction of apoptosis in oncogenically transformed cells.

Selective removal of oncogenically transformed cells by apoptosis induced via signalling by surrounding cells has been suggested to represent a natural anticarcinogenic process. To investigate its potential effect in detail, a mechanistic model of this process is proposed. The model is calibrated against in vitro data on apoptosis triggered in transformed cells by defined external inducers as well as through signalling by normal cells under coculture conditions. The model predicts that intercellular induction of apoptosis is capable of balancing the proliferation of oncogenically transformed cells and limiting the size of their populations over long times, even if their proliferation per se were unlimited. Experimental research is desired to verify whether the predicted stable population of transformed cells corresponds to a kind of dormancy during early-stage carcinogenesis (dormant preneoplastic lesions), and how this process relates to other anticarcinogenic mechanisms taking place under in vivo conditions.

MODELLING GLOBAL CARBON AND RADIOCARBON CYCLES.

Carbon cycle receives growing attention, in particular in connection with the climate change. Radiocarbon (14C) serves not only as the well-known basis of a dating technique but also as a tracer of the global carbon cycle, enabling one to assess the sizes of diverse compartments, fluxes between them and the related characteristic times. Mathematical modelling of the carbon cycle helps integrate the measurements, estimate the roles of underpinning processes and provide predictions, for instance on future CO2 concentrations in the atmosphere for various emission scenarios. We present a model based on a single-box atmosphere, ocean surface layer, one-dimensional diffusive ocean and two-box biota. We discuss its validation against measured data, predictions on future CO2 levels and interpretation of past events on the radiocarbon calibration curve.

EXPERIMENTAL IN VITRO DOSIMETRY OF 223RA AND 177LU.

Targeted alpha therapy with radionuclides undergoing multiple alpha-particle decays is a promising method of nuclear medicine. To study the effectiveness of alpha versus beta emitters, survival of DU145 prostate cancer cells exposed to 223Ra or 177Lu was assessed. Per decay, the cells were much more sensitive to the alpha than beta emitter. However, per unit dose the sensitivities would be comparable, contrary to the well-known evidence, if the decay energy were deposited within the sample completely and homogeneously. Measurements by Timepix detectors showed about three times higher counts of alpha particles above than below the sample. After the first alpha decay of 223Ra to 219Rn, this gas likely moves upwards and its subsequent three alpha decays occur in the upper part of the sample. Correct estimation of absorbed dose is a critical issue when analysing in vitro data and when translating their results to clinical applications.

Anatomy-dependent lung doses from 3D-conformal breast-cancer radiotherapy.

This study aims to identify key anatomic features that govern the individual variability of lung doses from breast-cancer radiotherapy. 3D conformal, intensity-modulated and hybrid techniques with 50.4 Gy whole-breast dose were planned for 128 patients. From their CT images, 17 anatomic measures were assessed and tested as predictors for lung dose-volume characteristics. Tangential techniques yielded mean ipsilateral lung doses in the range of 3-11 Gy. This inter-patient variability was explained to almost 40% by central lung distance, and to almost 60% if this measure was complemented by midplane lung width and maximum heart distance. Also the variability in further dose-volume metrics such as volume fractions receiving 5, 20 or 40 Gy could be largely explained by the anatomy. Multi-field intensity-modulated radiotherapy reduced high-exposed lung volumes, but resulted in higher mean ipsilateral lung doses and larger low-dose burden. Contralateral lung doses ranged from 0.3 to 1 Gy. The results highlight that there are large differences in lung doses among breast-cancer patients. Most of this inter-individual variability can be explained by a few anatomic features. The results will be implemented in a dedicated software tool to provide personalized estimates of long-term health risks related to breast-cancer radiotherapy. The results may also be used to identify favourable as well as problematic anatomies, and serve as a quick quantitative benchmark for individual treatment plans.

BORON-ENHANCED BIOLOGICAL EFFECTIVENESS OF PROTON IRRADIATION: STRATEGY TO ASSESS THE UNDERPINNING MECHANISM.

Proton radiotherapy for the treatment of cancer offers an excellent dose distribution. Cellular experiments have shown that in terms of biological effects, the sharp dose distribution is further amplified, by as much as 75%, in the presence of boron. It is a matter of debate whether the underlying physical processes involve the nuclear reaction of 11B with protons or 10B with secondary neutrons, both producing densely ionizing short-ranged particles. Likewise, potential roles of intercellular communication or boron acting as a radiosensitizer are not clear. We present an ongoing research project based on a multiscale approach to elucidate the mechanism by which boron enhances the effectiveness of proton irradiation in the Bragg peak. It combines experimental with simulation tools to study the physics of proton-boron interactions, and to analyze intra- and inter-cellular boron biology upon proton irradiation.

RADIATION DAMAGE TO DNA PLASMIDS IN THE PRESENCE OF BOROCAPTATES.

Boron derivatives have great potential in cancer diagnostics and treatment. Borocaptates are used in boron neutron capture therapy and potentially in proton boron fusion therapy. This work examines modulation effects of two borocaptate compounds on radiation-induced DNA damage. Aqueous solutions of pBR322 plasmid containing increasing concentrations of borocaptates were irradiated with 60Co gamma rays or 30 MeV protons. Induction of single and double DNA strand breaks was investigated using agarose gel electrophoresis. In this model system, representing DNA without the intervention of cellular repair mechanisms, the boron derivatives acted as antioxidants. Clinically relevant boron concentrations of 40 ppm reduced the DNA single strand breakage seven-fold. Possible mechanisms of the observed effect are discussed.

RADIOCARBON DATING OF CHARCOALS FROM HISTORICAL MORTARS FROM TÝROV AND PYSOLEC CASTLES.

Organic inclusions in lime binders provide useful samples for radiocarbon dating of historical objects. Two Czech castles Týrov and Pysolec from Late Middle Ages were explored, and tens of charcoals were found in their walls. The radiocarbon content of the charcoals was measured with accelerator mass spectrometry. The dating results showed that none of the charcoals were younger than the known historical ages (Týrov: 1260 - 1270, Pysolec: 1300 - 1340), but some were considerably older. Two charcoals from Pysolec castle dated to Palaeolithic, likely originating from fluvial sediments added as an aggregate to the mortar. When excluding these two charcoals, the others indicated most likely dates being 50-100 y older than the building dates of the castles. This systemic effect corresponds to the age of wood used for lime burning and shall be accounted for when dating mortars using charcoals.

Track structures, DNA targets and radiation effects in the biophysical Monte Carlo simulation code PARTRAC.

This review describes the PARTRAC suite of comprehensive Monte Carlo simulation tools for calculations of track structures of a variety of ionizing radiation qualities and their biological effects. A multi-scale target model characterizes essential structures of the whole genomic DNA within human fibroblasts and lymphocytes in atomic resolution. Calculation methods and essential results are recapitulated regarding the physical, physico-chemical and chemical stage of track structure development of radiation damage induction. Recent model extension towards DNA repair processes extends the time dimension by about 12 orders of magnitude and paves the way for superior predictions of radiation risks.

Potential Morbidity Reduction for Lung Stereotactic Body Radiation Therapy Using Respiratory Gating.

We investigated the potential of respiratory gating to mitigate the motion-caused misdosage in lung stereotactic body radiotherapy (SBRT). For fourteen patients with lung tumors, we investigated treatment plans for a gating window (GW) including three breathing phases around the maximum exhalation phase, GW40-60. For a subset of six patients, we also assessed a preceding three-phase GW20-40 and six-phase GW20-70. We analyzed the target volume, lung, esophagus, and heart doses. Using normal tissue complication probability (NTCP) models, we estimated radiation pneumonitis and esophagitis risks. Compared to plans without gating, GW40-60 significantly reduced doses to organs at risk without impairing the tumor doses. On average, the mean lung dose decreased by 0.6 Gy (p < 0.001), treated lung V20Gy by 2.4% (p = 0.003), esophageal dose to 5cc by 2.0 Gy (p = 0.003), and maximum heart dose by 3.2 Gy (p = 0.009). The model-estimated mean risks of 11% for pneumonitis and 12% for esophagitis without gating decreased upon GW40-60 to 7% and 9%, respectively. For the highest-risk patient, gating reduced the pneumonitis risk from 43% to 32%. Gating is most beneficial for patients with high-toxicity risks. Pre-treatment toxicity risk assessment may help optimize patient selection for gating, as well as GW selection for individual patients.

Analytical formulas representing track-structure simulations on DNA damage induced by protons and light ions at radiotherapy-relevant energies.

Track structure based simulations valuably complement experimental research on biological effects of ionizing radiation. They provide information at the highest level of detail on initial DNA damage induced by diverse types of radiation. Simulations with the biophysical Monte Carlo code PARTRAC have been used for testing working hypotheses on radiation action mechanisms, for benchmarking other damage codes and as input for modelling subsequent biological processes. To facilitate such applications and in particular to enable extending the simulations to mixed radiation field conditions, we present analytical formulas that capture PARTRAC simulation results on DNA single- and double-strand breaks and their clusters induced in cells irradiated by ions ranging from hydrogen to neon at energies from 0.5 GeV/u down to their stopping. These functions offer a means by which radiation transport codes at the macroscopic scale could easily be extended to predict biological effects, exploiting a large database of results from micro-/nanoscale simulations, without having to deal with the coupling of spatial scales and running full track-structure calculations.

Longitudinal atherosclerotic changes after radio(chemo)therapy of hypopharyngeal carcinoma.

BACKGROUND: Radiotherapy treatment of head and neck cancer affects local arteries and increases the risk of stroke. This study aimed at a closer characterization of this damage and its development in time with a longitudinal study set up. METHODS: Male patients treated between 2011 and 2016 for hypopharyngeal carcinoma were identified from the in-house clinical data base. They were included into the study if besides the planning CT at least one additional CT image was available from follow-up (13 patients) or at least two MRI scans (16 patients of which 2 were already included). All patients received radiotherapy, and chemotherapy was administered to 16 patients. The time from the beginning of radiotherapy to the last available image ranged from 2 months to 4.5 years. For six segments of the carotid arteries, the number and volume of atherosclerotic plaques were determined from the CT scans, and the intima media thickness from the MRI scans. Information on comorbid cardiovascular disease, hypertension and diabetes mellitus was retrieved from medical records. RESULTS: Total plaque volume rose from 0.25 cm3 before to 0.33 cm3 after therapy but this was not significant (p = 0.26). The mean number of plaques increased from 5.7 to 8.1 (p = 0.002), and the intima media thickened from 1.17 mm to 1.35 mm (p = 0.002). However, the mean intima media thickness practically did not change in patients with comorbid diabetes mellitus (p-value for homogeneity: 0.03). For patients without diabetes mellitus, dynamics of both plaque number and intima media thickness, was consistent with an increase until about one year after therapy and no further progression thereafter. CONCLUSION: Our study confirmed the thickening of artery walls and the increase in the number of plaques. Results imply that definitive radiation damage to the artery walls can be determined not earlier than about one year after radiotherapy and there is no substantial deterioration thereafter. Reasons for the absence of an observable intima media thickening in patients with diabetes are unclear.

MECHANISTIC MODELING PREDICTS ANTI-CARCINOGENIC RADIATION EFFECTS ON INTERCELLULAR SIGNALING IN VITRO TURN PRO-CARCINOGENIC IN VIVO.

Oncogenic transformed cells represent an in vitro system mimicking early-stage carcinogenesis. These precancerous cells are subject to a selective removal via apoptosis induced by neighbor cells. By modulating the underpinning intercellular signaling mediated by cytokines and reactive oxygen/nitrogen species, ionizing radiation enhances this removal of precancerous cells in vitro, at doses from a few mGy to a few Gy. However, epidemiological data demonstrate that radiation exposure induces cancer, at least above 100 mGy. Mechanistic modeling of the given anti-carcinogenic process explains this discrepancy: The model reproduces in vitro data on apoptosis and its enhancement by radiation. For in vivo-like conditions with signal lifetimes shorter and cell densities higher than in vitro, radiation is predicted to reduce this anti-carcinogenic mechanism. Early-stage lesions that would be turned dormant or completely removed may grow large and escape this control mechanism upon irradiation.

BIOPHYSICAL SIMULATION TOOL PARTRAC: MODELLING PROTON BEAMS AT THERAPY-RELEVANT ENERGIES.

The biophysical simulation tool PARTRAC has been primarily developed to model radiation physics, chemistry and biology on nanometre to micrometre scales. However, the tool can be applied in simulating radiation effects in an event-by-event manner over macroscopic volumes as well. Benchmark simulations are reported showing that PARTRAC does reproduce the macroscopic Bragg peaks of proton beams, although the penetration depths are underestimated by a few per cent for high-energy beams. PARTRAC also quantifies the increase in DNA damage and its complexity along the beam penetration depth. Enhanced biological effectiveness is predicted in particular within distal Bragg peak parts of therapeutic proton beams.

Track-structure simulations of energy deposition patterns to mitochondria and damage to their DNA.

PURPOSE: Mitochondria have been implicated in initiating and/or amplifying the biological effects of ionizing radiation not mediated via damage to nuclear DNA. To help elucidate the underlying mechanisms, energy deposition patterns to mitochondria and radiation damage to their DNA have been modelled. METHODS: Track-structure simulations have been performed with PARTRAC biophysical tool for 60Co γ-rays and 5 MeV α-particles. Energy deposition to the cell's mitochondria has been analyzed. A model of mitochondrial DNA reflecting experimental information on its structure has been developed and used to assess its radiation-induced damage. RESULTS: Energy deposition to mitochondria is highly inhomogeneous, especially at low doses. Although a dose-dependent fraction of mitochondria sees no energy deposition at all, the hit ones receive rather high amounts of energy. Nevertheless, only little damage to mitochondrial DNA occurs, even at large doses. CONCLUSION: Mitochondrial DNA does not represent a critical target for radiation effects. Likely, the key role of mitochondria in radiation-induced biological effects arises from the communication between mitochondria and/or with the nucleus. Through this signaling, initial modifications in a few heavily hit mitochondria seem to be amplified to a massive long-term effect manifested in the whole cell or even tissue.

WHAT ANATOMIC FEATURES GOVERN PERSONAL LONG-TERM HEALTH RISKS FROM BREAST CANCER RADIOTHERAPY?

Breast cancer radiotherapy may in the long term lead to radiation-induced secondary cancer or heart disease. These health risks hugely vary among patients, partially due to anatomy-driven differences in doses deposited to the heart, ipsilateral lung and contralateral breast. We identify four anatomic features that largely cover these dosimetric variations to enable personalized risk estimates. For three exemplary, very different risk scenarios, the given parameter set reproduces 63-74% of the individual risk variability for left-sided breast cancer patients. These anatomic features will be used in the PASSOS software to support decision processes in breast-cancer therapy.