Predicting drug pharmacokinetics and effect in vascularized tumors using computer simulation

In this paper, we investigate the pharmacokinetics and effect of doxorubicin and cisplatin in vascularized tumors through two-dimensional simulations. We take into account especially vascular and morphological heterogeneity as well as cellular and lesion-level pharmacokinetic determinants like P-glycoprotein (Pgp) efflux and cell density. To do this we construct a multi-compartment PKPD model calibrated from published experimental data and simulate 2-h bolus administrations followed by 18-h drug washout. Our results show that lesion-scale drug and nutrient distribution may significantly impact therapeutic efficacy and should be considered as carefully as genetic determinants modulating, for example, the production of multidrug-resistance protein or topoisomerase II. We visualize and rigorously quantify distributions of nutrient, drug, and resulting cell inhibition. A main result is the existence of significant heterogeneity in all three, yielding poor inhibition in a large fraction of the lesion, and commensurately increased serum drug concentration necessary for an average 50% inhibition throughout the lesion (the IC50 concentration). For doxorubicin the effect of hypoxia and hypoglycemia (“nutrient effect”) is isolated and shown to further increase cell inhibition heterogeneity and double the IC50, both undesirable. We also show how the therapeutic effectiveness of doxorubicin penetration therapy depends upon other determinants affecting drug distribution, such as cellular efflux and density, offering some insight into the conditions under which otherwise promising therapies may fail and, more importantly, when they will succeed. Cisplatin is used as a contrast to doxorubicin since both published experimental data and our simulations indicate its lesion distribution is more uniform than that of doxorubicin. Because of this some of the complexity in predicting its therapeutic efficacy is mitigated. Using this advantage, we show results suggesting that in vitro monolayer assays using this drug may more accurately predict in vivo performance than for drugs like doxorubicin. The nonlinear interaction among various determinants representing cell and lesion phenotype as well as therapeutic strategies is a unifying theme of our results. Throughout it can be appreciated that macroscopic environmental conditions, notably drug and nutrient distributions, give rise to considerable variation in lesion response, hence clinical resistance. Moreover, the synergy or antagonism of combined therapeutic strategies depends heavily upon this environment

A Computational Model for Predicting Nanoparticle Accumulation in Tumor Vasculature

Vascular targeting of malignant tissues with systemically injected nanoparticles (NPs) holds promise in molecular imaging and anti-angiogenic therapies. Here, a computational model is presented to predict the development of tumor neovasculature over time and the specific, vascular accumulation of blood-borne NPs. A multidimensional tumor-growth model is integrated with a mesoscale formulation for the NP adhesion to blood vessel walls. The fraction of injected NPs depositing within the diseased vasculature and their spatial distribution is computed as a function of tumor stage, from 0 to day 24 post-tumor inception. As the malignant mass grows in size, average blood flow and shear rates increase within the tumor neovasculature, reaching values comparable with those measured in healthy, pre-existing vessels already at 10 days. The NP vascular affinity, interpreted as the likelihood for a blood-borne NP to firmly adhere to the vessel walls, is a fundamental parameter in this analysis and depends on NP size and ligand density, and vascular receptor expression. For high vascular affinities, NPs tend to accumulate mostly at the inlet tumor vessels leaving the inner and outer vasculature depleted of NPs. For low vascular affinities, NPs distribute quite uniformly intra-tumorally but exhibit low accumulation doses. It is shown that an optimal vascular affinity can be identified providing the proper balance between accumulation dose and uniform spatial distribution of the NPs. This balance depends on the stage of tumor development (vascularity and endothelial receptor expression) and the NP properties (size, ligand density and ligand-receptor molecular affinity). Also, it is demonstrated that for insufficiently developed vascular networks, NPs are transported preferentially through the healthy, pre-existing vessels, thus bypassing the tumor mass. The computational tool described here can effectively select an optimal NP formulation presenting high accumulation doses and uniform spatial intra-tumor distributions as a function of the development stage of the malignancy

Mathematical oncology:how are the mathematical and physical sciences contributing to the war on breast cancer?

Mathematical modeling has recently been added as a tool in the fight against cancer. The field of mathematical oncology has received great attention and increased enormously, but over-optimistic estimations about its ability have created unrealistic expectations. We present a critical appraisal of the current state of mathematical models of cancer. Although the field is still expanding and useful clinical applications may occur in the future, managing over-expectation requires the proposal of alternative directions for mathematical modeling. Here, we propose two main avenues for this modeling: 1) the identification of the elementary biophysical laws of cancer development, and 2) the development of a multiscale mathematical theory as the framework for models predictive of tumor growth. Finally, we suggest how these new directions could contribute to addressing the current challenges of understanding breast cancer growth and metastasis