Abstract
Background:
α particle emitting bismuth (212Bi) as decay product of 212Pb-labeled pharmaceuticals has been effective in targeted α particle therapy (TAT). Estimating the contribution of 212Bi released from its chelator to the absorbed doses in nontarget tissues is challenging in TAT. Physiologically based pharmacokinetic (PBPK) modeling can help overcome this limitation. Therefore, a whole-body 212Bi–PBPK model was developed to describe the pharmacokinetics (PKs) of 212Bi in rats.
Materials and Methods:
The rat 212Bi–PBPK model was implemented using the modeling software SAAM II with data and parameter values from the literature. Besides other mechanisms, 212Bi interactions with red blood cells, high molecular weight plasma protein, and intracellular biological thiols are described. Important PK parameters were fitted to time–activity data. Absorbed dose coefficients (ADCs) were calculated for injecting 0.774 fmol of 212Bi.
Results:
212Bi uptake rates of liver, bone, small intestine, bone marrow, skin, and muscle were (0.86 ± 0.13), (3.85 ± 0.63), (0.27 ± 0.05), (1.44 ± 0.29), (0.04 ± 0.01), and (0.007 ± 0.007) per min with corresponding ADCs of 0.09, 0.03, 0.03, 0.07, 0.01, and 0.003 mGy/kBq, respectively. An ADC of 0.70 mGy/kBq was determined for kidneys.
Conclusions:
Kidneys are the dose-limiting organs in 212Bi-based TAT. The 212Bi–PBPK model is an effective tool to investigate the 212Bi biodistribution in murine models. Integrating the 212Bi–PBPK model into other murine and human PBPK models of α particle generators can help study the efficacy and safety of TAT.
Introduction
In vivo α particle generators of long-lived radioactive parents such as 212Pb have been recently used as in vivo sources of α particles through their α particle emitting daughters. 1 –3 However, this concept is limited by the release of 212Bi from the chelator that is subsequently redistributed within the body. The contribution of free 212Bi to absorbed doses in organs at risk depends on the stability of the 212Bi–chelator complex and the physiological uptake in different tissues. 4 Investigating the pharmacokinetics (PKs) of the free 212Bi is, therefore, vital to estimate the risk of radiation toxicity for accumulating organs, such as kidneys, due to the high biological effectiveness of the emitted α particles. 5
The systemic biokinetics of bismuth compounds as therapeutic agents were described in a number of preclinical and clinical studies. 6 –12 However, the reports of the PKs of 212Bi are scarce and limited to concentration–time curves because the 60.6 min half-life (t 1/2) of 212Bi is inconvenient for preparations and dissections. 13 Therefore, most of the studies were limited to either bismuth compounds or radioisotopes with long half-lives, such as 205Bi (t 1/2 = 14.5 d) and 206Bi (t 1/2 = 6.4 d). Consequently, the investigated time interval does not sufficiently cover the early time period after injection, which is relevant for 212Bi. Mathematical modeling can help to understand the PKs and the localization of 212Bi ions and estimate absorbed doses due to the uptake of 212Bi in normal tissues. 14 –16
In this study, a whole-body physiologically based pharmacokinetic (PBPK) model for 212Bi in rats was developed including all relevant physiological mechanisms with parameter values and biokinetic data from the literature. 13,17 The uptake rates of 212Bi by blood and tissues were estimated and tissue absorbed doses were subsequently calculated.
Materials and Methods
212Bi–PBPK model
The developed whole body PBPK model of 212Bi distribution in rats includes all relevant physiological mechanisms such as the distribution of 212Bi through blood flow, transcapillary transport, cellular uptake, and excretion with parameter values from the literature (Supplementary Appendix A1). 13,17 –26 212Bi binding to high molecular weight plasma proteins (HWPPs) and its uptake by red blood cells (RBCs) were also taken into account. 27 –29
The rat 212Bi–PBPK model comprises tissues, namely kidneys, liver, spleen, small intestine (SI) representing gastrointestinal track, bones, pancreas, brain, lung, heart, skin, bone marrow (BM), muscle, fat, and remainder of body. Each tissue is divided into vascular, interstitial, and cellular spaces (Fig. 1).

