Search PubMed⌕ Search

Biomedical subjects

F Ballarini

Publications and source records attributed to F Ballarini.

At least 19 recordsLinked to original sources

The FLUKA code: new developments and application to 1 GeV/n iron beams.

The modeling of ion transport and interactions in matter is a subject of growing interest, driven by the continuous increase of possible application fields. These include hadron therapy, dosimetry, and space missions, but there are also several issues involving fundamental research, accelerator physics, and cosmic ray physics, where a reliable description of heavy ion induced cascades is important. In the present work, the capabilities of the FLUKA code for ion beams will be briefly recalled and some recent developments presented. Applications of the code to the simulation of therapeutic carbon, nitrogen and oxygen ion beams, and of iron beams, which are of direct interest for space mission related experiments, will be also presented together with interesting consideration relative to the evaluation of dosimetric quantities. Both applications involve ion beams in the AGeV range.

Carbon↗

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 gamma-rays and light ions. MATERIALS AND METHODS: Human fibroblasts were exposed in vitro to different doses (in the range from 40 - 200 Gy) of (60)Co gamma-rays and 0.84 MeV protons (Linear Energy Transfer, LET, in tissue 28.5 keV/microm). 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 (60)Co gamma-rays). RESULTS: For both gamma-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) x 10(-9) Gy(-1) bp(-1) for gamma-rays and (7.1 +/- 0.6) x 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 gamma-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.

Cells, Cultured↗

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.

Algorithms↗

The FLUKA code for space applications: recent developments.

The FLUKA Monte Carlo transport code is widely used for fundamental research, radioprotection and dosimetry, hybrid nuclear energy system and cosmic ray calculations. The validity of its physical models has been benchmarked against a variety of experimental data over a wide range of energies, ranging from accelerator data to cosmic ray showers in the earth atmosphere. The code is presently undergoing several developments in order to better fit the needs of space applications. The generation of particle spectra according to up-to-date cosmic ray data as well as the effect of the solar and geomagnetic modulation have been implemented and already successfully applied to a variety of problems. The implementation of suitable models for heavy ion nuclear interactions has reached an operational stage. At medium/high energy FLUKA is using the DPMJET model. The major task of incorporating heavy ion interactions from a few GeV/n down to the threshold for inelastic collisions is also progressing and promising results have been obtained using a modified version of the RQMD-2.4 code. This interim solution is now fully operational, while waiting for the development of new models based on the FLUKA hadron-nucleus interaction code, a newly developed QMD code, and the implementation of the Boltzmann master equation theory for low energy ion interactions.

Computer Simulation↗

Role of shielding in modulating the effects of solar particle events: Monte Carlo calculation of absorbed dose and DNA complex lesions in different organs.

Distributions of absorbed dose and DNA clustered damage yields in various organs and tissues following the October 1989 solar particle event (SPE) were calculated by coupling the FLUKA Monte Carlo transport code with two anthropomorphic phantoms (a mathematical model and a voxel model), with the main aim of quantifying the role of the shielding features in modulating organ doses. The phantoms, which were assumed to be in deep space, were inserted into a shielding box of variable thickness and material and were irradiated with the proton spectra of the October 1989 event. Average numbers of DNA lesions per cell in different organs were calculated by adopting a technique already tested in previous works, consisting of integrating into "condensed-history" Monte Carlo transport codes--such as FLUKA--yields of radiobiological damage, either calculated with "event-by-event" track structure simulations, or taken from experimental works available in the literature. More specifically, the yields of "Complex Lesions" (or "CL", defined and calculated as a clustered DNA damage in a previous work) per unit dose and DNA mass (CL Gy-1 Da-1) due to the various beam components, including those derived from nuclear interactions with the shielding and the human body, were integrated in FLUKA. This provided spatial distributions of CL/cell yields in different organs, as well as distributions of absorbed doses. The contributions of primary protons and secondary hadrons were calculated separately, and the simulations were repeated for values of Al shielding thickness ranging between 1 and 20 g/cm2. Slight differences were found between the two phantom types. Skin and eye lenses were found to receive larger doses with respect to internal organs; however, shielding was more effective for skin and lenses. Secondary particles arising from nuclear interactions were found to have a minor role, although their relative contribution was found to be larger for the Complex Lesions than for the absorbed dose, due to their higher LET and thus higher biological effectiveness.

