Evidence of extranuclear cell sensitivity to alpha-particle radiation using a microdosimetric model. I. Presentation and validation of a microdosimetric model

Electron tracks simulation in water: Performance comparison between GPU CPU and the EUMED grid installation

PURPOSE

We explored different technologies to minimize simulation time of the Monte-Carlo method for track generation following the Geant4-DNA processes for electrons in water.

METHODS

A GPU software tool is developed for electron track simulations. A similar CPU version is also developed using the same collision models. CPU simulations were carried out on a single user desktop computer and on the computing grid France Grilles using 10 and 100 computing nodes. Computing time results for CPU, GPU, and grid simulations are compared with those using Geant4-DNA processes.

RESULTS

The CPU simulations better performs when the number of electrons is less than 10⁴ with 100 eV initial energy, this number decreases as the energy increases. The GPU simulations gives better results when the number of electrons is more than 10⁴ with initial energy of 100 eV, this number decreases to 10³ for electrons with 10KeV and increases back with higher energy. The use of the grid introduces an additional queuing time which slows down the overall simulation performance. Thus, the Grid gives better performance when the number of electrons is over 10⁵ with initial energy of 10KeV, and this number decreases as the energy increases.

CONCLUSIONS

The CPU is best suited for small numbers of primary incident electrons. The GPU is best suited when the number of primary incident particles occupies sufficient resources on GPU card in order to get an important computing power. The grid is best suited for simulations with high number of primary incident electrons with high initial energy.

Small-scale dosimetry for alpha particle 241Am source cell irradiation and estimation of γ-H2AX foci distribution in prostate cancer cell line PC3

Background

The development of new targeted alpha therapies motivates improving alpha particle dosimetry. For alpha particles, microscopic targets must be considered to estimate dosimetric quantities that can predict the biological response. As double-strand breaks (DSB) on DNA are the main cause of cell death by ionizing radiation, cell nuclei are relevant volumes necessary to consider as targets. Since a large variance is expected of alpha particle hits in individual cell nuclei irradiated by an uncollimated alpha-emitting source, the damage induced should have a similar distribution. The induction of DSB can be measured by immunofluorescent γ-H2AX staining. The cell γ-H2AX foci distribution and alpha particle hits distribution should be comparable and thereby verify the necessity to consider the relevant dosimetric volumes.

Methods

A Monte Carlo simulation model of an ²⁴¹Am source alpha particle irradiation setup was combined with two versions of realistic cell nuclei phantoms. These were generated from DAPI-stained PC3 cells imaged with fluorescent microscopy, one consisting of elliptical cylinders and the other of segmented mesh volumes. PC3 cells were irradiated with the ²⁴¹Am source for 4, 8 and 12 min, and after 30 min fixated and stained with immunofluorescent γ-H2AX marker. The detected radiation-induced foci (RIF) were compared to simulated RIF.

Results

The mesh volume phantom detected a higher mean of alpha particle hits and energy imparted (MeV) per cell nuclei than the elliptical cylinder phantom, but the mean specific energy (Gy) was very similar. The mesh volume phantom detected a slightly larger variance between individual cells, stemming from the more extreme and less continuous distribution of cell nuclei sizes represented in this phantom. The simulated RIF distribution from both phantoms was in good agreement with the detected RIF, although the detected distribution had a zero-inflated shape not seen in the simulated distributions. An estimate of undetected foci was used to correct the detected RIF distribution and improved the agreement with the simulations.

Conclusion

Two methods to generate cell nuclei phantoms for Monte Carlo dosimetry simulations were tested and generated similar results. The simulated and detected RIF distributions from alpha particle-irradiated PC3 cells were in good agreement, proposing the necessity to consider microscopic targets in alpha particle dosimetry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-022-00475-x.

Revisiting the Radiobiology of Targeted Alpha Therapy

