Proliferation and the advantage of longer-lived radionuclides in radioimmunotherapy

Generalized methods for predicting biological response to mixed radiation types and calculating equieffective doses (EQDX)

BACKGROUND

Predicting biological responses to mixed radiation types is of considerable importance when combining radiation therapies that use multiple radiation types and delivery regimens. These may include the use of both low- and high-linear energy transfer (LET) radiations. A number of theoretical models have been developed to address this issue. However, model predictions do not consistently match published experimental data for mixed radiation exposures. Furthermore, the models are often computationally intensive. Accordingly, there is a need for efficient analytical models that can predict responses to mixtures of low- and high-LET radiations. Additionally, a general formalism to calculate equieffective dose (EQDX) for mixed radiations is needed.

PURPOSE

To develop a computationally efficient analytical model that can predict responses to complex mixtures of low- and high-LET radiations as a function of either absorbed dose or EQDX.

METHODS

The Zaider-Rossi model (ZRM) was modified by replacing the geometric mean of the quadratic coefficients in the interaction term with the arithmetic mean. This modified ZRM model (mZRM) was then further generalized to any number of radiation types and its validity was tested against published experimental observations. Comparisons between the predictions of the ZRM and mZRM, and other models, were made using two and three radiation types. In addition, a generalized formalism for calculating EQDX for mixed radiations was developed within the context of mZRM and validated with published experimental results.

RESULTS

The predictions of biological responses to mixed-LET radiations calculated with the mZRM are in better agreement with experimental observations than ZRM, especially when high- and low-LET radiations are mixed. In these situations, the ZRM overestimated the surviving fraction. Furthermore, the EQDX calculated with mZRM are in better agreement with experimental observations.

CONCLUSION

The mZRM is a computationally efficient model that can be used to predict biological response to mixed radiations that have low- and high-LET characteristics. Importantly, interaction terms are retained in the calculation of EQDX for mixed radiation exposures within the mZRM framework. The mZRM has application in a wide range of radiation therapies, including radiopharmaceutical therapy.

Dosimetry of bone seeking beta emitters for bone pain palliation metastases

Amongst cancer patients, bone pain due to skeletal metastases is a major cause of morbidity. A number of beta-emitting radiopharmaceuticals have been used to provide internal radiotherapy of bone metastases and provide palliative pain relief. In this article we describe the different physical characteristics of the various beta emitting radionuclides which have been used in this clinical setting and the potential impact of differences in dose-rate on radiobiological outcomes. A detailed review of the biodistribution of these treatments, based on both in-vivo clinical investigations and post mortem autoradiography assessments is provided. These treatments result in physiological delivery of radiation doses to the target disease as well as to critical healthy organs. Particular attention is paid to the radiation doses received by normal bone tissue, bone marrow as well as metastatic bone disease. The underlying models of radiation transport within bone and bone marrow are reviewed alongside the practical steps that must be taken to acquire and analyse the information require for clinical dosimetry assessments. The role of whole body measurements, blood and faecal assays as well as both planar and tomographic gamma camera imaging are considered. In addition we review the rationale for allocating measured bone uptake between trabecular and cortical bone tissue. The difference between bone volume and bone surface seeking radiopharmaceuticals are also discussed. This review also extends to the development of preclinical models of bone metastases which may inform future dosimetric calculations. Finally, we also present a comprehensive review of the dosimetry of the established treatments ⁸⁹Strontium-chloride; ³²Phosphorus; ¹⁸⁸Rhenium-hydroxyethylidine disphosphonate; ¹⁸⁶Rhenium-1,1-hydroxyethylidene disphosphonate (¹⁸⁶Re-HEDP); ¹⁵³Samarium-ethylenediaminetetramethylene phosphonate; as well as the emerging treatments ¹⁸⁸Rhenium-zoledronic acid; ¹⁸⁸Rhenium-ibedronat; ¹⁷⁷Lutetium-zoledronic acid; and ¹⁷⁷Lutetium ethylenediaminetetramethylene phosphonate. This review highlights not only the inter treatment differences in the radiation absorbed doses delivered to metastatic disease by different radiopharmaceuticals but also the intra treatment differences which result in a large range of observed doses between patients.

Population exposure-response model of 131I in patients with benign thyroid disease

PURPOSE

The study aimed to explore the relationship of different exposure measures with ¹³¹I therapy response in patients with benign thyroid disease, estimate the variability in the response, investigate possible covariates, and discuss dosing implications of the results.

METHODS

