Effect of tumor size on drug delivery to lung tumors M. Soltani, M. Sefidgar, H. azmara C. Marcus, R. M. Subramaniam, A. Rahmim Abstract Drug delivery to solid tumors can be exressed hysically using biomechanical henomena such as convection, diffusion of drug in extracellular matrices, and drug extravasation from microvessels. Alying comutational methods to solve governing conservation equations clarifies the mechanisms of drug delivery from the injection site to a solid tumor. In this study, multile tumor geometries were obtained from PET/CT images. An advanced numerical method was used to solve fluid flow and solute transort equations simultaneously to investigate the effect of tumor size on drug delivery to lung tumors. Data from 20 atients with lung tumors were analyzed and the tumor geometrical information including size, shae, and asect ratios were classified. In order to investigate effect of tumor size, tumors with similar shaes but different sizes ranging from 1 to 28.6 cm 3 were selected and analyzed. A hyothetical tumor, similar to one of the analyzed tumors but scaled to reduce its size to 0.2 cm 3, was also analyzed. An ideal bolus injection was considered for the model. The effects of two transort mechanisms, namely convection and diffusion, were considered in this study. The results show because of size of considered lung tumor, the diffusion transort rate is higher than convection transort rate. ased on governing equations, the diffusion transort is only deended on concentration gradient and indeendent of size of tumor, therefore the redicted concentration rofile for considered tumors are similar. When size of tumor is decreased significantly, the drug concentration is also significantly increased. I. INTRODUCTION ung cancer is the most common cause of cancer-related death in men and women, and was resonsible for 1.56 million deaths annually, as of 2012 [1]. Chemotheray is one of the ways widely used for cancer treatment. ased on the findings from clinical alications, most cancer treatments with drugs fail to eliminate solid tumors comletely [2]. The comutational method can investigate why systemic administration cannot distribute drug uniformly in tumors[3]. The drug exchange between microvessels and extracellular Manuscrit received December 05, 2015. M. Soltani is with Division of Nuclear Medicine, Deartment of Radiology and Radiological Science, School of Medicine, Johns Hokins University, MD 21287, USA (e-mail: ma_soltani@yahoo.com), and Mechanical Engineering, KNT University, Tehran, Iran. M. Sefidgar is with Deartment of Mechanical Engineering, Pardis ranch Islamic Azad University, Pardis new city, Iran (e-mail: sefidgar@gmail.com). H. azmara is with Deartment of Mechanical Engineering, Pardis ranch Islamic Azad University, Pardis new city, Iran (e-mail: bazmara@gmail.com). C. Marcus is with the Deartment of Radiology, Johns Hokins University, altimore, MD, USA 21287 (e-mail: cmarcus7@jhmi.edu). M. Subramaniam is with the Deartments of Radiology, Oncology, Otolaryngology, Head and Neck Surgery, and Health Policy and Management, Johns Hokins University, altimore, MD, USA 21287 (e-mail: rsubram4@jhmi.edu). A. Rahmim is with the Deartments of Radiology, and Electrical & Comuter Engineering, Johns Hokins University, altimore, MD, USA 21287 (telehone: 410-502-8579, e-mail: arahmim1@jhmi.edu, webage: www.jhu.edu/rahmim). matrices, drug removal by lymhatic system, drug diffusion and convective transort in extracellular matrices should be included by mathematical simulation. Comutational fluid dynamics (CFD) can model the whole drug delivery rocess and clarify the mechanisms of drug delivery from the injection site to absortion by a solid tumor [4]. axter and Jain, based on the theoretical framework in their 1D mathematical method, found the effective factors on drug delivery such as microvessel ermeability, interstitial fluid ressure (IFP), and interstitial fluid velocity (IF) [5] [8]. Wang et al. [9] [11] develoed a simulation framework of drug delivery to tumors by considering the comlex 3D geometry. Wang and i [9] used modified MRI images for tumor geometry. They considered interstitial fluid flow with blood and lymhatic drainage in their model. Wang et al. [10] studied the effect of elevated interstitial ressure, convective flux, and blood drainage on the delivery of secified solute to brain tumors. The study of tissue transort roerty effect on drug delivery has been considered in a number of recent studies. Zhao et al. [12] used a 3D comutational model to redict the distribution of IF, IFP, and solute transort through a tumor. Arifin et al. [13], [14] studied the sensitivity of drug distribution to hysiochemical roerties in realistic models of brain tumors. A secific tumor catured by MRI was used by Pishko et al. [15] [17] for modeling drug distribution in tissue with satially-varying orosity and vascular ermeability. The sensitivity of solute distribution to tumor shae and size was not considered in the above-mentioned works. Develoed models of fluid flow in tumors have also focused on incororating satial and temoral variations in blood flow using microvascular network models [18] [22]. Soltani et al. [23] investigated the effect of tumor shae and size on drug delivery by modeling interstitial fluid flow and assuming that drug articles flow within the interstitial fluid. Sefidgar et al. [24] further develoed the mathematical model by adding the solute transort equation to reviously develoed models [25] [29]. New governing equations were subsequently investigated to find drug concentration in interstitial fluid (DCIF). Sherical and non-sherical tumors and their surrounding normal tissue were modeled assuming rigid orous media. We alied [30], [31] this numerical model to investigate drug distributions in clinically imaged ancreatic tumors. We also imlemented [32] [34] this mathematical method to multi scale models of tumor microvasculature for investigation of drug and tracer delivery distributions. Other arameters related to cancer fluid flow and solute transort were also investigated by our grou [35], [36]. In this study, multile tumor geometries as obtained from clinical PET/CT images of the lung were considered. An
advanced numerical method was used to simultaneously solve fluid flow and solute transort equations in order to investigate the effect of tumor shae and size on drug delivery to solid tumors. Data from n=20 lung cancer atients with nonresectable locoregional disease were analyzed, and geometrical information from the tumors including size were classified. To investigate effect of tumor size, tumors with similar shaes but different sizes, ranging from 1 to 28.6 cm 3, were selected and analyzed. A hyothetical tumor similar to one of the analyzed tumors, but scaled to reduce its size to 0.2 cm 3, was also analyzed. An ideal bolus injection was considered for the model. Solving fluid flow and solute transort equations simultaneously, the effects of lung tumor size on drug delivery to solid tumor are also investigated. II. METHODS A. Image Data Extraction PET/CT images from 20 atients with lung tumor were obtained and rocessed to extract the geometrical configuration of tumors. Tumors volume range was from 2.8 to 28.6 cm 3. The tumors used in this study are shown in Fig. 1. Fig. 1 A systematic flow chart of comutational technique. Mathematical Method: A systematic flow chart is illustrated in Fig. 2 to clarify the comutational techniques involved in this work. Next, the governing equation is introduced. Fig. 2 A systematic flow chart of comutational technique Tumor #2 olume=2.8cm 3 Tumor #3 olume=28.6cm 3 Tumor #8 olume=14.7cm 3 Tumor #10 olume=23.7cm 3 Tumor #11 olume=7.7cm 3 Tumor #18 olume=12.2cm 3 Tumor #20 olume=6.4cm 3 1) Mathematical Model of Interstitial Flow Since tumor tissue has characteristics the same as orous media, fluid flow behavior is defined by couling the fluid flow governing equations. Combination of Darcy s law and the mass balance s law results in [37] [39]: 2 κ P i = φ φ (1) P i : Interstitial fluid ressure, κ : hydraulic conductivity of the interstitium, φ : The source term, extravasation from microvessels, and φ : The drainage term, elimination by lymhatic system. In biological tissues, the fluid source is evaluated through Starling's law as follows[40] [42]: S P φ = ( P Pi σ s ( π b π i )) (2) P : lood ressure in microvessel, S : The surface area er unit volume of tissue for transort in the interstitium, π : Microvessel oncotic ressure, π : Interstitial oncotic ressure, i : The hydraulic conductivity of the microvessel wall, and s σ : Osmotic reflection coefficient. and the lymhatic system is related to the ressure difference between the interstitium and the lymhatic vessels and is considered only for normal tissues [4], [43]:
PS ( Pi P) Normaltissue φ ( r ) = (3) 0 Tumor tissue φ : The volumetric flow rate into the lymhatic, PS : The lymhatic filtration coefficient, and P : The hydrostatic ressure of the lymhatic. 2) Solute Transort: The interstitial transort of drug is governed by the convection diffusion equation; therefore, the general equation for the mass balance of solutes can be written as[44] [46]: C = D C vc + Φ Φ (4) ( eff ) ( i ) ( ) t C : The solute concentration based on tissue volume, Φ : The rate of solute transort er unit volume from microvessel into the interstitial sace, Φ : The rate of solute transort er unit volume from the interstitial sace into lymhatic vessels, and D eff : The effective diffusion tensor. For an isotroic and uniform diffusion in tissues, equation (4) can be written as: C 2 = D eff C ( v f C ) + ( Φ Φ ) (5) t The solute transort rate across the lymhatic vessels can be considered as [38], [47], [48]: φ C Normal Tissue Φ = (6) 0 Tumor Tissue The solute transort rate across the microvessel can be reresented by Patlak equation[18], [49]: PS Pe Φ = φ( 1 σ f ) CP + ( CP C) (7) Pe e 1 ( 1 σ ) φ f Pe = (8) PS φ : The fluid flow rate er unit volume of tissue across the microvessel wall, σ f : The filtration reflection coefficient, P : The microvessel ermeability coefficient, S/ : The microvessel surface area er unit volume of tissue, and C : Solute concentration in the lasma. The detail of solution and value of arameter used in above equation are mentioned in our revious work[24], [50]. C. Model arameterization The interstitial transort roerties for normal and tumor tissue are listed in Table 1. These values are used as baseline and some of them are investigated and changed in secified ranges for sensitivity analysis. The arameters of solute transort model taken from axter and Jain [43] are listed in Table 2. Although, the numerical model is alicable for any tye of drug, in resent study the roerties of Fragment antigen-binding (F(ab ) 2 ) as a samle is used. TAE 1. PARAMETER USED IN THIS SIMUATION [24] Parameter Descrition value P [m/pa s] Hydraulic conductivity of vessel 2.10 10-11 κ [m 2 Hydraulic conductivity of the /Pa s] interstitium 3.10 10-14 S/[m -1 ] Surface area er unit volume of tissue for transort in the 20000 interstitium P [Pa] Intravascular blood ressure 2075 π [Pa] Caillary oncotic ressure 2666 π i [Pa] Interstitial oncotic ressure 2000 σ Osmotic reflection coefficient 0.82 P [Pa] Hydrostatic ressure of the lymhatic 0 P S / ymhatic filtration coefficient [1/Pa s] 1 10-7 σ f Filtration reflection coefficient 0.9 D eff [m 2 /s] Effective diffusion tensor 2.0 10-12 P[m/s] Microvessel ermeability coefficient 17.3 10-11 III. RESUTS DCIF is simulated in case of the bolus injection in which the lasma concentration decreases with time exonentially 0 t ( C = C e τ ), in which τ is the drug half-life in lasma 0 (6.1hr). DCIF are non-dimensionalized by C. The Non-dimensionalized DCIF distribution is shown in Fig. 3 at 7 hr ost injection, for different sizes of tumor. The red color in Figure 2 shows the highest drug concentrations, which encomasses mostly the tumor boundary. Fig. 4 shows DCIF for tumor of a atient in which the size of the tumor was intentionally reduced to 0.2 cm 3. Fig. 5 shows the temoral variation of average drug concentration in each tumor. There is no significant difference in drug udate between tumors with different sizes. One imortant metric of disease develoment and resonse to cancer via drug delivery is volume of tumors. ased on our revious research [2, 5] when the tumor volume is in the order of >1 cm 3 (for all considered tumor in this study), the sensitivity of drug concentration to tumor size decreases. Only for a hyothetical tumor which has size less than 1 cm 3, DCIF is high and convection rate has effect on drug delivery. In the analyzed lung tumors, since the interstitial ressure is high in the tumor region, the convection rate vanishes and the diffusion rate reaches a constant value, and consequently the sensitivity of drug concentration to tumor size reduces.
Patient 2, Patient 20, Patient 3, Fig. 3 DCIF for different tyes of tumor at 7 hr ost injection for a bolus well-erfused injection Patient 8, Fig. 4 DCIF for the hyothetical tumor at 7 hr ost injection for a bolus well-erfused injection Patient 10, Patient 11, Fig. 5 Average drug concentration for different sizes of tumors shown in Fig. 1. Patient 18, I. DISCUSSION Drug transort in tissue is deended on two mechanisms, diffusion and convection. These mechanisms affect drug transort from vessels to the extracellular matrix, and transort within the extracellular matrix. The convection mechanism is deended on interstitial velocity for transort in tissue and on differences between intravascular and interstitial ressures for transort from vessels. Our revious studies [24], [35] indicated these arameters (velocity and ressure) as deended on shae and size of tumor. The diffusion mechanism is only
deendent on the concentration gradient and therefore is indeendent of shae and size. In large tumor (order of 10 cm 3 ), the interstitial velocity in tumors goes to zero and ressure differences are minimal [25], [38]. Therefore, the effect of convection transort is reduced and sometimes negligible in large tumors. In the considered lung tumors, because of the large sizes involved, drug delivery is only deended on diffusion mechanisms, and size of tumors does not affect drug delivery. For the hyothetical tumor, because of small size, the effect of convection transort is higher than in lager tumors, and as such, drug delivery in increasingly small tumors is greater than in large tumors.. CONCUSION In this study, we investigated drug delivery to different lung tumor size. Since, in large tumors, the convection mechanism is negligible and diffusion mechanism indeendent of tumor size is dominant, and the investigated tumors in size were of the order of 10 cm 3, drug concentrations are relatively similar. Only for smaller tumors with size lower than 1 cm 3, the convection mechanism imacts drug delivery, and drug concentration is subsequently increased more than 20%. ACKNOWEDGMENT We gratefully acknowledge suort by Siemens Medical Solutions. References [1] World Cancer Reort 2014, in World Health Organization, 2014,. Chater 1.1. ISN 9283204298. [2] S. lakeslee, Imenetrable tumors found to block even the newest cancer agents, The New York Times, 1989. [3] R. K. Jain, Delivery of Molecular and Cellular Medicine to Solid Tumors, Microcirculation, vol. 4, no. 1,. 1 23, 1997. [4] M. Wu, H.. Frieboes, S. R. McDougall, M. a J. Chalain,. Cristini, and J. owengrub, The effect of interstitial ressure on tumor growth: couling with the blood and lymhatic vascular systems., J. Theor. iol., vol. 320,. 131 151, Mar. 2013. [5]. T. axter and R. K. Jain, Transort of fluid and macromolecules in tumors. I. A microscoic model of the erivascular distribution., Microvasc. Res., vol. 41, no. 2,. 252 272, Mar. 1991. [6]. T. axter and R. K. Jain, Transort of fluid and macromolecules in tumors. (II) role of heterogeneous erfusion and lymhatics, Microvasc. Res., vol. 40,. 246 263, 1990. [7]. T. axter and R. K. Jain, Transort of Fluid and Macromolecules in Tumors III. Role of inding and Metabolism, Microvasc. Res., vol. 41,. 5 23, 1991. [8] R. K. Jain and. T. axter, Mechanisms of heterogeneous distribution of monoclonal antibodies and other macromolecules in tumors: significance of elevated interstitial ressure, Cancer. Res., vol. 48, no. 24 Pt 1,. 7022 7032, Dec. 1988. [9] C. H. Wang and J. i, Three-dimensional simulation of IgG delivery to tumors, Chem. Eng. Sci., vol. 53, no. 20,. 3579 3600, 1998. [10] C. H. Wang, J. i, C. S. Teo, and T. ee, The delivery of CNU to brain tumors., J. Control. Release, vol. 61, no. 1 2,. 21 41, 1999. [11] W. H. K. Tan, F. Wang, T. ee, and C. H. Wang, Comuter simulation of the delivery of etanidazole to brain tumor from PGA wafers: comarison between linear and double burst release systems., iotechnol. ioeng., vol. 82, no. 3,. 278 88, 2003. [12] J. Zhao, H. Salmon, and M. Sarntinoranont, Effect of heterogeneous vasculature on interstitial transort within a solid tumor, Microvasc. Res., vol. 73, no. 3,. 224 236, 2007. [13] D. Y. Arifin, K. Y. T. ee, and C.-H. Wang, Chemotheraeutic drug transort to brain tumor, J. Control. Release, vol. 137, no. 3,. 203 210, 2009. [14] D. Y. Arifin, C. Wang, and K. A. Smith, Patient-secific Chemotheraeutic Drug Delivery to rain Tumors brain tissue tumor ventricle, Mimics Innov. Award.,. 1 9, 2007. [15] G.. Pishko, G. W. Astary, T. H. Mareci, and M. Sarntinoranont, Sensitivity Analysis of an Image-ased Solid Tumor Comutational Model with Heterogeneous asculature and Porosity, Ann iomed Eng., vol. 39, no. 9,. 2360 2373, 2011. [16] K. N. Magdoom, G.. Pishko,. Rice, C. Pamo, D. W. Siemann, and M. Sarntinoranont, MRI-based comutational model of heterogeneous tracer transort following local infusion into a mouse hind limb tumor., PoS One, vol. 9, no. 3,. e89594, Jan. 2014. [17] G.. Pishko, G. W. Astary, J. Zhang, T. H. Mareci, and M. Sarntinoranont, Role of convection and diffusion on DCE-MRI arameters in low leakiness KHT sarcomas., Microvasc. Res., vol. 84, no. 3,. 306 13, Nov. 2012. [18]. P. Chauhan, T. Stylianooulos, J. D. Martin, O. Chen, W. S. Kamoun, M. G. awendi, D. Fukumura, and R. K. Jain, Normalization of tumour blood vessels imroves the delivery of nanomedicines in a size-deendent manner, Nat. Nanotechnol., vol. 7, no. 6,. 383 388, 2012. [19] T. Stylianooulos, K. Soteriou, D. Fukumura, and R. K. Jain, Cationic nanoarticles have suerior transvascular flux into solid tumors: insights from a mathematical model., Ann. iomed. Eng., vol. 41, no. 1,. 68 77, Jan. 2013. [20] T. Stylianooulos and R. K. Jain, Combining two strategies to imrove erfusion and drug delivery in solid tumors, PANS, vol. 110, no. 46,. 18632 18637, 2013. [21] M. Welter and H. Rieger, Interstitial fluid flow and drug delivery in vascularized tumors: a comutational model., PoS One, vol. 8, no. 8,. e70395, Jan. 2013. [22]. Q. Cai Y, Xu S, Wu J, Y. Cai, S. Xu, J. Wu, and Q. ong, Couled modelling of tumour angiogenesis, tumour growth and blood erfusion., J. Theor. iol., vol. 279, no. 1,. 90 101, Jun. 2011. [23] M. Soltani and P. Chen, Effect of tumor shae and size on drug delivery to solid tumors, J. iol. Eng., vol. 6, no. 4,. 4, 2012. [24] M. Sefidgar, M. Soltani, K. Raahemifar, H. azmara, S. Nayinian, and M. azargan, Effect of tumor shae, size, and tissue transort roerties on drug delivery to solid tumors, J. iol. Eng., vol. 8, no. 1,. 12, 2014. [25] M. Soltani and P. Chen, Numerical Modeling of Fluid Flow in Solid Tumors, PoS One, vol. 6, no. 6,. 1 15, 2011. [26] M. Soltani and P. Chen, Numerical Modeling of Interstitial Fluid Flow Couled with lood Flow through a Remodeled Solid Tumor Microvascular Network., PoS One, vol. 8, no. 6,. e67025, 2013. [27] M. Soltani, Numerical Modeling of Drug Delivery to Solid Tumor Microvasculature, PhD thesis, Chem. Eng. (Nanotechnology), Waterloo, Ontario, Canada, 2012. [28] M. Sefidgar, M. Soltani, H. azmara, M. Mousavi, M. azargan, and A. Elkamel, Interstitial Flow in Cancerous Tissue : Effect of Considering Remodeled Caillary Network, J. Tissue Sci. Eng., vol. 4, no. 3,. 1 8, 2014. [29] M. Sefidgar, K. Raahemifar, H. azmara, M. azargan, S. M.
Mousavi, and M. Soltani, Effect of remodeled tumor-induced caillary network on interstitial flow in cancerous tissue, in 2nd Middle East Conference on iomedical Engineering, 2014,. 212 215. [30] M. Soltani, H. azmara, M. Sefidgar, R. Subramaniam, and A. Rahmim, SU-D-201-04: Study On the Imact of Tumor Shae and Size On Drug Delivery to Pancreatic Tumors, Med. Phys., vol. 42, no. 6,. 3220, 2015. [31] M. Soltani, M. Sefidgar, H. azmara, R. Subramaniam, and A. Rahmim, Effect of tumor shae and size on drug delivery, J. Nucl. Med., vol. 56, no. sulement 3,. 1220, May 2015. [32] M. Sefidgar, M. Soltani, K. Raahemifar, M. Sadeghi, H. azmara, M. azargan, and M. Mousavi Naeenian, Numerical modeling of drug delivery in a dynamic solid tumor microvasculature, Microvasc. Res., vol. 99, no. 0,. 43 56, May 2015. [33] M. Soltani, M. Sefidgar, H. azmara, S. Sheikhbahaei, C. Marcus, S. Ashrafinia, R. Subramaniam, and A. M. Rahmim, WE-A-204-07: Satiotemoral Distribution of the FDG PET Tracer in Solid Tumors: Contributions of Diffusion and Convection Mechanisms, Med. Phys., vol. 42, no. 6,. 3660, 2015. [34] M. Soltani, M. Sefidgar, H. azmara, and A. Rahmim, Enhanced modeling of satiotemoral distribution of PET tracers in solid tumors and estimation of transort arameters, J. Nucl. Med., vol. 56, no. sulement 3,. 1221, May 2015. [35] M. Sefidgar, M. Soltani, K. Raahemifar, and H. azmara, Effect of Fluid Friction on Interstitial Fluid Flow Couled with lood Flow through Solid Tumor Microvascular Network, Comut. Math. Methods Med., vol. 2015, no. Article ID 673426, 2015. [36] H. azmara, M. Soltani, M. Sefidgar, and M. azargan, The ital Role of lood Flow-Induced Proliferation and Migration in Caillary Network Formation in a Multiscale Model of Angiogenesis, PoS One, vol. June 5,. 1 31, 2015. [37] K. afai, POROUS Media, Alication in iological systems and biotechnology. Taylor and Francis, 2011. [38] R. K. Jain, R. T. Tong, and.. Munn, Effect of vascular normalization by antiangiogenic theray on interstitial hyertension, eritumor edema, and lymhatic metastasis: insights from a mathematical model., Cancer Res., vol. 67, no. 6,. 2729 2735, 2007. [39] A. Scheidegger, The hysics of flow through orous media. University of Toronto Press., 1963. [40] E. H. Starling, On the absortion of uids from the connective tissue sace, J. hisilogy, vol. 19,. 312 326, 1986. [41] K. R. Olson, Transvascular and intravascular fluid transort in the rainbow trout: revisiting Starling s forces, the secondary circulation and interstitial comliance, J. Ex. iol., vol. 206, no. 3,. 457 467, Feb. 2003. [42] J. T. Tamsma, H. J. Keizer, and A. E. Meinders, Review Pathogenesis of malignant ascites : Starling s law of caillary hemodynamics revisited,. 1353 1357, 2001. [43]. T. axter and R. K. Jain, Transort of fluid and macromolecules in tumors. (I) role of interstitial ressure and convection, Microvasc. Res., vol. 37,. 77 104, 1989. [44] E. A. Swabb, J. Wei, and P. M. Gullino, Diffusion and Convection in Normal and Neolastic Tissues, Cancer Res., vol. 34, no. 10,. 2814 2822, Oct. 1974. [45] X. Zhao and M. N. Toksoz, Solute Transort In Heterogeneous Porous Media, Massachusetts Inst. Technol. Earth Resour. ab., 1994. [46]. Wu,. Gao, Y. Tian, and R. Muñoz-Carena, Analytical and exerimental analysis of solute transort in heterogeneous orous media, J. Environ. Sci. Heal. Part A, vol. 49, no. 3,. 338 343, Nov. 2013. [47] A. J. eu, D. A. erk, A. ymboussaki, K. Alitalo, and R. K. Jain, Absence of Functional ymhatics within a Murine Sarcoma: A Molecular and Functional Evaluation, Cancer Res., vol. 60, no. 16,. 4324 4327, Aug. 2000. [48] R. K. Jain and. T. Fenton, Intratumoral ymhatic essels: A Case of Mistaken Identity or Malfunction?, J. Natl. Cancer Inst., vol. 94, no. 6,. 417 421, Mar. 2002. [49] C. S. Patlak, D. A. Goldstein, and J. F. Hoffman, The flow of solute and solvent across a two-membrane system, J. Theor. iol., vol. 5, no. 3,. 426 442, Nov. 1963. [50] M. Soltani, M. Sefidgar, M. E. Casey, R.. Wahl, R. M. Subramaniam, and A. Rahmim, Comrehensive Modeling of the Satiotemoral Distribution of PET Tracer Utake in Solid Tumors based on the Convection-Diffusion-Reaction Equation, in 21th IEEE nuclear science symosium & medical imaging conference, 2014.