Elsevier

Journal of Biomechanics

Volume 59, 5 July 2017, Pages 80-89
Journal of Biomechanics

Perfusion kinetics in human brain tumor with DCE-MRI derived model and CFD analysis

https://doi.org/10.1016/j.jbiomech.2017.05.017Get rights and content

Abstract

Cancer is one of the leading causes of death all over the world. Among the strategies that are used for cancer treatment, the effectiveness of chemotherapy is often hindered by factors such as irregular and non-uniform uptake of drugs inside tumor. Thus, accurate prediction of drug transport and deposition inside tumor is crucial for increasing the effectiveness of chemotherapeutic treatment. In this study, a computational model of human brain tumor is developed that incorporates dynamic contrast enhanced-magnetic resonance imaging (DCE-MRI) data into a voxelized porous media model. The model takes into account realistic transport and perfusion kinetics parameters together with realistic heterogeneous tumor vasculature and accurate arterial input function (AIF), which makes it patient specific. The computational results for interstitial fluid pressure (IFP), interstitial fluid velocity (IFV) and tracer concentration show good agreement with the experimental results. The computational model can be extended further for predicting the deposition of chemotherapeutic drugs in tumor environment as well as selection of the best chemotherapeutic drug for a specific patient.

Introduction

Tumors are abnormal mass of tissues, which usually depend on other normal tissues for their nutritional needs. Tumors contain highly tortuous, fenestrated, discontinuous vessels and large vascular areas, making them highly heterogeneous (Jain, 1987). The pathophysiological state of tumor comprises of accumulated solid stress, abnormal and heterogeneous blood vessel networks, elevated IFP and a dense interstitial structure. This leads to transport barriers such as irregular perfusion, inadequate drug delivery, sporadic and uneven uptake. Further, these barriers cause heterogeneous extravasation of therapeutic agents into tumor tissues that may significantly limit the rate and extent of drug delivery to tumor (Jain et al., 2014, Narang and Varia, 2011). Since transport of macromolecular therapeutic agents in tumor microvasculature is vital for effective treatment of solid tumors, there is an imperative need to develop non-invasive approaches to understand and prevent various transport obstacles.

A number of mathematical models have been developed to study transport of drugs in tumor. Baxter and Jain used a homogeneous porous media model of both uniformly and non-uniformly perfused tumor to solve for interstitial fluid flow and solute transport parameters (Baxter and Jain, 1989, Baxter and Jain, 1990). Eikenberry used a tumor cord model to simulate the delivery of a chemotherapeutic drug to a solid tumor (Eikenberry, 2009). Pozrikidis modelled the tumor capillaries and evaluated blood and interstitial flow (Pozrikidis, 2010). Recently, computational fluid dynamics (CFD) approaches are increasingly being used to predict transport of macromolecules within tissues. For a homogenous tumor model, Soltani and Chen observed that IFP is less than the effective pressure below a certain tumor radius, and transport of chemotherapeutic drug to tumor site is easier (Soltani and Chen, 2011). Zhan et al. modelled heterogeneous vasculature of tumor and concluded that the drug accumulates faster in well vascularized regions and clears out from these regions quickly, resulting in less tumor cell killing (Zhan et al., 2014a). Pishko et al. modelled the tracer transport through a heterogeneous tumor tissue by using a non-voxelized approach (Pishko et al., 2011). Later, Magdoom et al. improved the model by using a voxelized approach that resulted in reduced computational time and increased accuracy (Magdoom et al., 2012).

However, existing studies have a number of limitations. Many of them assumed spherical shape and homogeneous vasculature of tumor, which is far from reality. Even studies employing heterogeneous vasculature of tumor have mainly focused on animal models, with Simple Tofts Kinetics Model used for the analysis of MRI data and the intravascular term ignored. Further, these models did not employ the patient specific AIF and AIF values were taken from literature (also called global AIF). This resulted in inaccurate determination of perfusion kinetic parameters and tracer concentration.

The objective of the present paper is to overcome these limitations by developing a computational model for the transport of contrast agent (tracer) in realistic human brain tumors. Realistic heterogeneous vasculature, permeability and porosity maps of the tumor are determined by analysing the DCE-MRI data. The computational model is employed to model tracer concentration along with IFP and IFV. Similar approach has been used by Magdoom et al. to model interstitial transport in mice. However, the current study introduces realistic heterogeneous human brain tumor and employs patient specific AIF to determine accurate perfusion kinetic parameters with the use of General Tracer Kinetic Model (GTKM). To the best of our knowledge, no study has reported the CFD analysis of realistic human brain tumors based on DCE-MRI data and patient specific AIF.

Section snippets

Methods

In this study, GTKM also called the Extended Tofts Model (Tofts, 1997, Tofts and Parker, 2013) has been used for analyzing DCE-MRI data. The obtained permeability and porosity maps were imported into the computational model developed in OpenFOAM, an open source CFD code, to solve for fluid flow parameters and tracer transport. GTKM model takes the blood volume fraction or intravascular term into account and is known to be more accurate for highly perfused tumors that is the case for humans due

Magnetic resonance imaging (MRI) protocol

