Multiscale approach predictions for biological outcomes in ion-beam cancer therapy

Effects of BPQ binding on the nonadiabatic dynamics of excited electrons in poly(dG)-poly(dC) DNA under proton irradiation

The interaction between DNA and small molecules is important for understanding the mechanisms of DNA-based multifunctional devices and has been extensively studied. However, there are few reports on such interactions in irradiation environments. Here, we investigate the nonadiabatic dynamic behaviors of excited electrons in double-stranded DNA bound to BPQ molecule upon proton irradiation, focusing on the energy deposition and electronic excitation dynamics as protons traverse the DNA along different channels. Our results reveal that the presence of BPQ significantly influences charge migration and DNA damage, with notable differences between CGvdW and CGcovalent adducts. The energy deposition process is highly dependent on charge density, and guanine exhibits higher excitation propensity than cytosine due to its structural characteristics. The BPQ molecule enhances DNA charge migration and promotes damage through secondary electron migration. These findings provide insights into the nonadiabatic dynamics of DNA under ionizing radiation and have implications for designing targeted electrophilic organics to improve radiotherapy efficacy.

Modelling Cellular Response to Ionizing Radiation: Mechanistic, Semi-Mechanistic, and Phenomenological Approaches - A Historical Perspective

Radiation research is a multidisciplinary field, and among its many branches, mathematical and computational modelers have played a significant role in advancing boundaries of knowledge. A fundamental contribution is modelling cellular response to ionizing radiation as that is the key to not only understanding how radiation can kill cancer cells, but also cause cancer and other health issues. The invention of microdosimetry in the 1950s by Harold Rossi paved the way for brilliant scientists to study the mechanism of radiation at cellular and sub-cellular scales. This paper reviews some snippets of ingenious mathematical and computational models published in microdosimetry symposium proceedings and publications of the radiation research community. Among these are simulations of radiation tracks at atomic and molecular levels using Monte Carlo methods, models of cell survival, quantification of the amount of energy required to create a single strand break, and models of DNA-damage-repair. These models can broadly be categorized into mechanistic, semi-mechanistic, and phenomenological approaches, and this review seeks to provide historical context of their development. We salute pioneers of the field and great teachers who supported and educated the younger members of the community and showed them how to build upon their work.

Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment

This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.

Dynamic interplay between thionine and DNA under carbon ion irradiation: a real-time first-principles study

Understanding the interactions between deoxyribonucleic acid (DNA) and photosensitizer under ion irradiation benefits the development of aptasensors, DNA biosensors and cancer diagnosis. Using real-time time-depended density functional theory, by simulating high-energy C ion passing through DNA with poly(dG)·poly(dC) sequence and that with embedded thionine (3,7-diamino-5-phenothiazinium, TH), we compared the electronic stopping power (ESP), evolution of the structure and charge, and absorption spectrum. TH inserting leads the increase in space charge density, a larger electron de-excitation and a larger ESP, but the speed corresponding to the maximum ESP is almost same. When C ion passes through TH-DNA, the structure of TH slightly changes and there still exists noncovalent interaction between TH and DNA, but the absorption coefficient depends on the electron occupied state of TH when the ion passes through. These results indicate that at low radiation doses, TH still can be a DNA detector, although its response wavelength and intensity have been slightly changed, and provide a theoretical reference to improve the possible application of phenothiazine dye in DNA biosensor under ion irradiation.

The Mayo Clinic Florida Microdosimetric Kinetic Model of Clonogenic Survival: Application to Various Repair-Competent Rodent and Human Cell Lines

The relative biological effectiveness (RBE) calculations used during the planning of ion therapy treatments are generally based on the microdosimetric kinetic model (MKM) and the local effect model (LEM). The Mayo Clinic Florida MKM (MCF MKM) was recently developed to overcome the limitations of previous MKMs in reproducing the biological data and to eliminate the need for ion-exposed in vitro data as input for the model calculations. Since we are considering to implement the MCF MKM in clinic, this article presents (a) an extensive benchmark of the MCF MKM predictions against corresponding in vitro clonogenic survival data for 4 rodent and 10 cell lines exposed to ions from 1H to 238U, and (b) a systematic comparison with published results of the latest version of the LEM (LEM IV). Additionally, we introduce a novel approach to derive an approximate value of the MCF MKM model parameters by knowing only the animal species and the mean number of chromosomes. The overall good agreement between MCF MKM predictions and in vitro data suggests the MCF MKM can be reliably used for the RBE calculations. In most cases, a reasonable agreement was found between the MCF MKM and the LEM IV.

