F. Ballarini

First steps towards systems radiation biology studies concerned with DNA and chromosome structure within living cells

For the understanding of radiation action on biological systems like cellular macromolecules (e.g., DNA in its higher structures) a synergistic approach of experiments and quantitative modelling of working hypotheses is necessary. Further on, the influence on calculated results of certain assumptions in such working hypotheses must critically be evaluated. In the present work, this issue is highlighted in two aspects for the case of DNA damage in single cells. First, yields of double-strand breaks and frequency distributions of DNA fragment lengths after ion irradiation were calculated using different assumptions on the DNA target model. Compared to a former target model now a moderate effect due to the inclusion of a spherical chromatin domain model has been found. Second, the influence of assumptions on particular geometric chromosome models on calculated chromosome aberration data is illustrated with two target-modelling approaches for this end point.

DNA DSB induced in human cells by charged particles and gamma rays: Experimental results and theoretical approaches

Purpose:To quantify the role played by radiation track structure and background fragments in modulating DNA fragmentation in human cells exposed to γ-rays and light ions. Materials and methods: Human fibroblasts were exposed in vitro to different doses (in the range from 40 – 200 Gy) of 60Co γ-rays and 0.84 MeV protons (Linear Energy Transfer, LET, in tissue 28.5 keV/μm). The resulting DNA fragments were scored under two electrophoretic conditions, in order to optimize separation in the size ranges 0.023 – 1.0 Mbp and 1.0 – 5.7 Mbp. In parallel, DNA fragmentation was simulated both with a phenomenological approach based on the “generalized broken-stick” model, and with a mechanistic approach based on the PARTRAC (acronym of PARticle TRACk) Monte Carlo code (1.32 MeV photons were used for the simulation of 60Co γ-rays). Results: For both γ-rays and protons, the experimental dose response in the range 0.023 – 5.7 Mbp could be approximated as a straight line, the slope of which provided a yield of (5.3 ± 0.4) • 10−9 Gy−1 bp−1 for γ-rays and (7.1 ± 0.6) • 10−9 Gy−1 bp−1 for protons, leading to a Relative Biological Effectiveness (RBE) of 1.3 ± 0.2. From both theoretical analyses it appeared that, while γ-ray data were consistent with double-strand breaks (DSB) random induction, protons at low doses showed significant deviation from randomness, implying enhanced production of small fragments in the low molecular weight part of the experimental range. The theoretical analysis of fragment production was then extended to ranges where data were not available, i.e. to fragments larger than 5.7 Mbp and smaller than 23 kbp. The main outcome was that small fragments (<23 kbp) are produced almost exclusively via non-random processes, since their number is considerably higher than that produced by a random insertion of DSB. Furthermore, for protons the number of these small fragments is a significant fraction (about 20%) of the total number of fragments; these fragments remain undetected in these experiments. Calculations for 3.3 MeV alpha particle irradiation (for which no experimental data were available) were performed to further investigate the role of fragments smaller than 23 kbp; in this case, besides the non-random character of their production, their number resulted to be at least as much as half of the total number of fragments. Conclusion: Comparison between experimental data and two different theoretical approaches provided further support to the hypothesis of an important role of track structure in modulating DNA damage. According to the theoretical approaches, non-randomness of fragment production was found for proton irradiation for the smaller fragments in the experimental size range and, in a significantly larger extent, for fragments of size less than 23 kbp, both for protons and alpha particles.

MRI-guided neutron capture therapy by use of a dual gadolinium/boron agent targeted at tumour cells through upregulated low-density lipoprotein transporters.

