Integrated intravital microscopy and mathematical modeling to optimize nanotherapeutics delivery to tumors

Spatiotemporal dynamics of tumor - CAR T-cell interaction following local administration in solid cancers

The success of chimeric antigen receptor (CAR) T-cell therapy in treating hematologic malignancies has generated widespread interest in translating this technology to solid cancers. However, issues like tumor infiltration, the immunosuppressive tumor microenvironment, and tumor heterogeneity limit its efficacy in the solid tumor setting. Recent experimental and clinical studies propose local administration directly into the tumor or at the tumor site to increase CAR T-cell infiltration and improve treatment outcomes. Characteristics of the types of solid tumors that may be the most receptive to this treatment approach remain unclear. In this work, we develop a spatiotemporal model for CAR T-cell treatment of solid tumors, and use numerical simulations to compare the effect of introducing CAR T cells via intratumoral injection versus intracavitary administration in diverse cancer types. We demonstrate that the model can recapitulate tumor and CAR T-cell data from imaging studies of local administration of CAR T cells in mouse models. Our results suggest that locally administered CAR T cells will be most successful against slowly proliferating, highly diffusive tumors, which have the lowest average tumor cell density. These findings affirm the clinical observation that CAR T cells will not perform equally across different types of solid tumors, and suggest that measuring tumor density may be helpful when considering the feasibility of CAR T-cell therapy and planning dosages for a particular patient. We additionally find that local delivery of CAR T cells can result in deep tumor responses, provided that the initial CAR T-cell dose does not contain a significant fraction of exhausted cells.

An Insight into Perfusion Anisotropy within Solid Murine Lung Cancer Tumors

Blood vessels are essential for maintaining tumor growth, progression, and metastasis, yet the tumor vasculature is under a constant state of remodeling. Since the tumor vasculature is an attractive therapeutic target, there is a need to predict the dynamic changes in intratumoral fluid pressure and velocity that occur across the tumor microenvironment (TME). The goal of this study was to obtain insight into perfusion anisotropy within lung tumors. To achieve this goal, we used the perfusion marker Hoechst 33342 and vascular endothelial marker CD31 to stain tumor sections from C57BL/6 mice harboring Lewis lung carcinoma tumors on their flank. Vasculature, capillary diameter, and permeability distribution were extracted at different time points along the tumor growth curve. A computational model was generated by applying a unique modeling approach based on the smeared physical fields (Kojic Transport Model, KTM). KTM predicts spatial and temporal changes in intratumoral pressure and fluid velocity within the growing tumor. Anisotropic perfusion occurs within two domains: capillary and extracellular space. Anisotropy in tumor structure causes the nonuniform distribution of pressure and fluid velocity. These results provide insights regarding local vascular distribution for optimal drug dosing and delivery to better predict distribution and duration of retention within the TME.

Anti-tumor effect of pH-sensitive drug-loaded nanoparticles optimized via an integrated computational/experimental approach

The acidic pH of tumor tissue has been used to trigger drug release from nanoparticles. However, dynamic interactions between tumor pH and vascularity present challenges to optimize therapy to particular microenvironment conditions. Despite recent development of pH-sensitive nanomaterials that can accurately quantify drug release from nanoparticles, tailoring release to maximize tumor response remains elusive. This study hypothesizes that a computational modeling-based platform that simulates the heterogeneously vascularized tumor microenvironment can enable evaluation of the complex intra-tumoral dynamics involving nanoparticle transport and pH-dependent drug release, and predict optimal nanoparticle parameters to maximize the response. To this end, SPNCD nanoparticles comprising superparamagnetic cores of iron oxide (Fe3O4) and a poly(lactide-co-glycolide acid) shell loaded with doxorubicin (DOX) were fabricated. Drug release was measured in vitro as a function of pH. A 2D model of vascularized tumor growth was calibrated to experimental data and used to evaluate SPNCD effect as a function of drug release rate and tissue vascular heterogeneity. Simulations show that pH-dependent drug release from SPNCD delays tumor regrowth more than DOX alone across all levels of vascular heterogeneity, and that SPNCD significantly inhibit tumor radius over time compared to systemic DOX. The minimum tumor radius forecast by the model was comparable to previous in vivo SPNCD inhibition data. Sensitivity analyses of the SPNCD pH-dependent drug release rate indicate that slower rates are more inhibitory than faster rates. We conclude that an integrated computational and experimental approach enables tailoring drug release by pH-responsive nanomaterials to maximize the tumor response.

Rational nanoparticle design: Optimization using insights from experiments and mathematical models

Polymeric nanoparticles are highly tunable drug delivery systems that show promise in targeting therapeutics to specific sites within the body. Rational nanoparticle design can make use of mathematical models to organize and extend experimental data, allowing for optimization of nanoparticles for particular drug delivery applications. While rational nanoparticle design is attractive from the standpoint of improving therapy and reducing unnecessary experiments, it has yet to be fully realized. The difficulty lies in the complexity of nanoparticle structure and behavior, which is added to the complexity of the physiological mechanisms involved in nanoparticle distribution throughout the body. In this review, we discuss the most important aspects of rational design of polymeric nanoparticles. Ultimately, we conclude that many experimental datasets are required to fully model polymeric nanoparticle behavior at multiple scales. Further, we suggest ways to consider the limitations and uncertainty of experimental data in creating nanoparticle design optimization schema, which we call quantitative nanoparticle design frameworks.

Real-time imaging of nanobubble ultrasound contrast agent flow, extravasation, and diffusion through an extracellular matrix using a microfluidic model

Lipid shell-stabilized nanoparticles with a perfluorocarbon gas-core, or nanobubbles, have recently attracted attention as a new contrast agent for molecular ultrasound imaging and image-guided therapy. Due to their small size (∼275 nm diameter) and flexible shell, nanobubbles have been shown to extravasate through hyperpermeable vasculature (e.g., in tumors). However, little is known about the dynamics and depth of extravasation of intact, acoustically active nanobubbles. Accordingly, in this work, we developed a microfluidic chip with a lumen and extracellular matrix (ECM) and imaging method that allows real-time imaging and characterization of the extravasation process with high-frequency ultrasound. The microfluidic device has a lumen and is surrounded by an extracellular matrix with tunable porosity. The combination of ultrasound imaging and the microfluidic chip advantageously produces real-time images of the entire length and depth of the matrix. This captures the matrix heterogeneity, offering advantages over other imaging techniques with smaller fields of view. Results from this study show that nanobubbles diffuse through a 1.3 μm pore size (2 mg mL-1) collagen I matrix 25× faster with a penetration depth that was 0.19 mm deeper than a 3.7 μm (4 mg mL-1) matrix. In the 3.7 μm pore size matrix, nanobubbles diffused 92× faster than large nanobubbles (∼875 nm diameter). Decorrelation time analysis was successfully used to differentiate flowing and extra-luminally diffusing nanobubbles. In this work, we show for the first time that combination of an ultrasound-capable microfluidic chip and real-time imaging provided valuable insight into spatiotemporal nanoparticle movement through a heterogeneous extracellular matrix. This work could help accurately predict parameters (e.g., injection dosage) that improve translation of nanoparticles from in vitro to in vivo environments.

Evaluation of Lung Cancer Patient Response to First-Line Chemotherapy by Integration of Tumor Core Biopsy Metabolomics with Multiscale Modeling

The standard of care for intermediate (Stage II) and advanced (Stages III and IV) non-small cell lung cancer (NSCLC) involves chemotherapy with taxane/platinum derivatives, with or without radiation. Ideally, patients would be screened a priori to allow non-responders to be initially treated with second-line therapies. This evaluation is non-trivial, however, since tumors behave as complex multiscale systems. To address this need, this study employs a multiscale modeling approach to evaluate first-line chemotherapy response of individual patient tumors based on metabolomic analysis of tumor core biopsies obtained during routine clinical evaluation. Model parameters were calculated for a patient cohort as a function of these metabolomic profiles, previously obtained from high-resolution 2DLC-MS/MS analysis. Evaluation metrics were defined to classify patients as Disease-Control (DC) [encompassing complete-response (CR), partial-response (PR), and stable-disease (SD)] and Progressive-Disease (PD) following first-line chemotherapy. Response was simulated for each patient and compared to actual response. The results show that patient classifications were significantly separated from each other, and also when grouped as DC vs. PD and as CR/PR vs. SD/PD, by fraction of initial tumor radius metric at 6 days post simulated bolus drug injection. This study shows that patient first-line chemotherapy response can in principle be evaluated from multiscale modeling integrated with tumor tissue metabolomic data, offering a first step towards individualized lung cancer treatment prognosis.

