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

<p>Abstract</p> <p>Background</p> <p>Doxorubicin is a common anticancer agent used in the treatment of a number of neoplasms, with the lifetime dose limited due to the potential for cardiotoxocity. This has motivated efforts to develop optimal dosage regimes that maximize anti-tumor activity while minimizing cardiac toxicity, which is correlated with peak plasma concentration. Doxorubicin is characterized by poor penetration from tumoral vessels into the tumor mass, due to the highly irregular tumor vasculature. I model the delivery of a soluble drug from the vasculature to a solid tumor using a tumor cord model and examine the penetration of doxorubicin under different dosage regimes and tumor microenvironments.</p> <p>Methods</p> <p>A coupled ODE-PDE model is employed where drug is transported from the vasculature into a tumor cord domain according to the principle of solute transport. Within the tumor cord, extracellular drug diffuses and saturable pharmacokinetics govern uptake and efflux by cancer cells. Cancer cell death is also determined as a function of peak intracellular drug concentration.</p> <p>Results</p> <p>The model predicts that transport to the tumor cord from the vasculature is dominated by diffusive transport of free drug during the initial plasma drug distribution phase. I characterize the effect of all parameters describing the tumor microenvironment on drug delivery, and large intercapillary distance is predicted to be a major barrier to drug delivery. Comparing continuous drug infusion with bolus injection shows that the optimum infusion time depends upon the drug dose, with bolus injection best for low-dose therapy but short infusions better for high doses. Simulations of multiple treatments suggest that additional treatments have similar efficacy in terms of cell mortality, but drug penetration is limited. Moreover, fractionating a single large dose into several smaller doses slightly improves anti-tumor efficacy.</p> <p>Conclusion</p> <p>Drug infusion time has a significant effect on the spatial profile of cell mortality within tumor cord systems. Therefore, extending infusion times (up to 2 hours) and fractionating large doses are two strategies that may preserve or increase anti-tumor activity and reduce cardiotoxicity by decreasing peak plasma concentration. However, even under optimal conditions, doxorubicin may have limited delivery into advanced solid tumors.</p

The Johnson-Segalman model with a diffusion term in Couette flow

We study the Johnson-Segalman (JS) model as a paradigm for some complex fluids which are observed to phase separate, or ``shear-band'' in flow. We analyze the behavior of this model in cylindrical Couette flow and demonstrate the history dependence inherent in the local JS model. We add a simple gradient term to the stress dynamics and demonstrate how this term breaks the degeneracy of the local model and prescribes a much smaller (discrete, rather than continuous) set of banded steady state solutions. We investigate some of the effects of the curvature of Couette flow on the observable steady state behavior and kinetics, and discuss some of the implications for metastability.Comment: 14 pp, to be published in Journal of Rheolog

Coexistence and Phase Separation in Sheared Complex Fluids

