Research Papers

Effect of Tumor Volume on Drug Delivery in Heterogeneous Vasculature of Human Brain Tumors

[+] Author and Article Information
Ajay Bhandari

Department of Mechanical Engineering,
Indian Institute of Technology,
Kanpur 208016, India
e-mail: ajayb@iitk.ac.in

Ankit Bansal

Department of Mechanical and
Industrial Engineering,
Indian Institute of Technology,
Roorkee 247677, India
e-mail: abansfme@iitr.ac.in

Rishav Jain

Department of Mechanical Engineering,
Indian Institute of Technology,
Kanpur 208016, India;
Engineering and Design Section,
Mitsubishi Heavy Industries
Compressor Corporation,
Hiroshima 7330036, Japan
e-mail: rishav_jain@compressor.mhi.co.jp

Anup Singh

Centre for Biomedical Engineering,
Indian Institute of Technology,
Delhi 110016, India;
Department of Biomedical Engineering,
All Indian Institute of Medical Sciences,
Delhi 110016, India
e-mail: anupsm@cbme.iitd.ac.in

Niraj Sinha

Department of Mechanical Engineering,
Indian Institute of Technology,
Kanpur 208016, India
e-mail: nsinha@iitk.ac.in

1Corresponding author.

Manuscript received July 10, 2018; final manuscript received November 29, 2018; published online January 18, 2019. Assoc. Editor: Ali Sadegh.

ASME J of Medical Diagnostics 2(2), 021004 (Jan 18, 2019) (10 pages) Paper No: JESMDT-18-1031; doi: 10.1115/1.4042195 History: Received July 10, 2018; Revised November 29, 2018

Drug distribution in tumors is strongly dependent on tumor biological properties such as tumor volume, vasculature, and porosity. An understanding of the drug distribution pattern in tumors can help in enhancing the effectiveness of anticancer treatment. A numerical model is employed to study the distribution of contrast agent in the heterogeneous vasculature of human brain tumors of different volumes. Dynamic contrast enhanced-magnetic resonance imaging (DCE-MRI) has been done for a number of patients with different tumor volumes. Leaky tracer kinetic model (LTKM) is employed to obtain perfusion parameters from the DCE-MRI data. These parameters are used as input in the computational fluid dynamics (CFD) model to predict interstitial fluid pressure (IFP), interstitial fluid velocity (IFV), and distribution of the contrast agent in different tumors. Numerical results demonstrate that the IFP is independent of tumor volume. On the other hand, the IFV increases as the tumor volume increases. Further, the concentration of contrast agent also increases with the tumor volume. The results obtained in this work are in line with the experimental DCE-MRI data. It is observed that large volume tumors tend to retain a higher concentration of contrast agent for a longer duration of time because of large extravasation flux and slow washout as compared to smaller tumors. These results may be qualitatively extrapolated to chemotherapeutic drug delivery, implying faster healing in large volume tumors. This study helps in understanding the effect of tumor volume on the treatment outcome for a wide range of human tumors.

Copyright © 2019 by ASME
Your Session has timed out. Please sign back in to continue.