Multiscale modeling of stochastic dynamics processes with MBN Explorer

Stochastic dynamics describes processes in complex systems having the probabilistic nature. They can involve very different dynamical systems and occur on very different temporal and spatial scale. This paper discusses the concept of stochastic dynamics and its implementation in the popular program MBN Explorer. Stochastic dynamics in MBN Explorer relies on the Monte Carlo approach and permits simulations of physical, chemical, and biological processes. The paper presents the basic theoretical concepts underlying stochastic dynamics implementation and provides several examples highlighting its applicability to different systems, such as diffusing proteins seeking an anchor point on a cell membrane, deposition of nanoparticles on a surface leading to structures with fractal morphologies, and oscillations of compounds in an autocatalytic reaction. The chosen examples illustrate the diversity of applications that can be modeled by means of stochastic dynamics with MBN Explorer.

Energy Deposition around Swift Carbon-Ion Tracks in Liquid Water

Energetic carbon ions are promising projectiles used for cancer radiotherapy. A thorough knowledge of how the energy of these ions is deposited in biological media (mainly composed of liquid water) is required. This can be attained by means of detailed computer simulations, both macroscopically (relevant for appropriately delivering the dose) and at the nanoscale (important for determining the inflicted radiobiological damage). The energy lost per unit path length (i.e., the so-called stopping power) of carbon ions is here theoretically calculated within the dielectric formalism from the excitation spectrum of liquid water obtained from two complementary approaches (one relying on an optical-data model and the other exclusively on ab initio calculations). In addition, the energy carried at the nanometre scale by the generated secondary electrons around the ion's path is simulated by means of a detailed Monte Carlo code. For this purpose, we use the ion and electron cross sections calculated by means of state-of-the art approaches suited to take into account the condensed-phase nature of the liquid water target. As a result of these simulations, the radial dose around the ion's path is obtained, as well as the distributions of clustered events in nanometric volumes similar to the dimensions of DNA convolutions, contributing to the biological damage for carbon ions in a wide energy range, covering from the plateau to the maximum of the Bragg peak.

Lethal DNA damage caused by ion-induced shock waves in cells

The elucidation of fundamental mechanisms underlying ion-induced radiation damage of biological systems is crucial for advancing radiotherapy with ion beams and for radiation protection in space. The study of ion-induced biodamage using the phenomenon-based multiscale approach (MSA) to the physics of radiation damage with ions has led to the prediction of nanoscale shock waves created by ions in a biological medium at the high linear energy transfer (LET). The high-LET regime corresponds to the keV and higher-energy losses by ions per nanometer, which is typical for ions heavier than carbon at the Bragg peak region in biological media. This paper reveals that the thermomechanical stress of the DNA molecule caused by the ion-induced shock wave becomes the dominant mechanism of complex DNA damage at the high-LET ion irradiation. Damage of the DNA molecule in water caused by a projectile-ion-induced shock wave is studied by means of reactive molecular dynamics simulations. Five projectile ions (carbon, oxygen, silicon, argon, and iron) at the Bragg peak energies are considered. For the chosen segment of the DNA molecule and the collision geometry, the number of DNA strand breaks is evaluated for each projectile ion as a function of the bond dissociation energy and the distance from the ion's path to the DNA strands. Simulations reveal that argon and especially iron ions induce the breakage of multiple bonds in a DNA double convolution containing 20 DNA base pairs. The DNA damage produced in segments of such size leads to complex irreparable lesions in a cell. This makes the shock-wave-induced thermomechanical stress the dominant mechanism of complex DNA damage at the high-LET ion irradiation. A detailed theory for evaluating the DNA damage caused by ions at high-LET is formulated and integrated into the MSA formalism. The theoretical analysis reveals that a single ion hitting a cell nucleus at high-LET is sufficient to produce highly complex, lethal damages to a cell by the shock-wave-induced thermomechanical stress. Accounting for the shock-wave-induced thermomechanical mechanism of DNA damage provides an explanation for the "overkill" effect observed experimentally in the dependence of cell survival probabilities on the radiation dose delivered with iron ions. This important observation provides strong experimental evidence of the ion-induced shock-wave effect and the related mechanism of radiation damage in cells.