We demonstrate how to construct dynamic phase diagrams for complex fluids that undergo transitions under flow, in which the conserved composition variable and the broken-symmetry order parameter (nematic, smectic, crystalline, etc.) are coupled to shear rate. Our construction relies on a selection criterion, the existence of a steady interface connecting two stable homogeneous states. We use the (generalized) Doi model of lyotropic nematic liquid crystals as a model system, but the method can be easily applied to other systems, provided non-local effects are included.Comment: 4 pages REVTEX, 5 figures using epsf macros. To appear in Physical Review E (Rapid Communications

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 Switching Mechanism in Doxorubicin Bioactivation Can Be Exploited to Control Doxorubicin Toxicity

Although doxorubicin toxicity in cancer cells is multifactorial, the enzymatic bioactivation of the drug can significantly contribute to its cytotoxicity. Previous research has identified most of the components that comprise the doxorubicin bioactivation network; however, adaptation of the network to changes in doxorubicin treatment or to patient-specific changes in network components is much less understood. To investigate the properties of the coupled reduction/oxidation reactions of the doxorubicin bioactivation network, we analyzed metabolic differences between two patient-derived acute lymphoblastic leukemia (ALL) cell lines exhibiting varied doxorubicin sensitivities. We developed computational models that accurately predicted doxorubicin bioactivation in both ALL cell lines at high and low doxorubicin concentrations. Oxygen-dependent redox cycling promoted superoxide accumulation while NADPH-dependent reductive conversion promoted semiquinone doxorubicin. This fundamental switch in control is observed between doxorubicin sensitive and insensitive ALL cells and between high and low doxorubicin concentrations. We demonstrate that pharmacological intervention strategies can be employed to either enhance or impede doxorubicin cytotoxicity in ALL cells due to the switching that occurs between oxygen-dependent superoxide generation and NADPH-dependent doxorubicin semiquinone formation

Quantifying the Proteolytic Release of Extracellular Matrix-Sequestered VEGF with a Computational Model

BACKGROUND: VEGF proteolysis by plasmin or matrix metalloproteinases (MMPs) is believed to play an important role in regulating vascular patterning in vivo by releasing VEGF from the extracellular matrix (ECM). However, a quantitative understanding of the kinetics of VEGF cleavage and the efficiency of cell-mediated VEGF release is currently lacking. To address these uncertainties, we develop a molecular-detailed quantitative model of VEGF proteolysis, used here in the context of an endothelial sprout. METHODOLOGY AND FINDINGS: To study a cell's ability to cleave VEGF, the model captures MMP secretion, VEGF-ECM binding, VEGF proteolysis from VEGF165 to VEGF114 (the expected MMP cleavage product of VEGF165) and VEGF receptor-mediated recapture. Using experimental data, we estimated the effective bimolecular rate constant of VEGF165 cleavage by plasmin to be 328 M(-1) s(-1) at 25 degrees C, which is relatively slow compared to typical MMP-ECM proteolysis reactions. While previous studies have implicated cellular proteolysis in growth factor processing, we show that single cells do not individually have the capacity to cleave VEGF to any appreciable extent (less than 0.1% conversion). In addition, we find that a tip cell's receptor system will not efficiently recapture the cleaved VEGF due to an inability of cleaved VEGF to associate with Neuropilin-1. CONCLUSIONS: Overall, VEGF165 cleavage in vivo is likely to be mediated by the combined effect of numerous cells, instead of behaving in a single-cell-directed, autocrine manner. We show that heparan sulfate proteoglycans (HSPGs) potentiate VEGF cleavage by increasing the VEGF clearance time in tissues. In addition, we find that the VEGF-HSPG complex is more sensitive to proteases than is soluble VEGF, which may imply its potential relevance in receptor signaling. Finally, according to our calculations, experimentally measured soluble protease levels are approximately two orders of magnitude lower than that needed to reconcile levels of VEGF cleavage seen in pathological situations

Effect of cell arrangement and interstitial volume fraction on the diffusivity of monoclonal antibodies in tissue.

We present theoretical calculations relating the effective diffusivity of monoclonal antibodies in tissue (Deff) to the actual diffusivity in the interstitium (Dint) and the interstitial volume fraction phi. Measured diffusivity values are effective values, deduced from concentration profiles with the tissue treated as a continuum. By using homogenization theory, the ratio Deff/Dint is calculated for a range of interstitial volume fractions from 10 to 65%. It is assumed that only diffusion in the interstitial spaces between cells contributes to the effective diffusivity. The geometries considered have cuboidal cells arranged periodically, with uniform gaps between cells. Deff/Dint is found to generally be between (2/3) phi and phi for these geometries. In general, the pathways for diffusion between cells are not straight. The effect of winding pathways on Deff/Dint is examined by varying the arrangement of the cells, and found to be slight. Also, the estimates of Deff/Dint are shown to be insensitive to typical nonuniformities in the widths of gaps between cells. From our calculations and from published experimental measurements of the effective diffusivity of an IgG polyclonal antibody both in water and in tumor tissue, we deduce that the diffusivity of this molecule in the interstitium is one-tenth to one-twentieth its diffusivity in water. We also conclude that exclusion of molecules from cells (an effect independent of molecular weight) contributes as much as interstitial hindrance to the reduction of effective diffusivity, for small interstitial volume fractions (around 20%). This suggests that the increase in the rate of delivery to tissues resulting from the use of smaller molecular-weight molecules (such as antibody fragments or bifunctional antibodies) may be less than expected