212Bi–PBPK model for rats. The model includes main physiological mechanisms to describe the PKs of 212Bi in rats based on its uptake by tissue cells and its association and dissociation with HWPPs and RBCs. Free 212Bi, HWPP- and RBC-bound 212Bi are distributed in the circulation to vascular spaces of body organs followed by free 212Bi diffusion and uptake by tissue cells. 212Bi binds to intracellular thiol molecules and accumulates intracellularly. The α particle emitter 212Bi in each compartment in the model decays with a half-life of 60.6 min. 4 BM, SI, and RB are BM, SI, and remainder of body. BM, bone marrow; HWPPs, high molecular weight plasma proteins; PBPK, physiologically based pharmacokinetic; PK, pharmacokinetic; RB, remainder of body; RBCs, red blood cells; SI, small intestine. Color images are available online.
The interactions of free 212Bi with HWPPs and RBCs were described by first-order kinetics. Both, HWPP- and RBC-bound 212Bi are confined to the corresponding subspaces in the vascular space of each organ. Free 212Bi diffuses to the interstitium and is passively transported into cells where it binds with high affinity to intracellular metallothionein (MT) metalloproteins and glutathione (GSH) ligands as reported by Hong et al. 30 Cellular uptake of free 212Bi was assumed to follow first-order kinetics using an effective value of the intracellular transfer rate of 212Bi due to its interactions with biological thiols. 27,30 Free 212Bi uptake rates by tissue cells were assumed to be proportional to the amount of MT in normal rats because 212Bi ions have higher affinity to MT than to GSH. 31 Cellular uptake rates of spleen, lung, and pancreas were assumed to be equal to that of SI cells. The uptake rates of heart, brain, fat, and remainder of body were assumed to be equal to that of muscle cells. This assumption was based on studies showing that muscle, skin, brain, and fat have low concentrations of accumulated bismuth. 6
HWPP- and RBC-bound 212Bi molecules in the vascular space of kidneys are not filtrated through the renal glomerular filtration barrier (GFB) because of their large molecular weight. In contrast, the small molecular weights of free 212Bi and MT- or GSH-bound 212Bi result in rapid filtration through the glomeruli and reabsorption by proximal tubular cells. 6 A fraction of filtered free 212Bi and MT- or GSH-bound 212Bi is excreted.
Model evaluation and parameter estimation
The 212Bi–PBPK model for rats was implemented using modeling software SAAM II version 2.3 (The Epsilon Group, USA). The SAAM II software supports systems of ordinary differential equations to design and simulate experiments and analyze biokinetic data. 32 The 212Bi–PBPK model was evaluated using nondecay corrected percentage of injected dose per gram (%ID/g) data for 212Bi in different tissues of unanesthetized rats (Supplementary Appendix A2). 13 The %ID/g of 212Bi in the plasma of the virtual rat, that is, the 212Bi–PBPK model, includes the plasma content of free 212Bi and HWPP-bound 212Bi. The injected amount of 212Bi ions in the virtual rat was 0.774 fmol of 88.8 kBq (2.4 μCi), which is the amount used in the experiment. 13 Values of PK parameters of the model were estimated by data fitting with the implemented Rosenbrock least-squares algorithm using the highest available threshold of the convergence criterion of 10−7. A relative data-weighting scheme with a standard deviation (SD) of 10% was assigned for all data. The goodness of fit was evaluated using visual inspection, coefficient of variation (CV), and the off-diagonal values of the parameters of the correlation matrix (−0.8 < CM <0.8 for most elements). 33
Dosimetry
Tissue absorbed dose coefficients (ADCs) due to distributed bound and free 212Bi in the virtual rat were calculated. The radioactive daughters of the branched decay of 212Bi, that is, 208Tl (t 1/2 = 3.1 min) and 212Po (t 1/2 = 0.3 μs), were assumed to decay at the site of production. The average energy emitted per decay of 212Bi was assumed to be equal to 7.8 MeV. 18
Results
The time–activity curves for different tissues, that is, kidneys, liver, bone, SI, BM, skin, muscle, plasma, and RBCs, showed acceptable fit (Supplementary Appendix A3) with R 2 ≥ 0.94. The estimated values of the PK parameters with the corresponding SD are given in Table 1. The cellular uptake rates of kidneys, liver, BM, bone, and SI were highest. For skin and muscle, low cellular uptake rates were observed. The calculated ADCs (mGy/kBq) in the tissues are given in Table 2. 26,34,35
Model Parameters Fitted to the Experimental Data
The CV of the fitted parameters was <50% for all tissues except muscles with CV ∼100%. All elements of the correlation matrix of the fitted parameters were <0.8.
CV, coefficient of variation; HWPP, high molecular weight plasma protein; RBCs, red blood cells.
Absorbed Dose Coefficients Resulting from the Distribution of 0.77 fmol of Free 212Bi Administered into the Plasma Vein Compartment of the Virtual Rat
ADC, absorbed dose coefficient.
Discussion and Conclusions
The PKs of the free α particle emitting 212Bi as a decay product of 212Pb in vivo α particle generators are a critical determinant of the efficacy and safety of targeted α particle therapy (TAT). A whole body 212Bi–PBPK model for rats was developed to determine the values of the PK parameters of the distributed 212Bi and the resulting absorbed doses in tissues.
The estimated first-order association rate of bismuth ions with HWPPs (0.20 min−1) was higher than that with RBCs (0.01/min) within a time frame of two half-lives of 212Bi. Sun et al. and Li et al. reported that bismuth ions bind with high affinity to transferrin and albumin in plasma. 36 –38 The ionic form of bismuth results in rapid binding of bismuth to plasma proteins. Because the GFB depends on the molecular size, the HWPP- and RBC-bound 212Bi would be less likely to enter the ultrafiltrate. 6 In addition, at low bismuth concentrations, Rao and Feldman pointed out that a smaller fraction of bismuth ions is available for the uptake by rat erythrocytes than at high concentrations. 28
Despite the variation in the systemic biodistribution of bismuth depending on the form of injection or route of administration, kidneys are considered as the main sequestering organ for bismuth. 6,11,39 –41 Bismuth ions are absorbed and sequestered by cysteine-rich intracellular ligands and proteins in the proximal epithelium, such as GSH and MT. 30,42,43 Renal uptake rate was fixed to 1/min based on the results of a local sensitivity analysis. Different uptake values of renal epithelial cells of 0.1, 1, and 10/min were fixed to check how they affect other estimated parameters. The estimated values of fitted parameters deviated <2% for the different rates of renal uptake.
The high uptake rates of 212Bi by bone and BM can be linked to the MT role in bone growth and hematopoiesis. 44 –48 The liver uptake rate was high as expected for a detoxifying organ and was followed by the cellular uptake rates for SI, skin, and muscle. SI, spleen, lung, and pancreas were assumed to have the same bismuth uptake rate. This assumption is based on the report of Russ et al. that spleen, SI, and lung show high concentration of bismuth. 6 Also, the results of Onosaka and Cherian were considered that SI has the same ability as pancreas to induce and store MT when exposed to Zn salts. 31 However, the uptakes of these tissues were proposed to be lower than that of the liver as the latter is a detoxifying organ. As gonads are reported to show low concentration of bismuth, the reproductive organs, that is, testes, are expected to have an uptake rate similar to that of muscle. The high estimated CV of the muscle uptake rate is due to the assumption that 212Bi uptake rates of most body tissues, that is, heart, brain, fat, and remainder of body, are equal to that of muscle. The order of the estimated uptake rates is consistent with the reported rank of bismuth concentrations in the tissues of rats. 6,8,40
Tissue sequestering of 212Bi may be related to the abundance of the intracellular MT and GSH relative to bismuth ions, the ability of 212Bi to induce the de novo synthesis of these biological thiols, and the short observation period of the simulation (150 min). The simulated retention of bismuth in tissue cells is compatible with the documented prolonged cell retention of free 212Bi. 49
Kidneys have the highest ADC (6.5 mGy/kBq) among other tissues. This is in line with the studies that showed that kidneys are the major dose limiting organs in TAT. 50,51 However, the estimated ADC values in normal tissues (Table 2) are expected to be larger than the cases in 212Bi-based TAT. All of the intravenously administered free 212Bi in the tumor-free rat will result in nontarget radiation toxicity. In contrast, in TAT, bismuth will be most probably released from bound or internalized pharmaceuticals at different locations within the body. Therefore, the estimated ADC values may demonstrate the worst-case scenario when all produced 212Bi are dissociated from the targeting pharmaceuticals in the circulation.
The estimated PK parameters of 212Bi can also describe the PKs of the α particle emitter 213Bi. The latter is a decay product of the in vivo α particle generator 225Ac with an even shorter half-life of 45.6 min. As PBPK modeling allows for the separation between radiopharmaceutical-specific and physiological parameters, it holds promise to integrate the estimated PK parameters of bismuth ions into other PBPK models for mice and humans. As a consequence, the 212Bi–PBPK model may be integrated into other human PBPK models of in vivo α particle generators, such as 212Pb- or 225Ac-labeled pharmaceuticals, to investigate the fate of free 212/213Bi radioisotopes as the main challenge encountered in TAT. The β− decay of the bound 212Pb produces different fractions of free 212Bi depending on the chelating systems used. Mirzadeh et al. reported that the 212Pb decay in the 212Pb–1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA) complex results in 36% dissociated 212Bi atoms due to the electronic excitation of the produced 212Bi daughter atoms. 52 The energy of the emitted gamma after 88% of the β− decays of the parent is transferred to the inner orbital electrons in ∼36% of the decays. The reorganization of valence electrons, induced by the cascade of conversion electrons, emits characteristic X-rays and Auger electrons. Subsequently, these electronic excitations generate a highly ionized Bi5+ cation, which dissociates from its DOTA chelator. 52 Using 1,4,7,10-tetrakis(carbamoylmethyl)-1,4,7,10-tetraazacyclododecane chelators in 212Pb in vivo α particle generators results in the dissociation of ∼16% of the produced 212Bi from the chelating system. 4 Also, the model enables for estimating the absorbed doses associated with nonspecific uptake of the distributed free 212Bi in TAT. In consequence, mathematical modeling can be an effective tool to estimate the radiation toxicity associated with the distributed free α particle emitters and assess the efficacy and safety of using in vivo α particle generators in TAT.
Footnotes
Authors' Contributions
N.R.R.Z. designed the study, developed, implemented, and tested the model, performed the simulations, and wrote the article. P.K. evaluated the model, wrote, and edited the article. G.W. discussed the model and edited the article. A.J.B. wrote and edited the article. G.G. designed the study, evaluated the implementation of the model and the simulations, wrote, and edited the article. All coauthors have reviewed and approved the article before submission.
Disclosure Statement
No competing financial interests exist.
Funding Information
The authors gratefully acknowledge the funding from the DAAD (German Academic Exchange Service, Research Grants, Doctoral Programs in Germany 2018/19-57381412) and the Research Campus M2 OLIE (German Federal Ministry of Education and Research (BMBF) within the Framework “Forschungscampus–public-private partnership for Innovations” funding code 13GW0389B).
Supplementary Material
Supplementary Appendix SA1
Supplementary Appendix SA2
Supplementary Appendix SA3
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