Targeted alpha therapy (TAT) using alpha particle-emitting radionuclides is in the spotlight after the approval of ²²³RaCl₂ for patients with metastatic castration-resistant prostate cancer and the development of several alpha emitter-based radiopharmaceuticals. It is acknowledged that alpha particles are highly cytotoxic because they produce complex DNA lesions. Hence, the nucleus is considered their critical target, and many studies did not report any effect in other subcellular compartments. Moreover, their physical features, including their range in tissues (<100 μm) and their linear energy transfer (50–230 keV/μm), are well-characterized. Theoretically, TAT is indicated for very small-volume, disseminated tumors (e.g., micrometastases, circulating tumor cells). Moreover, due to their high cytotoxicity, alpha particles should be preferred to beta particles and X-rays to overcome radiation resistance. However, clinical studies showed that TAT might be efficient also in quite large tumors, and biological effects have been observed also away from irradiated cells. These distant effects are called bystander effects when occurring at short distance (<1 mm), and systemic effects when occurring at much longer distance. Systemic effects implicate the immune system. These findings showed that cells can die without receiving any radiation dose, and that a more complex and integrated view of radiobiology is required. This includes the notion that the direct, bystander and systemic responses cannot be dissociated because DNA damage is intimately linked to bystander effects and immune response. Here, we provide a brief overview of the paradigms that need to be revisited.

Internal microdosimetry of alpha-emitting radionuclides

At the tissue level, energy deposition in cells is determined by the microdistribution of alpha-emitting radionuclides in relation to sensitive target cells. Furthermore, the highly localized energy deposition of alpha particle tracks and the limited range of alpha particles in tissue produce a highly inhomogeneous energy deposition in traversed cell nuclei. Thus, energy deposition in cell nuclei in a given tissue is characterized by the probability of alpha particle hits and, in the case of a hit, by the energy deposited there. In classical microdosimetry, the randomness of energy deposition in cellular sites is described by a stochastic quantity, the specific energy, which approximates the macroscopic dose for a sufficiently large number of energy deposition events. Typical examples of the alpha-emitting radionuclides in internal microdosimetry are radon progeny and plutonium in the lungs, plutonium and americium in bones, and radium in targeted radionuclide therapy. Several microdosimetric approaches have been proposed to relate specific energy distributions to radiobiological effects, such as hit-related concepts, LET and track length-based models, effect-specific interpretations of specific energy distributions, such as the dual radiation action theory or the hit-size effectiveness function, and finally track structure models. Since microdosimetry characterizes only the initial step of energy deposition, microdosimetric concepts are most successful in exposure situations where biological effects are dominated by energy deposition, but not by subsequently operating biological mechanisms. Indeed, the simulation of the combined action of physical and biological factors may eventually require the application of track structure models at the nanometer scale.

New beam monitoring tool for radiobiology experiments at the cyclotron ARRONAX

The ARRONAX cyclotron is able to deliver alpha particles at 68 MeV. In the frame of radiological research, a new method is studied to infer in situ the deposited dose: it is based on the online measurement of the bremsstrahlung (>1 keV) produced by the interaction of the incident particle with the medium. Experiments are made using bombarded poly(methyl methacrylate) (PMMA)-equivalent water targets in order to characterise this continuous X-ray spectrum. The intensity of the bremsstrahlung spectrum allows for the beam monitoring. A simulation code of the bremsstrahlung has been built, and a good agreement is found with the experimental spectra. With this simulation, it is possible to predict the sensibility of this method: it varies with the target thickness, showing a good sensibility for thin target (<1000 µm) and saturation for thicker ones. Bremsstrahlung spectrum also shows a sensibility on the target's chemical composition.

Introduction to radiobiology of targeted radionuclide therapy