DCE-MR imaging was performed on a 1.5-T GE scanner (Signa, Lx Echospeed plus; General Electric Medical Systems, Milwaukee, WI, USA). DCE-MR imaging was done for four patients. Written consent from each patient was obtained before MRI study. Out of those four patients, MR imaging of one patient was done for 14 min, whereas for others it was done for 4 min. DCE imaging was performed using a 3D-SPGR sequence (TR/TE = 5.0 ms/1.4 ms, flip angle = 15°, field of view (FOV) = 240 × 240 mm2, slice thickness = 6 mm,

Tissue transport model

The computational model for the transport of fluid and tracer through tissues is based on flow through porous media. Each point in this porous media transport model consists of vascular and tissue compartment. The governing equations of transport are mentioned in Eqs. (3), (4), (5). IFP is calculated by the continuity equation with added source and sink terms and is expressed as·V=KtransKavgtransLpSV(Peff-Pi)-Lp,lySlV(Pi-Pl)where Peff=Pv-σT(πv-πi)

The first term on right hand side of Eq. (3)

Computational method

To reduce the computation time, only the tumor part and the surrounding normal tissue were taken as the computational domain. A rectangular volume of size 32 × 24 × 72 mm3 was created and meshed in a structured Cartesian grid as shown in Fig. 1. The voxel size was taken in accordance with the MRI resolution (0.9375 × 0.9375 × 6 mm3), maintaining complete analogy between MRI data and CFD model. Tracer kinetic parameters (Ktrans,ve) and transport properties present in Eqs. (3), (4), (5) of the tissue were

Results and discussion

Fig. 2 shows the pre-contrast and post-contrast images of brain of one slice at two different times. Local AIF (Fig. 2(d)) has been used for the patient in this study as it leads to calculation of accurate kinetic parameters. Permeability and porosity maps obtained from the MRI data are shown in Fig. 3. The porosity values highly depend on the type and grade of tumors. A number of studies have reported that high grade tumors have higher values of porosity as compared to low grade tumors (

Conclusion

A computational model based on MR imaging and governing equations of fluid mechanics was developed to study fluid flow and transport of contrast agent in three-dimensional heterogeneous realistic human brain tumor. The heterogeneous vasculature was determined on voxel scale by processing DCE-MRI data. A good agreement was found between simulated contrast agent concentration and those obtained using the experimental MRI data. Also, simulated IFP and IFV values correlated well with the

Conflict of interest statement

The authors declare that they have no competing interest.

Acknowledgements

The authors thank Dr. R.K. Gupta for providing clinical data and Prof. R.K.S. Rathore for technical support in DCE-MRI data analysis. This research was supported by grants from IIT Kanpur and Science and Engineering Research Board (Grant number: YSS/2014/000092).

References (36)

  • M. Bradbury

    The Concept of a Blood Brain Barrier

    (1979)
  • Y. Boucher et al.

    Interstitial fluid pressure in intracranial tumors in patients and in rodents

    Br. J. Cancer

    (1997)
  • S. Eikenberry

    A tumor cord model for doxorubicin delivery and dose optimization in solid tumors

    Theor. Biol. Med. Model.

    (2009)
  • Y.M.F. Goh et al.

    Simulation of delivery of doxorubicin to hepatoma

    Pharm. Res.

    (2001)
  • R. Guttman et al.

    Interstital hypertension in head and neck tumors in patients: correlation with tumor size

    Can. Res.

    (1992)
  • R.K. Jain

    Transport of molecules in the tumor interstitium: a review

    Can. Res.

    (1987)
  • R.K. Jain et al.

    The role of mechanical forces in tumor growth and therapy

    Annu. Rev. Biomed. Eng.

    (2014)
  • N.A. Jessen et al.

    The glymphatic system: a beginner’s guide

    Neurochem. Res.

    (2015)
  • Cited by (30)

    • Interstitial fluid pressure as an emerging biomarker in solid tumors

      2022, Biochimica et Biophysica Acta - Reviews on Cancer
    • Simulation of magnetic nanoparticles crossing through a simplified blood-brain barrier model for Glioblastoma multiforme treatment

      2021, Computer Methods and Programs in Biomedicine
      Citation Excerpt :

      Drochon [36] presented an analytical solution for the blood flow in a cylindrical vessel with an external magnetic force, indicating that the magnetic field influences the fluid velocity profile and the wall displacements. Bhandari et al. [37] developed a computational model of the human brain tumor to investigate the perfusion kinetics using CFD methods, while Tehrani et al. [38] presented a novel numerical scheme for nanoparticles steering by performing simulations on both Y-shaped channels and realistic three-dimensional (3-D) model channels extracted from brain vessels. The magnetically targeted drug delivery to human parts has been investigated using computational methods.

    • Convection enhanced delivery of anti-angiogenic and cytotoxic agents in combination therapy against brain tumour

      2020, European Journal of Pharmaceutical Sciences
      Citation Excerpt :

      Parameters adopted in this work correspond to the averaged and representative values from literatures, as given in Tables 1 and 2. This assumption can be relaxed by using dynamic contrast enhanced MR images (Zhan et al., 2014; Bhandari et al., 2017), on which the signal intensity correlates to the concentration of imaging tracer. By applying the drug transvascular transport equation (Zhao et al., 2007), the microvascular density can be calculated from spatiotemporal profiles of the tracer. (

    View all citing articles on Scopus
    View full text