Update of the particle irradiation data ensemble (PIDE) for cell survival

The particle irradiation data ensemble (PIDE) is the largest database of cell survival data measured after exposure to ion beams and photon reference radiation. We report here on the updated version of the PIDE database and demonstrate how to investigate generic properties of radiation dose response using these sets of raw data. The database now contains information of over 1100 pairs of photon and ion dose response curves. It provides the originally published raw data of cell survival in addition to given linear quadratic (LQ) model parameters. If available, the raw data were used to derive LQ model parameters in the same way for all experiments. To demonstrate the extent of the database and the variability among experiments we focus on the dose response curves after ion and photon radiation separately in a first step. Furthermore, we discuss the capability and the limitations of the database for analyzing properties of low and high linear energy transfer (LET) radiation response based on multiple experiments. PIDE is freely available to the research community under www.gsi.de/bio-pide.

Relative Role of Physical Mechanisms on Complex Biodamage Induced by Carbon Irradiation

The effective use of swift ion beams in cancer treatment (known as hadrontherapy) as well as appropriate protection in manned space missions rely on the accurate understanding of the energy delivery to cells that damages their genetic information. The key ingredient characterizing the response of a medium to the perturbation induced by charged particles is its electronic excitation spectrum. By using linear-response time-dependent density functional theory, we obtained the energy and momentum transfer excitation spectrum (the energy-loss function, ELF) of liquid water (the main constituent of biological tissues), which was in excellent agreement with experimental data. The inelastic scattering cross sections obtained from this ELF, together with the elastic scattering cross sections derived by considering the condensed phase nature of the medium, were used to perform accurate Monte Carlo simulations of the energy deposited by swift carbon ions in liquid water and carried away by the generated secondary electrons, producing inelastic events such as ionization, excitation, and dissociative electron attachment (DEA). The latter are strongly correlated with cellular death, which is scored in sensitive volumes with the size of two DNA convolutions. The sizes of the clusters of damaging events for a wide range of carbon-ion energies, from those relevant to hadrontherapy up to those for cosmic radiation, predict with unprecedented statistical accuracy the nature and relative magnitude of the main inelastic processes contributing to radiation biodamage, confirming that ionization accounts for the vast majority of complex damage. DEA, typically regarded as a very relevant biodamage mechanism, surprisingly plays a minor role in carbon-ion induced clusters of harmful events.

5-Nitro-2,4-Dichloropyrimidine as an Universal Model for Low-Energy Electron Processes Relevant for Radiosensitization

We report experimental results of low-energy electron interactions with.

Modeling the effect of ion-induced shock waves and DNA breakage with the reactive CHARMM force field

Ion-induced DNA damage is an important effect underlying ion beam cancer therapy. This article introduces the methodology of modeling DNA damage induced by a shock wave caused by a projectile ion. Specifically it is demonstrated how single- and double strand breaks in a DNA molecule could be described by the reactive CHARMM (rCHARMM) force field implemented in the program MBN Explorer. The entire workflow of performing the shock wave simulations, including obtaining the crucial simulation parameters, is described in seven steps. Two exemplary analyses are provided for a case study simulation serving to: (a) quantify the shock wave propagation and (b) describe the dynamics of formation of DNA breaks. The article concludes by discussing the computational cost of the simulations and revealing the possible maximal computational time for different simulation set-ups.

Ionizing radiation and natural constituents of living cells: Low-energy electron interaction with coenzyme Q analogs

Resonance electron attachment to short-tail analogs of coenzyme Q₁₀ is investigated in the electron energy range 0 eV-14 eV under gas-phase conditions by means of dissociative electron attachment spectroscopy. Formation of long-lived (milliseconds) molecular negative ions is detected at 1.2 eV, but not at thermal energy. A huge increase in the electron detachment time as compared with the reference para-benzoquinone (40 µs) is ascribed to the presence of the isoprene side chains. Elimination of a neutral CH₃ radical is found to be the most intense decay detected on the microsecond time scale. The results give some insight into the timescale of electron-driven processes stimulated in living tissues by high-energy radiation and are of importance in prospective fields of radiobiology and medicine.

Ionizing Radiation and Complex DNA Damage: Quantifying the Radiobiological Damage Using Monte Carlo Simulations