During the last decades, new radionuclide-based targeted therapies have emerged as efficient tools for cancer treatment. Targeted radionuclide therapies (TRTs) are based on a multidisciplinary approach that involves the cooperation of specialists in several research fields. Among them, radiobiologists investigate the biological effects of ionizing radiation, specifically the molecular and cellular mechanisms involved in the radiation response. Most of the knowledge about radiation effects concerns external beam radiation therapy (EBRT) and radiobiology has then strongly contributed to the development of this therapeutic approach. Similarly, radiobiology and dosimetry are also assumed to be ways for improving TRT, in particular in the therapy of solid tumors, which are radioresistant. However, extrapolation of EBRT radiobiology to TRT is not straightforward. Indeed, the specific physical characteristics of TRT (heterogeneous and mixed irradiation, protracted exposure, and low absorbed dose rate) differ from those of conventional EBRT (homogeneous irradiation, short exposure, and high absorbed dose rate), and consequently the response of irradiated tissues might be different. Therefore, specific TRT radiobiology needs to be explored. Determining dose-effect correlation is also a prerequisite for rigorous preclinical radiobiology studies because dosimetry provides the necessary referential to all TRT situations. It is required too for developing patient-tailored TRT in the clinic in order to estimate the best dose for tumor control, while protecting the healthy tissues, thereby improving therapeutic efficacy. Finally, it will allow to determine the relative contribution of targeted effects (assumed to be dose-related) and non-targeted effects (assumed to be non-dose-related) of ionizing radiation. However, conversely to EBRT where it is routinely used, dosimetry is still challenging in TRT. Therefore, it constitutes with radiobiology, one of the main challenges of TRT in the future.

The cytoplasm as a radiation target: an in silico study of microbeam cell irradiation

We performed in silico microbeam cell irradiation modelling to quantitatively investigate ionisations resulting from soft x-ray and alpha particle microbeams targeting the cytoplasm of a realistic cell model. Our results on the spatial distribution of ionisations show that as x-rays are susceptible to scatter within a cell that can lead to ionisations in the nucleus, soft x-ray microbeams may not be suitable for investigating the DNA damage response to radiation targeting the cytoplasm alone. In contrast, ionisations from an ideal alpha microbeam are tightly confined to the cytoplasm, but a realistic alpha microbeam degrades upon interaction with components upstream of the cellular target. Thus it is difficult to completely rule out a contribution from alpha particle hits to the nucleus when investigating DNA damage response to cytoplasmic irradiation. We find that although the cytoplasm targeting efficiency of an alpha microbeam is better than that of a soft x-ray microbeam (the probability of stray alphas hitting the nucleus is 0.2% compared to 3.6% for x-rays), stray alphas produce more ionisations in the nucleus and thus have greater potential for initiating damage responses therein. Our results suggest that observed biological responses to cytoplasmic irradiation include a small component that can be attributed to stray ionisations in the nucleus resulting from the stochastic nature of particle interactions that cause out-of-beam scatter. This contribution is difficult to isolate experimentally, thus demonstrating the value of the in silico approach.

A review of the use and potential of the GATE Monte Carlo simulation code for radiation therapy and dosimetry applications

In this paper, the authors' review the applicability of the open-source GATE Monte Carlo simulation platform based on the GEANT4 toolkit for radiation therapy and dosimetry applications. The many applications of GATE for state-of-the-art radiotherapy simulations are described including external beam radiotherapy, brachytherapy, intraoperative radiotherapy, hadrontherapy, molecular radiotherapy, and in vivo dose monitoring. Investigations that have been performed using GEANT4 only are also mentioned to illustrate the potential of GATE. The very practical feature of GATE making it easy to model both a treatment and an imaging acquisition within the same framework is emphasized. The computational times associated with several applications are provided to illustrate the practical feasibility of the simulations using current computing facilities.

Dosimetry for radiopharmaceutical therapy