A population exposure-response analysis was performed using nonlinear mixed-effects modelling. Data from 95 adult patients with benign thyroid disease were analysed. Evaluated exposure parameters were: administered radioactivity dose (Aa) [MBq], total absorbed dose (ABD) [Gy], maximum of absorbed dose-rate (MXR) [Gy/h] and biologically effective dose (BED) [Gy]. The response was modelled as ordered categorical data: hyper-, eu- and hypothyroidism. The final model performance was evaluated by a visual predictive check.

RESULTS

The probability of the outcome following ¹³¹I therapy was best described by a proportional-odds model, including the log-linear model of ¹³¹I effect and the exponential model of the response-time relationship. All exposure measures were statistically significant with p<0.001, with BED and ABD being statistically better than the other two. Nevertheless, as BED resulted in the lowest AIC value, it was included in the final model. Accordingly, BED value of 289.7 Gy is associated with 80% probability of successful treatment outcome 12 months after ¹³¹I application in patients with median thyroid volume (32.28 mL). The target thyroid volume was a statistically significant covariate. The visual predictive check of the final model showed good model performance.

CONCLUSION

Our results imply that BED formalism could aid in therapy individualisation. The larger thyroid volume is associated with a lower probability of a successful outcome.

Monte Carlo single-cell dosimetry using Geant4-DNA: the effects of cell nucleus displacement and rotation on cellular S values

Investigation of biological effects of low-dose ionizing radiation at the (sub-) cellular level, which is referred to as microdosimetry, remains a major challenge of today's radiobiology research. Monte Carlo simulation of radiation tracks can provide a detailed description of the physical processes involved in dimensions as small as the critical substructures of the cell. Hereby, in the present study, microdosimetric calculations of cellular S values for mono-energetic electrons and six Auger-emitting radionuclides were performed in single-cell models of liquid water using Geant4-DNA. The effects of displacement and rotation of the nucleus within the cell on the cellular S values were studied in spherical and ellipsoidal geometries. It was found that for the examined electron energies and radionuclides, in the case of nucleus cross-absorption where the radioactivity is either localized in the cytoplasm of the cell or distributed on the cell surface, rotation of the nucleus within the cell affects cellular S values less than displacement of the nucleus. Especially, the considerable differences observed in S(nucleus ← cell surface) values between an eccentric and a concentric cell-nucleus configuration in spherical and ellipsoidal geometries (up to 63% and up to 44%, respectively) suggests that the approximation of concentricity should be used with caution, at least for localized irradiation of the cell membrane by an Auger-emitter in targeted radionuclide cancer therapy. The obtained results, which are based on a more realistic modeling of the cell than was done before, provide more accurate information about nuclear dose. This can be useful for theranostic applications.

Design and testing of a microcontroller that enables alpha particle irradiators to deliver complex dose rate patterns

There is increasing interest in using alpha particle emitting radionuclides for cancer therapy because of their unique cytotoxic properties which are advantageous for eradicating tumor cells. The high linear energy transfer (LET) of alpha particles produces a correspondingly high density of ionizations along their track. Alpha particle emitting radiopharmaceuticals deposit this energy in tissues over prolonged periods with complex dose rate patterns that depend on the physical half-life of the radionuclide, and the biological uptake and clearance half-times in tumor and normal tissues. We have previously shown that the dose rate increase half-time that arises as a consequence of these biokinetics can have a profound effect on the radiotoxicity of low-LET radiation. The microcontroller hardware and software described here offer a unique way to deliver these complex dose rate patterns with a broad-beam alpha particle irradiator, thereby enabling experiments to study the radiobiology of complex dose rate patterns of alpha particles. Complex dose rate patterns were created by precise manipulation of the timing of opening and closing of the electromechanical shutters of an α-particle irradiator. An Arduino Uno and custom circuitry was implemented to control the shutters. The software that controls the circuits and shutters has a user-friendly Graphic User Interface (GUI). Alpha particle detectors were used to validate the programmed dose rate profiles. Circuit diagrams and downloadable software are provided to facilitate adoption of this technology by other radiobiology laboratories.

Cellular Response to Exponentially Increasing and Decreasing Dose Rates: Implications for Treatment Planning in Targeted Radionuclide Therapy