Ionizing radiation is a common tool in medical procedures. Monte Carlo (MC) techniques are widely used when dosimetry is the matter of investigation. The scientific community has invested, over the last 20 years, a lot of effort into improving the knowledge of radiation biology. The present article aims to summarize the understanding of the field of DNA damage response (DDR) to ionizing radiation by providing an overview on MC simulation studies that try to explain several aspects of radiation biology. The need for accurate techniques for the quantification of DNA damage is crucial, as it becomes a clinical need to evaluate the outcome of various applications including both low- and high-energy radiation medical procedures. Understanding DNA repair processes would improve radiation therapy procedures. Monte Carlo simulations are a promising tool in radiobiology studies, as there are clear prospects for more advanced tools that could be used in multidisciplinary studies, in the fields of physics, medicine, biology and chemistry. Still, lot of effort is needed to evolve MC simulation tools and apply them in multiscale studies starting from small DNA segments and reaching a population of cells.

Optimization of treatment planning for hypoxic tumours and re-modulation of radiation intensity in heavy-ion radiotherapy

Aim

The purpose of this study is to optimize treatment planning in carbon ion radiotherapy, taking into account the effect of tumour hypoxia.

Background

In conventional hadron therapy, the goal is to create a homogenous dose in the tumour area and, thus, achieve a uniform cell survival level. Since the induction of a specific damage to cells is directly influenced by the level of hypoxia in the tissue, the varying oxygen pressure in the different regions of hypoxic tumours would disrupt the uniformity of the cell survival level.

Materials and methods

Using the Geant4 Monte Carlo Code, the physical dose profile and dose-averaged linear energy transfer were calculated in the tumour. Then, the oxygen enhancement ratio in different areas of the tumour were compared with different pressures.

Results

Modulations of radiation intensities as well as energies of ion beams were calculated, both considering and disregarding the effect of hypoxia, and the required dose profiles were compared with each other. Cell survival levels were also compared between the two methods. An equation was obtained for re-modulating the beams in the presence of hypoxia, and radiation weighting factors were extracted for the beam intensities.

Conclusion

The results show that taking the effect of hypoxia into account would cause the reduction of average doses delivered to the tumour tissues up to 1.54 times. In this regard, the required dose is reduced by 1.63 times in the healthy tissues before the tumour. This will result in an effective protection of healthy tissues around the tumour.

Reaction of Electrons with DNA: Radiation Damage to Radiosensitization

This review article provides a concise overview of electron involvement in DNA radiation damage. The review begins with the various states of radiation-produced electrons: Secondary electrons (SE), low energy electrons (LEE), electrons at near zero kinetic energy in water (quasi-free electrons, (e-qf)) electrons in the process of solvation in water (presolvated electrons, e-pre), and fully solvated electrons (e-aq). A current summary of the structure of e-aq, and its reactions with DNA-model systems is presented. Theoretical works on reduction potentials of DNA-bases were found to be in agreement with experiments. This review points out the proposed role of LEE-induced frank DNA-strand breaks in ion-beam irradiated DNA. The final section presents radiation-produced electron-mediated site-specific formation of oxidative neutral aminyl radicals from azidonucleosides and the evidence of radiosensitization provided by these aminyl radicals in azidonucleoside-incorporated breast cancer cells.

Renormalization of radiobiological response functions by energy loss fluctuations and complexities in chromosome aberration induction: deactivation theory for proton therapy from cells to tumor control

We employ a multi-scale mechanistic approach built upon our recent phenomenological/computational methodologies [R. Abolfath et al., Sci. Rep. 7, 8340 (2017)] to investigate radiation induced cell toxicities and deactivation mechanisms as a function of linear energy transfer in hadron therapy. Our theoretical model consists of a system of Markov chains in microscopic and macroscopic spatio-temporal landscapes, i.e., stochastic birth-death processes of cells in millimeter-scale colonies that incorporates a coarse-grained driving force to account for microscopic radiation induced damage. The coupling, hence the driving force in this process, stems from a nano-meter scale radiation induced DNA damage that incorporates the enzymatic end-joining repair and mis-repair mechanisms. We use this model for global fitting of the high-throughput and high accuracy clonogenic cell-survival data acquired under exposure of the therapeutic scanned proton beams, the experimental design that considers γ-H2AX as the biological endpoint and exhibits maximum observed achievable dose and LET, beyond which the majority of the cells undergo collective biological deactivation processes. An estimate to optimal dose and LET calculated from tumor control probability by extension to ~106 cells per mm-size voxels is presented. We attribute the increase in degree of complexity in chromosome aberration to variabilities in the observed biological responses as the beam linear energy transfer (LET) increases, and verify consistency of the predicted cell death probability with the in vitro cell survival assay of approximately 100 non-small cell lung cancer (NSCLC) cells. The present model provides an interesting interpretation to variabilities in α and β indices via perturbative expansion of the cell survival fraction (SF) in terms of specific and lineal energies, z and y, corresponding to continuous transitions in pair-wise to ternary, quaternary and more complex recombination of broken chromosomes from the entrance to the end of the range of proton beam.

