Tumour control probability derived from dose distribution in homogeneous and heterogeneous models: assuming similar pharmacokinetics, (125)Sn-(177)Lu is superior to (90)Y-(177)Lu in peptide receptor radiotherapy

A Review on Tumor Control Probability (TCP) and Preclinical Dosimetry in Targeted Radionuclide Therapy (TRT)

Targeted radionuclide therapy (TRT) uses radiopharmaceuticals to specifically irradiate tumor cells while sparing healthy tissue. Response to this treatment highly depends on the absorbed dose. Tumor control probability (TCP) models aim to predict the tumor response based on the absorbed dose by taking into account the different characteristics of TRT. For instance, TRT employs radiation with a high linear energy transfer (LET), which results in an increased effectiveness. Furthermore, a heterogeneous radiopharmaceutical distribution could result in a heterogeneous dose distribution at a tissue, cellular as well as subcellular level, which will generally reduce the tumor response. Finally, the dose rate in TRT is protracted, relatively low, and variable over time. This allows cells to repair more DNA damage, which may reduce the effectiveness of TRT. Within this review, an overview is given on how these characteristics can be included in TCP models, while some experimental findings are also discussed. Many parameters in TCP models are preclinically determined and TCP models also play a role in the preclinical stage of radiopharmaceutical development; however, this all depends critically on the calculated absorbed dose. Accordingly, an overview of the existing preclinical dosimetry methods is given, together with their limitation and applications. It can be concluded that although the theoretical extension of TCP models from external beam radiotherapy towards TRT has been established quite well, the experimental confirmation is lacking. Thus, requiring additional comprehensive studies at the sub-cellular, cellular, and organ level, which should be provided with accurate preclinical dosimetry.

Optimized cocktail of 90Y/177Lu for radionuclide therapy of neuroendocrine tumors of various sizes: a simulation study

BACKGROUND AND OBJECTIVES

There is significant interest and potential in the treatment of neuroendocrine tumors via peptide receptor radionuclide therapy (PRRT) using one or both of 90Y and 177Lu-labeled peptides. Given the presence of different tumor sizes in patients and differing radionuclide dose delivery properties, the present study aims to use Monte Carlo simulations to estimate S-values to spherical tumors of various sizes with 90Y and 177Lu separately and in combination. The goal is to determine ratios of 90Y to 177Lu that result in the largest absorbed doses per decay of the radionuclides and the most suitable dose profiles to treat tumors of specific sizes.

MATERIAL AND METHODS

Particle transfer calculations and simulations were performed using the Monte Carlo GATE simulation software. Spherical tumors of different sizes, ranging from 0.5 to 20 mm in radius, were designed. Activities of 177Lu and 90Y, individually and in combination, were homogeneously placed within the total volume of the tumors. We determined the S-values to the tumors, and to the external volume outside of the tumors (cross-dose) which was used to approximate background tissue. The dose profiles were obtained for each of the different tumor sizes, and the uniformity of dose within each tumor was calculated.

RESULTS

For all tumor sizes, the self-dose and cross-dose per decay from 90Y were higher than that from 177Lu. We observed that 177Lu had the most uniform dose distribution within tumors with radii less than 5 mm. For tumors greater than 5 mm in radius, a ratio of 25% 90Y to 75% 177Lu resulted in the most uniform doses. When the ratio of 177Lu to 90Y was smaller, the uniformity improved more with increasing tumor size. The cross-dose stayed approximately constant for tumors larger than 15 mm for all ratios of 177Lu to 90Y. Finally, as the size of the tumor increased, differences in the S-values between different ratios of 177Lu to 90Y decreased.

CONCLUSION

Our work showed that to achieve a more uniform dose distribution within the tumor, 177Lu alone is more effective for small tumors. For medium and large tumors, a ratio of 90Y to 177Lu with more or less 177Lu, respectively, is recommended.

Renal and Red Marrow Dosimetry in Peptide Receptor Radionuclide Therapy: 20 Years of History and Ahead

