Spatiotemporal distribution modeling of PET tracer uptake in solid tumors

Radiopharmaceutical transport in solid tumors via a 3-dimensional image-based spatiotemporal model

Lutetium-177 prostate-specific membrane antigen (177Lu-PSMA)-targeted radiopharmaceutical therapy is a clinically approved treatment for patients with metastatic castration-resistant prostate cancer (mCRPC). Even though common practice reluctantly follows "one size fits all" approach, medical community believes there is significant room for deeper understanding and personalization of radiopharmaceutical therapies. To pursue this aim, we present a 3-dimensional spatiotemporal radiopharmaceutical delivery model based on clinical imaging data to simulate pharmacokinetic of 177Lu-PSMA within the prostate tumors. The model includes interstitial flow, radiopharmaceutical transport in tissues, receptor cycles, association/dissociation with ligands, synthesis of PSMA receptors, receptor recycling, internalization of radiopharmaceuticals, and degradation of receptors and drugs. The model was studied for a range of values for injection amount (100-1000 nmol), receptor density (10-500 nmol•l-1), and recycling rate of receptors (10-4 to 10-1 min-1). Furthermore, injection type, different convection-diffusion-reaction mechanisms, characteristic time scales, and length scales are discussed. The study found that increasing receptor density, ligand amount, and labeled ligands improved radiopharmaceutical uptake in the tumor. A high receptor recycling rate (0.1 min-1) increased radiopharmaceutical concentration by promoting repeated binding to tumor cell receptors. Continuous infusion results in higher radiopharmaceutical concentrations within tumors compared to bolus administration. These insights are crucial for advancing targeted therapy for prostate cancer by understanding the mechanism of radiopharmaceutical distribution in tumors. Furthermore, measures of characteristic length and advection time scale were computed. The presented spatiotemporal tumor transport model can analyze different physiological parameters affecting 177Lu-PSMA delivery.

Spatiotemporal modeling of radiopharmaceutical transport in solid tumors: Application to 177Lu-PSMA therapy of prostate cancer

BACKGROUND AND OBJECTIVE

177Lu-labeled prostate-specific membrane antigen (PSMA) radiopharmaceutical therapy (RPT) represents a pivotal advancement in addressing prostate cancer. However, existing therapies, while promising, remain incompletely understood and optimized. Computational models offer potential insights into RPTs, aiding in clinical drug delivery enhancement. In this study, we investigate the impact of various physiological parameters on the delivery of 177Lu-PSMA-617 RPT using the convection-diffusion-reaction (CDR) model.

METHODS

Our investigation encompasses tumor geometry and surrounding tissue, characterized by well-defined boundaries and initial conditions. Utilizing the finite element method, we solve governing equations across a range of parameters: dissociation constant KD (1, 0.1, 0.01 [nM]), internalization rate (0.01-0.0001 [min-1]), diverse tumor shapes, and variable necrotic zone sizes. This model can provide an accurate analysis of radiopharmaceutical delivery from the injection site to the tumor cell, including drug transport in the vascular, interstitial, and intracellular spaces, and considering important parameters (e.g., drug extravasation from microvessels or to lymphatic vessels, the extracellular matrix, receptors, and intracellular space).

RESULTS

Our findings reveal significant enhancements in tumor-absorbed doses as KD decreases. This outcome can be attributed to the higher affinity of radiopharmaceuticals for PSMA receptors as KD diminishes, facilitating a more efficient binding and retention of the therapeutic agent within the tumor microenvironment. Additionally, tumor-absorbed doses for KD ∼ 1 [nM] show an upward trend with higher internalization rates. This observation can be rationalized by considering that a greater internalization rate would result in a higher proportion of radiopharmaceuticals being taken up by tumor cells after binding to receptors on the cell surface. Notably, tumor shape and necrotic zone size exhibit limited influence on tumor absorbed dose.

CONCLUSIONS

The present study employs the CDR model to explore the role of physiological parameters in shaping 177Lu-PSMA-617 RPT delivery. These findings provide insights for improving prostate cancer therapy by understanding radiopharmaceutical transport dynamics. This computational approach contributes to advancing our understanding of radiopharmaceutical delivery mechanisms and has implications for enhancing treatment efficacy.

Multi-Scale and Multi-Physics Models of the Transport of Therapeutic/Diagnostic Cancer Agents

The effectiveness of tumor treatment heavily relies on the successful delivery of anticancer drugs [...].