Mechanistic Modelling of Radiation Responses

Radiobiological modelling has been a key part of radiation biology and therapy for many decades, and many aspects of clinical practice are guided by tools such as the linear-quadratic model. However, most of the models in regular clinical use are abstract and empirical, and do not provide significant scope for mechanistic interpretation or making predictions in novel cell lines or therapies. In this review, we will discuss the key areas of ongoing mechanistic research in radiation biology, including physical, chemical, and biological steps, and review a range of mechanistic modelling approaches which are being applied in each area, highlighting the possible opportunities and challenges presented by these techniques.

First principles simulation of damage to solvated nucleotides due to shock waves

We present a first-principles molecular dynamics study of the effect of shock waves (SWs) propagating in a model biological medium. We find that the SW can cause chemical modifications through varied and complex mechanisms, in particular, phosphate-sugar and sugar-base bond breaks. In addition, the SW promotes the dissociation of water molecules, thus enhancing the ionic strength of the medium. Freed protons can hydrolyze base and sugar rings previously opened by the shock. However, many of these events are only temporary, and bonds reform rapidly. Irreversible damage is observed for pressures above 15-20 GPa. These results are important to gain a better understanding of the microscopic damage mechanisms underlying cosmic-ray irradiation in space and ion-beam cancer therapy.

A New Standard DNA Damage (SDD) Data Format

Our understanding of radiation-induced cellular damage has greatly improved over the past few decades. Despite this progress, there are still many obstacles to fully understand how radiation interacts with biologically relevant cellular components, such as DNA, to cause observable end points such as cell killing. Damage in DNA is identified as a major route of cell killing. One hurdle when modeling biological effects is the difficulty in directly comparing results generated by members of different research groups. Multiple Monte Carlo codes have been developed to simulate damage induction at the DNA scale, while at the same time various groups have developed models that describe DNA repair processes with varying levels of detail. These repair models are intrinsically linked to the damage model employed in their development, making it difficult to disentangle systematic effects in either part of the modeling chain. These modeling chains typically consist of track-structure Monte Carlo simulations of the physical interactions creating direct damages to DNA, followed by simulations of the production and initial reactions of chemical species causing so-called "indirect" damages. After the induction of DNA damage, DNA repair models combine the simulated damage patterns with biological models to determine the biological consequences of the damage. To date, the effect of the environment, such as molecular oxygen (normoxic vs. hypoxic), has been poorly considered. We propose a new standard DNA damage (SDD) data format to unify the interface between the simulation of damage induction in DNA and the biological modeling of DNA repair processes, and introduce the effect of the environment (molecular oxygen or other compounds) as a flexible parameter. Such a standard greatly facilitates inter-model comparisons, providing an ideal environment to tease out model assumptions and identify persistent, underlying mechanisms. Through inter-model comparisons, this unified standard has the potential to greatly advance our understanding of the underlying mechanisms of radiation-induced DNA damage and the resulting observable biological effects when radiation parameters and/or environmental conditions change.

Modelling of Cellular Survival Following Radiation-Induced DNA Double-Strand Breaks