Radiopharmaceutical therapy (RPT) involves the use of radionuclides that are either conjugated to tumor-targeting agents (e.g., nanoscale constructs, antibodies, peptides, and small molecules) or concentrated in tissue through natural physiological mechanisms that occur predominantly in neoplastic or otherwise targeted cells (e.g., Graves disease). The ability to collect pharmacokinetic data by imaging and use this to perform dosimetry calculations for treatment planning distinguishes RPT from other systemic treatment modalities such as chemotherapy, wherein imaging is not generally used. Treatment planning has not been widely adopted, in part, because early attempts to relate dosimetry to outcome were not successful. This was partially because a dosimetry methodology appropriate to risk evaluation rather than efficacy and toxicity was being applied to RPT. The weakest links in both diagnostic and therapeutic dosimetry are the accuracy of the input and the reliability of the radiobiological models used to convert dosimetric data to the relevant biologic end points. Dosimetry for RPT places a greater demand on both of these weak links. To date, most dosimetric studies have been retrospective, with a focus on tumor dose-response correlations rather than prospective treatment planning. In this regard, transarterial radioembolization also known as intra-arterial radiation therapy, which uses radiolabeled ((90)Y) microspheres of glass or resin to treat lesions in the liver holds much promise for more widespread dosimetric treatment planning. The recent interest in RPT with alpha-particle emitters has highlighted the need to adopt a dosimetry methodology that specifically accounts for the unique aspects of alpha particles. The short range of alpha-particle emitters means that in cases in which the distribution of activity is localized to specific functional components or cell types of an organ, the absorbed dose will be equally localized and dosimetric calculations on the scale of organs or even voxels (~5mm) are no longer sufficient. This limitation may be overcome by using preclinical models to implement macromodeling to micromodeling. In contrast to chemotherapy, RPT offers the possibility of evaluating radiopharmaceutical distributions, calculating tumor and normal tissue absorbed doses, and devising a treatment plan that is optimal for a specific patient or specific group of patients.

Redefining relative biological effectiveness in the context of the EQDX formalism: implications for alpha-particle emitter therapy

Alpha-particle radiopharmaceutical therapy (αRPT) is currently enjoying increasing attention as a viable alternative to chemotherapy for targeting of disseminated micrometastatic disease. In theory, αRPT can be personalized through pre-therapeutic imaging and dosimetry. However, in practice, given the particularities of α-particle emissions, a dosimetric methodology that accurately predicts the thresholds for organ toxicity has not been reported. This is in part due to the fact that the biological effects caused by α-particle radiation differ markedly from the effects caused by traditional external beam (photon or electron) radiation or β-particle emitting radiopharmaceuticals. The concept of relative biological effectiveness (RBE) is used to quantify the ratio of absorbed doses required to achieve a given biological response with alpha particles versus a reference radiation (typically a beta emitter or external beam radiation). However, as conventionally defined, the RBE varies as a function of absorbed dose and therefore a single RBE value is limited in its utility because it cannot be used to predict response over a wide range of absorbed doses. Therefore, efforts are underway to standardize bioeffect modeling for different fractionation schemes and dose rates for both nuclear medicine and external beam radiotherapy. Given the preponderant use of external beams of radiation compared to nuclear medicine in cancer therapy, the more clinically relevant quantity, the 2 Gy equieffective dose, EQD2(α/β), has recently been proposed by the ICRU. In concert with EQD2(α/β), we introduce a new, redefined RBE quantity, named RBE2(α/β), as the ratio of the two linear coefficients that characterize the α particle absorbed dose-response curve and the low-LET megavoltage photon 2 Gy fraction equieffective dose-response curve. The theoretical framework for the proposed new formalism is presented along with its application to experimental data obtained from irradiation of a breast cancer cell line. Radiobiological parameters are obtained using the linear quadratic model to fit cell survival data for MDA-MB-231 human breast cancer cells that were irradiated with either α particles or a single fraction of low-LET (137)Cs γ rays. From these, the linear coefficient for both the biologically effective dose (BED) and the EQD2(α/β) response lines were derived for fractionated irradiation. The standard RBE calculation, using the traditional single fraction reference radiation, gave RBE values that ranged from 2.4 for a surviving fraction of 0.82-6.0 for a surviving fraction of 0.02, while the dose-independent RBE2(4.6) value was 4.5 for all surviving fraction values. Furthermore, bioeffect modeling with RBE2(α/β) and EQD2(α/β) demonstrated the capacity to predict the surviving fraction of cells irradiated with acute and fractionated low-LET radiation, α particles and chronic exponentially decreasing dose rates of low-LET radiation. RBE2(α/β) is independent of absorbed dose for α-particle emitters and it provides a more logical framework for data reporting and conversion to equieffective dose than the conventional dose-dependent definition of RBE. Moreover, it provides a much needed foundation for the ongoing development of an α-particle dosimetry paradigm and will facilitate the use of tolerance dose data available from external beam radiation therapy, thereby helping to develop αRPT as a single modality as well as for combination therapies.