Siegel, R. , and Deepa Naishadham, A. J. , 2013, “ Cancer Statistics, 2013,” CA-Cancer J. Clin., 63(1), pp. 11–30. [CrossRef] [PubMed]
Jain, R. K. , 2005, “ Normalization of Tumor Vasculature: An Emerging Concept in Antiangiogenic Therapy,” Science, 307(5706), pp. 58–62. [CrossRef] [PubMed]
Jain, R. K. , 1996, “ Delivery of Molecular Medicine to Solid Tumors,” Science, 271(5252), pp. 1079–1080. [CrossRef] [PubMed]
Gray, L. H. , 1953, “ The Concentration of Oxygen Dissolved in Tissues at the Time of Irradiation as a Factor of Radiotherapy,” Br. J. Radiol., 26(312), pp. 638–648. [CrossRef] [PubMed]
Baxter, L. T. , and Jain, R. K. , 1989, “ Transport of Fluid and Macromolecules in Tumors—Part I: Role of Interstitial Pressure and Convection,” Microvasc. Res., 37(1), pp. 77–104. [CrossRef] [PubMed]
Baxter, L. T. , and Jain, R. K. , 1990, “ Transport of Fluid and Macromolecules in Tumors—Part II: Role of Heterogeneous Perfusion and Lymphatics,” Microvasc. Res., 40(2), pp. 246–263. [CrossRef] [PubMed]
Goh, Y. M. F. , Kong, H. L. , and Wang, C. H. , 2001, “ Simulation of the Delivery of Doxorubicin to Hepatoma,” Pharm. Res., 18(6), pp. 761–770. [CrossRef] [PubMed]
Arifin, D. Y. , Lee, K. Y. T. , and Wang, C. H. , 2009, “ Chemotherapeutic Drug Transport to Brain Tumor,” J. Control. Release, 137(3), pp. 203–210. [CrossRef] [PubMed]
Pishko, G. R. , Astary, G. W. , Mareci, T. H. , and Sarntinoranont, M. , 2011, “ Sensitivity Analysis of An Image-Based Solid Tumor Computational Model With Heterogeneous Vasculature and Porosity,” Ann. Biomed. Eng., 39(9), pp. 2360–2373. [CrossRef] [PubMed]
Magdoom, K. N. , Pishko, G. L. , and Kim, J. H. , 2015, “ Evaluation of a Voxelized Model Based on DCE-MRI for Tracer Transport in Tumor,” ASME J. Biomech. Eng., 134(9), p. 091004.
Soltani, M. , and Chen, P. , 2012, “ Effect of Tumor Shape and Size on Drug Delivery to Solid Tumors,” J. Biol. Eng., 6, pp. 1–15. [CrossRef] [PubMed]
Soltani, M. , and Chen, P. , 2014, “ Effect of Tumor Shape, Size, and Tissue Transport Properties on Drug Delivery to Solid Tumors,” J. Biol. Eng., 8(12), pp. 1–13. [PubMed]
Moghadam, M. C. , Deyranlou, A. , Sharifi, A. , and Niazmand, H. , 2015, “ Numerical Simulation of the Tumor Interstitial Fluid Transport: Consideration of Drug Delivery Mechanism,” Microvasc. Res., 101, pp. 62–71. [CrossRef] [PubMed]
Zhan, W. , Gedroyc, W. , and Xu, X. Y. , 2017, “ The Effect of Tumour Size on Drug Transport and Uptake in 3-D Tumour Models Reconstructed From Magnetic Resonance Images,” PLoS One, 12(2), p. e0172276
Tofts, P. S. , Brix, G. , Buckley, D. L. , L., Evelhoch , J., Henderson , E., Knopp , M. V., Larsson , H. B. W., Lee , T.-Y., Mayr , N. A., Parker , G. J. M., Port , R. E., Taylor , J. , and Weisskoff, R. M. , 1999, “ Estimating Kinetic Parameters From Dynamic Contrast-Enhanced T1-Weighted\tMRI of a Diffusable Tracer: Standardized Quantities and Symbols,” J. Magn. Reson. Imaging, 10(3), pp. 223–232. [CrossRef] [PubMed]
Sahoo, P. , Rathore, R. K. S. , Awasthi, R. , Roy, B. , Verma, S. , Rathore, D. , Behari, S. , Husain, M. , Husain, N. , Pandey, C. M. , Mohakud, S. , and Gupta, R. K. , 2012, “ Subcompartmentalization of Extracellular Extravascular Space (EES) Into Permeability and Leaky Space With Local Arterial Input Function (AIF) Results in Improved Discrimination Between High- and Low-Grade Glioma Using Dynamic Contrast-Enhanced (DCE) MRI,” J. Magn. Reson. Imaging, 38(3), pp. 677–688.
Singh, A. , Haris, M. , Rathore, D. , Purwar, A. , Sarma, M. , Bayu, G. , Husain, N. , Rathore, R. K. S. , and Gupta, R. K. , 2007, “ Quantification of Physiological and Hemodynamic Indices Using T(1) Dynamic Contrast-Enhanced MRI in Intracranial Mass Lesions,” J. Magn. Reson. Imaging, 26(4), pp. 871–880. [CrossRef] [PubMed]
Pintaske, J. , Martirosian, P. , Graf, H. , Erb, G. , Lodemann, K.-P. , Claussen, C. D. , and Schick, F. , 2006, “ Relaxivity of Gadopentetate Dimeglumine (Magnevist), Gadobutrol (Gadovist), and Gadobenate Dimeglumine (MultiHance) in Human Blood Plasma at 0.2, 1.5, and 3 Tesla,” Invest. Radiol., 41(3), pp. 213–221. [CrossRef] [PubMed]
Singh, A. , Rathore, R. K. S. , Haris, M. , Verma, S. K. , Husain, N. , and Gupta, R. K. , 2009, “ Improved Bolus Arrival Time and Arterial Input Function Estimation for Tracer Kinetic Analysis in DCE-MRI,” J. Magn. Reson. Imaging, 29(1), pp. 166–176. [CrossRef] [PubMed]
Tofts, P. S. , 1997, “ Modeling Tracer Kinetics in Dynamic GD-DTPA MR Imaging,” J. Magn. Reson. Imaging, 7(1), pp. 91–101. [CrossRef] [PubMed]
Bhandari, A. , Bansal, A. , Singh, A. , and Sinha, N. , 2017, “ Study of Perfusion Kinetics in Human Brain Tumor Using Leaky Tracer Kinetic Model of DCE-MRI Data and CFD,” Communications in Computer and Information Science, M. Fei , S. Ma , X. Li , X. Sun , L. Jia , and Z. Su , eds., Springer, Singapore, pp. 63–73.
Jessen, N. A. , Munk, A. S. F. , Lundgaard, I. , and Nedergaard, M. , 2015, “ The Glymphatic System: A Beginner's Guide,” Neurochem. Res., 40(12), pp. 2583–2599. [CrossRef] [PubMed]
Paulson, O. B. , Hertz, M. M. , Bolwig, T. G. , and Lassen, N. A. , 1977, “ Filtration and Diffusion of Water Across the Blood-Brain Barrier in Man,” Microvasc. Res., 13(1), pp. 113–123. [CrossRef] [PubMed]
Smith, J. H. , and Humphrey, J. A. C. , 2007, “ Interstitial Transport and Transvascular Fluid Exchange During Infusion Into Brain and Tumor Tissue,” Microvasc. Res., 73(1), pp. 58–73. [CrossRef] [PubMed]
Bradbury, M. , 1979, The Concept of Blood Brain Barrier, Wiley, Chichester, UK.
Wang, C. H. , Li, J. , Teo, C. S. , and Lee, T. , 1999, “ The Delivery of BCNU to Brain Tumors,” J. Control. Release, 61(1–2), pp. 21–41. [CrossRef] [PubMed]
Swabb, E. A. , Wei, J. , Gullino, P. M. , Swabb, E. A. , Wei, J. , and Cullino, P. M. , 1974, “ Diffusion and Convection in Normal and Neoplastic Tissues,” Cancer Res., 34(10), pp. 2814–2822. http://cancerres.aacrjournals.org/content/34/10/2814.short [PubMed]
Netti, P. A. , Berk, D. A. , Swartz, M. A. , Grodzinsky, A. J. , and Jain, R. K. , 2000, “ Role of Extracellular Matrix Assembly in Interstitial Transport in Solid Tumors,” Cancer Res., 60(9), pp. 2497–2503. https://www.ncbi.nlm.nih.gov/pubmed/10811131
Zhan, W. , Gedroyc, W. , and Xu, X. Y. , 2014, “ Mathematical Modelling of Drug Transport and Uptake in a Realistic Model of Solid Tumour,” Protein Pept. Lett., 21(11), pp. 1146–1156. [CrossRef] [PubMed]
Gerlowski, L. E. , and Jain, R. K. , 1986, “ Microvascular Permeability of Normal and Neoplastic Tissues,” Microvasc. Res., 31(3), pp. 288–305. [CrossRef] [PubMed]
Anderson, D. A. , Tannehill, J. C. , and Pletcher, R. H. , 1984, Computational Fluid Mechanics and Heat Transfer, Hemisphere, New York, pp. 671–674.
Abe, T. , Mizobuchi, Y. , Nakajima, K. , Otomi, Y. , Irahara, S. , Obama, Y. , Majigsuren, M. , Khashbat, D. , Kageji, T. , Nagahiro, S. , and Harada, M. , 2015, “ Diagnosis of Brain Tumors Using Dynamic Contrast-Enhanced Perfusion Imaging With a Short Acquisition Time,” Springerplus, 4, p. 88. [CrossRef] [PubMed]
Bhandari, A. , Bansal, A. , Singh, A. , and Sinha, N. , 2017, “ Perfusion Kinetics in Human Brain Tumor With DCE-MRI Derived Model and CFD Analysis,” J. Biomech., 59, pp. 80–89. [CrossRef] [PubMed]
Gutmann, R. , Leunig, M. , Feyh, J. , Size, T. , Goetz, A. E. , Messmer, K. , Kastenbauer, E. , and Jain, R. K. , 1992, “ Interstitial Hypertension in Head and Neck Tumors in Patients: Correlation With Tumor Size,” Cancer Res., 52(7), pp. 1993–1995. http://cancerres.aacrjournals.org/content/52/7/1993.short [PubMed]
Jain, R. K. , 2012, “ Delivery of Molecular and Cellular Medicine to Solid Tumors,” Adv. Drug Delivery Rev., 64(Suppl.), pp. 353–365. [CrossRef]
DiResta, G. R. , Lee, J. , Larson, S. M. , and Arbit, E. , 1993, “ Characterization of Neuroblastoma Xenograft in Rat Flank—Part I: Growth, Interstitial Fluid Pressure, and Interstitial Fluid Velocity Distribution Profiles,” Microvasc. Res., 46(2), pp. 158–177. [CrossRef] [PubMed]
Pishko, G. L. , Astary, G. W. , Zhang, J. , Mareci, T. H. , and Sarntinoranont, M. , 2012, “ Role of Convection and Diffusion on DCE-MRI Parameters in Low Leakiness KHT Sarcomas,” Microvasc. Res., 84(3), pp. 306–313. [CrossRef] [PubMed]
Bhandari, A. , Bansal, A. , Singh, A. , and Sinha, N. , 2017, “ Numerical Study of Transport of Anti-Cancer Drugs in Heterogeneous Vasculature of Human Brain Tumors Using DCE-MRI,” ASME J. Biomech. Eng., 140(5), p. 051010. [CrossRef]
Barnett, P. A. , Roman-Goldstein, S. , Ramsey, F. , McCormick, C. I. , Sexton, G. , Szumowski, J. , and Neuwelt, E. A. , 1995, “ Differential Permeability and Quantitative MR Imaging of a Human Lung Carcinoma Brain Xenograft in the Nude Rat,” Am. J. Pathol., 146(2), pp. 436–449. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1869863/ [PubMed]
Zhang, R. D. , Price, J. E. , Fujimaki, T. , Bucana, C. D. , and Fidler, I. J. , 1992, “ Differential Permeability of the Blood-Brain Barrier in Experimental Brain Metastases Produced by Human Neoplasms Implanted Into Nude Mice,” Am. J. Pathol., 141(5), pp. 1115–1124. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1886664/ [PubMed]
Greig, N. , Jones, H. , and Cavanagh, J. , 1983, “ Blood-Brain Barrier Integrity and Host Responses in Experimental Metastatic Brain Tumours,” Clin. Exp. Metastasis, 1(3), pp. 229–246. [CrossRef] [PubMed]
Jain, R. , Gutierrez, J. , Narang, J. , Scarpace, L. , Schultz, L. R. , Lemke, N. , Patel, S. C. , Mikkelsen, T. , and Rock, J. P. , 2011, “ In Vivo Correlation of Tumor Blood Volume and Permeability With Histologic and Molecular Angiogenic Markers in Gliomas,” Am. J. Neuroradiol., 32(2), pp. 388–394. [CrossRef]
Jain, R. K. , Tong, R. T. , and Munn, L. L. , 2007, “ Effect of Vascular Normalization by Antiangiogenic Therapy on Interstitial Hypertension, Peritumor Edema, and Lymphatic Metastasis: Insights From a Mathematical Model,” Cancer Res., 67(6), pp. 2729–2735. http://cancerres.aacrjournals.org/content/67/6/2729.short [PubMed]