The upregulation of low-density lipoprotein (LDL) transporters in tumour cells has been exploited to deliver a sufficient amount of gadolinium/boron/ligand (Gd/B/L) probes for neutron capture therapy, a binary chemio-radiotherapy for cancer treatment. The Gd/B/L probe consists of a carborane unit (ten B atoms) bearing an aliphatic chain on one side (to bind LDL particles), and a Gd(III)/1,4,7,10-tetraazacyclododecane monoamide complex on the other (for detection by magnetic resonance imaging (MRI)). Up to 190 Gd/B/L probes were loaded per LDL particle. The uptake from tumour cells was initially assessed on cell cultures of human hepatoma (HepG2), murine melanoma (B16), and human glioblastoma (U87). The MRI assessment of the amount of Gd/B/L taken up by tumour cells was validated by inductively coupled plasma-mass-spectrometric measurements of the Gd and B content. Measurements were undertaken in vivo on mice bearing tumours in which B16 tumour cells were inoculated at the base of the neck. From the acquisition of magnetic resonance images, it was established that after 4-6 hours from the administration of the Gd/B/L-LDL particles (0.1 and 1 mmol kg(-1) of Gd and (10)B, respectively) the amount of boron taken up in the tumour region is above the threshold required for successful NCT treatment. After neutron irradiation, tumour growth was followed for 20 days by MRI. The group of treated mice showed markedly lower tumour growth with respect to the control group.

A Model of Chromosome Aberration Induction: Applications to Space Research

Abstract Ballarini, F. and Ottolenghi, A. A Model of Chromosome Aberration Induction: Applications to Space Research. Radiat. Res. 164, 567–570 (2005). A mechanistic model and Monte Carlo code simulating chromosome aberration induction in human lymphocytes is presented. The model is based on the assumption that aberrations arise from clustered DNA lesions and that only the free ends of clustered lesions created in neighboring chromosome territories or in the same territory can join and produce exchanges. The lesions are distributed in the cell nucleus according to the radiation track structure. Interphase chromosome territories are modeled as compact intranuclear regions with volumes proportional to the chromosome DNA contents. Both Giemsa staining and FISH painting can be simulated, and background aberrations can be taken into account. The good agreement with in vitro data provides validation of the model in terms of both the assumptions adopted and the simulation techniques. As an application in the field of space research, the model predictions were compared with aberration yields measured among crew members of long-term missions on board Mir and ISS, assuming an average radiation quality factor of 2.4. The agreement obtained also validated the model for in vivo exposure scenarios and suggested possible applications to the prediction of other relevant aberrations, typically translocations.

Design, development and characterization of multi-functionalized gold nanoparticles for biodetection and targeted boron delivery in BNCT applications

The aim of this study is to optimize targeted boron delivery to cancer cells and its tracking down to the cellular level. To this end, we describe the design and synthesis of novel nanovectors that double as targeted boron delivery agents and fluorescent imaging probes. Gold nanoparticles were coated with multilayers of polyelectrolytes functionalized with the fluorescent dye (FITC), boronophenylalanine and folic acid. In vitro confocal fluorescence microscopy demonstrated significant uptake of the nanoparticles in cancer cells that are known to overexpress folate receptors.

DNA fragmentation induced in human fibroblasts by 56Fe ions: Experimental data and monte carlo simulations

Abstract Campa, A., Alloni, D., Antonelli, F., Ballarini, F., Belli, M., Dini, V., Esposito, G., Facoetti, A., Friedland, W., Furusawa, Y., Liotta, M., Ottolenghi, A., Paretzke, H. G., Simone, G., Sorrentino, E. and Tabocchini, M. A. DNA Fragmentation Induced in Human Fibroblasts by 56Fe Ions: Experimental Data and Monte Carlo Simulations. Radiat. Res. 171, 438–445 (2009). We studied the DNA fragmentation induced in human fibroblasts by iron-ion beams of two different energies: 115 MeV/nucleon and 414 MeV/nucleon. Experimental data were obtained in the fragment size range 1–5700 kbp; Monte Carlo simulations were performed with the PARTRAC code; data analysis was also performed through the Generalized Broken Stick (GBS) model. The comparison between experimental and simulated data for the number of fragments produced in two different size ranges, 1–23 kbp and 23–5700 kbp, gives a satisfactory agreement for both radiation qualities. The Monte Carlo simulations also allow the counting of fragments outside the experimental range: The number of fragments smaller than 1 kbp is large for both beams, although with a strong difference between the two cases. As a consequence, we can compute different RBEs depending on the size range considered for the fragment counting. The PARTRAC evaluation takes into account fragments of all sizes, while the evaluation from the experimental data considers only the fragments in the range of 1–5700 kbp. When the PARTRAC evaluation is restricted to this range, the agreement between experimental and computed RBE values is again good. When fragments smaller than 1 kbp are also considered, the RBE increases considerably, since γ rays produce a small number of such fragments. The analysis performed with the GBS model proved to be quite sensitive to showing, with a phenomenological single parameter, variations in double-strand break (DSB) correlation.