A mechanistic model of cellular survival following radiation-induced DNA double-strand breaks (DSBs) was proposed in this study. DSBs were assumed as the initial lesions in the DNA of the cell nucleus induced by ionizing radiation. The non-homologous end-joining (NHEJ) pathway was considered as the domain pathway of DSB repair in mammalian cells. The model was proposed to predict the relationship between radiation-induced DSBs in nucleus and probability of cell survival, which was quantitatively described by two input parameters and six fitting parameters. One input parameter was the average number of primary particles which caused DSB, the other input parameter was the average number of DSBs yielded by each primary particle that caused DSB. The fitting parameters were used to describe the biological characteristics of the irradiated cells. By determining the fitting parameters of the model with experimental data, the model is able to estimate surviving fractions for the same type of cells exposed to particles with different physical parameters. The model further revealed the mechanism of cell death induced by the DSB effect. Relative biological effectiveness (RBE) of charged particles at different survival could be calculated with the model, which would provide reference for clinical treatment.

Activation of the DNA-repair mechanism through NBS1 and MRE11 diffusion

The non-homologous end joining of a DNA double strand break is initiated by the MRE11-NBS1-RAD50 complex whose subunits are the first three proteins to arrive to the breakage site thereby making the recruitment time of MRE11, NBS1 and RAD50 essential for cell survival. In the present investigation, the nature of MRE11 and NBS1 transportation from the cytoplasm to the nucleus, hosting the damaged DNA strand, is hypothesized to be a passive diffusive process. The feasibility of such a mechanism is addressed through theoretical and computational approaches which permit establishing the characteristic recruitment time of MRE11 and NBS1 by the nucleus. A computational model of a cell is constructed from a set of biological parameters and the kinetic Monte Carlo algorithm is used to simulate the diffusing MRE11 and NBS1 particles as a random walk process. To accurately describe the experimented data, it is discovered that MRE11 and NBS1 should start diffusion from significantly different starting positions which suggests that diffusion might not be the only transport mechanism of repair protein recruitment to the DNA break.

The DNA repair mechanism is crucial for a cell to avoid apoptosis, and is a complicated process involving many different repair proteins. The mean of transportation of these repair proteins is largely unknown as their transportation mechanisms need clarification. We have focused on the transportation of some of the first proteins to arrive to the damaged DNA from the cytoplasm, namely the MRE11 and NBS1. Our hypothesis of diffusion as the transportation mechanism of MRE11 and NBS1 is tested against a theoretical model as well as simulation data and experiments, revealing a possibility for diffusion to be the method of recruiting MRE11 and NBS1 to the DNA break.

Proton relative biological effectiveness (RBE): a multiscale problem

Proton radiotherapy is undergoing rapid expansion both within the UK and internationally, but significant challenges still need to be overcome if maximum benefit is to be realised from this technique. One major limitation is the persistent uncertainty in proton relative biological effectiveness (RBE). While RBE values are needed to link proton radiotherapy to our existing experience with photon radiotherapy, RBE remains poorly understood and is typically incorporated as a constant dose scaling factor of 1.1 in clinical plans. This is in contrast to extensive experimental evidence indicating that RBE is a function of dose, tissue type, and proton linear energy transfer, among other parameters. In this article, we discuss the challenges associated with obtaining clinically relevant values for proton RBE through commonly-used assays, and highlight the wide range of other experimental end points which can inform our understanding of RBE. We propose that accurate and robust optimization of proton radiotherapy ultimately requires a multiscale understanding of RBE, integrating subcellular, cellular, and patient-level processes.

Toward A variable RBE for proton beam therapy

In the clinic, proton beam therapy (PBT) is based on the use of a generic relative biological effectiveness (RBE) of 1.1 compared to photons in human cancers and normal tissues. However, the experimental basis for this RBE lacks any significant number of representative tumor models and clinically relevant endpoints for dose-limiting organs at risk. It is now increasingly appreciated that much of the variations of treatment responses in cancers are due to inter-tumoral genomic heterogeneity. Indeed, recently it has been shown that defects in certain DNA repair pathways, which are found in subsets of many cancers, are associated with a RBE increase in vitro. However, there currently exist little in vivo or clinical data that confirm the existence of similarly increased RBE values in human cancers. Furthermore, evidence for variable RBE values for normal tissue toxicity has been sparse and conflicting to date. If we could predict variable RBE values in patients, we would be able to optimally use and personalize PBT. For example, predictive tumor biomarkers may facilitate selection of patients with proton-sensitive cancers previously ineligible for PBT. Dose de-escalation may be possible to reduce normal tissue toxicity, especially in pediatric patients. Knowledge of increased tumor RBE may allow us to develop biologically optimized therapies to enhance local control while RBE biomarkers for normal tissues could lead to a better understanding and prevention of unusual PBT-associated toxicity. Here, we will review experimental data on the repair of proton damage to DNA that impact both RBE values and biophysical modeling to predict RBE variations. Experimental approaches for studying proton sensitivity in vitro and in vivo will be reviewed as well and recent clinical findings discussed. Ultimately, therapeutically exploiting the understudied biological advantages of protons and developing approaches to limit treatment toxicity should fundamentally impact the clinical use of PBT.