A spatiotemporal multi-scale computational model for FDG PET imaging at different stages of tumor growth and angiogenesis

A deeper understanding of the tumor microenvironment (TME) and its role in metabolic activity at different stages of vascularized tumors can provide useful insights into cancer progression and better support clinical assessments. In this study, a robust and comprehensive multi-scale computational model for spatiotemporal transport of F-18 fluorodeoxyglucose (FDG) is developed to incorporate important aspects of the TME, spanning subcellular-, cellular-, and tissue-level scales. Our mathematical model includes biophysiological details, such as radiopharmaceutical transport within interstitial space via convection and diffusion mechanisms, radiopharmaceutical exchange between intracellular and extracellular matrices by glucose transporters, cellular uptake of radiopharmaceutical, as well as its intracellular phosphorylation by the enzyme. Further, to examine the effects of tumor size by varying microvascular densities (MVDs) on FDG dynamics, four different capillary networks are generated by angiogenesis modeling. Results demonstrate that as tumor grows, its MVD increases, and hence, the spatiotemporal distribution of total FDG uptake by tumor tissue changes towards a more homogenous distribution. In addition, spatiotemporal distributions in tumor with lower MVD have relatively smaller magnitudes, due to the lower diffusion rate of FDG as well as lower local intravenous FDG release. Since mean standardized uptake value (SUV_mean) differs at various stages of microvascular networks with different tumor sizes, it may be meaningful to normalize the measured values by tumor size and the MVD prior to routine clinical reporting. Overall, the present framework has the potential for more accurate investigation of biological phenomena within TME towards personalized medicine.

Synthetic 18F-FDG PET Image Generation Using a Combination of Biomathematical Modeling and Machine Learning

No previous works have attempted to combine generative adversarial network (GAN) architectures and the biomathematical modeling of positron emission tomography (PET) radiotracer uptake in tumors to generate extra training samples. Here, we developed a novel computational model to produce synthetic 18F-fluorodeoxyglucose (18F-FDG) PET images of solid tumors in different stages of progression and angiogenesis. First, a comprehensive biomathematical model is employed for creating tumor-induced angiogenesis, intravascular and extravascular fluid flow, as well as modeling of the transport phenomena and reaction processes of 18F-FDG in a tumor microenvironment. Then, a deep convolutional GAN (DCGAN) model is employed for producing synthetic PET images using 170 input images of 18F-FDG uptake in each of 10 different tumor microvascular networks. The interstitial fluid parameters and spatiotemporal distribution of 18F-FDG uptake in tumor and healthy tissues have been compared against previously published numerical and experimental studies, indicating the accuracy of the model. The structural similarity index measure (SSIM) and peak signal-to-noise ratio (PSNR) of the generated PET sample and the experimental one are 0.72 and 28.53, respectively. Our results demonstrate that a combination of biomathematical modeling and GAN-based augmentation models provides a robust framework for the non-invasive and accurate generation of synthetic PET images of solid tumors in different stages.

Computational modeling of PET tracer distribution in solid tumors integrating microvasculature

BACKGROUND

We present computational modeling of positron emission tomography radiotracer uptake with consideration of blood flow and interstitial fluid flow, performing spatiotemporally-coupled modeling of uptake and integrating the microvasculature. In our mathematical modeling, the uptake of fluorodeoxyglucose F-18 (FDG) was simulated based on the Convection-Diffusion-Reaction equation given its high accuracy and reliability in modeling of transport phenomena. In the proposed model, blood flow and interstitial flow are solved simultaneously to calculate interstitial pressure and velocity distribution inside cancer and normal tissues. As a result, the spatiotemporal distribution of the FDG tracer is calculated based on velocity and pressure distributions in both kinds of tissues.

RESULTS

Interstitial pressure has maximum value in the tumor region compared to surrounding tissue. In addition, interstitial fluid velocity is extremely low in the entire computational domain indicating that convection can be neglected without effecting results noticeably. Furthermore, our results illustrate that the total concentration of FDG in the tumor region is an order of magnitude larger than in surrounding normal tissue, due to lack of functional lymphatic drainage system and also highly-permeable microvessels in tumors. The magnitude of the free tracer and metabolized (phosphorylated) radiotracer concentrations followed very different trends over the entire time period, regardless of tissue type (tumor vs. normal).

CONCLUSION

Our spatiotemporally-coupled modeling provides helpful tools towards improved understanding and quantification of in vivo preclinical and clinical studies.

Evaluation of solid tumor response to sequential treatment cycles via a new computational hybrid approach