Anti-prostate-specific membrane antigen liposomes loaded with 225Ac for potential targeted antivascular α-particle therapy of cancer

This study evaluates targeted liposomes loaded with the α-particle generator (225)Ac to selectively kill prostate-specific membrane antigen (PSMA)-expressing cells with the aim to assess their potential for targeted antivascular radiotherapy.

METHODS

In this study, PEGylated liposomes were loaded with (225)Ac and labeled with the mouse antihuman PSMA J591 antibody or with the A10 PSMA aptamer. The targeting selectivity, extent of internalization, and killing efficacy of liposomes were evaluated on monolayers of prostate cancer cells intrinsically expressing PSMA (human LNCaP and rat Mat-Lu cells) and on monolayers of HUVEC induced to express PSMA (induced HUVEC).

RESULTS

The loading efficiency of (225)Ac into preformed liposomes ranged from 58.0% ± 4.6% to 85.6% ± 11.7% of introduced radioactivity. The conjugation reactions resulted in approximately 17 ± 2 J591 antibodies and 9 ± 2 A10 aptamers per liposome. The average size of liposomes, 107 ± 2 nm in diameter, was not affected by conjugation or loading. LNCaP cells exhibit 2:1:0.5 relative PSMA expression, compared with MatLu and induced HUVEC, respectively, based on flow cytometry detecting association of the J591 antibody. J591-labeled liposomes display higher levels of total specific binding to all cell lines than A10 aptamer-labeled liposomes. Specific cell association of targeted liposomes increases with incubation time. Cytotoxicity studies demonstrate that radiolabeled J591-labeled liposomes are most cytotoxic, with median lethal dose values, after 24 h of incubation, equal to 1.96 (5.3 × 10(-5)), 2.92 × 10(2) (7.9 × 10(-3)), and 2.33 × 10(1) Bq/mL (6.3 × 10(-4) μCi/mL) for LNCaP, Mat-Lu, and induced HUVEC, respectively, which are comparable to the values for the radiolabeled J591 antibody. For A10 aptamer-labeled liposomes, the corresponding values are 3.70 × 10(1) (1.0 × 10(-3)), 1.85 × 10(3) (5.0 × 10(-2)), and 4.07 × 10(3) Bq/mL (1.1 × 10(-1) μCi/mL), respectively.

CONCLUSION

Our studies demonstrate that anti-PSMA-targeted liposomes loaded with (225)Ac selectively bind, become internalized, and kill PSMA-expressing cells including endothelial cells induced to express PSMA. These findings-combined with the unique ability of liposomes to be easily tuned, in terms of size and surface modification, for optimizing biodistributions-suggest the potential of PSMA-targeting liposomes encapsulating α-particle emitters for selective antivascular α radiotherapy.

Comparison between internalizing anti-HER2 mAbs and non-internalizing anti-CEA mAbs in alpha-radioimmunotherapy of small volume peritoneal carcinomatosis using 212Pb

Background and Purpose

We assessed the contribution of antibody internalization in the efficacy and toxicity of intraperitoneal α-radioimmunotherapy (RIT) of small volume carcinomatosis using ²¹²Pb-labeled monoclonal antibodies (mAbs) that target HER2 (internalizing) or CEA (non-internalizing) receptors.

Materials and Methods

Athymic nude mice bearing 2–3 mm intraperitoneal tumor xenografts were intraperitoneally injected with similar activities (370, 740 and 1480 kBq; 37 MBq/mg) of ²¹²Pb-labeled 35A7 (anti-CEA), trastuzumab (anti-HER2) or PX (non-specific) mAbs, or with equivalent amounts of unlabeled mAbs, or with NaCl. Tumor volume was monitored by bioluminescence and survival was reported. Hematologic toxicity and body weight were assessed. Biodistribution of ²¹²Pb-labeled mAbs and absorbed dose-effect relationships using MIRD formalism were established.