The treatment of cancer using targeted radionuclide therapy is of interest to nuclear medicine and radiation oncology because of its potential for killing tumor cells while minimizing dose-limiting toxicities to normal tissue. The ionizing radiations emitted by radiopharmaceuticals deliver radiation absorbed doses over protracted periods of time with continuously varying dose rates. As targeted radionuclide therapy becomes a more prominent part of cancer therapy, accurate models for estimating the biologically effective dose (BED) or equieffective dose (EQD2α/β) will become essential for treatment planning. This study examines the radiobiological impact of the dose rate increase half-time during the uptake phase of the radiopharmaceutical. MDA-MB-231 human breast cancer cells and V79 Chinese hamster lung fibroblasts were irradiated chronically with 662 keV γ rays delivered with time-varying dose rates that are clinically relevant. The temporal dose-rate patterns were: 1. acute, 2. exponential decrease with a half-time of 64 h (Td = 64 h), 3. initial exponential increase to a maximum (half time Ti = 2, 8 or 24 h) followed by exponential decrease (Td = 64 h). Cell survival assays were conducted and surviving fractions were determined. There was a marked reduction in biological effect when Ti was increased. Cell survival data were tested against existing dose-response models to assess their capacity to predict response. Currently accepted models that are used in radiation oncology overestimated BED and EQD2α/β at low-dose rates and underestimated them at high-dose rates. This appears to be caused by an adaptive response arising as a consequence of the initial low-dose-rate phase of exposure. An adaptive response function was derived that yields more accurate BED and EQD2α/β values over the spectrum of dose rates and absorbed doses delivered. Our experimental data demonstrate a marked increase in cell survival when the dose-rate-increase half-time is increased, thereby suggesting an adaptive response arising as a consequence of this phase of exposure. We have modified conventional radiobiological models used in the clinic for brachytherapy and external beams of radiation to account for this phenomenon and facilitate their use for treatment planning in targeted radionuclide therapy.

Neutron activation of In(iii) complexes with thiosemicarbazones leads to the production of potential radiopharmaceuticals for the treatment of breast cancer

^114mIn(iii) complexes with 2-acetylpyridine-derived thiosemicarbazones show potent cytotoxic activity.

The role of nuclear medicine in modern therapy of cancer

Nuclear medicine is a multidisciplinary field that develops and uses instrumentation and tracers (radiopharmaceuticals) to study physiological processes and noninvasively diagnose, stage, and treat diseases. Particularly, it offers a unique means to study cancer biology in vivo and to optimize cancer therapy for individual patients. A tracer is either a radionuclide alone, such as iodine-131 or a radiolabel in a carrier molecule such as (18)F in fluorodeoxyglucose ((18)F-FDG), or other feasible radionuclide attached to a drug, a protein, or a peptide, which when introduced into the body, would accumulate in the tissue of interest. Nuclear medicine imaging, including single-photon emission computer tomography and positron emission tomography, can provide important quantitative and functional information about normal tissues or disease conditions, in contrast to conventional, anatomical imaging techniques such as ultrasound, computed tomography, or magnetic resonance imaging. For treatment, tumor-targeting agents, conjugated with therapeutic radionuclides, may be used to deposit lethal radiation at tumor sites. This review outlines the role of nuclear medicine in modern cancer therapy.

Effect of different cell cluster models on the radiobiological output for (211)At-radioimmunotherapy

The cell cluster modeling is a widely used method to estimate the small-scale dosimetry and provides the implication for a clinic. This work evaluated the effect of different regular cluster models on the radiobiological outputs for (211)At-radioimmunotherapy. The cell activity threshold was estimated using a tumor control probability of 0.90. Basically, regular models show similar features with cluster configuration and cell dimension variation. However, their individual results such as the cumulated activity threshold per cell and the prescription dose per volume should not be substituted reciprocally. The tissue composed of smaller cells or midcell packing will need a little more high prescription dose per volume. The radiation sensitivity parameters in a linear-quadratic model are critical to decide the radiobiological response with dose. The cumulated cell activity threshold increases exponentially with α decreasing, and its influence on the big cell dimension is more than on the small one. The different subsources affect radioresistant organs or tissues more remarkably than radiosensitive ones, especially the cells with large cytoplasm. The heterogeneous activity of Gaussian distribution will decrease the therapeutical effectiveness for the nucleus source, but its influence on the cytoplasm and cell surface sources is a little uncertain, as their real mean value is always higher than its set mean value by assuming the cell activity uptakes from zero. Careful usage of underdose with heterogeneous activity distribution should be practiced in clinics. The deteriorated heterogeneous distribution will salvage the potential subversive and lead to the failure of tumor local control. Some cells with no or little activity that are located on the edge or vertex of cube or corner models will have the ability to survive, as there is a lack of a part of the cross-fire dose effect, and so more attention should be paid in selecting the dosage. Although this work focuses on the clinic implication of (211)At in α-radioimmunotherapy, these cell cluster models can be generalized to other radionuclides.

A radiobiological investigation on dose and dose rate for permanent implant brachytherapy of breast using 125I or 103Pd sources

PURPOSE