The development of dosimetry and studies in peptide receptor radionuclide therapy (PRRT) over the past two decades are reviewed. Differences in kidney and bone marrow toxicity reported between ⁹⁰Y, ¹⁷⁷Lu and external beam radiotherapy (EBRT) are discussed with regard to the physical properties of these beta emitter radionuclides. The impact of these properties on the response to small and large tumors is also considered. Capacities of the imaging modalities to assess the dosimetry to target tissues are evaluated. Studies published in the past two years that confirm a red marrow uptake in ¹⁷⁷Lu-DOTATATE therapy, as already observed 20 years ago in ⁸⁶Y-DOTATOC PET studies, are analyzed in light of the recent developments in the transferrin transport mechanism. The review enlightens the importance (i) of using state-of-the-art imaging modalities, (ii) of individualizing the activity to be injected with regard to the huge tissue uptake variability observed between patients, (iii) of challenging the currently used but inappropriate blood-based red marrow dosimetry and (iv) of considering individual tandem therapy. Last, a smart individually optimized tandem therapy taking benefit of the bi-orthogonal toxicity-response pattern of ¹⁷⁷Lu-DOTATATE and of ⁹⁰Y-DOTATOC is proposed.

68Ga-DOTATOC PET/CT-Based Radiomic Analysis and PRRT Outcome: A Preliminary Evaluation Based on an Exploratory Radiomic Analysis on Two Patients

Aim: This work aims to evaluate whether the radiomic features extracted by 68Ga-DOTATOC-PET/CT of two patients are associated with the response to peptide receptor radionuclide therapy (PRRT) in patients affected by neuroendocrine tumor (NET).

Methods: This is a pilot report in two NET patients who experienced a discordant response to PRRT (responder vs. non-responder) according to RECIST1.1. The patients presented with liver metastasis from the rectum and pancreas G3-NET, respectively. Whole-body total-lesion somatostatin receptor-expression (TLSREwb-50) and somatostatin receptor-expressing tumor volume (SRETV wb-50) were obtained in pre- and post-PRRT PET/CT. Radiomic analysis was performed, extracting 38 radiomic features (RFs) from the patients' lesions. The Mann–Whitney test was used to compare RFs in the responder patient vs. the non-responder patient. Pearson correlation and principal component analysis (PCA) were used to evaluate the correlation and independence of the different RFs.

Results: TLSREwb-50 and SRETVwb-50 modifications correlate with RECIST1.1 response. A total of 28 RFs extracted on pre-therapy PET/CT showed significant differences between the two patients in the Mann–Whitney test (p < 0.05). A total of seven second-order features, with poor correlation with SUVmax and PET volume, were identified by the Pearson correlation matrix. Finally, the first two PCA principal components explain 83.8% of total variance.

Conclusion: TLSREwb-50 and SRETVwb-50 are parameters that might be used to predict and to assess the PET response to PRRT. RFs might have a role in defining inter-patient heterogeneity and in the prediction of therapy response. It is important to implement future studies with larger and more homogeneous patient populations to confirm the efficacy of these biomarkers.

Targeting Micrometastases: The Effect of Heterogeneous Radionuclide Distribution on Tumor Control Probability

The spatial distribution of radiopharmaceuticals that emit short-range high linear-energy-transfer electrons greatly affects the absorbed dose and their biological effectiveness. The purpose of this study was to investigate the effect of heterogeneous radionuclide distribution on tumor control probability (TCP) in a micrometastases model. Methods: Cancer cell lines; MDA-MB-468, SQ20B and 231-H2N were grown as spheroids to represent micrometastases. The intracellular distribution of a representative radiopeptide (111In-labelled epidermal growth factor, EGF) and radioimmunotherapeutic (111In-labelled Trastuzumab) was determined in cell internalization experiments. The intratumoral distribution was evaluated by microautoradiography of spheroids. γH2AX staining was performed on spheroid sections to correlate DNA damage with radionuclide distribution. Experimental surviving fractions (SFexp ) were obtained using clonogenic assays. A random closed-packed algorithm, which models the random packing behavior of cells and reflects variation in the radii of cells and nuclei, was used to simulate 3-D spheroids. Calculated survival fractions (SFcal ) were generated using an iterative modelling method based on Monte Carlo determined absorbed dose with the PENELOPE code and were compared to (SFexp ). Radiobiological parameters deduced from experimental results and MC simulations were used to predict the TCP for a 3-D spheroid model. Results: Calculated SFs were in good agreement with experimental data, particularly when an increased value for relative biological effectiveness (RBE) was applied to self-dose deposited by sources located in the nucleus and when radiobiological parameters were adjusted to account for dose protraction. Only in MDA-MB-468 spheroids treated with 111In-EGF was a TCP>0.5 achieved, indicating that for this cell type the radiopeptide would be curative when targeting micrometastases. This is attributed to the relative radiosensitivity of MDA-MB-468 cells, high nuclear uptake of the radiopeptide and uniform distribution of radioactivity throughout the spheroid. Conclusion: It is imperative to include biological endpoints when evaluating the distribution of radionuclides in models emulating micrometastatic disease. The spatial distribution of radioactivity is a clear determinant of biological effect and TCP as demonstrated in this study.

