Marie Davídková

Track structure modeling in liquid water: A review of the Geant4-DNA very low energy extension of the Geant4 Monte Carlo simulation toolkit

Understanding the fundamental mechanisms involved in the induction of biological damage by ionizing radiation remains a major challenge of todays radiobiology research. The Monte Carlo simulation of physical, physicochemical and chemical processes involved may provide a powerful tool for the simulation of early damage induction. The Geant4-DNA extension of the general purpose Monte Carlo Geant4 simulation toolkit aims to provide the scientific community with an open source access platform for the mechanistic simulation of such early damage. This paper presents the most recent review of the Geant4-DNA extension, as available to Geant4 users since June 2015 (release 10.2 Beta). In particular, the review includes the description of new physical models for the description of electron elastic and inelastic interactions in liquid water, as well as new examples dedicated to the simulation of physicochemical and chemical stages of water radiolysis. Several implementations of geometrical models of biological targets are presented as well, and the list of Geant4-DNA examples is described.

Radiation-induced oxidative damage to the DNA-binding domain of the lactose repressor

Understanding the cellular effects of radiation-induced oxidation requires the unravelling of key molecular events, particularly damage to proteins with important cellular functions. The Escherichia coli lactose operon is a classical model of gene regulation systems. Its functional mechanism involves the specific binding of a protein, the repressor, to a specific DNA sequence, the operator. We have shown previously that upon irradiation with gamma-rays in solution, the repressor loses its ability to bind the operator. Water radiolysis generates hydroxyl radicals (OH* radicals) which attack the protein. Damage of the repressor DNA-binding domain, called the headpiece, is most likely to be responsible of this loss of function. Using CD, fluorescence spectroscopy and a combination of proteolytic cleavage with MS, we have examined the state of the irradiated headpiece. CD measurements revealed a dose-dependent conformational change involving metastable intermediate states. Fluorescence measurements showed a gradual degradation of tyrosine residues. MS was used to count the number of oxidations in different regions of the headpiece and to narrow down the parts of the sequence bearing oxidized residues. By calculating the relative probabilities of reaction of each amino acid with OH. radicals, we can predict the most probable oxidation targets. By comparing the experimental results with the predictions we conclude that Tyr7, Tyr12, Tyr17, Met42 and Tyr47 are the most likely hotspots of oxidation. The loss of repressor function is thus correlated with chemical modifications and conformational changes of the headpiece.

Investigating the Implications of a Variable RBE on Proton Dose Fractionation Across a Clinical Pencil Beam Scanned Spread-Out Bragg Peak

Purpose To investigate the clinical implications of a variable relative biological effectiveness (RBE) on proton dose fractionation. Using acute exposures, the current clinical adoption of a generic, constant cell killing RBE has been shown to underestimate the effect of the sharp increase in linear energy transfer (LET) in the distal regions of the spread-out Bragg peak (SOBP). However, experimental data for the impact of dose fractionation in such scenarios are still limited. Methods and Materials Human fibroblasts (AG01522) at 4 key depth positions on a clinical SOBP of maximum energy 219.65 MeV were subjected to various fractionation regimens with an interfraction period of 24 hours at Proton Therapy Center in Prague, Czech Republic. Cell killing RBE variations were measured using standard clonogenic assays and were further validated using Monte Carlo simulations and parameterized using a linear quadratic formalism. Results Significant variations in the cell killing RBE for fractionated exposures along the proton dose profile were observed. RBE increased sharply toward the distal position, corresponding to a reduction in cell sparing effectiveness of fractionated proton exposures at higher LET. The effect was more pronounced at smaller doses per fraction. Experimental survival fractions were adequately predicted using a linear quadratic formalism assuming full repair between fractions. Data were also used to validate a parameterized variable RBE model based on linear α parameter response with LET that showed considerable deviations from clinically predicted isoeffective fractionation regimens. Conclusions The RBE-weighted absorbed dose calculated using the clinically adopted generic RBE of 1.1 significantly underestimates the biological effective dose from variable RBE, particularly in fractionation regimens with low doses per fraction. Coupled with an increase in effective range in fractionated exposures, our study provides an RBE dataset that can be used by the modeling community for the optimization of fractionated proton therapy.

Heterochromatinization associated with cell differentiation as a model to study DNA double strand break induction and repair in the context of higher-order chromatin structure

Cell differentiation is associated with extensive gene silencing, heterochromatinization and potentially decreasing need for repairing DNA double-strand breaks (DSBs). Differentiation stages of blood cells thus represent an excellent model to study DSB induction, repair and misrepair in the context of changing higher-order chromatin structure. We show that immature granulocytes form γH2AX and 53BP1 foci, contrary to the mature cells; however, these foci colocalize only rarely and DSB repair is inefficient. Moreover, specific chromatin structure of granulocytes probably influences DSB induction.

Function of chromatin structure and dynamics in DNA damage, repair and misrepair: γ-rays and protons in action.