Astronauts↗

Models of chromosome aberration induction: an example based on radiation track structure.

A few examples of models of chromosome aberration induction are summarised and discussed on the basis of the three main theories of aberration formation, that is "breakage-and-reunion", "exchange" and "one-hit". A model and code developed at the Universities of Milan and Pavia is then presented in detail. The model provides dose-response curves for different aberration types (dicentrics, translocations, rings, complex exchanges and deletions) induced in human lymphocytes by gamma rays, protons and alpha particles of different energies, both as monochromatic fields and as mixed fields. The main assumptions are that only clustered - and thus severe - DNA breaks ("Complex Lesions", CL) can participate in the production of aberrations, and that only break free ends in neighbouring chromosome territories can interact and form exchanges. The yields of CLs induced by the various radiation types of interest are taken from a previous modelling work. These lesions are distributed within a sphere representing the cell nucleus according to the radiation track structure, e.g. randomly for gamma rays and along straight lines for light ions. Interphase chromosome territories are explicitly simulated and configurations are obtained in which each chromosome occupies an intranuclear domain with volume proportional to its DNA content. In order to allow direct comparisons with experimental data, small fragments can be neglected since usually they cannot be detected in experiments. The presence of a background level of aberrations is also taken into account. The results of the simulations are in good agreement with experimental dose-response curves available in the literature, that provides a validation of the model both in terms of the adopted assumptions and in terms of the simulation techniques. To address the question of "true" incompleteness, simulations were also run in which all fragments were assumed to be visible.

Animals↗

Chromosome aberrations as biomarkers of radiation exposure: modelling basic mechanisms.

The space radiation environment is a mixed field consisting of different particles having different energies, including high charge and energy (HZE) ions. Conventional measurements of absorbed doses may not be sufficient to completely characterise the radiation field and perform reliable estimates of health risks. Biological dosimetry, based on the observation of specific radiation-induced endpoints (typically chromosome aberrations), can be a helpful approach in case of monitored exposure to space radiation or other mixed fields, as well as in case of accidental exposure. Furthermore, various ratios of aberrations (e.g. dicentric chromosomes to centric rings and complex exchanges to simple exchanges) have been suggested as possible fingerprints of radiation quality, although all of them have been subjected to some criticisms. In this context a mechanistic model and a Monte Carlo code for the simulation of chromosome aberration induction were developed. The model, able to provide dose-responses for different aberrations (e.g. dicentrics, rings, fragments, translocations, insertions and other complex exchanges), was further developed to assess the dependence of various ratios of aberrations on radiation quality. The predictions of the model were compared with available data, whose experimental conditions were faithfully reproduced. Particular attention was devoted to the scoring criteria adopted in different laboratories and to possible biases introduced by interphase death and mitotic delay. This latter aspect was investigated by taking into account both metaphase data and data obtained with Premature Chromosome Condensation (PCC).

Alpha Particles↗

Modelling study on the protective role of OH radical scavengers and DNA higher-order structures in induction of single- and double-strand break by gamma-radiation.