The present report addresses the question of what could be the appropriate dose and dose rate for 125I and 103PD permanent seed implants for breast cancer as monotherapy for early stage breast cancer. This is addressed by employing a radiobiological methodology, which is based on the linear quadratic model, to identify a biologically effective dose (BED) to the prescription point of the brachytherapy implant, which would produce equivalent cell killing (or same cell survival) when compared to a specified external radiotherapy scheme.

METHODS

In the present analysis, the tumor and normal tissue BED ratios of brachytherapy and external radiotherapy are examined for different combinations of tumor proliferation constant (K), alpha/beta ratios, initial dose rate (R0), and reference external radiotherapy scheme (50 or 60 Gy in 2 Gy per fraction). The results of the radiobiological analysis are compared against other reports and clinical protocols in order to examine possible opportunities of improvement.

RESULTS

The analysis indicates that physical doses of approximately 100-110 Gy delivered with an initial dose rate of around 0.05 Gyh(-1) and 78-80 Gy delivered at 0.135 Gyh(-1) for 125I and 103Pd permanent implants, respectively, are equivalent to 50 Gy external beam radiotherapy (EBRT) in 2 Gy per fraction. Similarly, for physical doses of approximately 115-127 Gy delivered with an initia dose rate of around 0.059 Gyh(-1) and 92 Gy delivered at 0.157 Gyh(-1) for 125I and 103Pd, respectively, are equivalent to 60 Gy EBRT in 2 Gy per fraction. It is shown that the initial dose rate required to produce isoeffective tumor response with 50 or 60 Gy EBRT in 2 Gy per fraction increases as the repopulation factor K increases, even though repopulation is also considered in EBRT. Also, the initial dose rate increases as the value of the alpha/beta ratio decreases. The impact of the different alpha/beta ratios on the ratio of the tumor BEDs is significantly large for both the 125I and 103Pd implants with the deviation between the alpha/beta = 10.0 Gy ratios and those using the 4.0 and 3.5 Gy values ranging between 18% and 22% in most of the cases.

CONCLUSIONS

For the cases of 125I and 103Pd, the equivalent physical doses to 50 Gy EBRT in 2 Gy per fraction are associated with an overdosage of the involved normal tissue in the range of 4%-16% and an underdosage by 10%-15% for a BED for normal tissue, using an alpha/beta value of 3.0 Gy (BEDNT,3 Gy) of 100 Gy. These values are lower by 10%-20% than the published value of 124 Gy for 125I and by about 13% when compared to the published isoeffective dose of 90 Gy for 103Pd. Similarly, the equivalent physical doses to 60 Gy EBRT in 2 Gy per fraction are associated with an overdosage of the involved normal tissue by 10%-20% and an underdosage by 4%-10% for BEDNT,3 Gy of 110 Gy.

Extension of the biological effective dose to the MIRD schema and possible implications in radionuclide therapy dosimetry

In dosimetry-based treatment planning protocols, patients with rapid clearance of the radiopharmaceutical require a larger amount of initial activity than those with slow clearance to match the absorbed dose to the critical organ. As a result, the dose-rate to the critical organ is higher in patients with rapid clearance and may cause unexpected toxicity compared to patients with slow clearance. In order to account for the biological impact of different dose-rates, radiobiological modeling is beginning to be applied to the analysis of radionuclide therapy patient data. To date, the formalism used for these analyses is based on kinetics derived from activity in a single organ, the target. This does not include the influence of other source organs to the dose and dose-rate to the target organ. As a result, only self-dose irradiation in the target organ contributes to the dose-rate. In this work, the biological effective dose (BED) formalism has been extended to include the effect of multiple source organ contributions to the net dose-rate in a target organ. The generalized BED derivation has been based on the Medical Internal Radionuclide Dose Committee (MIRD) schema assuming multiple source organs following exponential effective clearance of the radionuclide. A BED-based approach to determine the largest safe dose to critical organs has also been developed. The extended BED formalism is applied to red marrow dosimetry, as well as kidney dosimetry considering the cortex and the medulla separately, since both those organs are commonly dose limiting in radionuclide therapy. The analysis shows that because the red marrow is an early responding tissue (high alpha/beta), it is less susceptible to unexpected toxicity arising from rapid clearance of high levels of administered activity in the marrow or in the remainder of the body. In kidney dosimetry, the study demonstrates a complex interplay between clearance of activity in the cortex and the medulla, as well as the initial activity ratio and the S value ratio between the two. In some scenarios, projected BED based on both the cortex and the medulla is a more appropriate constraint on the administered activity than the BED based on the cortex only. Furthermore, different fractionated regimens were considered to reduce renal toxicity. The MIRD-based BED formalism is expected to be useful for patient-specific adjustments of activity and to facilitate the investigation of dose-toxicity correlations with respect to dose-rate and tissue repair mechanism.