Nuclear Models in FLUKA: Present Capabilities, Open Problems, and Future Improvements

The nuclear reaction models embedded in the FLUKA code cover hadron, ion, photon and neutrino induced nuclear interactions from energies as low as few tens of MeV up to several tens of TeV. A short description of the main physics ingredients in the FLUKA nuclear models is given, with emphasis on the intermediate energy range and on “exotic” reactions. The treatment of electromagnetic dissociation as recently implemented in FLUKA is described. Examples of performances are presented for illustrative situations covering some of the most typical FLUKA applications.

The application of FLUKA to dosimetry and radiation therapy

The FLUKA Monte Carlo code has been evolving over the last several decades and is now widely used for radiation shielding calculations. In order to facilitate the use of FLUKA in dosimetry and therapy applications, supporting software has been developed to allow the direct conversion of the output files from standard CT-scans directly into a voxel geometry for transport within FLUKA. Since the CT-scan information essentially contains only the electron density information over the scanned volume, one needs the specific compositions for each voxel individually. We present here the results of a simple algorithm to assign tissues in the human body to one of four categories: soft-tissue, hard-bone, trabecular-bone and porous-lung. In addition, we explore the problem of the pathlength distributions in porous media such as trabecular bone. A mechanism will be implemented within FLUKA to allow for variable multipal fixed density materials to accommodate the pathlength distributions discovered.

Boron uptake measurements in a rat model for Boron Neutron Capture Therapy of lung tumours.

Lung carcinoma is the leading cause of cancer mortality in the Western countries. Despite the introduction over the last few years of new therapeutic agents, survival from lung cancer has shown no discernible improvement in the last 20 years. For these reasons any efforts to find and validate new effective therapeutic procedures for lung cancer are very timely. The selective boron uptake in the tumour with respect to healthy tissues makes Boron Neutron Capture Therapy a potentially advantageous option in the treatment of tumours that affect whole vital organs, and that are surgically inoperable. To study the possibility of applying BNCT to the treatment of diffuse pulmonary tumours, an animal model for boron uptake measurements in lung metastases was developed. Both healthy and tumour-bearing rats were infused with Boronophenylalanine (BPA) and sacrificed at different time intervals after drug administration. The lungs were extracted, and prepared for boron analysis by neutron autoradiography and α-spectroscopy. The boron concentrations in tumour and normal lung were plotted as a function of the time elapsed after BPA administration. The concentration in tumour is almost constant within the error bars for all the time intervals of the experiment (1-8 h), while the curve in normal lung decreases after 4 h from BPA infusion. At 4 h, the ratio of boron concentration in tumour to boron concentration in healthy lung is higher than 3, and it stays above this level up to 8 h. Also the images of boron distribution in the samples, obtained by neutron autoradiography, show a selective absorption in the metastases.

Calculations of dose distributions in the lungs of a rat model irradiated in the thermal column of the TRIGA reactor in Pavia

To test the possibility to apply boron neutron capture therapy (BNCT) to lung tumors, some rats are planned to be irradiated in the thermal column of the TRIGA reactor of the University of Pavia. Before the irradiation, lung metastases will be induced in BDIX rats, which will be subsequently infused with boronophenylalanine (BPA). During the irradiation, the rats will be positioned in a box designed to shield the whole animal except the thorax area. In order to optimize the irradiation set-up and to design a suitable shielding box, a set of calculations were performed with the MCNP Monte Carlo transport code. A rat model was constructed using the MCNP geometry capabilities and was positioned in a box with walls filled with lithium carbonate. A window was opened in front of the lung region. Different shapes of the holder and of the window were tested and analyzed in terms of the dose distribution obtained in the lungs and of the dose absorbed by the radiosensitive organs in the rat. The best configuration of the holder ensures an almost uniform thermal neutron flux inside the lungs (Phi(max)/Phi(min)=1.5), an irradiation time about 10 min long, to deliver at least 40 Gy(w) to the tumor, a mean lung dose of 5.9+/-0.4 Gy(w), and doses absorbed by all the other healthy tissues below the tolerance limits.