According to their physical characteristics, protons and ion beams promise a revolution in cancer radiotherapy. Curing protocols however reflect rather the empirical knowledge than experimental data on DNA repair. This especially holds for the spatio-temporal organization of repair processes in the context of higher-order chromatin structure-the problematics addressed in this work. The consequences for the mechanism of chromosomal translocations are compared for gamma rays and proton beams.

A High-Power Laser-Driven Source of Sub-nanosecond Soft X-Ray Pulses for Single-Shot Radiobiology Experiments

Abstract Davídková, M., Juha, L., Bittner, M., Koptyaev, S., Krása, J., Pfeifer, M., Štísová, V., Bartnik, A., Fiedorowicz, H., Mikolajczyk, J., Ryc, L., Pína, L., Horváth, M., Babánková, D., Cihelka, J. and Civiš, S. A High-Power Laser-Driven Source of Sub-nanosecond Soft X-Ray Pulses for Single-Shot Radiobiology Experiments. Radiat. Res. 168, 382–387 (2007). A large-scale, double-stream gas puff target has been illuminated by sub-kJ, near-infrared (NIR) focused laser pulses at the PALS facility (Prague Asterix Laser System) to produce high-energy pulses of soft X rays from hot, dense plasma. The double-puff arrangement ensures high gas density and conversion efficiency from NIR to X rays approaching that typical for solid targets. In addition, its major advantage over solid targets is that it is free of debris and has substantially suppressed charged-particle emission. The X-ray emission characteristics of the source were determined for a range of gases that included krypton, xenon, N2, CO and N2-CO. A demonstrated application of the xenon-based source is a single-shot damage induction to plasmid DNA. The yields of single-strand breaks (SSBs) and double-strand breaks (DSBs) were determined as a function of energy fluence adjusted by varying distance of sample from the source and thickness of aluminum filters.

Radiation-induced tetramer-to-dimer transition of Escherichia coli lactose repressor

The wild type lactose repressor of Escherichia coli is a tetrameric protein formed by two identical dimers. They are associated via a C-terminal 4-helix bundle (called tetramerization domain) whose stability is ensured by the interaction of leucine zipper motifs. Upon in vitro gamma-irradiation the repressor losses its ability to bind the operator DNA sequence due to damage of its DNA-binding domains. Using an engineered dimeric repressor for comparison, we show here that irradiation induces also the change of repressor oligomerisation state from tetramer to dimer. The splitting of the tetramer into dimers can result from the oxidation of the leucine residues of the tetramerization domain.

Lethal events in V79 cells irradiated by low-energy protons and correlations with distribution patterns of energy deposition, radical concentration and DNA damage.

Published survival data of V79 cells irradiated by 0.5-5.0 MeV (7-40 keV/microm) protons have been analyzed with a detailed radiobiological model to estimate the per-track yields of lethal lesions. Their correlations with distribution patterns of deposited energy, radical concentrations and with the yields of specific classes of DNA damage have been studied. The observed correlations indicate a potential interpretation of DNA damage lethal for the cell and the initial physical and chemical processes leading to such damage.

Contribution of indirect effects to clustered damage in DNA irradiated with protons.

Protons are the dominant particles both in galactic cosmic rays and in solar particle events and, furthermore, proton irradiation becomes increasingly used in tumour treatment. It is believed that complex DNA damage is the determining factor for the consequent cellular response to radiation. DNA plasmid pBR322 was irradiated at U120-M cyclotron with 30 MeV protons and treated with two Escherichia coli base excision repair enzymes. The yields of SSBs and DSBs were analysed using agarose gel electrophoresis. DNA has been irradiated in the presence of hydroxyl radical scavenger (coumarin-3-carboxylic acid) in order to distinguish between direct and indirect damage of the biological target. Pure scavenger solution was used as a probe for measurement of induced OH· radical yields. Experimental OH· radical yield kinetics was compared with predictions computed by two theoretical models-RADAMOL and Geant4-DNA. Both approaches use Geant4-DNA for description of physical stages of radiation action, and then each of them applies a distinct model for description of the pre-chemical and chemical stage.

Impact of oxygen concentration on yields of DNA damages caused by ionizing radiation

Local hypoxia-induced radiation resistance is one of the major problems in current radiation therapy of solid tumors. Concentration of cellular oxygen is not uniform within tumor volume and changes during the course of radiotherapy. The aim of this work is to evaluate the influence of oxygen concentration on the water radiolysis process and subsequently on the yields of primary DNA damages caused by ionizing radiation. Monte Carlo based computer modeling approach was used for this purpose. Detailed track structures provided by TRIOL code and molecular structure of 100 bp DNA oligomer generated using Amber 8 molecular dynamics package were used as input for stochastic model RADAMOL. Diffusion and chemical reactions of radicals produced in the water radiolysis process were followed in a step-by-step simulation. Both unscavengeable and scavengeable DNA damage was taken into account. It was shown that concentration of dilluted oxygen in water significantly influences yields of eaq-, H and O- and consequently modifies yields of primary DNA damage. However, this itself cannot explain the oxygen radiosensitizing effect in living cells. The influence of oxygen on DNA damage chemical repair and fixation should be included in further research.