Results

Transient hematological toxicity, as revealed by white blood cells and platelets numbering, was reported in mice treated with the highest activities of ²¹²Pb-labeled mAbs. The median survival (MS) was significantly higher in mice injected with 1.48 MBq of ²¹²Pb-35A7 (non-internalizing mAbs) (MS = 94 days) than in animals treated with the same activity of ²¹²Pb-PX mAbs or with NaCl (MS = 18 days). MS was even not reached after 130 days when follow-up was discontinued in mice treated with 1.48 MBq of ²¹²Pb-trastuzumab. The later efficacy was unexpected since final absorbed dose resulting from injection of 1.48 MBq, was higher for ²¹²Pb-35A7 (35.5 Gy) than for ²¹²Pb-trastuzumab (27.6 Gy). These results also highlight the lack of absorbed dose-effect relationship when mean absorbed dose was calculated using MIRD formalism and the requirement to perform small-scale dosimetry.

Conclusions

These data indicate that it might be an advantage of using internalizing anti-HER2 compared with non-internalizing anti-CEA ²¹²Pb-labeled mAbs in the therapy of small volume xenograft tumors. They support clinical investigations of ²¹²Pb-mAbs RIT as an adjuvant treatment after cytoreductive surgery in patients with peritoneal carcinomatosis.

Ex vivo activity quantification in micrometastases at the cellular scale using the α-camera technique

Targeted α-therapy (TAT) appears to be an ideal therapeutic technique for eliminating malignant circulating, minimal residual, or micrometastatic cells. These types of malignancies are typically infraclinical, complicating the evaluation of potential treatments. This study presents a method of ex vivo activity quantification with an α-camera device, allowing measurement of the activity taken up by tumor cells in biologic structures a few tens of microns.

METHODS

We examined micrometastases from a murine model of ovarian carcinoma after injection of a radioimmunoconjugate labeled with (211)At for TAT. At different time points, biologic samples were excised and cryosectioned. The activity level and the number of tumor cells were determined by combined information from 2 adjacent sections: one exposed to the α-camera and the other stained with hematoxylin and eosin. The time-activity curves for tumor cell clusters, comprising fewer than 10 cells, were derived for 2 different injected activities (6 and 1 MBq).

RESULTS

High uptake and good retention of the radioimmunoconjugate were observed at the surface of tumor cells. Dosimetric calculations based on the measured time-integrated activity indicated that for an injected activity of 1 MBq, isolated tumor cells received at least 12 Gy. In larger micrometastases (≤ 100 μm in diameter), the activity uptake per cell was lower, possibly because of hindered penetration of radiolabeled antibodies; however, the mean absorbed dose delivered to tumor cells was above 30 Gy, due to cross-fire irradiation.

CONCLUSION

Using the α-camera, we developed a method of ex vivo activity quantification at the cellular scale, which was further applied to characterize the behavior of a radiolabeled antibody administered in vivo against ovarian carcinoma. This study demonstrated a reliable measurement of activity. This method of activity quantification, based on experimentally measured data, is expected to improve the relevance of small-scale dosimetry studies and thus to accelerate the optimization of TAT.

Comparison of 211At-PRIT and 211At-RIT of ovarian microtumors in a nude mouse model

Abstract Purpose: Pretargeted radioimmunotherapy (PRIT) against intraperitoneal (i.p.) ovarian microtumors using avidin-conjugated monoclonal antibody MX35 (avidin-MX35) and (211)At-labeled, biotinylated, succinylated poly-l-lysine ((211)At-B-PLsuc) was compared with conventional radioimmunotherapy (RIT) using (211)At-labeled MX35 in a nude mouse model.

METHODS