A radiobiological model of metastatic burden reduction for molecular radiotherapy: application to patients with bone metastases

Skeletal tumour burden is a biomarker of prognosis and survival in cancer patients. This study proposes a novel method based on the linear quadratic model to predict the reduction in metastatic tumour burden as a function of the absorbed doses delivered from molecular radiotherapy treatments. The range of absorbed doses necessary to eradicate all the bone lesions and to reduce the metastatic burden was investigated in a cohort of 22 patients with bone metastases from castration-resistant prostate cancer. A metastatic burden reduction curve was generated for each patient, which predicts the reduction in metastatic burden as a function of the patient mean absorbed dose, defined as the mean of all the lesion absorbed doses in any given patient. In the patient cohort studied, the median of the patient mean absorbed dose predicted to reduce the metastatic burden by 50% was 89 Gy (interquartile range: 83-105 Gy), whilst a median of 183 Gy (interquartile range: 107-247 Gy) was found necessary to eradicate all metastases in a given patient. The absorbed dose required to eradicate all the lesions was strongly correlated with the variability of the absorbed doses delivered to multiple lesions in a given patient (r  =  0.98, P  <  0.0001). The metastatic burden reduction curves showed a potential large reduction in metastatic burden for a small increase in absorbed dose in 91% of patients. The results indicate the range of absorbed doses required to potentially obtain a significant survival benefit. The metastatic burden reduction method provides a simple tool that could be used in routine clinical practice for patient selection and to indicate the required administered activity to achieve a predicted patient mean absorbed dose and reduction in metastatic tumour burden.

A dosimetric model for the heterogeneous delivery of radioactive nanoparticles In vivo: a feasibility study

ᅟ

Accurate and quantitative dosimetry for internal radiation therapy can be especially challenging, given the heterogeneity of patient anatomy, tumor anatomy, and source deposition. Internal radiotherapy sources such as nanoparticles and monoclonal antibodies require high resolution imaging to accurately model the heterogeneous distribution of these sources in the tumor. The resolution of nuclear imaging modalities is not high enough to measure the heterogeneity of intratumoral nanoparticle deposition or intratumoral regions, and mathematical models do not represent the actual heterogeneous dose or dose response. To help answer questions at the interface of tumor dosimetry and tumor biology, we have modeled the actual 3-dimensional dose distribution of heterogeneously delivered radioactive nanoparticles in a tumor after systemic injection.

Methods

24 h after systemic injection of dually fluorescent and radioactive nanoparticles into a tumor-bearing mouse, the tumor was cut into 342 adjacent sections and imaged to quantify the source distribution in each section. The images were stacked to generate a 3D model of source distribution, and a novel MATLAB code was employed to calculate the dose to cells on a middle section in the tumor using a low step size dose kernel.

Results

The average dose calculated by this novel 3D model compared closely with standard ways of calculating average dose, and showed a positive correlation with experimentally determined cytotoxicity in vivo. The high resolution images allowed us to determine that the dose required to initiate radiation-induced H2AX phosphorylation was approximately one Gray. The nanoparticle distribution was further used to model the dose distribution of two other radionuclides.

Conclusions

The ability of this model to quantify the absorbed dose and dose response in different intratumoral regions allows one to investigate how source deposition in different tumor areas can affect dose and cytotoxicity, as well as how characteristics of the tumor microenvironment, such as hypoxia or high stromal areas, may affect the potency of a given dose.

Electronic supplementary material

The online version of this article (doi:10.1186/s13014-017-0794-z) contains supplementary material, which is available to authorized users.

Influence of tumour size on the efficacy of targeted alpha therapy with (213)Bi-[DOTA(0),Tyr(3)]-octreotate

Background

Targeted alpha therapy has been postulated to have great potential for the treatment of small clusters of tumour cells as well as small metastases. ²¹³Bismuth, an α-emitter with a half-life of 46 min, has shown to be effective in preclinical as well as in clinical applications. In this study, we evaluated whether ²¹³Bi-[DOTA⁰, Tyr³]-octreotate (²¹³Bi-DOTATATE), a ²¹³Bi-labelled somatostatin analogue with high affinity for somatostatin receptor subtype 2 (SSTR₂), is suitable for the treatment of larger neuroendocrine tumours overexpressing SSTR₂ in comparison to its effectiveness for smaller tumours. We performed a preclinical targeted radionuclide therapy study with ²¹³Bi-DOTATATE in animals bearing tumours of different sizes (50 and 200 mm³) using two tumour models: H69 (human small cell lung carcinoma) and CA20948 (rat pancreatic tumour).

