Explore the vast range of titles available.

Tumor associated angiogenesis is the development of new blood vessels in response to proteins secreted by tumor cells. These new blood vessels allow tumors to continue to grow beyond what the pre-existing vasculature could support. Here, we construct a mathematical model to simulate tumor angiogenesis by considering each endothelial cell as an agent, and allowing the vascular endothelial growth factor (VEGF) and nutrient fields to impact the dynamics and phenotypic transitions of each tumor and endothelial cell. The phenotypes of the endothelial cells (i.e., tip, stalk, and phalanx cells) are selected by the local VEGF field, and govern the migration and growth of vessel sprouts at the cellular level. Over time, these vessels grow and migrate to the tumor, forming anastomotic loops to supply nutrients, while interacting with the tumor through mechanical forces and the consumption of VEGF. The model is able to capture collapsing and breaking of vessels caused by tumor-endothelial cell interactions. This is accomplished through modeling the physical interaction between the vasculature and the tumor, resulting in vessel occlusion and tumor heterogeneity over time due to the stages of response in angiogenesis. Key parameters are identified through a sensitivity analysis based on the Sobol method, establishing which parameters should be the focus of subsequent experimental efforts. During the avascular phase (i.e., before angiogenesis is triggered), the nutrient consumption rate, followed by the rate of nutrient diffusion, yield the greatest influence on the number and distribution of tumor cells. Similarly, the consumption and diffusion of VEGF yield the greatest influence on the endothelial and tumor cell numbers during angiogenesis. In summary, we present a hybrid mathematical approach that characterizes vascular changes via an agent-based model, while treating nutrient and VEGF changes through a continuum model. The model describes the physical interaction between a tumor and the surrounding blood vessels, explicitly allowing the forces of the growing tumor to influence the nutrient delivery of the vasculature.

While acquired chemoresistance is recognized as a key challenge to treating many types of cancer, the dynamics with which drug sensitivity changes after exposure are poorly characterized. Most chemotherapeutic regimens call for repeated dosing at regular intervals, and if drug sensitivity changes on a similar time scale then the treatment interval could be optimized to improve treatment performance. Theoretical work suggests that such optimal schedules exist, but experimental confirmation has been obstructed by the difficulty of deconvolving the simultaneous processes of death, adaptation, and regrowth taking place in cancer cell populations. Here we present a method of optimizing drug schedules in vitro through iterative application of experimentally calibrated models, and demonstrate its ability to characterize dynamic changes in sensitivity to the chemotherapeutic doxorubicin in three breast cancer cell lines subjected to treatment schedules varying in concentration, interval between pulse treatments, and number of sequential pulse treatments. Cell populations are monitored longitudinally through automated imaging for 600–800 hours, and this data is used to calibrate a family of cancer growth models, each consisting of a system of ordinary differential equations, derived from the bi-exponential model which characterizes resistant and sensitive subpopulations. We identify a model incorporating both a period of growth arrest in surviving cells and a delay in the death of chemosensitive cells which outperforms the original bi-exponential growth model in Akaike Information Criterion based model selection, and use the calibrated model to quantify the performance of each drug schedule. We find that the inter-treatment interval is a key variable in determining the performance of sequential dosing schedules and identify an optimal retreatment time for each cell line which extends regrowth time by 40%-239%, demonstrating that the time scale of changes in chemosensitivity following doxorubicin exposure allows optimization of drug scheduling by varying this inter-treatment interval.

The goal of this study is to calibrate a multiscale model of tumor angiogenesis with time-resolved data to allow for systematic testing of mathematical predictions of vascular sprouting. The multi-scale model consists of an agent-based description of tumor and endothelial cell dynamics coupled to a continuum model of vascular endothelial growth factor concentration. First, we calibrate ordinary differential equation models to time-resolved protein concentration data to estimate the rates of secretion and consumption of vascular endothelial growth factor by endothelial and tumor cells, respectively. These parameters are then input into the multiscale tumor angiogenesis model, and the remaining model parameters are then calibrated to time resolved confocal microscopy images obtained within a 3D vascularized microfluidic platform. The microfluidic platform mimics a functional blood vessel with a surrounding collagen matrix seeded with inflammatory breast cancer cells, which induce tumor angiogenesis. Once the multi-scale model is fully parameterized, we forecast the spatiotemporal distribution of vascular sprouts at future time points and directly compare the predictions to experimentally measured data. We assess the ability of our model to globally recapitulate angiogenic vasculature density, resulting in an average relative calibration error of 17.7% ± 6.3% and an average prediction error of 20.2% ± 4% and 21.7% ± 3.6% using one and four calibrated parameters, respectively. We then assess the model’s ability to predict local vessel morphology (individualized vessel structure as opposed to global vascular density), initialized with the first time point and calibrated with two intermediate time points. In this study, we have rigorously calibrated a mechanism-based, multiscale, mathematical model of angiogenic sprouting to multimodal experimental data to make specific, testable predictions.