Mice were inoculated i.p. with 1×10(7) NIH:OVCAR-3 cells. After 3 weeks, they received PRIT (1.0 or 1.5 MBq), RIT (0.9 MBq), or no treatment. Concurrently, 10 additional animals were sacrificed and examined to determine disease progression at the start of therapy. Treated animals were analyzed with regard to presence of tumors and ascites (tumor-free fraction; TFF), 8 weeks after therapy.

RESULTS

Tumor status at baseline was advanced: 70% of sacrificed animals exhibited ascites. The TFFs were 0.35 (PRIT 1.0 MBq), 0.45 (PRIT 1.5 MBq), and 0.45 (RIT). The 1.5-MBq PRIT group exhibited lower incidence of ascites and fewer tumors >1 mm than RIT-treated animals.

CONCLUSIONS

PRIT was as effective as RIT with regard to TFF; however, the size distribution of tumors and presence of ascites indicated that 1.5-MBq PRIT was more efficient. Despite advanced disease in many animals at the time of treatment, PRIT demonstrated good potential to treat disseminated ovarian cancer.

Microdosimetry for targeted alpha therapy of cancer

Targeted alpha therapy (TAT) has the advantage of delivering therapeutic doses to individual cancer cells while reducing the dose to normal tissues. TAT applications relate to hematologic malignancies and now extend to solid tumors. Results from several clinical trials have shown efficacy with limited toxicity. However, the dosimetry for the labeled alpha particle is challenging because of the heterogeneous antigen expression among cancer cells and the nature of short-range, high-LET alpha radiation. This paper demonstrates that it is inappropriate to investigate the therapeutic efficacy of TAT by macrodosimetry. The objective of this work is to review the microdosimetry of TAT as a function of the cell geometry, source-target configuration, cell sensitivity, and biological factors. A detailed knowledge of each of these parameters is required for accurate microdosimetric calculations.

Microdosimetry of alpha particles for simple and 3D voxelised geometries using MCNPX and Geant4 Monte Carlo codes

Microdosimetry using Monte Carlo simulation is a suitable technique to describe the stochastic nature of energy deposition by alpha particle at cellular level. Because of its short range, the energy imparted by this particle to the targets is highly non-uniform. Thus, to achieve accurate dosimetric results, the modelling of the geometry should be as realistic as possible. The objectives of the present study were to validate the use of the MCNPX and Geant4 Monte Carlo codes for microdosimetric studies using simple and three-dimensional voxelised geometry and to study their limit of validity in this last case. To that aim, the specific energy (z) deposited in the cell nucleus, the single-hit density of specific energy f(1)(z) and the mean-specific energy were calculated. Results show a good agreement when compared with the literature using simple geometry. The maximum percentage difference found is <6 %. For voxelised phantom, the study of the voxel size highlighted that the shape of the curve f(1)(z) obtained with MCNPX for <1 µm voxel size presents a significant difference with the shape of non-voxelised geometry. When using Geant4, little differences are observed whatever the voxel size is. Below 1 µm, the use of Geant4 is required. However, the calculation time is 10 times higher with Geant4 than MCNPX code in the same conditions.

A nephron-based model of the kidneys for macro-to-micro α-particle dosimetry

Targeted α-particle therapy is a promising treatment modality for cancer. Due to the short path-length of α-particles, the potential efficacy and toxicity of these agents is best evaluated by microscale dosimetry calculations instead of whole-organ, absorbed fraction-based dosimetry. Yet time-integrated activity (TIA), the necessary input for dosimetry, can still only be quantified reliably at the organ or macroscopic level. We describe a nephron- and cellular-based kidney dosimetry model for α-particle radiopharmaceutical therapy, more suited to the short range and high linear energy transfer of α-particle emitters, which takes as input kidney or cortex TIA and through a macro to micro model-based methodology assigns TIA to micro-level kidney substructures. We apply a geometrical model to provide nephron-level S-values for a range of isotopes allowing for pre-clinical and clinical applications according to the medical internal radiation dosimetry (MIRD) schema. We assume that the relationship between whole-organ TIA and TIA apportioned to microscale substructures as measured in an appropriate pre-clinical mammalian model also applies to the human. In both, the pre-clinical and the human model, microscale substructures are described as a collection of simple geometrical shapes akin to those used in the Cristy-Eckerman phantoms for normal organs. Anatomical parameters are taken from the literature for a human model, while murine parameters are measured ex vivo. The murine histological slides also provide the data for volume of occupancy of the different compartments of the nephron in the kidney: glomerulus versus proximal tubule versus distal tubule. Monte Carlo simulations are run with activity placed in the different nephron compartments for several α-particle emitters currently under investigation in radiopharmaceutical therapy. The S-values were calculated for the α-emitters and their descendants between the different nephron compartments for both the human and murine models. The renal cortex and medulla S-values were also calculated and the results compared to traditional absorbed fraction calculations. The nephron model enables a more optimal implementation of treatment and is a critical step in understanding toxicity for human translation of targeted α-particle therapy. The S-values established here will enable a MIRD-type application of α-particle dosimetry for α-emitters, i.e. measuring the TIA in the kidney (or renal cortex) will provide meaningful and accurate nephron-level dosimetry.