PURPOSE: To quantify the protective effects of (non-histonic) OH-radical scavengers and DNA higher-order structures in induction of single- (ssbs) and double-strand breaks (dsbs) by gamma-rays. MATERIALS AND METHODS: Spatial distributions of energy depositions by gamma-rays in liquid water were modelled with the track structure modules of the biophysical simulation code PARTRAC. Such distributions were superimposed on different DNA structure models (e.g. linear DNA, SV40 'minichromosomes' and compact chromatin), and direct energy depositions in the sugar-phosphate were considered as potential (direct) ssbs. The diffusion and interaction of the main chemical species produced in liquid water radiolysis were explicitly simulated, and reactions of *OH with the sugar-phosphate were considered as potential (indirect) ssbs. Two ssb on opposite DNA strands within 10 base pairs were considered as one dsb. Yields of ssb and dsb Gy(-1) Dalton(-1) in different DNA target structures were calculated as a function of the *OH mean lifetime, whose inverse value was taken as representative of the scavenging capacity of the DNA environment. RESULTS AND CONCLUSIONS: A further validation of the models implemented in the PARTRAC code has been provided, thus allowing a better understanding of the mechanisms underlying DNA damage. More specifically, the protection due to *OH scavengers was separately quantified with respect to that due to histones and chromatin folding, which could be 'switched off' in the simulations. As expected, for a given value of the environment scavenging capacity, linear DNA was more susceptible to strand breakage than SV40 minichromosomes, which in turn showed higher damage yields with respect to cellular DNA due to the larger accessibility offered to *OH. Furthermore, by increasing the scavenging capacity, the break yields decreased in all structures and tended to coincide with direct damage yields. Very good agreement was found with available experimental data. Comparisons with data on 'nucleoid' DNA (i.e. unfolded and histone-depleted DNA) also suggested that the experimental procedures used to obtain such structures might lower the environment scavenging capacity owing to the loss of cellular scavengers.

DNA↗

Estimating mixed field effects: an application supporting the lack of a non-linear component for chromosome aberration induction by neutrons.

The action of neutron fields on biological structures was investigated on the basis of chromosome aberration induction in human cells. Available experimental data on aberration induction by neutrons and their interaction products were reviewed. Present criteria adopted in neutron radiation protection were discussed. The linear coefficient alpha and the quadratic coefficient beta describing dose-response curves for dicentric chromosomes induced by neutrons of different energies were calculated via integration of experimental data on dicentric induction by photons and charged particles into the Monte Carlo transport code FLUKA. The predicted values of the linear coefficients for neutron beams of different energies showed good agreement with the corresponding experimental values, whereas the data themselves indicated that the neutron quadratic coefficient cannot be obtained by 'averaging' the beta values of recoil ions and other nuclear reaction products. This supports the hypothesis that neutron induced aberrations increase substantially linearly with dose, a question that has been object of debate for a long time and is still open.

Chromosome Aberrations↗

Nuclear architecture and radiation induced chromosome aberrations: models and simulations.

Knowledge of radiation track structure and its interaction with biological targets is a fundamental starting point in understanding the mechanisms underlying the induction of biological damage. In this context Monte Carlo codes are a powerful tool of investigation, allowing one to simulate both track structure and the features of the target(s) of interest at different scales, from nanometres (linear dimensions of DNA) to micrometres (linear dimensions of human cell nuclei and interphase chromosome territories). In the light of recent experimental findings on nuclear architecture, different approaches in modelling chromosome structure and aberration induction are discussed. In particular, a model is presented in which chromosome territories were explicitly described as subnuclear regions and aberration induction was modelled by coupling the structure of the target with that of the radiation track. Comparisons between model predictions and experimental results from the literature are also reported.

Cell Nucleus↗

Met protein and hepatocyte growth factor (HGF) in papillary carcinoma of the thyroid: evidence for a pathogenetic role in tumourigenesis.

In the last 10 years, evidence has accumulated that overexpression of Met protein is a distinguishing feature of almost every case of well-differentiated papillary carcinoma. Increased expression of the protein is probably due to enhanced transcription of the MET gene and/or to post-transcriptional mechanisms. So far, alterations of the MET gene have not been recognized, but evidence has been provided that activated RAS and RET can cause accumulation of MET RNA. Thus, the possibility exists that dysregulation of MET is the final result of different molecular pathways capable of inducing thyroid cell transformation; RET rearrangements might account for some of the cases, but the demonstration that the majority of papillary carcinomas do not have recognized alterations of the RET gene strongly suggests that MET gene dysregulation can also be achieved through other molecular pathways. Dysregulation of MET causes marked accumulation of Met protein in tumour cells that is promptly detected by immunohistochemistry. Thus, overexpression of Met protein might represent an immunohistochemical marker of papillary carcinoma, potentially helpful in problematic cases, but caution is required; moderate expression of Met protein is observed in non-neoplastic thyroid diseases, such as Graves' and Hashimoto's thyroiditis, and reagents active on paraffin sections may have a low affinity and/or low specificity for Met protein, leading to artifactual staining. Met protein-positive papillary carcinoma cells may produce hepatocyte growth factor (HGF) and may activate HGF through the urokinase-type plasminogen activator (uPA) bound to urokinase-type plasminogen activator receptor (uPA-R). Thus, papillary carcinoma cells possess the molecular machinery necessary for a productive HGF/Met interaction. In vitro studies have demonstrated that HGF enhances the motility and invasiveness of tumour cells and induces the synthesis and release of chemokines active in the recruitment of dendritic cells. These observations provide a rational basis for the understanding of two distinguishing features of papillary carcinoma. First, the tumour is often characterized by early metastatic spread to regional lymph nodes and by multifocal involvement of the gland, which suggests highly invasive behaviour. Second, a prominent peritumoural inflammatory reaction is often observed, which suggests cross-talk between tumour cells and the immune system.