Using Parallel Coordinates in Optimization of Nano-Particle Drug Delivery

Nanoparticle drug delivery better targets neoplastic lesions than free drugs and thus has emerged as a safer form of cancer therapy. Nanoparticle design variables are important determinants of efficacy as they influence the drug biodistribution and pharmacokinetics. Previously, we determined optimal designs through mechanistic modeling and optimization. However, the numerical nature of the tumor model and numerous candidate nanoparticle designs hinder hypothesis generation and treatment personalization. In this paper, we utilize the parallel coordinates technique to visualize high-dimensional optimal solutions and extract correlations between nanoparticle design and treatment outcomes. We found that at optimality, two major design variables are dependent, and thus the optimization problem can be reduced. In addition, we obtained an analytical relationship between optimal nanoparticle sizes and optimal distribution, which could facilitate the utilization of tumors models in preclinical studies. Our approach has simplified the results of the previously integrated modeling and optimization framework developed for nanotherapy and enhanced the interpretation and utilization of findings. Integrated mathematical frameworks are increasing in the medical field, and our method can be applied outside nanotherapy to facilitate the clinical translation of computational methods.

Modeling of Tumor Growth with Input from Patient-Specific Metabolomic Data

Advances in omic technologies have provided insight into cancer progression and treatment response. However, the nonlinear characteristics of cancer growth present a challenge to bridge from the molecular- to the tissue-scale, as tumor behavior cannot be encapsulated by the sum of the individual molecular details gleaned experimentally. Mathematical modeling and computational simulation have been traditionally employed to facilitate analysis of nonlinear systems. In this study, for the first time tumor metabolomic data are linked via mathematical modeling to the tumor tissue-scale behavior, showing the capability to mechanistically simulate cancer progression personalized to omic information obtainable from patient tumor core biopsy analysis. Generally, a higher degree of metabolic dysregulation has been correlated with more aggressive tumor behavior. Accordingly, key parameters influenced by metabolomic data in this model include tumor proliferation, vascularization, aggressiveness, lactic acid production, monocyte infiltration and macrophage polarization, and drug effect. The model enables evaluating interactions of interest between these parameters which drive tumor growth based on the metabolomic data. The results show that the model can group patients consistently with the clinically observed outcomes of response/non-response to chemotherapy. This modeling approach provides a first step towards evaluation of tumor growth based on tumor-specific metabolomic data.

Multispectral intravital microscopy for simultaneous bright-field and fluorescence imaging of the microvasculature

Intravital video microscopy permits the observation of microcirculatory blood flow. This often requires fluorescent probes to visualize structures and dynamic processes that cannot be observed with conventional bright-field microscopy. Conventional light microscopes do not allow for simultaneous bright-field and fluorescent imaging. Moreover, in conventional microscopes, only one type of fluorescent label can be observed. This study introduces multispectral intravital video microscopy, which combines bright-field and fluorescence microscopy in a standard light microscope. The technique enables simultaneous real-time observation of fluorescently-labeled structures in relation to their direct physical surroundings. The advancement provides context for the orientation, movement, and function of labeled structures in the microcirculation.

Modeling of tumor response to macrophage and T lymphocyte interactions in the liver metastatic microenvironment

The dynamic interactions between macrophages and T-lymphocytes in the tumor microenvironment exert both antagonistic and synergistic functions affecting tumor growth. Extensive experimental effort has been expended to investigate immunotherapeutic strategies targeting macrophage polarization as well as T-cell activation with the goal to promote tumor cell killing and cancer elimination. However, these interactions remain poorly understood, and cancer immunotherapeutic strategies are often disappointing. The complex system encompassing innate and adaptive immune cell activity in response to tumor growth could benefit from a systems perspective built upon mathematical modeling. This study develops a modeling system to help evaluate the effects of macrophage and T-lymphocyte interactions on tumor growth. The system enables simulating the combined cytotoxic and tumor-promoting interactions of these two immune cell populations in a vascularized organ microenvironment, such as in liver metastases. A hypothetical immunotherapeutic strategy is simulated to increase the number of tumor-suppressive (M1-phenotype) vs. tumor-promoting (M2-phenotype) macrophages to gauge their effects on CD8+ T-cells and CD4+ T-helper cells, which in turn affect the macrophage functions. The results highlight the dynamic interactions between macrophages and T-lymphocytes in the tumor microenvironment and show that with the chosen set of parameter values, the overall cytotoxic effect from macrophages and T-lymphocytes obtained by driving the M1:M2 ratio higher could saturate and fail to achieve tumor regression. Further expansion of this modeling platform to include additional tumor-immune cell interactions, coupled with parameters representing particular tumor characteristics, could enable systematic evaluation of immunotherapeutic strategies tailored to patient-tumor specific conditions, including metastatic disease.

Tumor ablation due to inhomogeneous anisotropic diffusion in generic three-dimensional topologies

In recent decades computer-aided technologies have become prevalent in medicine, however, cancer drugs are often only tested on in vitro cell lines from biopsies. We derive a full three-dimensional model of inhomogeneous -anisotropic diffusion in a tumor region coupled to a binary population model, which simulates in vivo scenarios faster than traditional cell-line tests. The diffusion tensors are acquired using diffusion tensor magnetic resonance imaging from a patient diagnosed with glioblastoma multiform. Then we numerically simulate the full model with finite element methods and produce drug concentration heat maps, apoptosis hotspots, and dose-response curves. Finally, predictions are made about optimal injection locations and volumes, which are presented in a form that can be employed by doctors and oncologists.

A computational diffusion model to study antibody transport within reconstructed tumor microenvironments

Background

Antibodies revolutionized cancer treatment over the past decades. Despite their successfully application, there are still challenges to overcome to improve efficacy, such as the heterogeneous distribution of antibodies within tumors. Tumor microenvironment features, such as the distribution of tumor and other cell types and the composition of the extracellular matrix may work together to hinder antibodies from reaching the target tumor cells. To understand these interactions, we propose a framework combining in vitro and in silico models. We took advantage of in vitro cancer models previously developed by our group, consisting of tumor cells and fibroblasts co-cultured in 3D within alginate capsules, for reconstruction of tumor microenvironment features.

Results

In this work, an experimental-computational framework of antibody transport within alginate capsules was established, assuming a purely diffusive transport, combined with an exponential saturation effect that mimics the saturation of binding sites on the cell surface. Our tumor microenvironment in vitro models were challenged with a fluorescent antibody and its transport recorded using light sheet fluorescence microscopy. Diffusion and saturation parameters of the computational model were adjusted to reproduce the experimental antibody distribution, with root mean square error under 5%. This computational framework is flexible and can simulate different random distributions of tumor microenvironment elements (fibroblasts, cancer cells and collagen fibers) within the capsule. The random distribution algorithm can be tuned to follow the general patterns observed in the experimental models.

Conclusions

We present a computational and microscopy framework to track and simulate antibody transport within the tumor microenvironment that complements the previously established in vitro models platform. This framework paves the way to the development of a valuable tool to study the influence of different components of the tumor microenvironment on antibody transport.

Modeling of Nanotherapy Response as a Function of the Tumor Microenvironment: Focus on Liver Metastasis

The tumor microenvironment (TME) presents a challenging barrier for effective nanotherapy-mediated drug delivery to solid tumors. In particular for tumors less vascularized than the surrounding normal tissue, as in liver metastases, the structure of the organ itself conjures with cancer-specific behavior to impair drug transport and uptake by cancer cells. Cells and elements in the TME of hypovascularized tumors play a key role in the process of delivery and retention of anti-cancer therapeutics by nanocarriers. This brief review describes the drug transport challenges and how they are being addressed with advanced in vitro 3D tissue models as well as with in silico mathematical modeling. This modeling complements network-oriented techniques, which seek to interpret intra-cellular relevant pathways and signal transduction within cells and with their surrounding microenvironment. With a concerted effort integrating experimental observations with computational analyses spanning from the molecular- to the tissue-scale, the goal of effective nanotherapy customized to patient tumor-specific conditions may be finally realized.