The development of an in silico approach that evaluates and identifies appropriate treatment protocols for individuals could help grow personalized treatment and increase cancer patient lifespans. With this motivation, the present study introduces a novel approach for sequential treatment cycles based on simultaneously examining drug delivery, tumor growth, and chemotherapy efficacy. This model incorporates the physical conditions of tumor geometry, including tumor, capillary network, and normal tissue assuming real circumstances, as well as the intravascular and interstitial fluid flow, drug concentration, chemotherapy efficacy, and tumor recurrence. Three treatment approaches-maximum tolerated dose (MTD), metronomic chemotherapy (MC), and chemo-switching (CS)-as well as different chemotherapy schedules are investigated on a real tumor geometry extracted from image. Additionally, a sensitivity analysis of effective parameters of drug is carried out to evaluate the potential of using different other drugs in cancer treatment. The main findings are: (i) CS, MC, and MTD have the best performance in reducing tumor cells, respectively; (ii) multiple doses raise the efficacy of drugs that have slower clearance, higher diffusivity, and lower to medium binding affinities; (iii) the suggested approach to eradicating tumors is to reduce their cells to a predetermined rate through chemotherapy and then apply adjunct therapy.

Effects of hypoxia and nanocarrier size on pH-responsive nano-delivery system to solid tumors

One of the special features of solid tumors is the acidity of the tumor microenvironment, which is mainly due to the presence of hypoxic regions. Therefore, pH-responsive drug delivery systems have recently been highly welcomed. In the present study, a comprehensive mathematical model is presented based on extravascular drug release paradigm. Accordingly, drug delivery system using pH-responsive nanocarriers is taken into account to examine the impacts of hypoxic regions as well as the size of nanocarriers for cancerous cell-death. The extent of hypoxic regions is controlled by vascular density. This means that regions with very low vascular density represent regions of hypoxia. Using this mathematical model, it is possible to simulate the extracellular and intracellular concentrations of drug by considering the association/disassociation of the free drug to the cell-surface receptors and cellular uptake. Results show that nanocarriers with smaller sizes are more effective due to higher accumulation in the tumor tissue interstitium. The small size of the nanocarriers also allows them to penetrate deeper, so they can expose a larger portion of the tumor to the drug. Additionally, the presence of hypoxic regions in tumor reduces the fraction of killed cancer cells due to reduced penetration depth. The proposed model can be considered for optimizing and developing pH-sensitive delivery systems to reduce both cost and time of the process.

Mathematical simulation and prediction of tumor volume using RBF artificial neural network at different circumstances in the tumor microenvironment

Uncontrolled proliferation of cells in a tissue caused by genetic mutations inside a cell is referred to as a tumor. A tumor which grows rapidly encounters a barrier when it grows to a certain size in presence of preexisting vasculature. This is the time when it has to find a way to go on the growth. The tumor starts to secrete tumor angiogenic factors (TAFs) and stimulate preexisting vessels to grow new sprouts. These new sprouts will find their way to the tumor in the extracellular matrix (ECM) by the gradient of TAF. As these new capillaries anastomose and reach tumor, fresh oxygen is available for the tumor and it will reinitiate the growth. Number of initial sprouts, distance of initial tumor cells from the vessel(s) and initial density of the tumor at the time of sprout formation are questions which are to be investigated. In the present study, the aim is to find the response of tumor cells and vessels to the reciprocal effects of each other in different circumstances in the tissue. Together with a mathematical formulation, a radial basis function (RBF) neural network is established to predict the number of tumor cells at different circumstances including size and distance of initial tumors from the parent vessel. A final formulation is given for the final number of tumor cells as a function of initial tumor size and distance between a parent vessel and a tumor. Results of this simulation demonstrate that, increasing the distance between a tumor and a parent vessel decreases the number of final tumor cells. Specially, this decrement becomes faster beyond a certain distance. Moreover, initial tumors in bigger domains must become much bigger before inducing angiogenesis which makes it harder for them to survive.

Enhanced Drug Delivery to Solid Tumors via Drug-Loaded Nanocarriers: An Image-Based Computational Framework

Objective

Nano-sized drug delivery systems (NSDDSs) offer a promising therapeutic technology with sufficient biocompatibility, stability, and drug-loading rates towards efficient drug delivery to solid tumors. We aim to apply a multi-scale computational model for evaluating drug delivery to predict treatment efficacy.

Methodology