Hydroperoxyl radical and formic acid formation from common DNA stabilizers upon low energy electron attachment

2-Amino-2-(hydroxymethyl)-1,3-propanediol (TRIS) and ethylenediaminetetraacetic acid (EDTA) are key components of biological buffers and are frequently used as DNA stabilizers in irradiation studies. Such surface or liquid phase studies are done with the aim to understand the fundamental mechanisms of DNA radiation damage and to improve cancer radiotherapy. When ionizing radiation is used, abundant secondary electrons are formed during the irradiation process, which are able to attach to the molecular compounds present on the surface. In the present study we experimentally investigate low energy electron attachment to TRIS and methyliminodiacetic acid (MIDA), an analogue of EDTA, supported by quantum chemical calculations. The most prominent dissociation channel for TRIS is through hydroperoxyl radical formation, whereas the dissociation of MIDA results in the formation of formic and acetic acid. These compounds are well-known to cause DNA modifications, like strand breaks. The present results indicate that buffer compounds may not have an exclusive protecting effect on DNA as suggested previously.

Nanoparticle radio-enhancement: principles, progress and application to cancer treatment

Enhancement of radiation effects by high-atomic number nanoparticles (NPs) has been increasingly studied for its potential to improve radiotherapeutic efficacy. The underlying principle of NP radio-enhancement is the potential to release copious electrons into a nanoscale volume, thereby amplifying radiation-induced biological damage. While the vast majority of studies to date have focused on gold nanoparticles with photon radiation, an increasing number of experimental, theoretical and simulation studies have explored opportunities offered by other NPs (e.g. gadolinium, platinum, iron oxide, hafnium) and other therapeutic radiation sources such as ion beams. It is thus of interest to the research community to consolidate findings from the different studies and summarise progress to date, as well as to identify strategies that offer promising opportunities for clinical translation. This is the purpose of this Topical Review.

A model for relative biological effectiveness of therapeutic proton beams based on a global fit of cell survival data

We introduce an approach for global fitting of the recently published high-throughput and high accuracy clonogenic cell-survival data for therapeutic scanned proton beams. Our fitting procedure accounts for the correlation between the cell-survival, the absorbed (physical) dose and the proton linear energy transfer (LET). The fitting polynomials and constraints have been constructed upon generalization of the microdosimetric kinetic model (gMKM) adapted to account for the low energy and high lineal-energy spectrum of the beam where the current radiobiological models may underestimate the reported relative biological effectiveness (RBE). The parameters (α, β) of the linear-quadratic (LQ) model calculated by the presented method reveal a smooth transition from low to high LETs which is an advantage of the current method over methods previously employed to fit the same clonogenic data. Finally, the presented approach provides insight into underlying microscopic mechanisms which, with future study, may help to elucidate radiobiological responses along the Bragg curve and resolve discrepancies between experimental data and current RBE models.

Comprehensive track-structure based evaluation of DNA damage by light ions from radiotherapy-relevant energies down to stopping

Track structures and resulting DNA damage in human cells have been simulated for hydrogen, helium, carbon, nitrogen, oxygen and neon ions with 0.25–256 MeV/u energy. The needed ion interaction cross sections have been scaled from those of hydrogen; Barkas scaling formula has been refined, extending its applicability down to about 10 keV/u, and validated against established stopping power data. Linear energy transfer (LET) has been scored from energy deposits in a cell nucleus; for very low-energy ions, it has been defined locally within thin slabs. The simulations show that protons and helium ions induce more DNA damage than heavier ions do at the same LET. With increasing LET, less DNA strand breaks are formed per unit dose, but due to their clustering the yields of double-strand breaks (DSB) increase, up to saturation around 300 keV/μm. Also individual DSB tend to cluster; DSB clusters peak around 500 keV/μm, while DSB multiplicities per cluster steadily increase with LET. Remarkably similar to patterns known from cell survival studies, LET-dependencies with pronounced maxima around 100–200 keV/μm occur on nanometre scale for sites that contain one or more DSB, and on micrometre scale for megabasepair-sized DNA fragments.