High-grade gliomas are an aggressive and invasive malignancy which are susceptible to treatment resistance due to heterogeneity in intratumoral properties such as cell proliferation and density and perfusion. Non-invasive imaging approaches can measure these properties, which can then be used to calibrate patient-specific mathematical models of tumor growth and response. We employed multiparametric magnetic resonance imaging (MRI) to identify tumor extent (via contrast-enhanced T 1 - weighted, and T 2 -FLAIR) and capture intratumoral heterogeneity in cell density (via diffusion-weighted imaging) to calibrate a family of mathematical models of chemoradiation response in nine patients with unresected or partially resected disease. The calibrated model parameters were used to forecast spatially-mapped individual tumor response at future imaging visits. We then employed the Akaike information criteria to select the most parsimonious member from the family, a novel two-species model describing the enhancing and non-enhancing components of the tumor. Using this model, we achieved low error in predictions of the enhancing volume (median: − 2.5%, interquartile range: 10.0%) and a strong correlation in total cell count (Kendall correlation coefficient 0.79) at 3-months post-treatment. These preliminary results demonstrate the plausibility of using multiparametric MRI data to inform spatially-informative, biologically-based predictive models of tumor response in the setting of clinical high-grade gliomas.

Hybrid multiscale agent-based models (ABMs) are unique in their ability to simulate individual cell interactions and microenvironmental dynamics. Unfortunately, the high computational cost of modeling individual cells, the inherent stochasticity of cell dynamics, and numerous model parameters are fundamental limitations of applying such models to predict tumor dynamics. To overcome these challenges, we have developed a coarse-grained two-scale ABM (cgABM) with a reduced parameter space that allows for an accurate and efficient calibration using a set of time-resolved microscopy measurements of cancer cells grown with different initial conditions. The multiscale model consists of a reaction-diffusion type model capturing the spatio-temporal evolution of glucose and growth factors in the tumor microenvironment (at tissue scale), coupled with a lattice-free ABM to simulate individual cell dynamics (at cellular scale). The experimental data consists of BT474 human breast carcinoma cells initialized with different glucose concentrations and tumor cell confluences. The confluence of live and dead cells was measured every three hours over four days. Given this model, we perform a time-dependent global sensitivity analysis to identify the relative importance of the model parameters. The subsequent cgABM is calibrated within a Bayesian framework to the experimental data to estimate model parameters, which are then used to predict the temporal evolution of the living and dead cell populations. To this end, a moment-based Bayesian inference is proposed to account for the stochasticity of the cgABM while quantifying uncertainties due to limited temporal observational data. The cgABM reduces the computational time of ABM simulations by 93% to 97% while staying within a 3% difference in prediction compared to ABM. Additionally, the cgABM can reliably predict the temporal evolution of breast cancer cells observed by the microscopy data with an average error and standard deviation for live and dead cells being 7.61±2.01 and 5.78±1.13, respectively.

We present the development and validation of a mathematical model that predicts how glucose dynamics influence metabolism and therefore tumor cell growth. Glucose, the starting material for glycolysis, has a fundamental influence on tumor cell growth. We employed time-resolved microscopy to track the temporal change of the number of live and dead tumor cells under different initial glucose concentrations and seeding densities. We then constructed a family of mathematical models (where cell death was accounted for differently in each member of the family) to describe overall tumor cell growth in response to the initial glucose and confluence conditions. The Akaikie Information Criteria was then employed to identify the most parsimonious model. The selected model was then trained on 75% of the data to calibrate the system and identify trends in model parameters as a function of initial glucose concentration and confluence. The calibrated parameters were applied to the remaining 25% of the data to predict the temporal dynamics given the known initial glucose concentration and confluence, and tested against the corresponding experimental measurements. With the selected model, we achieved an accuracy (defined as the fraction of measured data that fell within the 95% confidence intervals of the predicted growth curves) of 77.2 ± 6.3% and 87.2 ± 5.1% for live BT-474 and MDA-MB-231 cells, respectively.

This protocol describes a complete data acquisition, analysis and computational forecasting pipeline for employing quantitative MRI data to predict the response of locally advanced breast cancer to neoadjuvant therapy in a community-based care setting. The methodology has previously been successfully applied to a heterogeneous patient population. The protocol details how to acquire the necessary images followed by registration, segmentation, quantitative perfusion and diffusion analysis, model calibration, and prediction. The data collection portion of the protocol requires ~25 min of scanning, postprocessing requires 2–3 h, and the model calibration and prediction components require ~10 h per patient depending on tumor size. The response of individual breast cancer patients to neoadjuvant therapy is forecast by application of a biophysical, reaction–diffusion mathematical model to these data. Successful application of the protocol results in coregistered MRI data from at least two scan visits that quantifies an individual tumor’s size, cellularity and vascular properties. This enables a spatially resolved prediction of how a particular patient’s tumor will respond to therapy. Expertise in image acquisition and analysis, as well as the numerical solution of partial differential equations, is required to carry out this protocol. Quantitative MRI data acquired from patients with locally advanced breast cancer are used to calibrate a biophysical, reaction–diffusion mathematical model to predict response to neoadjuvant therapy on an individual patient basis.