Three strategies for drug delivery, namely conventional chemotherapy (one-stage), as well as chemotherapy through two- and three-stage NSDDSs, were simulated and compared. A geometric model of the tumor and the capillary network was obtained by processing a real image. Subsequently, equations related to intravascular and interstitial flows as well as drug transport in tissue were solved by considering real conditions as well as details such as drug binding to cells and cellular uptake. Finally, the role of periodic treatments was investigated considering tumor recurrence between treatments. The impact of different parameters, nanoparticle (NP) size, binding affinity of drug, and the kinetics of release rate, were additionally investigated to determine their therapeutic efficacy.

Results

Using NPs considerably increases the fraction of killed cells (FKCs) inside the tumor compared to conventional chemotherapy. Tumoral FKCs for two-stage DDS with smaller NP size (20nm) is higher than that of larger NPs (100nm), in all investigate release rates. Slower and continuous release of the chemotherapeutic agents from NPs have better treatment outcomes in comparison with faster release rate. In three-stage DDS, for intermediate and higher binding affinities, it is desirable for the secondary particle to be released at a faster rate, and the drug with slower rate. In lower binding affinities, high release rates have better performance. Results also demonstrate that after 5 treatments with three-stage DDS, 99.6% of tumor cells (TCs) are killed, while two-stage DDS and conventional chemotherapy kill 95.6% and 88.5% of tumor cells in the same period, respectively.

Conclusion

The presented framework has the potential to enable decision making for new drugs via computational modeling of treatment responses and has the potential to aid oncologists with personalized treatment plans towards more optimal treatment outcomes.

Controlled anti-cancer drug release through advanced nano-drug delivery systems: Static and dynamic targeting strategies

Advances in nanomedicine, including early cancer detection, targeted drug delivery, and personalized approaches to cancer treatment are on the rise. For example, targeted drug delivery systems can improve intracellular delivery because of their multifunctionality. Novel endogenous-based and exogenous-based stimulus-responsive drug delivery systems have been proposed to prevent the cancer progression with proper drug delivery. To control effective dose loading and sustained release, targeted permeability and individual variability can now be described in more-complex ways, such as by combining internal and external stimuli. Despite these advances in release control, certain challenges remain and are identified in this research, which emphasizes the control of drug release and applications of nanoparticle-based drug delivery systems. Using a multiscale and multidisciplinary approach, this study investigates and analyzes drug delivery and release strategies in the nanoparticle-based treatment of cancer, both mathematically and clinically.

Drug delivery to solid tumors with heterogeneous microvascular networks: Novel insights from image-based numerical modeling

The present study examines chemotherapy by incorporating multi-scale mathematical modeling to predict drug delivery and its effects. This approach leads to a more-realistic physiological tumor model than is possible with previous approaches, as it obtains the capillary network geometry from an image, and also considers the tumor's necrotic core, drug binding, and cellular uptake. Modeling of the fluid flow and drug transport is then performed in the extracellular matrix. The results demonstrate a 10% drop in the fraction of killed cancer cells 69% rather than the 79% reported earlier for a tumor of similar geometry a more-accurate value. This study examines how tumor-related parameters including the necrotic core size and tumor size, and also drug-related parameters drug dosage, binding affinity of drug, and drug degradation can affect the delivery of the drug to solid tumors. Results indicate that concentration of drug are high in the tumor, low in normal tissue, and remarkably low in the necrotic core. Results also offer a treatment of tumors with smaller necrotic core. Tumor size, which implies the tumor progression, has a considerable impact on treatment outcomes, so to be more effective, treatment should be applied at a specific size of tumor. It is demonstrated that binding affinity of drugs to cell-surface receptors and drug dosage have significant impact on treatment efficacy, so they should be regulated based on a balanced quantification between maximum treatment efficacy and minimum side effects. On the other hand, considering the effects of drug degradation in the model has not significant effect on treatment efficacy. The findings of the present study provide insight into the mechanism of drug delivery to solid tumors based on analyzing the effective parameters and modeling how their behavior in the tumor microenvironment affects treatment efficacy.

Evaluation of inverse methods for estimation of mechanical parameters in solid tumors