Three dimensional approach to investigating biological effects along energetic ion beam pathways

Heavy ion beams have many exciting applications, including radiotherapy of deep-seated tumors and simulation tests of space irradiation for astronauts. These beams often use a feature that concentrates the energy deposition largely along the end of the energy pathway, leading to different distributions of biological effects along the axial direction. Currently, there is relatively little information regarding the radial directional difference of biological effects along the heavy ion paths. This study utilized a filter membrane that was quantatively applied with cells to demonstrate a 3D distribution model of irradiation on biological effects in living organisms. Some results have indicated that there is excitatory effect on the non-irradiated regions with energetic ions, which may give new insights into the distribution of biological effects along the paths of heavy ion beams with mid-high energy.

Production of doubly-charged highly reactive species from the long-chain amino acid GABA initiated by Ar9+ionization

We present a combined experimental and theoretical study of the fragmentation of multiply-charged γ-aminobutyric acid molecules (GABA^z+,z= 2, 3) in the gas phase.

Gold nanoparticles for cancer radiotherapy: a review

Radiotherapy is currently used in around 50% of cancer treatments and relies on the deposition of energy directly into tumour tissue. Although it is generally effective, some of the deposited energy can adversely affect healthy tissue outside the tumour volume, especially in the case of photon radiation (gamma and X-rays). Improved radiotherapy outcomes can be achieved by employing ion beams due to the characteristic energy deposition curve which culminates in a localised, high radiation dose (in form of a Bragg peak). In addition to ion radiotherapy, novel sensitisers, such as nanoparticles, have shown to locally increase the damaging effect of both photon and ion radiation, when both are applied to the tumour area. Amongst the available nanoparticle systems, gold nanoparticles have become particularly popular due to several advantages: biocompatibility, well-established methods for synthesis in a wide range of sizes, and the possibility of coating of their surface with a large number of different molecules to provide partial control of, for example, surface charge or interaction with serum proteins. This gives a full range of options for design parameter combinations, in which the optimal choice is not always clear, partially due to a lack of understanding of many processes that take place upon irradiation of such complicated systems. In this review, we summarise the mechanisms of action of radiation therapy with photons and ions in the presence and absence of nanoparticles, as well as the influence of some of the core and coating design parameters of nanoparticles on their radiosensitisation capabilities.

Measurement of complex DNA damage induction and repair in human cellular systems after exposure to ionizing radiations of varying linear energy transfer (LET)

Detrimental effects of ionizing radiation (IR) are correlated to the varying efficiency of IR to induce complex DNA damage. A double strand break (DSB) can be considered the simpler form of complex DNA damage. These types of damage can consist of DSBs, single strand breaks (SSBs) and/or non-DSB lesions such as base damages and apurinic/apyrimidinic (AP; abasic) sites in different combinations. Enthralling theoretical (Monte Carlo simulations) and experimental evidence suggests an increase in the complexity of DNA damage and therefore repair resistance with linear energy transfer (LET). In this study, we have measured the induction and processing of DSB and non-DSB oxidative clusters using adaptations of immunofluorescence. Specifically, we applied foci colocalization approaches as the most current methodologies for the in situ detection of clustered DNA lesions in a variety of human normal (FEP18-11-T1) and cancerous cell lines of varying repair efficiency (MCF7, HepG2, A549, MO59K/J) and radiation qualities of increasing LET, that is γ-, X-rays 0.3-1 keV/μm, α-particles 116 keV/μm and ³⁶Ar ions 270 keV/μm. Using γ-H2AX or 53BP1 foci staining as DSB probes, we calculated a DSB apparent rate of 5-16 DSBs/cell/Gy decreasing with LET. A similar trend was measured for non-DSB oxidized base lesions detected using antibodies against the human repair enzymes 8-oxoguanine-DNA glycosylase (OGG1) or AP endonuclease (APE1), that is damage foci as probes for oxidized purines or abasic sites, respectively. In addition, using colocalization parameters previously introduced by our groups, we detected an increasing clustering of damage for DSBs and non-DSBs. We also make correlations of damage complexity with the repair efficiency of each cell line and we discuss the biological importance of these new findings with regard to the severity of IR due to the complex nature of its DNA damage.