Background Intra-and inter-tumoral heterogeneity in growth dynamics and vascularity influence tumor response to radiation therapy. Quantitative imaging techniques capture these dynamics non-invasively, and these data can initialize and constrain predictive models of response on an individual basis. Methods We have developed a family of 10 biologically-based mathematical models describing the spatiotemporal dynamics of tumor volume fraction, blood volume fraction, and response to radiation therapy. To evaluate this family of models, rats ( n  = 13) with C6 gliomas were imaged with magnetic resonance imaging (MRI) three times before, and four times following a single fraction of 20 Gy or 40 Gy whole brain irradiation. The first five 3D time series data of tumor volume fraction, estimated from diffusion-weighted (DW-) MRI, and blood volume fraction, estimated from dynamic contrast-enhanced (DCE-) MRI, were used to calibrate tumor-specific model parameters. The most parsimonious and well calibrated of the 10 models, selected using the Akaike information criterion, was then utilized to predict future growth and response at the final two imaging time points. Model predictions were compared at the global level (percent error in tumor volume, and Dice coefficient) as well as at the local or voxel level (concordance correlation coefficient). Result The selected model resulted in < 12% error in tumor volume predictions, strong spatial agreement between predicted and observed tumor volumes (Dice coefficient > 0.74), and high level of agreement at the voxel level between the predicted and observed tumor volume fraction and blood volume fraction (concordance correlation coefficient > 0.77 and > 0.65, respectively). Conclusions This study demonstrates that serial quantitative MRI data collected before and following radiation therapy can be used to accurately predict tumor and vasculature response with a biologically-based mathematical model that is calibrated on an individual basis. To the best of our knowledge, this is the first effort to characterize the tumor and vasculature response to radiation therapy temporally and spatially using imaging-driven mathematical models.

The goal of this study is to experimentally and computationally investigate combination trastuzumab-paclitaxel therapies and identify potential synergistic effects due to sequencing of the therapies with in vitro imaging and mathematical modeling. Longitudinal alterations in cell confluence are reported for an in vitro model of BT474 HER2+ breast cancer cells following various dosages and timings of paclitaxel and trastuzumab combination regimens. Results of combination drug regimens are evaluated for drug interaction relationships based on order, timing, and quantity of dose of the drugs. Altering the order of treatments, with the same total therapeutic dose, provided significant changes in overall cell confluence (p < 0.001). Two mathematical models are introduced that are constrained by the in vitro data to simulate the tumor cell response to the individual therapies. A collective model merging the two individual drug response models was designed to investigate the potential mechanisms of synergy for paclitaxel-trastuzumab combinations. This collective model shows increased synergy for regimens where trastuzumab is administered prior to paclitaxel and suggests trastuzumab accelerates the cytotoxic effects of paclitaxel. The synergy derived from the model is found to be in agreement with the combination index, where both indicate a spectrum of additive and synergistic interactions between the two drugs dependent on their dose order. The combined in vitro results and development of a mathematical model of drug synergy has potential to evaluate and improve standard-of-care combination therapies in cancer.

This study identifies physiological tumor habitats from quantitative magnetic resonance imaging (MRI) data and evaluates their alterations in response to therapy. Two models of breast cancer (BT-474 and MDA-MB-231) were imaged longitudinally with diffusion-weighted MRI and dynamic contrast-enhanced MRI to quantify tumor cellularity and vascularity, respectively, during treatment with trastuzumab or albumin-bound paclitaxel. Tumors were stained for anti-CD31, anti-Ki-67, and H&E. Imaging and histology data were clustered to identify tumor habitats and percent tumor volume (MRI) or area (histology) of each habitat was quantified. Histological habitats were correlated with MRI habitats. Clustering of both the MRI and histology data yielded three clusters: high-vascularity high-cellularity (HV-HC), low-vascularity high-cellularity (LV-HC), and low-vascularity low-cellularity (LV-LC). At day 4, BT-474 tumors treated with trastuzumab showed a decrease in LV-HC (p = 0.03) and increase in HV-HC (p = 0.03) percent tumor volume compared to control. MDA-MB-231 tumors treated with low-dose albumin-bound paclitaxel showed a longitudinal decrease in LV-HC percent tumor volume at day 3 (p = 0.01). Positive correlations were found between histological and imaging-derived habitats: HV-HC (BT-474: p = 0.03), LV-HC (MDA-MB-231: p = 0.04), LV-LC (BT-474: p = 0.04; MDA-MB-231: p < 0.01). Physiologically distinct tumor habitats associated with therapeutic response were identified with MRI and histology data in preclinical models of breast cancer.