Biomarkers, Tumor↗

Modelling chromosomal aberration induction by ionising radiation: the influence of interphase chromosome architecture.

Several advances have been achieved in the knowledge of nuclear architecture and functions during the last decade, thus allowing the identification of interphase chromosome territories and sub-chromosomal domains (e.g. arm and band domains). This is an important step in the study of radiation-induced chromosome aberrations; indeed, the coupling between track-structure simulations and reliable descriptions of the geometrical properties of the target is one of the main tasks in modelling aberration induction by radiation, since it allows one to clarify the role of the initial positioning of two DNA lesions in determining their interaction probability. In the present paper, the main recent findings on nuclear and chromosomal architecture are summarised. A few examples of models based on different descriptions of interphase chromosome organisation (random-walk models, domain models and static models) are presented, focussing on how the approach adopted in modelling the target nuclei and chromosomes can influence the simulation of chromosomal aberration yields. Each model is discussed by taking into account available experimental data on chromosome aberration induction and/or interphase chromatin organisation. Preliminary results from a mechanistic model based on a coupling between radiation track-structure features and explicitly-modelled, non-overlapping chromosome territories are presented.

Animals↗

A Monte Carlo code for a direct estimation of radiation risk.

An example of pragmatic approach for predicting mixed field effects is presented. The method was initially applied adopting the following, commonly used, assumptions: a) radiation risk (typically cancer) is correlated with chromosome aberration induction; b) radiation-induced chromosome-exchange yield can be well described by a linear-quadratic dependence on particle fluences (mostly linear with high-LET radiation), with parameters depending on particle types and energies. Information on monochromatic field radiobiological effects was integrated in a condensed-history Monte Carlo transport code (FLUKA), able to simulate nuclear interactions. The integrated code provides the chromosome aberration yield (and thus an estimation of radiation risk) in each voxel of any irradiated volume, given any external mixed-field irradiation; in the present work, the method was tested for neutron irradiation of a water phantom. FLUKA was then coupled with a geometrical human phantom provided with different radiation shielding, in order to apply this approach to estimate radiation risk in manned space missions.

Chromosome Aberrations↗

Mechanistic bases for modelling space radiation risk and planning radiation protection of astronauts.

The approaches generally adopted for planning radiation protection in ground-based facilities cannot be applied straightforward for astronaut protection in space. Indeed in such extreme conditions, modelling methods and shielding design must be based on a detailed mechanistic knowledge of the peculiar astronauts irradiation conditions. Great help can derive from mechanistic modelling, generally aimed to better understand the intermediate steps leading from the initial energy depositions to different biological endpoints, up to organ and organism level. In the present work, criteria will be illustrated for using mechanistic approaches in developing practical tools for astronauts radioprotection, once the external field and the interaction cross sections with the spacecraft walls are known; particular attention will be given to the treatment of mixed fields. Techniques for integrating into condensed-history codes stochastic information provided by event-by-event simulations will be presented.

Aerospace Medicine↗

Stochastic aspects and uncertainties in the prechemical and chemical stages of electron tracks in liquid water: a quantitative analysis based on Monte Carlo simulations.