To treat cancer, knowledge of mechanical parameters can be essential. This study proposes a new approach for estimating hydraulic conductivity (k) and hydraulic conductivity ratio (α) of a living tissue, based on inverse methods, allowing tissue parameter estimation using only a limited set of measurements. First, two population-based algorithms (Levenberg-Marquardt (LM) method and conjugate gradient (CG) method) and two gradient-based algorithms (genetic algorithm (GA) and particle swarm optimization (PSO) algorithm) are considered, and a comparative study between these different inverse methods is performed to determine which methods have a good performance in terms of convergence rate and stability. CG method is shown to perform well in the case of noise-free input data; however, in the case of noisy input data, it fails to converge. The other three methods (LM, GA, and PSO) converge with estimation errors <10% in both noise-free and noisy input data, suggesting their utility to tackle this problem. In the second part, the effectiveness and good accuracy of these robust algorithms (LM, GA, and PSO) are validated with experimental results. The hydraulic conductivity and hydraulic conductivity ratio of a specific living tumor tissue are then estimated for mammary adenocarcinoma (R3230AC). Moreover, assuming measurement of only one-point interstitial pressure inside the tumor, the effect of the location of this one-point on estimation accuracy of hydraulic conductivity is investigated. We show that estimation errors for points measured near the surface and center of the tumor are greater than at other points.

Spatio-temporal investigation of doxorubicin in a 3D heterogeneous tumor microenvironment

Doxorubicin (Adriamycin) is a type of chemotherapy drugs using to treat diseases such as breast cancer, bladder cancer, Kaposi's sarcoma, and lymphoma. Additionally, it can be first prescribed to reduce tumor size. The ratio of killed cells is varied depending on the clinical dosage regimen. Hence, the exact dosage of the drug must be administered to prevent the toxicity that could impair the immune system or leading to heart failure. In the present study, a 3D heterogeneous geometry with a solid tumor and healthy tissue is modeled for the drug delivery investigation. At the first stage, the physical properties of tumor microenvironment are obtained. Then, a five-compartmental model is used to evaluate the free, bound and internalized drug via the convection-diffusion-reaction (CDR) equation. Results are shown that a percent increase of 37.5% and 47.1% for the 75 mg m^-2 to 50 mg m^-2 in the AUC of bound drug and free drug concentration, respectively. The free and bound drugs have the same trend in time showing an apex at the earliest time of injection and then drops to the lowest amount about 9 hours after treatment. Moreover, the internalized drug has a different trend in time. It increases and reaches a constant amount of drug concentration in the cells. Besides, the fraction of surviving cells is also evaluated for both tumor and healthy tissues showing a 88.62% and 97.76% of surviving cells with 50 mg m^-2 of doxorubicin after the treatment, respectively. This model is developed to predict the heterogenous distribution of doxorubicin in three different drug concentrations for patient-specific drug treatment. This model could be used for different drugs to show the rate of perfusion and the ability to kill cancerous cells regarding their different doses and toxicity effects.

Numerical modeling of high-intensity focused ultrasound-mediated intraperitoneal delivery of thermosensitive liposomal doxorubicin for cancer chemotherapy

Although intraperitoneal chemotherapy (IPC) has been suggested as a promising method for the management of peritoneal dissemination (PD) of ovarian or colorectal cancers, the actual clinical use of this method has been restricted due to such problems as poor drug penetration into the tumor and high side effects. It is, therefore, necessary to develop new strategies to improve the efficacy of this approach. In the present work, a new strategy is proposed based on intraperitoneal (IP) injection of thermosensitive liposomal doxorubicin (TSL-Dox) with triggered release by mild hyperthermia induced by high intensity focused ultrasound (HIFU). A computational model is developed to evaluate the proposed drug delivery system. Results show an order of magnitude increase in drug penetration depth into the tumor compared to the conventional IP delivery. Furthermore, the effects of thermal conditions applied to the tumor, TSL size, tumor vessel permeability, and tumor size are investigated. Results indicate an improved efficiency of the drug delivery by expanding the heated region, yet, it increases the risk of unintentional TSL drug load release in the peritoneal cavity. Results also indicate that smaller TSLs have better treatment outcome. However, there is a significant reduction in treatment efficacy for TSLs with sizes smaller than the vessel wall pore size. Thus, tuning the size of TSL should be based on the tumor microvascular permeability. The simulation results suggest that the TSL-Dox delivery system in smaller tumors is far advantageous than larger ones. Results of our model can be used as guidelines for future preclinical studies.

Image-based spatio-temporal model of drug delivery in a heterogeneous vasculature of a solid tumor - Computational approach