The use of BED and EUD concepts in heterogeneous radioactivity distributions on a multicellular scale for targeted radionuclide therapy

There is evidence that nonuniform activity distributions within tumors might cause targeted radionuclide therapy (TRT) to fail. The aim of this study was to investigate the effects of the temporal and spatial behavior of the radioactivity in TRT, focusing on heterogeneous radiopharmaceutical distributions at a multicellular scale. Various activity distributions at the multicellular level from three radionuclides ((32)P, (90)Y, and (131)I) were simulated in cubic matrices (1- and 3-mm side). The in-house software package DOVE was used to calculate dose-rate maps, and survival fractions were calculated taking into account an up-take and a clearance phase. The effect from nonuniform activity distributions was analyzed in terms of dose volume histograms (DVHs), biologically effective dose (BED), and the effective uniform dose (EUD). The fraction of the absorbed dose that is "wasted," without producing a biological effect to the treatment, reaches 60% in the highly nonuniform distributions. For (32)P and (90)Y, the loss of therapeutic effectiveness was shown to be less than for (131)I. However, (90)Y, owing to its shorter physical half-life, presented lower mean BED values in almost every geometry, compared to (32)P and (131)I, and thus was less effective. (131)I, among all geometries, appeared to be more effective in more homogeneous activity distributions and in the 1-mm volume of interest, whereas it was the least effective radionuclide in the more heterogeneous activity distributions. (32)P presented the highest values of EUD, compared to (90)Y and (131)I. The EUD is a unique value that facilitates comparisons between different activity distributions in terms of treatment outcome. This study showed that as the degree of the heterogeneity in the dose distributions increases, the therapy effectiveness worsens. Nonuniform absorbed dose distributions can create a situation in which a fraction of cells are underirradiated, while another fraction of cells is "over-killed."

Theoretical study of the influence of a heterogeneous activity distribution on intratumoral absorbed dose distribution

Radioimmunotherapy of hematopoeitic cancers and micrometastases has been shown to have significant therapeutic benefit. The treatment of solid tumors with radionuclide therapy has been less successful. Previous investigations of intratumoral activity distribution and studies on intratumoral drug delivery suggest that a probable reason for the disappointing results in solid tumor treatment is nonuniform intratumoral distribution coupled with restricted intratumoral drug penetrance, thus inhibiting antineoplastic agents from reaching the tumor's center. This paper describes a nonuniform intratumoral activity distribution identified by limited radiolabeled tracer diffusion from tumor surface to tumor center. This activity was simulated using techniques that allowed the absorbed dose distributions to be estimated using different intratumoral diffusion capabilities and calculated for tumors of varying diameters. The influences of these absorbed dose distributions on solid tumor radionuclide therapy are also discussed. The absorbed dose distribution was calculated using the dose point kernel method that provided for the application of a three-dimensional (3D) convolution between a dose rate kernel function and an activity distribution function. These functions were incorporated into 3D matrices with voxels measuring 0.10 x 0.10 x 0.10 mm3. At this point fast Fourier transform (FFT) and multiplication in frequency domain followed by inverse FFT (iFFT) were used to effect this phase of the dose calculation process. The absorbed dose distribution for tumors of 1, 3, 5, 10, and 15 mm in diameter were studied. Using the therapeutic radionuclides of 131I, 186Re, 188Re, and 90Y, the total average dose, center dose, and surface dose for each of the different tumor diameters were reported. The absorbed dose in the nearby normal tissue was also evaluated. When the tumor diameters exceed 15 mm, a much lower tumor center dose is delivered compared with tumors between 3 and 5 mm in diameter. Based on these findings, the use of higher beta-energy radionuclides, such as 188Re and 90Y is more effective in delivering a higher absorbed dose to the tumor center at tumor diameters around 10 mm.

The Application of Dose-Rate Volume Histograms and Survival Fractions to Multicellular Dosimetry

The distribution of therapeutic radiopharmaceuticals in volumes smaller than those that can be fully resolved by the imaging system, such as by PET and SPECT scanners, is usually assumed to be homogeneous. The aim of this study was to investigate the implications of such an assumption at a scale that can be defined as multicellular for heterogeneous activity localizations of (32)P, (90)Y, and (131)I. Dose-rate distributions from heterogeneous radioactivity uptakes have been calculated in cubic volumes of 1, 3, and 4 mm using the in-house software package DOVE. These have been studied by the use of dose-rate volume histograms, and the influence of the heterogeneous dose distribution on the treatment outcome has been analyzed by the calculation of Integral Survival Fractions. The results showed that the effect of the heterogeneous localization of the compound can be overridden by the amount of radioactivity administered. However, significant variations in the survival probability distributions have been observed, depending on the amount of initial activity considered, the activity configuration, the radionuclide, and the time over which the energy was deposited. It has been shown, for example, that the ability of longer-range beta emitters, such as (32)P and (90)Y, to invalidate heterogeneous dose-rate distributions may be negated by the decay rate of the radioactivity.