Multi-objective optimization of tumor response to drug release from vasculature-bound nanoparticles

The pharmacokinetics of nanoparticle-borne drugs targeting tumors depends critically on nanoparticle design. Empirical approaches to evaluate such designs in order to maximize treatment efficacy are time- and cost-intensive. We have recently proposed the use of computational modeling of nanoparticle-mediated drug delivery targeting tumor vasculature coupled with numerical optimization to pursue optimal nanoparticle targeting and tumor uptake. Here, we build upon these studies to evaluate the effect of tumor size on optimal nanoparticle design by considering a cohort of heterogeneously-sized tumor lesions, as would be clinically expected. The results indicate that smaller nanoparticles yield higher tumor targeting and lesion regression for larger-sized tumors. We then augment the nanoparticle design optimization problem by considering drug diffusivity, which yields a two-fold tumor size decrease compared to optimizing nanoparticles without this consideration. We quantify the tradeoff between tumor targeting and size decrease using bi-objective optimization, and generate five Pareto-optimal nanoparticle designs. The results provide a spectrum of treatment outcomes - considering tumor targeting vs. antitumor effect - with the goal to enable therapy customization based on clinical need. This approach could be extended to other nanoparticle-based cancer therapies, and support the development of personalized nanomedicine in the longer term.

Nonlinear response to cancer nanotherapy due to macrophage interactions revealed by mathematical modeling and evaluated in a murine model via CRISPR-modulated macrophage polarization

Tumor-associated macrophages (TAMs) have been shown to both aid and hinder tumor growth, with patient outcomes potentially hinging on the proportion of M1, pro-inflammatory/growth-inhibiting, to M2, growth-supporting, phenotypes. Strategies to stimulate tumor regression by promoting polarization to M1 are a novel approach that harnesses the immune system to enhance therapeutic outcomes, including chemotherapy. We recently found that nanotherapy with mesoporous particles loaded with albumin-bound paclitaxel (MSV-nab-PTX) promotes macrophage polarization towards M1 in breast cancer liver metastases (BCLM). However, it remains unclear to what extent tumor regression can be maximized based on modulation of the macrophage phenotype, especially for poorly perfused tumors such as BCLM. Here, for the first time, a CRISPR system is employed to permanently modulate macrophage polarization in a controlled in vitro setting. This enables the design of 3D co-culture experiments mimicking the BCLM hypovascularized environment with various ratios of polarized macrophages. We implement a mathematical framework to evaluate nanoparticle-mediated chemotherapy in conjunction with TAM polarization. The response is predicted to be not linearly dependent on the M1:M2 ratio. To investigate this phenomenon, the response is simulated via the model for a variety of M1:M2 ratios. The modeling indicates that polarization to an all-M1 population may be less effective than a combination of both M1 and M2. Experimental results with the CRISPR system confirm this model-driven hypothesis. Altogether, this study indicates that response to nanoparticle-mediated chemotherapy targeting poorly perfused tumors may benefit from a fine-tuned M1:M2 ratio that maintains both phenotypes in the tumor microenvironment during treatment.

Image-guided mathematical modeling for pharmacological evaluation of nanomaterials and monoclonal antibodies

While plasma concentration kinetics has traditionally been the predictor of drug pharmacological effects, it can occasionally fail to represent kinetics at the site of action, particularly for solid tumors. This is especially true in the case of delivery of therapeutic macromolecules (drug-loaded nanomaterials or monoclonal antibodies), which can experience challenges to effective delivery due to particle size-dependent diffusion barriers at the target site. As a result, disparity between therapeutic plasma kinetics and kinetics at the site of action may exist, highlighting the importance of target site concentration kinetics in determining the pharmacodynamic effects of macromolecular therapeutic agents. Assessment of concentration kinetics at the target site has been facilitated by non-invasive in vivo imaging modalities. This allows for visualization and quantification of the whole-body disposition behavior of therapeutics that is essential for a comprehensive understanding of their pharmacokinetics and pharmacodynamics. Quantitative non-invasive imaging can also help guide the development and parameterization of mathematical models for descriptive and predictive purposes. Here, we present a review of the application of state-of-the-art imaging modalities for quantitative pharmacological evaluation of therapeutic nanoparticles and monoclonal antibodies, with a focus on their integration with mathematical models, and identify challenges and opportunities. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease Diagnostic Tools > in vivo Nanodiagnostics and Imaging Nanotechnology Approaches to Biology > Nanoscale Systems in Biology.

Modeling of Combination Chemotherapy and Immunotherapy for Lung Cancer

Cancer has traditionally been studied from a basic science perspective, focusing on the underlying biology, physiology, and biochemistry. Engineering has supplemented this effort via the development of technology, e.g., microscopy. In recent times, engineering and the physical sciences have positioned themselves as approaches on par with the traditional basic sciences to tackle the study of cancer. Mathematical modeling and computational simulation have become key elements of this engineering-focused effort, evaluating the growth of tumors and their response to therapy as problems that could benefit from a systems analysis perspective. Building upon previous work in this field, here is developed a modeling framework to help evaluate the response of tumors to the combination of chemotherapy and immunotherapy, focusing on non-small cell lung cancer (NSCLC). With system parameters set with patient tumor-specific parameters, the longer term goal of this work is to advance personalized cancer treatment.

Mesenchymal stem cells used as carrier cells of oncolytic adenovirus results in enhanced oncolytic virotherapy

Mesenchymal stem cells (MSCs) loaded with oncolytic viruses are presently being investigated as a new modality of advanced/metastatic tumors treatment and enhancement of virotherapy. MSCs can, however, either promote or suppress tumor growth. To address the critical question of how MSCs loaded with oncolytic viruses affect virotherapy outcomes and tumor growth patterns in a tumor microenvironment, we developed and analyzed an integrated mathematical-experimental model. We used the model to describe both the growth dynamics in our experiments of firefly luciferase-expressing Hep3B tumor xenografts and the effects of the immune response during the MSCs-based virotherapy. We further employed it to explore the conceptual clinical feasibility, particularly, in evaluating the relative significance of potential immune promotive/suppressive mechanisms induced by MSCs loaded with oncolytic viruses. We were able to delineate conditions which may significantly contribute to the success or failure of MSC-based virotherapy as well as generate new hypotheses. In fact, one of the most impactful outcomes shown by this investigation, not inferred from the experiments alone, was the initially counter-intuitive fact that using tumor-promoting MSCs as carriers is not only helpful but necessary in achieving tumor control. Considering the fact that it is still currently a controversial debate whether MSCs exert a pro- or anti-tumor action, mathematical models such as this one help to quantitatively predict the consequences of using MSCs for delivering virotherapeutic agents in vivo. Taken together, our results show that MSC-mediated systemic delivery of oncolytic viruses is a promising strategy for achieving synergistic anti-tumor efficacy with improved safety profiles.

Improving nanotherapy delivery and action through image-guided systems pharmacology

Despite recent advances in the translation of therapeutic nanoparticles (TNPs) into the clinic, the field continues to face challenges in predictably and selectively delivering nanomaterials for the treatment of solid cancers. The concept of enhanced permeability and retention (EPR) has been coined as a convenient but simplistic descriptor of high TNP accumulation in some tumors. However, in practice EPR represents a number of physiological variables rather than a single one (including dysfunctional vasculature, compromised lymphatics and recruited host cells, among other aspects of the tumor microenvironment) — each of which can be highly heterogenous within a given tumor, patient and across patients. Therefore, a clear need exists to dissect the specific biophysical factors underlying the EPR effect, to formulate better TNP designs, and to identify patients with high-EPR tumors who are likely to respond to TNP. The overall pharmacology of TNP is governed by an interconnected set of spatially defined and dynamic processes that benefit from a systems-level quantitative approach, and insights into the physiology have profited from the marriage between in vivo imaging and quantitative systems pharmacology (QSP) methodologies. In this article, we review recent developments pertinent to image-guided systems pharmacology of nanomedicines in oncology. We first discuss recent developments of quantitative imaging technologies that enable analysis of nanomaterial pharmacology at multiple spatiotemporal scales, and then examine reports that have adopted these imaging technologies to guide QSP approaches. In particular, we focus on studies that have integrated multi-scale imaging with computational modeling to derive insights about the EPR effect, as well as studies that have used modeling to guide the manipulation of the EPR effect and other aspects of the tumor microenvironment for improving TNP action. We anticipate that the synergistic combination of imaging with systems-level computational methods for effective clinical translation of TNPs will only grow in relevance as technologies increase in resolution, multiplexing capability, and in the ability to examine heterogeneous behaviors at the single-cell level.

