Background Targeted radionuclide therapy (TRT) is gaining importance. For TRT to be also used as adjuvant therapy or for treating minimal residual disease, there is a need to increase the radiation dose to small tumours. The aim of this in silico study was to compare the performances of 161Tb (a medium-energy β− emitter with additional Auger and conversion electron emissions) and 177Lu for irradiating single tumour cells and micrometastases, with various distributions of the radionuclide. Methods We used the Monte Carlo track-structure (MCTS) code CELLDOSE to compute the radiation doses delivered by 161Tb and 177Lu to single cells (14 μm cell diameter with 10 μm nucleus diameter) and to a tumour cluster consisting of a central cell surrounded by two layers of cells (18 neighbours). We focused the analysis on the absorbed dose to the nucleus of the single tumoral cell and to the nuclei of the cells in the cluster. For both radionuclides, the simulations were run assuming that 1 MeV was released per μm3 (1436 MeV/cell). We considered various distributions of the radionuclides: either at the cell surface, intracytoplasmic or intranuclear. Results For the single cell, the dose to the nucleus was substantially higher with 161Tb compared to 177Lu, regardless of the radionuclide distribution: 5.0 Gy vs. 1.9 Gy in the case of cell surface distribution; 8.3 Gy vs. 3.0 Gy for intracytoplasmic distribution; and 38.6 Gy vs. 10.7 Gy for intranuclear location. With the addition of the neighbouring cells, the radiation doses increased, but remained consistently higher for 161Tb compared to 177Lu. For example, the dose to the nucleus of the central cell of the cluster was 15.1 Gy for 161Tb and 7.2 Gy for 177Lu in the case of cell surface distribution of the radionuclide, 17.9 Gy for 161Tb and 8.3 Gy for 177Lu for intracytoplasmic distribution and 47.8 Gy for 161Tb and 15.7 Gy for 177Lu in the case of intranuclear location. Conclusion 161Tb should be a better candidate than 177Lu for irradiating single tumour cells and micrometastases, regardless of the radionuclide distribution.

Whether it is in radiobiology to identify DNA lesions or in medicine to adapt the radiotherapeutic protocols, a detailed understanding of the radiation-induced interactions in living matter is required. Monte Carlo track-structure codes have been successfully developed to describe these interactions and predict the radiation-induced energy deposits at the nanoscale level in the medium of interest. In this work, the quantum-mechanically based Monte Carlo track-structure code TILDA-V has been used to compute the slowing-down of protons in water and DNA. Stopping power and range are then reported and compared with existing data. Then, a first application of TILDA-V to cellular irradiations is also reported in order to highlight the absolute necessity of taking into account a realistic description of the cellular environment in microdosimetry.