A theoretical radiobiological assessment of the influence of radionuclide half-life on tumor response in targeted radiotherapy when a constant kidney toxicity is maintained

The potential of targeted radionuclide therapy may be limited if the antibody affinity to the tumor is relatively low and if significant normal tissue damage occurs while the tumor is sterilized. One way to increase the efficiency of the antibody-radionuclide complex might be to use knowledge of the radiobiological processes to select a near-optimal radionuclide half-life. In this paper, the role of physical half-life in targeted radiotherapy optimization is investigated using the linear quadratic (LQ) radiobiological model in conjunction with a range of radiobiological parameters relevant to the tumor. Five radionuclides ((211)At, (90)Y, (131)I, (86)Rb, and (114m)In) were selected, providing a half-life range from 0.3-49.5 days. The dose-limiting organ was assumed to be the kidney, with a simple fractional link between the initial (extrapolated) dose-rate to the tumor and the initial dose-rate to the kidney. The results suggest that short-lived radionuclides (half-life in the range of 1-10 days) have an advantage over medium- and long-lived radionuclides. Furthermore, for very rapid tumor uptake (uptake half-time of a few hours), very short-lived radionuclides (half-life of less than 1 day) could be efficiently employed. Ultimately, however, treatment outcome (in terms of tumor cell kill) is limited by the antibody affinity to the tumor.

Use of the linear-quadratic radiobiological model for quantifying kidney response in targeted radiotherapy

This paper reviews the generalized application of the linear quadratic (LQ) model of radiobiological effect to targeted radiotherapy. Special attention is given to formulations for normal tissue responses and these are applied, in particular, to the kidney. Because it is derived from self-consistent bio-physical principles, the LQ model currently remains the standard formalism for assessing biological responses for the whole range of radiotherapy treatments. A central feature of the model is the derivation of biologically effective doses (BEDs), which may be used to quantify the impact of a treatment on both tumors and normal tissues. BEDs are routinely derived for conventional external-beam treatments. The likely limits of targeted radiotherapy may, thus, be assessed by comparing the expected normal-tissue BEDs for such treatments with those known to be just tolerable in conventional therapy. The main parameters required in the model are defined, and data are provided which demonstrate the tentative link between targeted radiotherapy doses and those used in conventional radiotherapy. The extension of the LQ method to targeted radiotherapy involves using parameters for which the numerical values may not be accurately known at present. This places a restriction on the overall predictive accuracy of the model and the necessary caveats are, therefore, outlined.

Dosimetric model for intraperitoneal targeted liposomal radioimmunotherapy of ovarian cancer micrometastases

A simple model has been developed to investigate the dosimetry of micrometastases in the peritoneal cavity during intraperitoneal targeted liposomal radioimmunotherapy. The model is applied to free-floating tumours with radii between 0.005 cm and 0.1 cm. Tumour dose is assumed to come from two sources: free liposomes in solution in the peritoneal cavity and liposomes bound to the surface of the micrometastases. It is assumed that liposomes do not penetrate beyond the surface of the tumours and that the total amount of surface antigen does not change over the course of treatment. Integrated tumour doses are expressed as a function of biological parameters that describe the rates at which liposomes bind to and unbind from the tumour surface, the rate at which liposomes escape from the peritoneal cavity and the tumour surface antigen density. Integrated doses are translated into time-dependent tumour control probabilities (TCPs). The results of the work are illustrated in the context of a therapy in which liposomes labelled with Re-188 are targeted at ovarian cancer cells that express the surface antigen CA-125. The time required to produce a TCP of 95% is used to investigate the importance of the various parameters. The relative contributions of surface-bound radioactivity and unbound radioactivity are used to assess the conditions required for a targeted approach to provide an improvement over a non-targeted approach during intraperitoneal radiation therapy. Using Re-188 as the radionuclide, the model suggests that, for microscopic tumours, the relative importance of the surface-bound radioactivity increases with tumour size. This is evidenced by the requirement for larger antigen densities on smaller tumours to affect an improvement in the time required to produce a TCP of 95%. This is because for the smallest tumours considered, the unbound radioactivity is often capable of exerting a tumouricidal effect before the targeting agent has time to accumulate significantly on the tumour surface.