A bone marrow toxicity model for ²²³Ra alpha-emitter radiopharmaceutical therapy

Ra-223, an α-particle emitting bone-seeking radionuclide, has recently been used in clinical trials for osseous metastases of prostate cancer. We investigated the relationship between absorbed fraction-based red marrow dosimetry and cell level-dosimetry using a model that accounts for the expected localization of this agent relative to marrow cavity architecture. We show that cell level-based dosimetry is essential to understanding potential marrow toxicity. The GEANT4 software package was used to create simple spheres representing marrow cavities. Ra-223 was positioned on the trabecular bone surface or in the endosteal layer and simulated for decay, along with the descendants. The interior of the sphere was divided into cell-size voxels and the energy was collected in each voxel and interpreted as dose cell histograms. The average absorbed dose values and absorbed fractions were also calculated in order to compare those results with previously published values. The absorbed dose was predominantly deposited near the trabecular surface. The dose cell histogram results were used to plot the percentage of cells that received a potentially toxic absorbed dose (2 or 4 Gy) as a function of the average absorbed dose over the marrow cavity. The results show (1) a heterogeneous distribution of cellular absorbed dose, strongly dependent on the position of the cell within the marrow cavity; and (2) that increasing the average marrow cavity absorbed dose, or equivalently, increasing the administered activity resulted in only a small increase in potential marrow toxicity (i.e. the number of cells receiving more than 4 or 2 Gy), for a range of average marrow cavity absorbed doses from 1 to 20 Gy. The results from the trabecular model differ markedly from a standard absorbed fraction method while presenting comparable average dose values. These suggest that increasing the amount of radioactivity may not substantially increase the risk of toxicity, a result unavailable to the absorbed fraction method of dose calculation.

Evidence of extranuclear cell sensitivity to alpha-particle radiation using a microdosimetric model. II. Application of the microdosimetric model to experimental results

A microdosimetric model was used to analyze the results of experimental studies on cells of two lymphoid cell lines (T2 and Ada) irradiated with (213)Bi-radiolabeled antibodies. These antibodies targeted MHC/peptide complexes. The density of target antigen could be modulated by varying the concentration of the peptide loaded onto the cells. This offered the possibility of changing the ratio of specific (from cell-bound antibody) to non-specific (from antibody present in the supernatant) irradiation. For both cell lines, survival plotted as a function of the mean absorbed dose was a decreasing exponential. For the T2 cells, the microdosimetric sensitivity calculated for the whole cell was equal whether the irradiation was non-specific (z(0) = 0.12 +/- 0.02 Gy) or specific (z(0) = 0.12 +/- 0.09 Gy). Similar results were obtained for Ada cells. These results constitute a biological validation of the microdosimetric model. For both cells, the measured cell mortality was greater than the percentage of hit cells calculated with the model at low mean absorbed doses. This observation thus suggests bystander effects. It poses the question of the relevance of the mean absorbed dose to the cell nuclei. A new concept in cellular dosimetry taking into account cytoplasm or membrane irradiation and bystander modeling appears to be needed.