Methods

Pharmacokinetics was determined for calculation of dosimetry in organs and tumours. H69- or CA20948-xenografted mice with tumour volumes of approximately 120 mm³ were euthanized at 10, 30, 60 and 120 min post injection of a single dose of ²¹³Bi-DOTATATE (1.5–4.8 MBq). To investigate the therapeutic efficacy of ²¹³Bi-DOTATATE, xenografted H69 and CA20948 tumour-bearing mice with tumour sizes of 50 and 200 mm³ were administered daily with a therapeutic dose of ²¹³Bi-DOTATATE (0.3 nmol, 2–4 MBq) for three consecutive days. The animals were followed for 90 days after treatment. At day 90, mice were injected with 25 MBq ^99mTc-DMSA and imaged by SPECT/CT to investigate possible renal dysfunction due to ²¹³Bi-DOTATATE treatment.

Results

Higher tumour uptakes were found in CA20948 tumour-bearing animals compared to those in H69 tumour-bearing mice with the highest tumour uptake of 19.6 ± 6.6 %IA/g in CA20948 tumour-bearing animals, while for H69 tumour-bearing mice, the highest tumour uptake was found to be 9.8 ± 2.4 %IA/g. Nevertheless, as the anti-tumour effect was more pronounced in H69 tumour-bearing mice, the survival rate was higher. Furthermore, in the small tumour groups, no regrowth of tumour was found in two H69 tumour-bearing mice and in one of the CA20948 tumour-bearing mice. No renal dysfunction was observed in ²¹³Bi-DOTATATE-treated mice after the doses were applied.

Conclusions

²¹³Bi-DOTATATE demonstrated a great therapeutic effect in both small and larger tumour lesions. Higher probability for stable disease was found in animals with small tumours. ²¹³Bi-DOTATATE was effective in different neuroendocrine (H69 and CA20948) tumour models with overexpression of SSTR₂ in mice.

Electronic supplementary material

The online version of this article (doi:10.1186/s13550-016-0162-2) contains supplementary material, which is available to authorized users.

Translational radionanomedicine: a clinical perspective

Many nanomaterials were developed for the anticipated in vivo theranostic use exploiting their unique characteristics as a multifunctional platform. Nevertheless, only a few nanomaterials are under investigation for human use, most of which have not entered clinical trials yet. Radionanomedicine, a convergent discipline of radiotracer technology and use of nanomaterials in vivo, can facilitate clinical nanomedicine because of its advantages of radionuclide imaging and internal radiation therapy. In this review, we focuse on how radionanomedicine would impact profoundly on clinical translation of nanomaterial theranostics. Up-to-date advances and future challenges are critically reviewed regarding the issues of how to radiolabel and engineer radionanomaterials, in vivo behavior tracing of radionanomaterials and then the desired clinical radiation dosimetry. Radiolabeled extracellular vesicles were further discussed as endogenous nanomaterials radiolabeled for possible clinical use.

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.

Planning combined treatments of external beam radiation therapy and molecular radiotherapy

Molecular radiotherapy (MRT) with radiolabeled molecules has being constantly evolving, leading to notable results in cancer treatment. In some cases, the absorbed doses delivered to tumors by MRT are sufficient to obtain complete responses; in other cases, instead, to be effective, MRT needs to be combined with other therapeutic approaches. Recently, several studies proposed the combination of MRT with external beam radiation therapy (EBRT). Some describe the theoretical basis within radiobiological models, others report the results of clinical phase I-II studies aimed to assess the feasibility and tolerability. The latter includes the treatment of various tumors, such as meningiomas, paragangliomas, non-Hodgkin's lymphomas, bone, brain, hepatic, and breast lesions. The underlying principle of combined MRT and EBRT is the possibility of exploiting the full potential of each modality, given the different organs at risk. Target tissues can indeed receive a higher irradiation, while respecting the threshold limits of more than one critical tissue. Nevertheless, clinical trials are empirical and optimization is still a theoretical issue. This article describes the state of the art of combined MRT and EBRT regarding the rationale and the results of clinical studies, with special focus on the possibility of treatment improvement.

Objective assessment of image quality VI: imaging in radiation therapy