A new physical module for the biophysical simulation code PARTRAC has recently been developed, based on newly derived electron inelastic-scattering cross-sections in liquid water. In the present work, two modules of PARTRAC describing the production, diffusion and interaction of chemical species were developed with the specific purpose of quantifying the role of the uncertainties in the parameters controlling the early stages of liquid water radiolysis. A set of values for such parameters was identified, and time-dependent yields and frequency distributions of chemical species produced by electrons of different energies were calculated. The calculated yields were in good agreement with available data and simulations, thus confirming the reliability of the code. As the primary-electron energy decreases down to 1 keV, the *OH decay kinetics were found to get faster, reflecting variations in the spatial distribution of the initial energy depositions. In agreement with analogous works, an opposite trend was found for energies of a few hundred eV, due to the very small number of species involved. The spreading effects shown at long times by *OH frequency distributions following 1 keV irradiation were found to be essentially due to stochastic aspects of the chemical stage, whereas for 1 MeV tracks the physical and pre-chemical stages also were found to play a significant role. Relevant differences in the calculated e(aq) -yields were found by coupling the physics of PARTRAC with descriptions of the pre-chemical and chemical stages adopted in different models. This indicates a strict interrelation of the various stages, and thus a strong dependence of the parameter values on the assumptions made for the preceding and subsequent stages of the process. Although equally acceptable results can be obtained starting from different assumptions, it is necessary to keep control of such uncertainties, since they can significantly influence the modeling of radical attack on DNA and, more generally, radiobiological damage estimation. This study confirms the need for new, independently derived data on specific steps of water radiolysis, to be included in comprehensive biophysical simulation codes.

Algorithms↗

Modelling radiation-induced biological lesions: from initial energy depositions to chromosome aberrations.

The development of biophysical models of chromosome aberration induction has undergone considerable improvements in the past few years. This is mainly due to the development of new experimental techniques, such as fluorescence in situ hybridization (FISH) and premature chromosome condensation (PCC), and to a better knowledge of track structure characteristics (both in the physical and chemical stages). In particular, track structure simulations, providing a detailed description of the spatial distribution of energy depositions and relevant DNA lesions, represent a useful starting point for the development of 'ab initio' models. Various aspects of the processes determining the induction and the formation kinetics of chromosome aberrations are still under debate, concerning in particular the target description (interphase chromosome organization), the characterization of relevant DNA lesions, the possibility of inducing exchanges starting from single radiation-induced lesions, the rejoining mechanisms (proximity effects and possible induction of incomplete exchanges, i.e. one-way exchanges) and the influence of specific scoring criteria adopted both in experiments and models. Starting from Lea's breakage-and-reunion theory and Revell's exchange theory, an overview is given of various models recently reported in the literature. The assumptions adopted by the authors concerning the various processes involved in aberration formation are analysed in detail, in order to clarify the different approaches adopted in treating the open questions outlined above.

Algorithms↗

Chromosome aberrations induced by light ions: Monte Carlo simulations based on a mechanistic model.

PURPOSE: To investigate the mechanisms underlying the induction of chromosome aberrations by ionizing radiation, focusing attention on DNA damage severity, interphase chromosome geometry and the distribution of DNA strand breaks. METHODS: An ab initio biophysical model of aberration induction in human lymphocytes specific for light ions was developed, based on the assumption that 'complex lesions' (clustered DNA breaks) produce aberrations, whereas less severe breaks are repaired. It was assumed that interphase chromosomes are spatially localized and that chromosome break free-ends rejoin pairwise randomly; the unrejoining of a certain fraction of free-ends was assumed to be possible, and small fragments were neglected in order to reproduce experimental conditions. The yield of different aberrations was calculated and compared with some data obtained using Giemsa or FISH techniques. RESULTS: Dose-response curves for dicentrics and centric rings (Giemsa) and for reciprocal, complex and incomplete exchanges (FISH) were simulated; the ratio between complex and reciprocal exchanges was also calculated as a function of particle type and LET. The results showed agreement with data from lymphocyte irradiation with light ions. CONCLUSIONS: The results suggest that clustered DNA breaks are a critical damage type for aberration induction and that interphase chromosome localization plays an important role. Moreover, the effect of a given particle type is related both to the number of induced complex lesions and to their spatial distribution.

Alpha Particles↗