Pharmacokinetic/Pharmacodynamics Modeling of Drug-Loaded PLGA Nanoparticles Targeting Heterogeneously Vascularized Tumor Tissue

PURPOSE

Nanoparticle-mediated drug delivery and efficacy for cancer applications depends on systemic as well as local microenvironment characteristics. Here, a novel coupling of a nanoparticle (NP) kinetic model with a drug pharmacokinetic/pharmacodynamics model evaluates efficacy of cisplatin-loaded poly lactic-co-glycolic acid (PLGA) NPs in heterogeneously vascularized tumor tissue.

METHODS

Tumor lesions are modeled with various levels of vascular heterogeneity, as would be encountered with different types of tumors. The magnitude of the extracellular to cytosolic NP transport is varied to assess tumor-dependent cellular uptake. NP aggregation is simulated to evaluate its effects on drug distribution and tumor response.

RESULTS

Cisplatin-loaded PLGA NPs are most effective in decreasing tumor size in the case of high vascular-induced heterogeneity, a high NP cytosolic transfer coefficient, and no NP aggregation. Depending on the level of tissue heterogeneity, NP cytosolic transfer and drug half-life, NP aggregation yielding only extracellular drug release could be more effective than unaggregated NPs uptaken by cells and releasing drug both extra- and intra-cellularly.

CONCLUSIONS

Model-based customization of PLGA NP and drug design parameters, including cellular uptake and aggregation, tailored to patient tumor tissue characteristics such as proportion of viable tissue and vascular heterogeneity, could help optimize the NP-mediated tumor drug response.

Computational/experimental evaluation of liver metastasis post hepatic injury: interactions with macrophages and transitional ECM

The complex interactions between subclinical changes to hepatic extracellular matrix (ECM) in response to injury and tumor-associated macrophage microenvironmental cues facilitating metastatic cell seeding remain poorly understood. This study implements a combined computational modeling and experimental approach to evaluate tumor growth following hepatic injury, focusing on ECM remodeling and interactions with local macrophages. Experiments were performed to determine ECM density and macrophage-associated cytokine levels. Effects of ECM remodeling along with macrophage polarization on tumor growth were evaluated via computational modeling. For primary or metastatic cells in co-culture with macrophages, TNF-α levels were 5× higher with M1 vs. M2 macrophages. Metastatic cell co-culture exhibited 10× higher TNF-α induction than with primary tumor cells. Although TGFβ1 induction was similar between both co-cultures, levels were slightly higher with primary cells in the presence of M1. Simulated metastatic tumors exhibited decreased growth compared to primary tumors, due to high local M1-induced cytotoxicity, even in a highly vascularized microenvironment. Experimental analysis combined with computational modeling may provide insight into interactions between ECM remodeling, macrophage polarization, and liver tumor growth.

Nanoparticle Optimization for Enhanced Targeted Anticancer Drug Delivery

Nanoparticle (NP)-based drug delivery is a promising method to increase the therapeutic index of anticancer agents with low median toxic dose. The delivery efficiency, corresponding to the fraction of the injected NPs that adhere to the tumor site, depends on NP size a and aspect ratio AR. Values for these variables are currently chosen empirically, which may not result in optimal targeted drug delivery. This study applies rigorous optimization to the design of NPs. A preliminary investigation revealed that delivery efficiency increases monotonically with a and AR. However, maximizing a and AR results in nonuniform drug distribution, which impairs tumor regression. Therefore, a multiobjective optimization (MO) problem is formulated to quantify the trade-off between NPs accumulation and distribution. The MO is solved using the derivative-free mesh adaptive direct search algorithm. Theoretically, the Pareto-optimal set consists of an infinite number of mathematically equivalent solutions to the MO problem. However, interesting design solutions can be identified subjectively, e.g., the ellipsoid with a major axis of 720 nm and an aspect ratio of 7.45, as the solution closest to the utopia point. The MO problem formulation is then extended to optimize NP biochemical properties: ligand-receptor binding affinity and ligand density. Optimizing physical and chemical properties simultaneously results in optimal designs with reduced NP sizes and thus enhanced cellular uptake. The presented study provides an insight into NP structures that have potential for producing desirable drug delivery.

Efficacy of Surface-Modified PLGA Nanoparticles as a Function of Cervical Cancer Type

PURPOSE

Hypovascularization of cervical tumors, coupled with intrinsic and acquired drug resistance, has contributed to marginal therapeutic outcomes by hindering chemotherapeutic transport and efficacy. Recently, the heterogeneous penetration and distribution of cell penetrating peptide (CPP, here MPG) and polyethylene glycol (PEG) modified poly(lactic-co-glycolic acid) (PLGA) nanoparticles (NPs) were evaluated as a function of tumor type and morphology in cervical cancer spheroids modeling hypovascularized tumor nodules. Building upon this work, this study investigates the efficacy imparted by surface-modified Doxorubicin-loaded NPs transported into hypovascularized tissue.

METHODS

NP efficacy was measured in HeLa, CaSki, and SiHa cells. NP internalization and association, and associated cell viability, were determined in monolayer and spheroid models.

RESULTS

MPG and PEG-NP co-treatment was most efficacious in HeLa cells, while PEG NPs were most efficacious in CaSki cells. NP surface-modifications were unable to improve efficacy, relative to unmodified NPs, in SiHa cells.

CONCLUSIONS

The results highlight the dependence of efficacy on tumor type and the associated microenvironment. The results further relate previous NP transport studies to efficacy, as a function of surface-modification and cell type. Longer-term, this information may help guide the design of NP-mediated strategies to maximize efficacy based on patient-specific cervical tumor origin and characteristics.

Design Optimization of Tumor Vasculature-Bound Nanoparticles

Nanotherapy may constitute a promising approach to target tumors with anticancer drugs while minimizing systemic toxicity. Computational modeling can enable rapid evaluation of nanoparticle (NP) designs and numerical optimization. Here, an optimization study was performed using an existing tumor model to find NP size and ligand density that maximize tumoral NP accumulation while minimizing tumor size. Optimal NP avidity lies at lower bound of feasible values, suggesting reduced ligand density to prolong NP circulation. For the given set of tumor parameters, optimal NP diameters were 288 nm to maximize NP accumulation and 334 nm to minimize tumor diameter, leading to uniform NP distribution and adequate drug load. Results further show higher dependence of NP biodistribution on the NP design than on tumor morphological parameters. A parametric study with respect to drug potency was performed. The lower the potency of the drug, the bigger the difference is between the maximizer of NP accumulation and the minimizer of tumor size, indicating the existence of a specific drug potency that minimizes the differential between the two optimal solutions. This study shows the feasibility of applying optimization to NP designs to achieve efficacious cancer nanotherapy, and offers a first step towards a quantitative tool to support clinical decision making.

Modeling three-dimensional invasive solid tumor growth in heterogeneous microenvironment under chemotherapy

A systematic understanding of the evolution and growth dynamics of invasive solid tumors in response to different chemotherapy strategies is crucial for the development of individually optimized oncotherapy. Here, we develop a hybrid three-dimensional (3D) computational model that integrates pharmacokinetic model, continuum diffusion-reaction model and discrete cell automaton model to investigate 3D invasive solid tumor growth in heterogeneous microenvironment under chemotherapy. Specifically, we consider the effects of heterogeneous environment on drug diffusion, tumor growth, invasion and the drug-tumor interaction on individual cell level. We employ the hybrid model to investigate the evolution and growth dynamics of avascular invasive solid tumors under different chemotherapy strategies. Our simulations indicate that constant dosing is generally more effective in suppressing primary tumor growth than periodic dosing, due to the resulting continuous high drug concentration. In highly heterogeneous microenvironment, the malignancy of the tumor is significantly enhanced, leading to inefficiency of chemotherapies. The effects of geometrically-confined microenvironment and non-uniform drug dosing are also investigated. Our computational model, when supplemented with sufficient clinical data, could eventually lead to the development of efficient in silico tools for prognosis and treatment strategy optimization.