Grahic Jump Location
Fig. 2

(a) Rectangular computational domain in CFD and (b) segmented single slice enclosing tumor (shaded) and normal tissue (all dimensions in mm)

Grahic Jump Location
Fig. 3

Mesh independence study of one slice of tumor shown in Fig. 2: (a) IFP, (b) IFV, and (c) contrast agent concentration at 60 s

Grahic Jump Location
Fig. 4

(a) Precontrast and (b) Postcontrast image of MRI slice (outline shows manual tracing of tumor). Units of x-axis and y-axis are in mm.

Grahic Jump Location
Fig. 5

Perfusion parameters (a) permeability (Ktrans) (s−1) and (b) porosity of tumor shown in Fig. 2. Units of x-axis and y-axis are in mm.

Grahic Jump Location
Fig. 6

Contour maps of (a) IFP and (b) IFV of tumor shown in Fig. 2. Units of x-axis and y-axis are in mm. Unit of color scale bar for IFP is in Pa and for IFV is in m/s.

Grahic Jump Location
Fig. 7

Variation of IFP and IFV with tumor volume

Grahic Jump Location
Fig. 8

Contour maps of the time-averaged concentration of contrast agent of a particular slice of tumor (a) experimental and (b) simulated. Unit of the color scale bar for concentration is mmol/L. Units of x and y axis are same as shown in Fig. 2. (c) Averaged percentage error between experimental and simulated results.

Grahic Jump Location
Fig. 9

Extravasation flux pattern with time at one tumor voxel shown in the zoomed contour plot

Grahic Jump Location
Fig. 10

Variation of contrast agent concentration with tumor volume

Grahic Jump Location
Fig. 11

Variation of permeability (Ktrans) and extravasation flux with tumor volume

Grahic Jump Location
Fig. 12

Variation of plasma volume fraction with tumor volume

Grahic Jump Location
Fig. 13

Variation of the total contrast agent concentration for all the voxels per unit volume with time

Grahic Jump Location
Fig. 14

Variation of convective flux and diffusive flux with tumor volume



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In