The solute transport distribution in a tumor is an important criterion in the evaluation of the cancer treatment efficacy. The fraction of killed cells after each treatment can quantify the therapeutic effect and plays as a helpful tool to evaluate the chemotherapy treatment schedules. In the present study, an image-based spatio-temporal computational model of a solid tumor is provided for calculation of interstitial fluid flow and solute transport. Current model incorporates heterogeneous microvasculature for angiogenesis instead of synthetic mathematical modeling. In this modeling process, a comprehensive model according to Convection-Diffusion-Reaction (CDR) equations is employed due to its high accuracy for simulating the binding and the uptake of the drug by tumor cells. Based on the velocity and the pressure distribution, transient distribution of the different drug concentrations (free, bound, and internalized) is calculated. Then, the fraction of killed cells is obtained according to the internalized concentration. Results indicate the dependence of the drug distribution on both time and space, as well as the microvasculature density. Free and bound drug concentration have the same trend over time, whereas, internalized and total drug concentration increases over time and reaches a constant value. The highest amount of concentration occurred in the tumor region due to the higher permeability of the blood vessels. Moreover, the fraction of killed cells is approximately 78.87% and 24.94% after treatment with doxorubicin for cancerous and normal tissues, respectively. In general, the presented methodology may be applied in the field of personalized medicine to optimize patient-specific treatments. Also, such image-based modeling of solid tumors can be used in laboratories that working on drug delivery and evaluating new drugs before using them for any in vivo or clinical studies.

Modeling of FMISO [F18] nanoparticle PET tracer in normal-cancerous tissue based on real clinical image

Hypoxia as one of the principal properties of tumor cells is a reaction to the deprivation of oxygen. The location of tumor cells could be identified by assessment of oxygen and nutrient level in human body. Positron emission tomography (PET) is a well-known non-invasive method that is able to measure hypoxia based on the FMISO (Fluoromisonidazole) tracer dynamic. This paper aims to study the PET tracer concentration through convection-diffusion-reaction equations in a real human capillary-like network. A non-uniform oxygen pressure along the capillary path and convection mechanism for FMISO transport are taken into account to accurately model the characteristics of the tracer. To this end, a multi-scale model consists of laminar blood flow through the capillary network, interstitial pressure, oxygen pressure, FMISO diffusion and FMISO convection transport in the extravascular region is developed. The present model considers both normal and tumor tissue regions in computational domain. The accuracy of numerical model is verified with the experimental results available in the literature. The convection and diffusion types of transport mechanism are employed in order to calculate the concentration of FMISO in the normal and tumor sub-domain. The influences of intravascular oxygen pressure, FMISO transport mechanisms, capillary density and different types of tissue on the FMISO concentration have been investigated. According to result (Table 4) the convection mechanism of FMISO molecules transportation is negligible, but it causes more accuracy of the proposed model. The approach of present study can be employed in order to investigate the effects of various parameters, such as tumor shape, on the dynamic behavior of different PET tracers, such as FDG, can be extended to different case study problems, such as drug delivery.

Numerical modeling of nanodrug distribution in tumors with heterogeneous vasculature

The distribution and accumulation of nanoparticle dosage in a tumor are important in evaluating the effectiveness of cancer treatment. The cell survival rate can quantify the therapeutic effect, and the survival rates after multiple treatments are helpful to evaluate the efficacy of a chemotherapy plan. We developed a mathematical tumor model based on the governing equations describing the fluid flow and particle transport to investigate the drug transportation in a tumor and computed the resulting cumulative concentrations. The cell survival rate was calculated based on the cumulative concentration. The model was applied to a subcutaneous tumor with heterogeneous vascular distributions. Various sized dextrans and doxorubicin were respectively chosen as the nanodrug carrier and the traditional chemotherapeutic agent for comparison. The results showed that: 1) the largest nanoparticle drug in the current simulations yielded the highest cumulative concentration in the well vascular region, but second lowest in the surrounding normal tissues, which implies it has the best therapeutic effect to tumor and at the same time little harmful to normal tissue; 2) on the contrary, molecular chemotherapeutic agent produced the second lowest cumulative concentration in the well vascular tumor region, but highest in the surrounding normal tissue; 3) all drugs have very small cumulative concentrations in the tumor necrotic region, where drug transport is solely through diffusion. This might mean that it is hard to kill tumor stem cells hiding in it. The current model indicated that the effectiveness of the anti-tumor drug delivery was determined by the interplay of the vascular density and nanoparticle size, which governs the drug transport properties. The use of nanoparticles as anti-tumor drug carriers is generally a better choice than molecular chemotherapeutic agent because of its high treatment efficiency on tumor cells and less damage to normal tissues.