A model-based approach for the optimization of radioimmunotherapy through antibody design and radionuclide selection

BACKGROUND

The effectiveness of radioimmunotherapy (RIT) is known to depend on at least six factors: total absorbed dose and pattern of delivery, radiosensitivity, rate of repair of sublethal damage, ongoing proliferation during treatment, tumor heterogeneity, and tumor size. The purpose of this study was to develop a mathematic model that would relate the absorbed dose and its pattern of delivery to tumor response by incorporating information on each factor. This model was used to optimize therapeutic efficacy in mice by matching the antibody and radionuclide characteristics while ensuring recoverable marrow toxicity.

METHODS

Pharmacokinetic data were acquired in mice for a range of antibodies that varied in molecular weight, specificity, affinity, and avidity, and for a range of tumor sizes. This information was combined with the properties of iodine-131, rhenium-86, and yttrium-90 to determine the pattern of dose delivery. Tumor response was characterized in terms of radiosensitivity, rate of repair, and proliferation. Values for these parameters were obtained from in vitro assays and were incorporated into a response model based on the linear-quadratic model. Storage phosphor plate technology was used to acquire images of antibody distribution in tumor sections. These were registered with corresponding images showing tumor morphology, which were subsequently used to delineate regions that were distinct in terms of their response to radiation: oxygenated, radiosensitive areas that contained viable cells and hypoxic areas containing resistant viable cells and necrotic cells. Beta point dose kernels were then used to estimate the absorbed dose distribution in these regions.

RESULTS

Therapy in normoxic areas was more effective than in hypoxic areas. The multivalent, tumor-specific antibodies, with intermediate clearance rates, delivered the highest absorbed dose to viable tumor cells. Antibody affinity and avidity facilitated the prolonged retention in radiosensitive areas of tumor, where most of the dose was deposited. The effectiveness of therapy could be enhanced further by matching the radionuclide with the antibody and tumor size.

CONCLUSIONS

The model presented in this article allows the interaction between important radiobiologic parameters to be assessed and provides a tool for optimizing therapy in animal models and in patients.

Optimizing radioimmunotherapy by matching dose distribution with tumor structure using 3D reconstructions of serial images

The biological effect of radioimmunotherapy (RIT) is most commonly assessed in terms of the absorbed radiation dose. In tumor, conventional dosimetry methods assume a uniform radionuclide and calculate a mean dose throughout the tumor. However, the vasculature of solid tumors tends to be highly irregular and the systemic delivery of antibodies is therefore heterogeneous. Tumor-specific antibodies preferentially localize in the viable, radiosensitive parts of the tumor whereas non-specific antibodies can penetrate into the necrosis where the dose is wasted. As a result, the observed biological effect can be very different to the predicted effect from conventional dose estimates. The purpose of this study is to assess the potential for optimizing the biological effect of RIT by matching the dose-distribution with tumor structure through the selection of appropriate antibodies and radionuclides. Storage phosphor plate technology was used to acquire images of the antibody distribution in serial tumor sections. Images of the distributions of a trivalent (TFM), bivalent (A5B7-IgG), monovalent (MFE-23) and a non-specific antibody (MOPC) were obtained. These images were registered with corresponding images showing tumor morphology. Serial images were reconstructed to form 3D maps of the antibody distribution and tumor structure. Convolution of the image of antibody distribution with beta dose point kernals generated dose-rate distributions for 14C, 131I and 90Y. These were statistically compared with the tumor structure. The highest correlation was obtained for the multivalent antibodies combined with 131I, due to specific retention in viable areas of tumor coupled with the fact that much of the dose was deposted locally. With decreasing avidity the correlation also decreased and with the non-specific antibody this correlation was negative, indicating higher concentrations in the necrotic regions. In conclusion, the dose distribution can be optimized in tumor by selecting the appropriate antibodies and radionuclides. This has the potential to lead to a considerable enhancement of the efficacy of RIT in the clinic.

Theoretical estimation of absorbed dose to organs in radioimmunotherapy using radionuclides with multiple unstable daughters