Earlier work on objective assessment of image quality (OAIQ) focused largely on estimation or classification tasks in which the desired outcome of imaging is accurate diagnosis. This paper develops a general framework for assessing imaging quality on the basis of therapeutic outcomes rather than diagnostic performance. By analogy to receiver operating characteristic (ROC) curves and their variants as used in diagnostic OAIQ, the method proposed here utilizes the therapy operating characteristic or TOC curves, which are plots of the probability of tumor control versus the probability of normal-tissue complications as the overall dose level of a radiotherapy treatment is varied. The proposed figure of merit is the area under the TOC curve, denoted AUTOC. This paper reviews an earlier exposition of the theory of TOC and AUTOC, which was specific to the assessment of image-segmentation algorithms, and extends it to other applications of imaging in external-beam radiation treatment as well as in treatment with internal radioactive sources. For each application, a methodology for computing the TOC is presented. A key difference between ROC and TOC is that the latter can be defined for a single patient rather than a population of patients.

Radiobiologic optimization of combination radiopharmaceutical therapy applied to myeloablative treatment of non-Hodgkin lymphoma

Combination treatment is a hallmark of cancer therapy. Although the rationale for combination radiopharmaceutical therapy was described in the mid-1990s, such treatment strategies have only been implemented clinically recently and without a rigorous methodology for treatment optimization. Radiobiologic and quantitative imaging-based dosimetry tools are now available that enable rational implementation of combined targeted radiopharmaceutical therapy. Optimal implementation should simultaneously account for radiobiologic normal-organ tolerance while optimizing the ratio of 2 different radiopharmaceuticals required to maximize tumor control. We have developed such a methodology and applied it to hypothetical myeloablative treatment of non-Hodgkin lymphoma (NHL) patients using (131)I-tositumomab and (90)Y-ibritumomab tiuxetan.

METHODS

The range of potential administered activities (AAs) is limited by the normal-organ maximum-tolerated biologic effective doses (MTBEDs) arising from the combined radiopharmaceuticals. Dose-limiting normal organs are expected to be the lungs for (131)I-tositumomab and the liver for (90)Y-ibritumomab tiuxetan in myeloablative NHL treatment regimens. By plotting the limiting normal-organ constraints as a function of the AAs and calculating tumor biologic effective dose (BED) along the normal-organ MTBED limits, we obtained the optimal combination of activities. The model was tested using previously acquired patient normal-organ and tumor kinetic data and MTBED values taken from the literature.

RESULTS

The average AA value based solely on normal-organ constraints was 19.0 ± 8.2 GBq (range, 3.9-36.9 GBq) for (131)I-tositumomab and 2.77 ± 1.64 GBq (range, 0.42-7.54 GBq) for (90)Y-ibritumomab tiuxetan. Tumor BED optimization results were calculated and plotted as a function of AA for 5 different cases, established using patient normal-organ kinetics for the 2 radiopharmaceuticals. Results included AA ranges that would deliver 95% of the maximum tumor BED, allowing for informed inclusion of clinical considerations, such as a maximum-allowable (131)I administration.

CONCLUSION

A rational approach for combination radiopharmaceutical treatment has been developed within the framework of a proven 3-dimensional (3D) personalized dosimetry software, 3D-RD, and applied to the myeloablative treatment of NHL. We anticipate that combined radioisotope therapy will ultimately supplant single radioisotope therapy, much as combination chemotherapy has substantially replaced single-agent chemotherapy.

Kidney dosimetry in ¹⁷⁷Lu and ⁹⁰Y peptide receptor radionuclide therapy: influence of image timing, time-activity integration method, and risk factors

Kidney dosimetry in ¹⁷⁷Lu and ⁹⁰Y PRRT requires 3 to 6 whole-body/SPECT scans to extrapolate the peptide kinetics, and it is considered time and resource consuming. We investigated the most adequate timing for imaging and time-activity interpolating curve, as well as the performance of a simplified dosimetry, by means of just 1-2 scans. Finally the influence of risk factors and of the peptide (DOTATOC versus DOTATATE) is considered. 28 patients treated at first cycle with ¹⁷⁷Lu DOTATATE and 30 with ¹⁷⁷Lu DOTATOC underwent SPECT scans at 2 and 6 hours, 1, 2, and 3 days after the radiopharmaceutical injection. Dose was calculated with our simplified method, as well as the ones most used in the clinic, that is, trapezoids, monoexponential, and biexponential functions. The same was done skipping the 6 h and the 3 d points. We found that data should be collected until 100 h for ¹⁷⁷Lu therapy and 70 h for ⁹⁰Y therapy, otherwise the dose calculation is strongly influenced by the curve interpolating the data and should be carefully chosen. Risk factors (hypertension, diabetes) cause a rather statistically significant 20% increase in dose (t-test, P < 0.10), with DOTATATE affecting an increase of 25% compared to DOTATOC (t-test, P < 0.05).