Evaluation of Drug-Loaded Gold Nanoparticle Cytotoxicity as a Function of Tumor Vasculature-Induced Tissue Heterogeneity

The inherent heterogeneity of tumor tissue presents a major challenge to nanoparticle-mediated drug delivery. This heterogeneity spans from the molecular (genomic, proteomic, metabolomic) to the cellular (cell types, adhesion, migration) and to the tissue (vasculature, extra-cellular matrix) scales. In particular, tumor vasculature forms abnormally, inducing proliferative, hypoxic, and necrotic tumor tissue regions. As the vasculature is the main conduit for nanotherapy transport into tumors, vasculature-induced tissue heterogeneity can cause local inadequate delivery and concentration, leading to subpar response. Further, hypoxic tissue, although viable, would be immune to the effects of cell-cycle specific drugs. In order to enable a more systematic evaluation of such effects, here we employ computational modeling to study the therapeutic response as a function of vasculature-induced tumor tissue heterogeneity. Using data with three-layered gold nanoparticles loaded with cisplatin, nanotherapy is simulated interacting with different levels of tissue heterogeneity, and the treatment response is measured in terms of tumor regression. The results quantify the influence that varying levels of tumor vascular density coupled with the drug strength have on nanoparticle uptake and washout, and the associated tissue response. The drug strength affects the proportion of proliferating, hypoxic, and necrotic tissue fractions, which in turn dynamically affect and are affected by the vascular density. Higher drug strengths may be able to achieve stronger tumor regression but only if the intra-tumoral vascular density is above a certain threshold that affords sufficient transport. This study establishes an initial step towards a more systematic methodology to assess the effect of vasculature-induced tumor tissue heterogeneity on the response to nanotherapy.

Computational fluid dynamics with imaging of cleared tissue and of in vivo perfusion predicts drug uptake and treatment responses in tumours

Understanding the uptake of a drug by diseased tissue, and the drug's subsequent spatiotemporal distribution, are central factors in the development of effective targeted therapies. However, the interaction between the pathophysiology of diseased tissue and individual therapeutic agents can be complex, and can vary across tissue types and across subjects. Here, we show that the combination of mathematical modelling, high-resolution optical imaging of intact and optically cleared tumour tissue from animal models, and in vivo imaging of vascular perfusion predicts the heterogeneous uptake, by large tissue samples, of specific therapeutic agents, as well as their spatiotemporal distribution. In particular, by using murine models of colorectal cancer and glioma, we report and validate predictions of steady-state blood flow and intravascular and interstitial fluid pressure in tumours, of the spatially heterogeneous uptake of chelated gadolinium by tumours, and of the effect of a vascular disrupting agent on tumour vasculature.

Pharmacokinetics variability: Why nanoparticles are not just magic-bullets in oncology