The toxicity and clinical utility of long-lived alpha emitters such as Ac-225 and Ra-223 will depend upon the fate of alpha-particle emitting unstable intermediates generated after decay of the conjugated parent. For example, decay of Ac-225 to a stable element yields four alpha particles and seven radionuclides. Each of these progeny has its own free-state biodistribution and characteristic half-life. Therefore, their inclusion for a more accurate prediction of absorbed dose and potential toxicity requires a formalism that takes these factors into consideration as well. To facilitate the incorporation of such intermediates into the dose calculation, a previously developed methodology (model 1) has been extended. Two new models (models 2 and 3) for allocation of daughter products are introduced and are compared with the previously developed model. Model 1 restricts the transport to a function that yields either the place of origin or the place(s) of biodistribution depending on the half-life of the parent radionuclide. Model 2 includes the transient time within the bloodstream and model 3 incorporates additional binding at or within the tumor. This means that model 2 also allows for radionuclide decay and further daughter production while moving from one location to the next and that model 3 relaxes the constraint that the residence time within the tumor is solely based on the half-life of the parent. The models are used to estimate normal organ absorbed doses for the following parent radionuclides: Ac-225, Pb-212, At-211, Ra-223, and Bi-213. Model simulations are for a 0.1 g rapidly accessible tumor and a 10 g solid tumor. Additionally, the effects of varying radiolabled carrier molecule purity and amount of carrier molecules, as well as tumor cell antigen saturation are examined. The results indicate that there is a distinct advantage in using parent radionuclides such as Ac-225 or Ra-223, each having a half-life of more than 10 days and yielding four alpha particles per parent decay, in that lower doses to normal organs result for a given tumor dose in comparison to those radionuclides yielding fewer alpha particles. In model 2, which accounts for transit time through the blood, a dose of 20 Gy to a rapidly accessible 0.1 g tumor will result in a liver and kidney dose of 1.7 and 0.9 Gy, respectively from Ac-225. An equivalent dose to tumor from Ra-223 would yield a maximum normal organ dose of 0.4 and 0.3 Gy to bone and small intestines, respectively; the corresponding absorbed dose to small intestines from Pb-212 and Bi-213 is 2.2 and 3.0 Gy, respectively.

Physical and chemical properties of radionuclide therapy

As more radionuclide therapies move from laboratory feasibility studies into clinical reality, it becomes increasingly important for the labeling chemistry to produce consistently a stable radiopharmaceutical that remains intact under the challenge of human catabolism. Similarly, once proof of principle is established to bring a radionuclide conjugate into clinical therapy trials, dosimetric estimates should be made to select the appropriate radionuclide properties, which are based on animal-specific or patient-specific pharmacokinetics and match a set of specific clinical endpoints. These properties may include the radionuclide physical half-life, radiolabeled conjugate biological uptake and clearance, product-specific activity, range and type of emissions, and resultant effects on tumor and normal tissue cellular survival. The immunologist and labeling chemist have now produced a variety of strategies that have potential to increase the therapeutic ratio (tumor-to-normal tissue dose ratio). The advent of normal tissue clearing agents, fragmented or chimerized carriers to improve targeting, and the method of bispecific or two-step and three-step targeting agents has increased the need for realistic modeling of the carrier in vivo to guide prospectively the competitive development of these radiopharmaceuticals. In this article, examples have been taken from the literature to elucidate the benchmark of success that careful experimental design has fostered to bring these agents into clinical practice by creative and logical methodologies.

Effects of cellular repair and proliferation on targeted radionuclide therapy: a modelling study

A previous targeted radionuclide therapy modelling study has been extended to include the radiobiological effects of cellular repair and proliferation. Dose distributions have been converted into biologically effective dose (BED) distributions using a previously published formulation. With suitable estimated parameters, corrected tumour control probability (TCP) values were derived. The dependence of BED on the physical half-life of the radionuclide was also modelled. Results indicate that the TCP is greater when a shorter physical half-life is employed.

114mIn, a candidate for radionuclide therapy: low-energy cyclotron production and labeling of DTPA-D-phe-octreotide

A method for production of carrier-free (114m)In (half-life 49.5 days), which is a potential radionuclide for radionuclide therapy of slowly growing tumors, is presented. A target consisting of five enriched cadmium ((114)Cd) foils, each 50 microm thick, was irradiated by protons (from 12.6-6.5 MeV) giving a target yield of 0. 8 MBq/microAh. A simple and cost-efficient thermal diffusion method was used for the separation. The irradiated target foils were heated for 2 h at 306 degrees C and then etched in 0.05 M HCl. The obtained cadmium/indium solution was purified using a cation ion-exchange resin (AG 1 x 8, Bio-Rad Laboratories, Hercules, CA USA). An overall yield of approximately 60% was obtained, whereas the loss of the target material was <1% per separation cycle. The (114m)In production gave (114m)In with high specific radioactivity and was successfully used to label diethylenetriamine pentaacetic acid (DTPA)-D-Phe-octreotide. Furthermore, no difference in biodistribution between [(114m)In]- and [(111)In]-DTPA-D-Phe(1)-octreotide in tumor-bearing nude mice was seen. The high radionuclide uptake in the tumors indicates a good receptor binding of the labeled octreotide.