Developing nanoparticles to improve the specificity of anticancer agents towards tumor tissue and to better control drug delivery is a rising strategy in oncology. An increasing number of forms (e.g., conjugated nanoparticles, liposomes, immunoliposomes…) are now available on the shelves and numerous other scaffolds (e.g., dendrimeres, nanospheres, squalenes …) are currently at various stages of development. However, as of today most nanoparticles made available remain lipidic carriers. Pharmacokinetic variability is a major, yet largely underestimated issue with liposomal nanoparticles. A wide variety of causes (e.g., tumor type and disease staging, comorbidities, patient's immune system) can explain this variability, which can in return negatively impact pharmacodynamic endpoints such as poor efficacy or severe toxicities. This review aims to cover the main causes for erratic pharmacokinetics observed with most nanoparticles, especially liposomes used in oncology. Should the main causes of such variability be identified, specific studies in non-clinical or clinical development stages could be undertaken using dedicated models (i.e., mechanistic or semi-mechanistic mathematical models such as PBPK approaches) to better describe nanoparticles pharmacokinetics and decipher PK/PD relationships. In addition, identifying relevant biomarkers or parameters likely to impact nanoparticles pharmacokinetics would allow for either the modification of their characteristics to reduce the influence of the expected variability during development phases or the development of biomarker-based adaptive dosing strategies to maintain an optimal efficacy/toxicity balance. Broadly, we call for the development of comprehensive distribution studies and state-of-the-art modeling support to better understand and anticipate nanoparticle pharmacokinetics in oncology.

Nanoparticle-mediated drug delivery to treat infections in the female reproductive tract: evaluation of experimental systems and the potential for mathematical modeling

A variety of drug-delivery platforms have been employed to deliver therapeutic agents across cervicovaginal mucus (CVM) and the vaginal mucosa, offering the capability to increase the longevity and retention of active agents to treat infections of the female reproductive tract (FRT). Nanoparticles (NPs) have been shown to improve retention, diffusion, and cell-specific targeting via specific surface modifications, relative to other delivery platforms. In particular, polymeric NPs represent a promising option that has shown improved distribution through the CVM. These NPs are typically fabricated from nontoxic, non-inflammatory, US Food and Drug Administration-approved polymers that improve biocompatibility. This review summarizes recent experimental studies that have evaluated NP transport in the FRT, and highlights research areas that more thoroughly and efficiently inform polymeric NP design, including mathematical modeling. An overview of the in vitro, ex vivo, and in vivo NP studies conducted to date – whereby transport parameters are determined, extrapolated, and validated – is presented first. The impact of different NP design features on transport through the FRT is summarized, and gaps that exist due to the limitations of iterative experimentation alone are identified. The potential of mathematical modeling to complement the characterization and evaluation of diffusion and transport of delivery vehicles and active agents through the CVM and mucosa is discussed. Lastly, potential advancements combining experimental and mathematical knowledge are suggested to inform next-generation NP designs, such that infections in the FRT may be more effectively treated.

Pharmacokinetic/pharmacodynamic modeling of combination-chemotherapy for lung cancer

Chemotherapy for non-small cell lung cancer (NSCLC) typically involves a doublet regimen for a number of cycles. For any particular patient, a course of treatment is usually chosen from a large number of combinational protocols with drugs in concomitant or sequential administration. In spite of newer drugs and protocols, half of patients with early disease will live less than five years and 95% of those with advanced disease survive for less than one year. Here, we apply mathematical modeling to simulate tumor response to multiple drug regimens, with the capability to assess maximum tolerated dose (MTD) as well as metronomic drug administration. We couple pharmacokinetic-pharmacodynamic intracellular multi-compartment models with a model of vascularized tumor growth, setting input parameters from in vitro data, and using the models to project potential response in vivo. This represents an initial step towards the development of a comprehensive virtual system to evaluate tumor response to combinatorial drug regimens, with the goal to more efficiently identify optimal course of treatment with patient tumor-specific data. We evaluate cisplatin and gemcitabine with clinically-relevant dosages, and simulate four treatment NSCLC scenarios combining MTD and metronomic therapy. This work thus establishes a framework for systematic evaluation of tumor response to combination chemotherapy. The results with the chosen parameter set indicate that although a metronomic regimen may provide advantage over MTD, the combination of these regimens may not necessarily offer improved response. Future model evaluation of chemotherapy possibilities may help to assess their potential value to obtain sustained NSCLC regression for particular patients, with the ultimate goal of optimizing multiple-drug chemotherapy regimens in clinical practice.

Mathematical modeling of tumor-associated macrophage interactions with the cancer microenvironment

Background

Immuno-oncotherapy has emerged as a promising means to target cancer. In particular, therapeutic manipulation of tumor-associated macrophages holds promise due to their various and sometimes opposing roles in tumor progression. It is established that M1-type macrophages suppress tumor progression while M2-types support it. Recently, Tie2-expressing macrophages (TEM) have been identified as a distinct sub-population influencing tumor angiogenesis and vascular remodeling as well as monocyte differentiation.

Methods

This study develops a modeling framework to evaluate macrophage interactions with the tumor microenvironment, enabling assessment of how these interactions may affect tumor progression. M1, M2, and Tie2 expressing variants are integrated into a model of tumor growth representing a metastatic lesion in a highly vascularized organ, such as the liver. Behaviors simulated include M1 release of nitric oxide (NO), M2 release of growth-promoting factors, and TEM facilitation of angiogenesis via Angiopoietin-2 and promotion of monocyte differentiation into M2 via IL-10.

Results

The results show that M2 presence leads to larger tumor growth regardless of TEM effects, implying that immunotherapeutic strategies that lead to TEM ablation may fail to restrain growth when the M2 represents a sizeable population. As TEM pro-tumor effects are less pronounced and on a longer time scale than M1-driven tumor inhibition, a more nuanced approach to influence monocyte differentiation taking into account the tumor state (e.g., under chemotherapy) may be desirable.

Conclusions

The results highlight the dynamic interaction of macrophages within a growing tumor, and, further, establish the initial feasibility of a mathematical framework that could longer term help to optimize cancer immunotherapy.

Biotransport kinetics and intratumoral biodistribution of malonodiserinolamide-derivatized [60]fullerene in a murine model of breast adenocarcinoma

[60]Fullerene is a highly versatile nanoparticle (NP) platform for drug delivery to sites of pathology owing to its small size and both ease and versatility of chemical functionalization, facilitating multisite drug conjugation, drug targeting, and modulation of its physicochemical properties. The prominent and well-characterized role of the enhanced permeation and retention (EPR) effect in facilitating NP delivery to tumors motivated us to explore vascular transport kinetics of a water-soluble [60]fullerene derivatives using intravital microscopy in an immune competent murine model of breast adenocarcinoma. Herein, we present a novel local and global image analysis of vascular transport kinetics at the level of individual tumor blood vessels on the micron scale and across whole images, respectively. Similar to larger nanomaterials, [60]fullerenes displayed rapid extravasation from tumor vasculature, distinct from that in normal microvasculature. Temporal heterogeneity in fullerene delivery to tumors was observed, demonstrating the issue of nonuniform delivery beyond spatial dimensions. Trends in local region analysis of fullerene biokinetics by fluorescence quantification were in agreement with global image analysis. Further analysis of intratumoral vascular clearance rates suggested a possible enhanced penetration and retention effect of the fullerene compared to a 70 kDa vascular tracer. Overall, this study demonstrates the feasibility of tracking and quantifying the delivery kinetics and intratumoral biodistribution of fullerene-based drug delivery platforms, consistent with the EPR effect on short timescales and passive transport to tumors.

Development of Halofluorochromic Polymer Nanoassemblies for the Potential Detection of Liver Metastatic Colorectal Cancer Tumors Using Experimental and Computational Approaches

PURPOSE

To develop polymer nanoassemblies (PNAs) modified with halofluorochromic dyes to allow for the detection of liver metastatic colorectal cancer (CRC) to improve therapeutic outcomes.

METHODS

We combine experimental and computational approaches to evaluate macroscopic and microscopic PNA distributions in patient-derived xenograft primary and orthotropic liver metastatic CRC tumors. Halofluorochromic and non-halofluorochromic PNAs (hfPNAs and n-hfPNAs) were prepared from poly(ethylene glycol), fluorescent dyes (Nile blue, Alexa546, and IR820), and hydrophobic groups (palmitate), all of which were covalently tethered to a cationic polymer scaffold [poly(ethylene imine) or poly(lysine)] forming particles with an average diameter < 30 nm.

RESULTS

Dye-conjugated PNAs showed no aggregation under opsonizing conditions for 24 h and displayed low tissue diffusion and cellular uptake. Both hfPNAs and n-hfPNAs accumulated in primary and liver metastatic CRC tumors within 12 h post intravenous injection. In comparison to n-hfPNAs, hfPNAs fluoresced strongly only in the acidic tumor microenvironment (pH < 7.0) and distinguished small metastatic CRC tumors from healthy liver stroma. Computational simulations revealed that PNAs would steadily accumulate mainly in acidic (hypoxic) interstitium of metastatic tumors, independently of the vascularization degree of the tissue surrounding the lesions.

CONCLUSION

The combined experimental and computational data confirms that hfPNAs detecting acidic tumor tissue can be used to identify small liver metastatic CRC tumors with improved accuracy.

Macrophage Polarization Contributes to the Anti-Tumoral Efficacy of Mesoporous Nanovectors Loaded with Albumin-Bound Paclitaxel

Therapies targeted to the immune system, such as immunotherapy, are currently shaping a new, rapidly developing branch of promising cancer treatments, offering the potential to change the prognosis of previously non-responding patients. Macrophages comprise the most abundant population of immune cells in the tumor microenvironment (TME) and can undergo differentiation into functional phenotypes depending on the local tissue environment. Based on these functional phenotypes, tumor-associated macrophages (TAMs) can either aid tumor progression (M2 phenotype) or inhibit it (M1 phenotype). Presence of M2 macrophages and a high ratio of M2/M1 macrophages in the TME are clinically associated with poor prognosis in many types of cancers. Herein, we evaluate the effect of macrophage phenotype on the transport and anti-cancer efficacy of albumin-bound paclitaxel (nAb-PTX) loaded into porous silicon multistage nanovectors (MSV). Studies in a coculture of breast cancer cells (3D-spheroid) with macrophages and in vivo models were conducted to evaluate the therapeutic efficacy of MSV-nAb-PTX as a function of macrophage phenotype. Association with MSV increased drug accumulation within the macrophages and the tumor spheroids, shifting the inflammation state of the TME toward the pro-inflammatory, anti-tumorigenic milieu. Additionally, the treatment increased macrophage motility toward cancer cells, promoting the active transport of therapeutic nanovectors into the tumor lesion. Consequently, apoptosis of cancer cells was increased and proliferation decreased in the MSV-nAb-PTX-treated group as compared to controls. The results also confirmed that the tested system shifts the macrophage differentiation toward an M1 phenotype, possessing an anti-proliferative effect toward the breast cancer cells. These factors were further incorporated into a mathematical model to help analyze the synergistic effect of the macrophage polarization state on the efficacy of MSV-nAb-PTX in alleviating hypovascularized tumor lesions. In conclusion, the ability of MSV-nAb-PTX to polarize TAM to the M1 phenotype, causing (1) enhanced penetration of the drug-carrying macrophages to the center of the tumor lesion and (2) increased toxicity to tumor cells may explain the increased anti-cancer efficacy of the system in comparison to nAb-PTX and other controls.

Pancreatic adenocarcinoma response to chemotherapy enhanced with non-invasive radio frequency evaluated via an integrated experimental/computational approach

Although chemotherapy combined with radiofrequency exposure has shown promise in cancer treatment by coupling drug cytotoxicity with thermal ablation or thermally-induced cytotoxicity, limited access of the drug to tumor loci in hypo-vascularized lesions has hampered clinical application. We recently showed that high-intensity short-wave capacitively coupled radiofrequency (RF) electric-fields may reach inaccessible targets in vivo. This non-invasive RF combined with gemcitabine (Gem) chemotherapy enhanced drug uptake and effect in pancreatic adenocarcinoma (PDAC), notorious for having poor response and limited therapeutic options, but without inducing thermal injury. We hypothesize that the enhanced cytotoxicity derives from RF-facilitated drug transport in the tumor microenvironment. We propose an integrated experimental/computational approach to evaluate chemotherapeutic response combined with RF-induced phenotypic changes in tissue with impaired transport. Results show that RF facilitates diffusive transport in 3D cell cultures representing hypo-vascularized lesions, enhancing drug uptake and effect. Computational modeling evaluates drug vascular extravasation and diffusive transport as key RF-modulated parameters, with transport being dominant. Assessment of hypothetical schedules following current clinical protocol for Stage-IV PDAC suggests that unresponsive lesions may be growth-restrained when exposed to Gem plus RF. Comparison of these projections to experiments in vivo indicates that synergy may result from RF-induced cell phenotypic changes enhancing drug transport and cytotoxicity, thus providing a potential baseline for clinically-focused evaluation.

Imaging the pharmacology of nanomaterials by intravital microscopy: Toward understanding their biological behavior

Therapeutic nanoparticles (NPs) can deliver cytotoxic chemotherapeutics and other drugs more safely and efficiently to patients; furthermore, selective delivery to target tissues can theoretically be accomplished actively through coating NPs with molecular ligands, and passively through exploiting physiological "enhanced permeability and retention" features. However, clinical trial results have been mixed in showing improved efficacy with drug nanoencapsulation, largely due to heterogeneous NP accumulation at target sites across patients. Thus, a clear need exists to better understand why many NP strategies fail in vivo and not result in significantly improved tumor uptake or therapeutic response. Multicolor in vivo confocal fluorescence imaging (intravital microscopy; IVM) enables integrated pharmacokinetic and pharmacodynamic (PK/PD) measurement at the single-cell level, and has helped answer key questions regarding the biological mechanisms of in vivo NP behavior. This review summarizes progress to date and also describes useful technical strategies for successful IVM experimentation.

Multimodality imaging and mathematical modelling of drug delivery to glioblastomas

Patients diagnosed with glioblastoma, an aggressive brain tumour, have a poor prognosis, with a median overall survival of less than 15 months. Vasculature within these tumours is typically abnormal, with increased tortuosity, dilation and disorganization, and they typically exhibit a disrupted blood-brain barrier (BBB). Although it has been hypothesized that the 'normalization' of the vasculature resulting from anti-angiogenic therapies could improve drug delivery through improved blood flow, there is also evidence that suggests that the restoration of BBB integrity might limit the delivery of therapeutic agents and hence their effectiveness. In this paper, we apply mathematical models of blood flow, vascular permeability and diffusion within the tumour microenvironment to investigate the effect of these competing factors on drug delivery. Preliminary results from the modelling indicate that all three physiological parameters investigated-flow rate, vessel permeability and tissue diffusion coefficient-interact nonlinearly to produce the observed average drug concentration in the microenvironment.

Enhanced performance of macrophage-encapsulated nanoparticle albumin-bound-paclitaxel in hypo-perfused cancer lesions

Hypovascularization in tumors such as liver metastases originating from breast and other organs correlates with poor chemotherapeutic response and higher mortality. Poor prognosis is linked to impaired transport of both low- and high-molecular weight drugs into the lesions and to high washout rate. Nanoparticle albumin-bound-paclitaxel (nAb-PTX) has demonstrated benefits in clinical trials when compared to paclitaxel and docetaxel. However, its therapeutic efficacy for breast cancer liver metastasis is disappointing. As macrophages are the most abundant cells in the liver tumor microenvironment, we design a multistage system employing macrophages to deliver drugs into hypovascularized metastatic lesions, and perform in vitro, in vivo, and in silico evaluation. The system encapsulates nAb-PTX into nanoporous biocompatible and biodegradable multistage vectors (MSV), thus promoting nAb-PTX retention in macrophages. We develop a 3D in vitro model to simulate clinically observed hypo-perfused tumor lesions surrounded by macrophages. This model enables evaluation of nAb-PTX and MSV-nab PTX efficacy as a function of transport barriers. Addition of macrophages to this system significantly increases MSV-nAb-PTX efficacy, revealing the role of macrophages in drug transport. In the in vivo model, a significant increase in macrophage number, as compared to unaffected liver, is observed in mice, confirming the in vitro findings. Further, a mathematical model linking drug release and retention from macrophages is implemented to project MSV-nAb-PTX efficacy in a clinical setting. Based on macrophage presence detected via liver tumor imaging and biopsy, the proposed experimental/computational approach could enable prediction of MSV-nab PTX performance to treat metastatic cancer in the liver.

A Computational/Experimental Assessment of Antitumor Activity of Polymer Nanoassemblies for pH-Controlled Drug Delivery to Primary and Metastatic Tumors

PURPOSE

Polymer nanoassemblies (PNAs) with drug release fine-tuned to occur in acidic tumor regions (pH < 7) while sparing normal tissues (pH = 7.4) were previously shown to hold promise as nanoparticle drug carriers to effectively suppress tumor growth with reduced systemic toxicity. However, therapeutic benefits of pH-controlled drug delivery remain elusive due to complex interactions between the drug carriers, tumor cells with varying drug sensitivity, and the tumor microenvironment.

METHODS

We implement a combined computational and experimental approach to evaluate the in vivo antitumor activity of acid-sensitive PNAs controlling drug release in pH 5 ~ 7.4 at different rates [PNA1 (fastest) > PNA2 > PNA3 (slowest)].

RESULTS

Computational simulations projecting the transport, drug release, and antitumor activity of PNAs in primary and metastatic tumor models of colorectal cancer correspond well with experimental observations in vivo. The simulations also reveal that all PNAs could reach peak drug concentrations in tumors at 11 h post injection, while PNAs with slower drug release (PNA2 and PNA3) reduced tumor size more effectively than fast drug releasing PNA1 (24.5 and 20.3 vs 7.5%, respectively, as fraction of untreated control).

CONCLUSION

A combined computational/experimental approach may help to evaluate pH-controlled drug delivery targeting aggressive tumors that have substantial acidity.

Theory and Experimental Validation of a Spatio-temporal Model of Chemotherapy Transport to Enhance Tumor Cell Kill

It has been hypothesized that continuously releasing drug molecules into the tumor over an extended period of time may significantly improve the chemotherapeutic efficacy by overcoming physical transport limitations of conventional bolus drug treatment. In this paper, we present a generalized space- and time-dependent mathematical model of drug transport and drug-cell interactions to quantitatively formulate this hypothesis. Model parameters describe: perfusion and tissue architecture (blood volume fraction and blood vessel radius); diffusion penetration distance of drug (i.e., a function of tissue compactness and drug uptake rates by tumor cells); and cell death rates (as function of history of drug uptake). We performed preliminary testing and validation of the mathematical model using in vivo experiments with different drug delivery methods on a breast cancer mouse model. Experimental data demonstrated a 3-fold increase in response using nano-vectored drug vs. free drug delivery, in excellent quantitative agreement with the model predictions. Our model results implicate that therapeutically targeting blood volume fraction, e.g., through vascular normalization, would achieve a better outcome due to enhanced drug delivery.

Mathematical Based Calculation of Drug Penetration Depth in Solid Tumors

Cancer is a class of diseases characterized by out-of-control cells' growth which affect cells and make them damaged. Many treatment options for cancer exist. Chemotherapy as an important treatment option is the use of drugs to treat cancer. The anticancer drug travels to the tumor and then diffuses in it through capillaries. The diffusion of drugs in the solid tumor is limited by penetration depth which is different in case of different drugs and cancers. The computation of this depth is important as it helps physicians to investigate about treatment of infected tissue. Although many efforts have been made on studying and measuring drug penetration depth, less works have been done on computing this length from a mathematical point of view. In this paper, first we propose phase lagging model for diffusion of drug in the tumor. Then, using this model on one side and considering the classic diffusion on the other side, we compute the drug penetration depth in the solid tumor. This computed value of drug penetration depth is corroborated by comparison with the values measured by experiments.

Intravital Microscopy Imaging Approaches for Image-Guided Drug Delivery Systems

Rapid technical advances in the field of non-linear microscopy have made intravital microscopy a vital pre-clinical tool for research and development of imaging-guided drug delivery systems. The ability to dynamically monitor the fate of macromolecules in live animals provides invaluable information regarding properties of drug carriers (size, charge, and surface coating), physiological, and pathological processes that exist between point-of-injection and the projected of site of delivery, all of which influence delivery and effectiveness of drug delivery systems. In this Review, we highlight how integrating intravital microscopy imaging with experimental designs (in vitro analyses and mathematical modeling) can provide unique information critical in the design of novel disease-relevant drug delivery platforms with improved diagnostic and therapeutic indexes. The Review will provide the reader an overview of the various applications for which intravital microscopy has been used to monitor the delivery of diagnostic and therapeutic agents and discuss some of their potential clinical applications.

Enhanced uptake and transport of PLGA-modified nanoparticles in cervical cancer

BACKGROUND

Uncoordinated cellular proliferation and dysregulated angiogenesis in solid tumors are coupled with inadequate tissue, blood, and lymphatic vascularization. Consequently, tumors are often characterized by hypoxic regions with limited access to vascular-borne substances. In particular, systemically administered nanoparticles (NPs) targeting tumor cells and relying on vascular access to reach tumor tissue can suffer from limited therapeutic efficacy due to inhomogeneous intra-tumor distribution and insufficient cellular internalization of NPs. To circumvent these challenges, NP surfaces can be modified to facilitate tumor interstitial transport and cellular uptake.

RESULTS

We create poly(lactic-co-glycolic) acid NPs modified with MPG, polyethylene glycol (PEG), MPG/PEG, and Vimentin (VIM), and evaluate their cellular uptake in 2D (monolayer) cell culture of human cervical carcinoma (HeLa). We compare NP performance by evaluating uptake by non-cancerous vaginal (VK2) cells. We further assess NP interstitial transport in hypo-vascularized lesions by evaluating the effect of the various modifications on NP penetration in 3D cell culture of the HeLa cells. Results show that after 24 h incubation with HeLa cells in monolayer, MPG, MPG/PEG, PEG, and VIM NPs were internalized at 66×, 24×, 30×, and 15× that of unmodified NPs, respectively. In contrast, incubation with VK2 cells in monolayer showed that MPG , MPG/PEG , PEG , and VIM NPs internalized at 6.3×, 4.3×, 12.4×, and 3.0× that of unmodified NPs, respectively. Uptake was significantly enhanced in tumorigenic vs. normal cells, with internalization of MPG NPs by HeLa cells being twice that of PEG NPs by VK2 cells. After 24 h incubation in HeLa 3D cell culture, MPG and MPG/PEGNPs were internalized 2× and 3× compared to PEG and VIM NPs, respectively. Whereas MPG NPs were internalized mostly in the cell culture periphery (1.2×, 1.4×, and 2.7× that of PEG, MPG/PEG, and VIM NPs, respectively), PEG NPs at 250 μm penetrated 2× farther into the tissue culture than MPG NPs. For all NP types, cellular internalization was severely hindered in 3D compared to monolayer.

CONCLUSIONS

Although MPG surface modification enhances internalization and uptake in hypo-vascularized cervical tissue culture, coating with PEG reduces this internalization while enhancing penetration. A delivery strategy combining NPs with either modification may balance cellular internalization vs. tissue penetration in hypo-vascularized cervical cancer lesions.

An interdisciplinary computational/experimental approach to evaluate drug-loaded gold nanoparticle tumor cytotoxicity

AIM

Clinical translation of cancer nanotherapy has largely failed due to the infeasibility of optimizing the complex interaction of nano/drug/tumor/patient parameters. We develop an interdisciplinary approach modeling diffusive transport of drug-loaded gold nanoparticles in heterogeneously-vascularized tumors.

MATERIALS & METHODS

Evaluated lung cancer cytotoxicity to paclitaxel/cisplatin using novel two-layer (hexadecanethiol/phosphatidylcholine) and three-layer (with high-density-lipoprotein) nanoparticles. Computer simulations calibrated to in-vitro data simulated nanotherapy of heterogeneously-vascularized tumors.

RESULTS

Evaluation of free-drug cytotoxicity between monolayer/spheroid cultures demonstrates a substantial differential, with increased resistance conferred by diffusive transport. Nanoparticles had significantly higher efficacy than free-drug. Simulations of nanotherapy demonstrate 9.5% (cisplatin) and 41.3% (paclitaxel) tumor radius decrease.

CONCLUSION

Interdisciplinary approach evaluating gold nanoparticle cytotoxicity and diffusive transport may provide insight into cancer nanotherapy.

Redirecting Transport of Nanoparticle Albumin-Bound Paclitaxel to Macrophages Enhances Therapeutic Efficacy against Liver Metastases

Current treatments for liver metastases arising from primary breast and lung cancers are minimally effective. One reason for this unfavorable outcome is that liver metastases are poorly vascularized, limiting the ability to deliver therapeutics from the systemic circulation to lesions. Seeking to enhance transport of agents into the tumor microenvironment, we designed a system in which nanoparticle albumin-bound paclitaxel (nAb-PTX) is loaded into a nanoporous solid multistage nanovector (MSV) to enable the passage of the drug through the tumor vessel wall and enhance its interaction with liver macrophages. MSV enablement increased nAb-PTX efficacy and survival in mouse models of breast and lung liver metastasis. MSV-nAb-PTX also augmented the accumulation of paclitaxel and MSV in the liver, specifically in macrophages, whereas paclitaxel levels in the blood were unchanged after administering MSV-nAb-PTX or nAb-PTX. In vitro studies demonstrated that macrophages treated with MSV-nAb-PTX remained viable and were able to internalize, retain, and release significantly higher quantities of paclitaxel compared with treatment with nAb-PTX. The cytotoxic potency of the released paclitaxel was also confirmed in tumor cells cultured with the supernatants of macrophage treated with MSV-nAB-PTX. Collectively, our findings showed how redirecting nAb-PTX to liver macrophages within the tumor microenvironment can elicit a greater therapeutic response in patients with metastatic liver cancer, without increasing systemic side effects.

The Tumor Microenvironment as a Barrier to Cancer Nanotherapy

Although extensive research effort and resources have been dedicated to the development of nanotherapeutics to treat cancer, few formulations have reached clinical application. A major reason is that the large number of parameters available to tune nanotherapy characteristics coupled with the variability in tumor tissue precludes evaluation of complex interactions through experimentation alone. In order to optimize the nanotechnology design and gain further insight into these phenomena, mathematical modeling and computational simulation have been applied to complement empirical work. In this chapter, we discuss modeling work related to nanotherapy and the tumor microenvironment. We first summarize the biology underlying the dysregulated tumor microenvironment, followed by a description of major nano-scale parameters. We then present an overview of the mathematical modeling of cancer nanotherapy, including evaluation of nanotherapy in multi-dimensional tumor tissue, coupling of nanotherapy with vascular flow, modeling of nanotherapy in combination with in vivo imaging, modeling of nanoparticle transport based on in vitro data, modeling of vasculature-bound nanoparticles, evaluation of nanotherapy using pharmacokinetic modeling, and modeling of nano-based hyperthermia. We conclude that an even tighter interdisciplinary effort between biological, material, and physical scientists is needed in order to eventually overcome the tumor microenvironment barrier to successful nanotherapy.

Simulation of the Protein-Shedding Kinetics of a Fully Vascularized Tumor

Circulating biomarkers are of significant interest for cancer detection and treatment personalization. However, the biophysical processes that determine how proteins are shed from cancer cells or their microenvironment, diffuse through tissue, enter blood vasculature, and persist in circulation remain poorly understood. Since approaches primarily focused on experimental evaluation are incapable of measuring the shedding and persistence for every possible marker candidate, we propose an interdisciplinary computational/experimental approach that includes computational modeling of tumor tissue heterogeneity. The model implements protein production, transport, and shedding based on tumor vascularization, cell proliferation, hypoxia, and necrosis, thus quantitatively relating the tumor and circulating proteomes. The results highlight the dynamics of shedding as a function of protein diffusivity and production. Linking the simulated tumor parameters to clinical tumor and vascularization measurements could potentially enable this approach to reveal the tumor-specific conditions based on the protein detected in circulation and thus help to more accurately manage cancer diagnosis and treatment.

Computational Modeling of Tumor Response to Drug Release from Vasculature-Bound Nanoparticles

Systemically injected nanoparticle (NPs) targeting tumor vasculature offer a venue for anti-angiogenic therapies as well as cancer detection and imaging. Clinical application has been limited, however, due to the challenge of elucidating the complex interplay of nanotechnology, drug, and tumor parameters. A critical factor representing the likelihood of endothelial adhesion is the NP vascular affinity, a function of vascular receptor expression and NP size and surface-bound ligand density. We propose a theoretical framework to simulate the tumor response to vasculature-bound drug-loaded NPs and examine the interplay between NP distribution and accumulation as a function of NP vascular affinity, size, and drug loading and release characteristics. The results show that uniform spatial distribution coupled with high vascular affinity is achievable for smaller NPs but not for larger sizes. Consequently, small (100 nm) NPs with high vascular affinity are predicted to be more effective than larger (1000 nm) NPs with similar affinity, even though small NPs have lower drug loading and local drug release compared to the larger NPs. Medium vascular affinity coupled with medium or larger sized NPs is also effective due to a more uniform distribution with higher drug loading and release. Low vascular affinity hampered treatment efficacy regardless of NP size, with larger NPs additionally impeded by heterogeneous distribution and drug release. The results further show that increased drug diffusivity mainly benefits heterogeneously distributed NPs, and would negatively affect efficacy otherwise due to increased wash-out. This model system enables evaluation of efficacy for vascular-targeted drug-loaded NPs as a function of critical NP, drug, and tumor parameters.