Multiscale modelling and nonlinear simulation of vascular tumour growth

Personalized predictions of Glioblastoma infiltration: Mathematical models, Physics-Informed Neural Networks and multimodal scans

Predicting the infiltration of Glioblastoma (GBM) from medical MRI scans is crucial for understanding tumor growth dynamics and designing personalized radiotherapy treatment plans. Mathematical models of GBM growth can complement the data in the prediction of spatial distributions of tumor cells. However, this requires estimating patient-specific parameters of the model from clinical data, which is a challenging inverse problem due to limited temporal data and the limited time between imaging and diagnosis. This work proposes a method that uses Physics-Informed Neural Networks (PINNs) to estimate patient-specific parameters of a reaction-diffusion partial differential equation (PDE) model of GBM growth from a single 3D structural MRI snapshot. PINNs embed both the data and the PDE into a loss function, thus integrating theory and data. Key innovations include the identification and estimation of characteristic non-dimensional parameters, a pre-training step that utilizes the non-dimensional parameters and a fine-tuning step to determine the patient specific parameters. Additionally, the diffuse-domain method is employed to handle the complex brain geometry within the PINN framework. The method is validated on both synthetic and patient datasets, showing promise for personalized GBM treatment through parametric inference within clinically relevant timeframes.

Super-resolution ultrasound imaging of ischaemia flow: An in silico study

Super-resolution ultrasound (SRU) is a new ultrasound imaging mode that promises to facilitate the detection of microvascular disease by providing new vascular bio-markers that are directly linked to microvascular pathophysiology, thereby augmenting current knowledge and potentially enabling new treatment. Such a capability can be developed through thorough understanding as articulated by means of mathematical models. In this study, a 2D numerical flow model is adopted for generating flow adaptation in response to ischaemia, in order to determine the ability of SRU to register the resulting flow perturbations. The flow model results demonstrate that variations in flow behaviour in response to locally induced ischaemia can be significant throughout the entire vascular bed. Measured velocities have variations that are dependent on the location of ischaemia, with median values ranging between 2-7 mms-1. Moreover, the distinction between healthy and ischaemic networks are recorded accurately in the SRU results showing excellent agreement between SRU maps and the model. Up to 7-fold spatial resolution improvement to conventional contrast ultrasound was achieved in microbubble localisation while the detection precision and recall was consistently above 98%. The microbubble tracking precision was of a similar accuracy, whereas the recall was reduced (77%) under varying ischaemic impacted flow. Further, regions with velocities up to 30 mms-1 are in excellent agreement with SRU maps, while at regions that include a proportion of higher velocities, the median velocity values are within 1.28%-3.32% of the ground-truth. In conclusion, SRU is a highly promising methodology for the direct measurement of microvascular flow dynamics and may provide a valuable tool for the understanding and subsequent modelling of behaviour in the vascular bed.

Comparison of contrast-enhanced ultrasound imaging (CEUS) and super-resolution ultrasound (SRU) for the quantification of ischaemia flow redistribution: a theoretical study

The study of microcirculation can reveal important information related to pathology. Focusing on alterations that are represented by an obstruction of blood flow in microcirculatory regions may provide an insight into vascular biomarkers. The current in silico study assesses the capability of contrast enhanced ultrasound (CEUS) and super-resolution ultrasound imaging (SRU) flow-quantification to study occlusive actions in a microvascular bed, particularly the ability to characterise known and model induced flow behaviours. The aim is to investigate theoretical limits with the use of CEUS and SRU in order to propose realistic biomarker targets relevant for clinical diagnosis. Results from CEUS flow parameters display limitations congruent with prior investigations. Conventional resolution limits lead to signals dominated by large vessels, making discrimination of microvasculature specific signals difficult. Additionally, some occlusions lead to weakened parametric correlation against flow rate in the remainder of the network. Loss of correlation is dependent on the degree to which flow is redistributed, with comparatively minor redistribution correlating in accordance with ground truth measurements for change in mean transit time,dMTT(CEUS,R = 0.85; GT,R = 0.82) and change in peak intensity,dIp(CEUS,R = 0.87; GT,R = 0.96). Major redistributions, however, result in a loss of correlation, demonstrating that the effectiveness of time-intensity curve parameters is influenced by the site of occlusion. Conversely, results from SRU processing provides accurate depiction of the anatomy and dynamics present in the vascular bed, that extends to individual microvessels. Correspondence between model vessel structure displayed in SRU maps with the ground truth was>91%for cases of minor and major flow redistributions. In conclusion, SRU appears to be a highly promising technology in the quantification of subtle flow phenomena due ischaemia induced vascular flow redistribution.

AMBER: A Modular Model for Tumor Growth, Vasculature and Radiation Response

Computational models of tumor growth are valuable for simulating the dynamics of cancer progression and treatment responses. In particular, agent-based models (ABMs) tracking individual agents and their interactions are useful for their flexibility and ability to model complex behaviors. However, ABMs have often been confined to small domains or, when scaled up, have neglected crucial aspects like vasculature. Additionally, the integration into tumor ABMs of precise radiation dose calculations using gold-standard Monte Carlo (MC) methods, crucial in contemporary radiotherapy, has been lacking. Here, we introduce AMBER, an Agent-based fraMework for radioBiological Effects in Radiotherapy that computationally models tumor growth and radiation responses. AMBER is based on a voxelized geometry, enabling realistic simulations at relevant pre-clinical scales by tracking temporally discrete states stepwise. Its hybrid approach, combining traditional ABM techniques with continuous spatiotemporal fields of key microenvironmental factors such as oxygen and vascular endothelial growth factor, facilitates the generation of realistic tortuous vascular trees. Moreover, AMBER is integrated with TOPAS, an MC-based particle transport algorithm that simulates heterogeneous radiation doses. The impact of radiation on tumor dynamics considers the microenvironmental factors that alter radiosensitivity, such as oxygen availability, providing a full coupling between the biological and physical aspects. Our results show that simulations with AMBER yield accurate tumor evolution and radiation treatment outcomes, consistent with established volumetric growth laws and radiobiological understanding. Thus, AMBER emerges as a promising tool for replicating essential features of tumor growth and radiation response, offering a modular design for future expansions to incorporate specific biological traits.

Jitter solution in parameter identification based on cross-time scale fusion algorithm of lithium-ion batteries

Accurate state-of-charge (SOC) estimation is the core index of battery management system (BMS). When the battery equivalent circuit model (ECM) identifies the parameters under complex operating conditions, there is more jitter or even divergence, which will affect the estimation accuracy of battery SOC. To solve this problem, this paper proposes a new algorithm, namely the cross time scale fusion (CTSF) algorithm. Firstly, the cross-time scales Δ t 1 and Δ t 2 are determined, the number of cross-time cycles is calculated according to the total amount of complex operating condition data N. Then the ECM parameters are identified in Δ t 1 by using forgetting factor recursive least square (FFRLS), and the battery SOC is estimated in Δ t 2 based on the identified parameters, finally the battery parameters are identified and the SOC is estimated by cycling in the cross-time. The experimental results show that, no matter at the same temperature in different conditions or at different temperatures in the same condition, The proposed algorithm not only effectively solves the ECM parameter identification jitter problem, but also improves the accuracy of SOC estimation, the Mean Absolute Error (MAE) minimum of SOC result is 1.42% for different operating conditions at the same temperature and 0.25% for different temperatures at the same operating conditions, respectively.

Selected aspects of avascular tumor growth reproduced by a hybrid model of cell dynamics and chemical kinetics

We here propose a hybrid computational framework to reproduce and analyze aspects of the avascular progression of a generic solid tumor. Our method first employs an individual-based approach to represent the population of tumor cells, which are distinguished in viable and necrotic agents. The active part of the disease is in turn differentiated according to a set of metabolic states. We then describe the spatio-temporal evolution of the concentration of oxygen and of tumor-secreted proteolytic enzymes using partial differential equations (PDEs). A differential equation finally governs the local degradation of the extracellular matrix (ECM) by the malignant mass. Numerical realizations of the model are run to reproduce tumor growth and invasion in a number scenarios that differ for cell properties (adhesiveness, duplication potential, proteolytic activity) and/or environmental conditions (level of tissue oxygenation and matrix density pattern). In particular, our simulations suggest that tumor aggressiveness, in terms of invasive depth and extension of necrotic tissue, can be reduced by (i) stable cell-cell contact interactions, (ii) poor tendency of malignant agents to chemotactically move upon oxygen gradients, and (iii) presence of an overdense matrix, if coupled by a disrupted proteolytic activity of the disease.

A mixture-like model for tumor-immune system interactions

We introduce a mathematical model based on mixture theory intended to describe the tumor-immune system interactions within the tumor microenvironment. The equations account for the geometry of the tumor expansion, and the displacement of the immune cells, driven by diffusion and chemotactic mechanisms. They also take into account the constraints in terms of nutrient and oxygen supply. The numerical investigations analyze the impact of the different modeling assumptions and parameters. Depending on the parameters, the model can reproduce elimination, equilibrium or escape phases and it identifies a critical role of oxygen/nutrient supply in shaping the tumor growth. In addition, antitumor immune cells are key factors in controlling tumor growth, maintaining an equilibrium while protumor cells favor escape and tumor expansion.

"Patchiness" in mechanical stiffness across a tumor as an early-stage marker for malignancy

Mechanical phenotyping of tumors, either at an individual cell level or tumor cell population level is gaining traction as a diagnostic tool. However, the extent of diagnostic and prognostic information that can be gained through these measurements is still unclear. In this work, we focus on the heterogeneity in mechanical properties of cells obtained from a single source such as a tissue or tumor as a potential novel biomarker. We believe that this heterogeneity is a conventionally overlooked source of information in mechanical phenotyping data. We use mechanics-based in-silico models of cell-cell interactions and cell population dynamics within 3D environments to probe how heterogeneity in cell mechanics drives tissue and tumor dynamics. Our simulations show that the initial heterogeneity in the mechanical properties of individual cells and the arrangement of these heterogenous sub-populations within the environment can dictate overall cell population dynamics and cause a shift towards the growth of malignant cell phenotypes within healthy tissue environments. The overall heterogeneity in the cellular mechanotype and their spatial distributions is quantified by a "patchiness" index, which is the ratio of the global to local heterogeneity in cell populations. We observe that there exists a threshold value of the patchiness index beyond which an overall healthy population of cells will show a steady shift towards a more malignant phenotype. Based on these results, we propose that the "patchiness" of a tumor or tissue sample, can be an early indicator for malignant transformation and cancer occurrence in benign tumors or healthy tissues. Additionally, we suggest that tissue patchiness, measured either by biochemical or biophysical markers, can become an important metric in predicting tissue health and disease likelihood just as landscape patchiness is an important metric in ecology.

Computational modeling of angiogenesis: The importance of cell rearrangements during vascular growth

Angiogenesis is the process wherein endothelial cells (ECs) form sprouts that elongate from the pre-existing vasculature to create new vascular networks. In addition to its essential role in normal development, angiogenesis plays a vital role in pathologies such as cancer, diabetes and atherosclerosis. Mathematical and computational modeling has contributed to unraveling its complexity. Many existing theoretical models of angiogenic sprouting are based on the "snail-trail" hypothesis. This framework assumes that leading ECs positioned at sprout tips migrate toward low-oxygen regions while other ECs in the sprout passively follow the leaders' trails and proliferate to maintain sprout integrity. However, experimental results indicate that, contrary to the snail-trail assumption, ECs exchange positions within developing vessels, and the elongation of sprouts is primarily driven by directed migration of ECs. The functional role of cell rearrangements remains unclear. This review of the theoretical modeling of angiogenesis is the first to focus on the phenomenon of cell mixing during early sprouting. We start by describing the biological processes that occur during early angiogenesis, such as phenotype specification, cell rearrangements and cell interactions with the microenvironment. Next, we provide an overview of various theoretical approaches that have been employed to model angiogenesis, with particular emphasis on recent in silico models that account for the phenomenon of cell mixing. Finally, we discuss when cell mixing should be incorporated into theoretical models and what essential modeling components such models should include in order to investigate its functional role. This article is categorized under: Cardiovascular Diseases > Computational Models Cancer > Computational Models.

Enhanced perfusion following exposure to radiotherapy: A theoretical investigation

Tumour angiogenesis leads to the formation of blood vessels that are structurally and spatially heterogeneous. Poor blood perfusion, in conjunction with increased hypoxia and oxygen heterogeneity, impairs a tumour's response to radiotherapy. The optimal strategy for enhancing tumour perfusion remains unclear, preventing its regular deployment in combination therapies. In this work, we first identify vascular architectural features that correlate with enhanced perfusion following radiotherapy, using in vivo imaging data from vascular tumours. Then, we present a novel computational model to determine the relationship between these architectural features and blood perfusion in silico. If perfusion is defined to be the proportion of vessels that support blood flow, we find that vascular networks with small mean diameters and large numbers of angiogenic sprouts show the largest increases in perfusion post-irradiation for both biological and synthetic tumours. We also identify cases where perfusion increases due to the pruning of hypoperfused vessels, rather than blood being rerouted. These results indicate the importance of considering network composition when determining the optimal irradiation strategy. In the future, we aim to use our findings to identify tumours that are good candidates for perfusion enhancement and to improve the efficacy of combination therapies.

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.

Numerical investigation of the pharmacokinetics and pharmacodynamics of the chemotherapeutic drug in avascular and vascular stages of a brain tumor

One of the most commonly used approaches for treating solid tumors is the systemic delivery of chemotherapeutic drugs. However, our understanding of the factors influencing treatment efficacy through this method is still limited. This study presents a comprehensive and realistic mathematical model that incorporates the dynamics of tumor growth, capillary network extension, and drug delivery in a coupled and simultaneous manner. The model covers two stages of tumor growth: avascular and vascular. For tumor growth, a continuum model is employed using the phase field interface-capturing method. The neo-vascularization process is modeled using a hybrid discrete-continuum approach. Additionally, a multi-scale model is used to describe the pharmacokinetics of doxorubicin, considering various agents. The study investigates the effect of haptotaxis and reveals that a higher haptotaxis coefficient leads to faster tumor growth (up to 2.6 times) and a quicker progression to angiogenesis. The impact of tumor-related and drug-related parameters is also examined, including tumor size, tumor sensitivity to the drug, chemotherapy initialization, treatment cycle duration, drug affinity to cells, and drug dose. The findings indicate that chemotherapy is more effective during the angiogenesis stage when active loops have formed. Other clinical methods such as radiotherapy and surgery may be more appropriate during the avascular stage or the transition period between angiogenesis initialization and loop formation. The penetration depth of the drug decreases by approximately 50% with an increase in the drug binding rate to surface-cell receptors. As a result, high-associate-rate drugs are preferred for chemotherapy after active loops have formed, while low-associate-rate drugs are suitable for earlier stages.

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.

A proof-of-concept pipeline to guide evaluation of tumor tissue perfusion by dynamic contrast-agent imaging: Direct simulation and inverse tracer-kinetic procedures

Dynamic contrast-enhanced (DCE) perfusion imaging has shown great potential to non-invasively assess cancer development and its treatment by their characteristic tissue signatures. Different tracer kinetics models are being applied to estimate tissue and tumor perfusion parameters from DCE perfusion imaging. The goal of this work is to provide an in silico model-based pipeline to evaluate how these DCE imaging parameters may relate to the true tissue parameters. As histology data provides detailed microstructural but not functional parameters, this work can also help to better interpret such data. To this aim in silico vasculatures are constructed and the spread of contrast agent in the tissue is simulated. As a proof of principle we show the evaluation procedure of two tracer kinetic models from in silico contrast-agent perfusion data after a bolus injection. Representative microvascular arterial and venous trees are constructed in silico. Blood flow is computed in the different vessels. Contrast-agent input in the feeding artery, intra-vascular transport, intra-extravascular exchange and diffusion within the interstitial space are modeled. From this spatiotemporal model, intensity maps are computed leading to in silico dynamic perfusion images. Various tumor vascularizations (architecture and function) are studied and show spatiotemporal contrast imaging dynamics characteristic of in vivo tumor morphotypes. The Brix II also called 2CXM, and extended Tofts tracer-kinetics models common in DCE imaging are then applied to recover perfusion parameters that are compared with the ground truth parameters of the in silico spatiotemporal models. The results show that tumor features can be well identified for a certain permeability range. The simulation results in this work indicate that taking into account space explicitly to estimate perfusion parameters may lead to significant improvements in the perfusion interpretation of the current tracer-kinetics models.

On the role of mechanical signals on sprouting angiogenesis through computer modeling approaches

Sprouting angiogenesis, the formation of new vessels from preexisting vasculature, is an essential process in the regeneration of new tissues as well as in the development of some diseases like cancer. Although early studies identified chemical signaling as the main driver of this process, many recent studies have shown a strong role of mechanical signals in the formation of new capillaries. Different types of mechanical signals (e.g., external forces, cell traction forces, and blood flow-induced shear forces) have been shown to play distinct roles in the process; however, their interplay remains still largely unknown. During the last decades, mathematical and computational modeling approaches have been developed to investigate and better understand the mechanisms behind mechanically driven angiogenesis. In this manuscript, we review computational models of angiogenesis with a focus on models investigating the role of mechanics on the process. Our aim is not to provide a detailed review on model methodology but to describe what we have learnt from these models. We classify models according to the mechanical signals being investigated and describe how models have looked into their role on the angiogenic process. We show that a better understanding of the mechanobiology of the angiogenic process will require the development of computer models that incorporate the interactions between the multiple mechanical signals and their effect on cellular responses, since they all seem to play a key in sprout patterning. In the end, we describe some of the remaining challenges of computational modeling of angiogenesis and discuss potential avenues for future research.

The impact of tumor associated macrophages on tumor biology under the lens of mathematical modelling: A review

In this article, we review the role of mathematical modelling to elucidate the impact of tumor-associated macrophages (TAMs) in tumor progression and therapy design. We first outline the biology of TAMs, and its current application in tumor therapies, and their experimental methods that provide insights into tumor cell-macrophage interactions. We then focus on the mechanistic mathematical models describing the role of macrophages as drug carriers, the impact of macrophage polarized activation on tumor growth, and the role of tumor microenvironment (TME) parameters on the tumor-macrophage interactions. This review aims to identify the synergies between biological and mathematical approaches that allow us to translate knowledge on fundamental TAMs biology in addressing current clinical challenges.

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.

Inverse problem approaches for mutation laws in heterogeneous tumours with local and nonlocal dynamics

Cancer cell mutations occur when cells undergo multiple cell divisions, and these mutations can be spontaneous or environmentally-induced. The mechanisms that promote and sustain these mutations are still not fully understood. This study deals with the identification (or reconstruction) of the usually unknown cancer cell mutation law, which lead to the transformation of a primary tumour cell population into a secondary, more aggressive cell population. We focus on local and nonlocal mathematical models for cell dynamics and movement, and identify these mutation laws from macroscopic tumour snapshot data collected at some later stage in the tumour evolution. In a local cancer invasion model, we first reconstruct the mutation law when we assume that the mutations depend only on the surrounding cancer cells (i.e., the ECM plays no role in mutations). Second, we assume that the mutations depend on the ECM only, and we reconstruct the mutation law in this case. Third, we reconstruct the mutation when we assume that there is no prior knowledge about the mutations. Finally, for the nonlocal cancer invasion model, we reconstruct the mutation law that depends on the cancer cells and on the ECM. For these numerical reconstructions, our approximations are based on the finite difference method combined with the finite elements method. As the inverse problem is ill-posed, we use the Tikhonov regularisation technique in order to regularise the solution. Stability of the solution is examined by adding additive noise into the measurements.

A multiscale cell-based model of tumor growth for chemotherapy assessment and tumor-targeted therapy through a 3D computational approach

OBJECTIVES

Computational modeling of biological systems is a powerful tool to clarify diverse processes contributing to cancer. The aim is to clarify the complex biochemical and mechanical interactions between cells, the relevance of intracellular signaling pathways in tumor progression and related events to the cancer treatments, which are largely ignored in previous studies.

MATERIALS AND METHODS

A three-dimensional multiscale cell-based model is developed, covering multiple time and spatial scales, including intracellular, cellular, and extracellular processes. The model generates a realistic representation of the processes involved from an implementation of the signaling transduction network.

RESULTS

Considering a benign tumor development, results are in good agreement with the experimental ones, which identify three different phases in tumor growth. Simulating tumor vascular growth, results predict a highly vascularized tumor morphology in a lobulated form, a consequence of cells' motile behavior. A novel systematic study of chemotherapy intervention, in combination with targeted therapy, is presented to address the capability of the model to evaluate typical clinical protocols. The model also performs a dose comparison study in order to optimize treatment efficacy and surveys the effect of chemotherapy initiation delays and different regimens.

CONCLUSIONS

Results not only provide detailed insights into tumor progression, but also support suggestions for clinical implementation. This is a major step toward the goal of predicting the effects of not only traditional chemotherapy but also tumor-targeted therapies.

Biologically-Based Mathematical Modeling of Tumor Vasculature and Angiogenesis via Time-Resolved Imaging Data

Tumor-associated vasculature is responsible for the delivery of nutrients, removal of waste, and allowing growth beyond 2-3 mm3. Additionally, the vascular network, which is changing in both space and time, fundamentally influences tumor response to both systemic and radiation therapy. Thus, a robust understanding of vascular dynamics is necessary to accurately predict tumor growth, as well as establish optimal treatment protocols to achieve optimal tumor control. Such a goal requires the intimate integration of both theory and experiment. Quantitative and time-resolved imaging methods have emerged as technologies able to visualize and characterize tumor vascular properties before and during therapy at the tissue and cell scale. Parallel to, but separate from those developments, mathematical modeling techniques have been developed to enable in silico investigations into theoretical tumor and vascular dynamics. In particular, recent efforts have sought to integrate both theory and experiment to enable data-driven mathematical modeling. Such mathematical models are calibrated by data obtained from individual tumor-vascular systems to predict future vascular growth, delivery of systemic agents, and response to radiotherapy. In this review, we discuss experimental techniques for visualizing and quantifying vascular dynamics including magnetic resonance imaging, microfluidic devices, and confocal microscopy. We then focus on the integration of these experimental measures with biologically based mathematical models to generate testable predictions.

A multi-scale agent-based model for avascular tumour growth

In this paper, we have developed a multi-scale, lattice-free, agent based model of avascular tumour growth in epithelial tissue. The model integrates different events to identify the underlying diversity within intracellular, cellular, and extracellular layer dynamics. The model considers every cell as an agent. A cellular agent may proliferate, spawns two identical daughter agents, or it may be transformed into other phenotypes during its life time depending on its internal proteins' activity as well as its external microenvironment. In this context, a simplified age-structured cell cycle model is adopted from the existing literature. The model considers that the intracellular events are regulated by p27 gene expression. In this model, p27 protein controls the overall tumour growth dynamics. Moreover, p27 is controlled by the external oxygen and nutrients that are modelled with the reaction-diffusion equations. The model also considers several biophysical forces which directly effect on the tumour growth dynamics. This modelling framework offers biologically realistic outcomes and also covers important criteria of the hallmarks of cancer which include oxygen and nutrient consumptions, micro-environmental heterogeneity, tumour cell proliferation by avoiding growth suppressor signals, replication of tumour cells at an abnormally faster rate, and resistance of apoptosis. The avascular tumour growth model is validated with immunohistochemistry and histopathology data. The outcome of the proposed model is very close to the range of the patient data, which concludes that the model is capable enough to mimic these complex biophysical phenomena.

Phenotypic variation modulates the growth dynamics and response to radiotherapy of solid tumours under normoxia and hypoxia

In cancer, treatment failure and disease recurrence have been associated with small subpopulations of cancer cells with a stem-like phenotype. In this paper, we develop and investigate a phenotype-structured model of solid tumour growth in which cells are structured by a stemness level, which varies continuously between stem-like and terminally differentiated behaviours. Cell evolution is driven by proliferation and death, as well as advection and diffusion with respect to the stemness structure variable. Here, the magnitude and sign of the advective flux are allowed to vary with the oxygen level. We use the model to investigate how the environment, in particular oxygen levels, affects the tumour's population dynamics and composition, and its response to radiotherapy. We use a combination of numerical and analytical techniques to quantify how under physiological oxygen levels the cells evolve to a differentiated phenotype and under low oxygen level (i.e., hypoxia) they de-differentiate. Under normoxia, the proportion of cancer stem cells is typically negligible and the tumour may ultimately become extinct whereas under hypoxia cancer stem cells comprise a dominant proportion of the tumour volume, enhancing radio-resistance and favouring the tumour's long-term survival. We then investigate how such phenotypic heterogeneity impacts the tumour's response to treatment with radiotherapy under normoxia and hypoxia. Of particular interest is establishing how the presence of radio-resistant cancer stem cells can facilitate a tumour's regrowth following radiotherapy. We also use the model to show how radiation-induced changes in tumour oxygen levels can give rise to complex re-growth dynamics. For example, transient periods of hypoxia induced by damage to tumour blood vessels may rescue the cancer cell population from extinction and drive secondary regrowth.

Multiscale modeling of tumor growth and angiogenesis: Evaluation of tumor-targeted therapy

The dynamics of tumor growth and associated events cover multiple time and spatial scales, generally including extracellular, cellular and intracellular modifications. The main goal of this study is to model the biological and physical behavior of tumor evolution in presence of normal healthy tissue, considering a variety of events involved in the process. These include hyper and hypoactivation of signaling pathways during tumor growth, vessels' growth, intratumoral vascularization and competition of cancer cells with healthy host tissue. The work addresses two distinctive phases in tumor development-the avascular and vascular phases-and in each stage two cases are considered-with and without normal healthy cells. The tumor growth rate increases considerably as closed vessel loops (anastomoses) form around the tumor cells resulting from tumor induced vascularization. When taking into account the host tissue around the tumor, the results show that competition between normal cells and cancer cells leads to the formation of a hypoxic tumor core within a relatively short period of time. Moreover, a dense intratumoral vascular network is formed throughout the entire lesion as a sign of a high malignancy grade, which is consistent with reported experimental data for several types of solid carcinomas. In comparison with other mathematical models of tumor development, in this work we introduce a multiscale simulation that models the cellular interactions and cell behavior as a consequence of the activation of oncogenes and deactivation of gene signaling pathways within each cell. Simulating a therapy that blocks relevant signaling pathways results in the prevention of further tumor growth and leads to an expressive decrease in its size (82% in the simulation).

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.

Measuring Nanoparticle Penetration Through Bio-Mimetic Gels

Background

In cancer nanomedicine, drugs are transported by nanocarriers through a biological system to produce a therapeutic effect. The efficacy of the treatment is affected by the ability of the nanocarriers to overcome biological transport barriers to reach their target. In this work, we focus on the process of nanocarrier penetration through tumour tissue after extravasation. Visualising the dynamics of nanocarriers in tissue is difficult in vivo, and in vitro assays often do not capture the spatial and physical constraints relevant to model tissue penetration.

Methods

We propose a new simple, low-cost method to observe the transport dynamics of nanoparticles through a tissue-mimetic microfluidic chip. After loading a chip with triplicate conditions of gel type and loading with microparticles, microscopic analysis allows for tracking of fluorescent nanoparticles as they move through hydrogels (Matrigel and Collagen I) with and without cell-sized microparticles. A bespoke image-processing codebase written in MATLAB allows for statistical analysis of this tracking, and time-dependent dynamics can be determined.

Results

To demonstrate the method, we show size-dependence of transport mechanics can be observed, with diffusion of fluorescein dye throughout the channel in 8 h, while 20 nm carboxylate FluoSphere diffusion was hindered through both Collagen I and Matrigel™. Statistical measurements of the results are generated through the software package and show the significance of both size and presence of microparticles on penetration depth.

Conclusion

This provides an easy-to-understand output for the end user to measure nanoparticle tissue penetration, enabling the first steps towards future automated experimentation of transport dynamics for rational nanocarrier design.

Cell and extracellular matrix growth theory and its implications for tumorigenesis

Cells associated with an abnormal (cancerous) growth exchange flows, morph freely and grow hand-in-glove with their immediate environment, the extracellular matrix (ECM). The cell structure experiences two mass flows in counterflow. Flowing into the structure are nutrients and flowing out is refuse from the metabolically active biomass within. The physical effect of the evolution of the cell and extracellular structure is more flow and mixing in that space, that is, more mixing than in the absence of a biological growth in that space. The objective of the present theory is to predict the increase in the size of the cell cluster as a function of its structure, and also to predict the critical cluster sizes that mark the transitions from one distinct cluster configuration to the next. This amounts to predicting the timing and the main features of the transitions from single cell to clusters with two, four, eight and more cells, including larger clusters with cells organized on its outer surface. The predicted evolution of the size and configuration of the cell cluster is validated successfully by comparison with measurements from several independent studies of cancerous and non-cancerous growth patterns.

Mathematical simulation of tumour angiogenesis: angiopoietin balance is a key factor in vessel growth and regression

Excessive tumour growth results in a hypoxic environment around cancer cells, thus inducing tumour angiogenesis, which refers to the generation of new blood vessels from pre-existing vessels. This mechanism is biologically and physically complex, with various mathematical simulation models proposing to reproduce its formation. However, although temporary vessel regression is clinically known, few models succeed in reproducing this phenomenon. Here, we developed a three-dimensional simulation model encompassing both angiogenesis and tumour growth, specifically including angiopoietin. Angiopoietin regulates both adhesion and migration between vascular endothelial cells and wall cells, thus inhibiting the cell-to-cell adhesion required for angiogenesis initiation. Simulation results showed a regression, i.e. transient decrease, in the overall length of new vessels during vascular network formation. Using our model, we also evaluated the efficacy of administering the drug bevacizumab. The results highlighted differences in treatment efficacy: (1) earlier administration showed higher efficacy in inhibiting tumour growth, and (2) efficacy depended on the treatment interval even with the administration of the same dose. After thorough validation in the future, these results will contribute to the design of angiogenesis treatment protocols.

A multiscale model of complex endothelial cell dynamics in early angiogenesis

We introduce a hybrid two-dimensional multiscale model of angiogenesis, the process by which endothelial cells (ECs) migrate from a pre-existing vascular bed in response to local environmental cues and cell-cell interactions, to create a new vascular network. Recent experimental studies have highlighted a central role of cell rearrangements in the formation of angiogenic networks. Our model accounts for this phenomenon via the heterogeneous response of ECs to their microenvironment. These cell rearrangements, in turn, dynamically remodel the local environment. The model reproduces characteristic features of angiogenic sprouting that include branching, chemotactic sensitivity, the brush border effect, and cell mixing. These properties, rather than being hardwired into the model, emerge naturally from the gene expression patterns of individual cells. After calibrating and validating our model against experimental data, we use it to predict how the structure of the vascular network changes as the baseline gene expression levels of the VEGF-Delta-Notch pathway, and the composition of the extracellular environment, vary. In order to investigate the impact of cell rearrangements on the vascular network structure, we introduce the mixing measure, a scalar metric that quantifies cell mixing as the vascular network grows. We calculate the mixing measure for the simulated vascular networks generated by ECs of different lineages (wild type cells and mutant cells with impaired expression of a specific receptor). Our results show that the time evolution of the mixing measure is directly correlated to the generic features of the vascular branching pattern, thus, supporting the hypothesis that cell rearrangements play an essential role in sprouting angiogenesis. Furthermore, we predict that lower cell rearrangement leads to an imbalance between branching and sprout elongation. Since the computation of this statistic requires only individual cell trajectories, it can be computed for networks generated in biological experiments, making it a potential biomarker for pathological angiogenesis.

Angiogenesis, the process by which new blood vessels are formed by sprouting from the pre-existing vascular bed, plays a key role in both physiological and pathological processes, including tumour growth. The structure of a growing vascular network is determined by the coordinated behaviour of endothelial cells in response to various signalling cues. Recent experimental studies have highlighted the importance of cell rearrangements as a driver for sprout elongation. However, the functional role of this phenomenon remains unclear. We formulate a new multiscale model of angiogenesis which, by accounting explicitly for the complex dynamics of endothelial cells within growing angiogenic sprouts, is able to reproduce generic features of angiogenic structures (branching, chemotactic sensitivity, cell mixing, etc.) as emergent properties of its dynamics. We validate our model against experimental data and then use it to quantify the phenomenon of cell mixing in vascular networks generated by endothelial cells of different lineages. Our results show that there is a direct correlation between the time evolution of cell mixing in a growing vascular network and its branching structure, thus paving the way for understanding the functional role of cell rearrangements in angiogenesis.

Digitally Fabricated and Naturally Augmented In Vitro Tissues

Human in vitro tissues are extracorporeal 3D cultures of human cells embedded in biomaterials, commonly hydrogels, which recapitulate the heterogeneous, multiscale, and architectural environment of the human body. Contemporary strategies used in 3D tissue and organ engineering integrate the use of automated digital manufacturing methods, such as 3D printing, bioprinting, and biofabrication. Human tissues and organs, and their intra- and interphysiological interplay, are particularly intricate. For this reason, attentiveness is rising to intersect materials science, medicine, and biology with arts and informatics. This report presents advances in computational modeling of bioink polymerization and its compatibility with bioprinting, the use of digital design and fabrication in the development of fluidic culture devices, and the employment of generative algorithms for modeling the natural and biological augmentation of in vitro tissues. As a future direction, the use of serially linked in vitro tissues as human body-mimicking systems and their application in drug pharmacokinetics and metabolism, disease modeling, and diagnostics are discussed.

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.

Indirect Contributions to Tumor Dynamics in the First Stage of the Avascular Phase

A continuum model for tumor invasion in a two-dimensional spatial domain based on the interaction of the urokinase plasminogen activation system with a model for cancer cell dynamics is proposed. The arising system of partial differential equations is numerically solved using the finite element method. We simulated a portion of biological tissue imposing no flux boundary conditions. We monitored the cancer cell dynamics, as well the degradation of an extra cellular matrix representative, vitronectin, and the evolution of a specific degrading enzyme, plasmin, inside the biological tissue. The computations were parameterized as a function of the indirect cell proliferation induced by a plasminogen activator inhibitor binding to vitronectin and of the indirect plasmin deactivation due to the plasminogen activator inhibitor binding to the urokinase plasminogen activator. Their role during the cancer dynamical evolution was identified, together with a possible marker helping the mapping of the cancer invasive front. Our results indicate that indirect cancer cell proliferation biases the speed of the tumor invasive front as well as the heterogeneity of the cancer cell clustering and networking, as it ultimately acts on the proteolytic activity supporting cancer formation. Because of the initial conditions imposed, the numerical solutions of the model show a symmetrical dynamical evolution of heterogeneities inside the simulated domain. Moreover, an increase of up to about 12% in the invasion speed was observed, increasing the rate of indirect cancer cell proliferation, while increasing the plasmin deactivation rate inhibits heterogeneities and networking. As cancer cell proliferation causes vitronectin consumption and plasmin formation, the intensities of the concentration maps of both vitronectin and plasmin are superimposable to the cancer cell concentration maps. The qualitative imprinting that cancer cells leave on the extra cellular matrix during the time evolution as well their activity area is identified, framing the numerical results in the context of a methodology aimed at diagnostic and therapeutic improvement.

Tumor Growth Dynamic Modeling in Oncology Drug Development and Regulatory Approval: Past, Present, and Future Opportunities

Model-informed drug development (MIDD) approaches have rapidly advanced in drug development in recent years. Additionally, the Prescription Drug User Fee Act (PDUFA) VI has specific commitments to further enhance MIDD. Tumor growth dynamic (TGD) modeling, as one of the commonly utilized MIDD approaches in oncology, fulfills the purposes to accelerate the drug development, to support new drug and biologics license applications, and to guide the market access. Increasing knowledge of TGD modeling methodologies, encouraging applications in clinical setting for patients' survival, and complementing assessment of regulatory review for submissions, together fueled promising potentials for imminent enhancement of TGD in oncology. This review is to comprehensively summarize the history of TGD, and present case examples of the recent advance of TGD modeling (mixture model and joint model), as well as the TGD impact on regulatory decisions, thus illustrating challenges and opportunities. Additionally, this review presents the future perspectives for TGD approach.

Blackboard to Bedside: A Mathematical Modeling Bottom-Up Approach Toward Personalized Cancer Treatments

Cancers present with high variability across patients and tumors; thus, cancer care, in terms of disease prevention, detection, and control, can highly benefit from a personalized approach. For a comprehensive personalized oncology practice, this personalization should ideally consider data gathered from various information levels, which range from the macroscale population level down to the microscale tumor level, without omission of the central patient level. Appropriate data mined from each of these levels can significantly contribute in devising personalized treatment plans tailored to the individual patient and tumor. Mathematical models of solid tumors, combined with patient-specific tumor profiles, present a unique opportunity to personalize cancer treatments after detection using a bottom-up approach. Here, we discuss how information harvested from mathematical models and from corresponding in silico experiments can be implemented in preclinical and clinical applications. To conceptually illustrate the power of these models, one such model is presented, and various pertinent tumor and treatment scenarios are demonstrated in silico. The presented model, specifically a multiscale, hybrid cellular automaton, has been fully validated in vitro using multiple cell-line-specific data. We discuss various insights provided by this model and other models like it and their role in designing predictive tools that are both patient, and tumor specific. After refinement and parametrization with appropriate data, such in silico tools have the potential to be used in a clinical setting to aid in treatment protocols and decision making.

A Mathematical Dissection of the Adaptation of Cell Populations to Fluctuating Oxygen Levels

The disordered network of blood vessels that arises from tumour angiogenesis results in variations in the delivery of oxygen into the tumour tissue. This brings about regions of chronic hypoxia (i.e. sustained low oxygen levels) and regions with alternating periods of low and relatively higher oxygen levels, and makes it necessary for cancer cells to adapt to fluctuating environmental conditions. We use a phenotype-structured model to dissect the evolutionary dynamics of cell populations exposed to fluctuating oxygen levels. In this model, the phenotypic state of every cell is described by a continuous variable that provides a simple representation of its metabolic phenotype, ranging from fully oxidative to fully glycolytic, and cells are grouped into two competing populations that undergo heritable, spontaneous phenotypic variations at different rates. Model simulations indicate that, depending on the rate at which oxygen is consumed by the cells, dynamic nonlinear interactions between cells and oxygen can stimulate chronic hypoxia and cycling hypoxia. Moreover, the model supports the idea that under chronic-hypoxic conditions lower rates of phenotypic variation lead to a competitive advantage, whereas higher rates of phenotypic variation can confer a competitive advantage under cycling-hypoxic conditions. In the latter case, the numerical results obtained show that bet-hedging evolutionary strategies, whereby cells switch between oxidative and glycolytic phenotypes, can spontaneously emerge. We explain how these results can shed light on the evolutionary process that may underpin the emergence of phenotypic heterogeneity in vascularised tumours.

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.

Mathematical Modeling Shows That the Response of a Solid Tumor to Antiangiogenic Therapy Depends on the Type of Growth

It has been hypothesized that solid tumors with invasive type of growth should possess intrinsic resistance to antiangiogenic therapy, which is aimed at cessation of the formation of new blood vessels and subsequent shortage of nutrient inflow to the tumor. In order to investigate this effect, a continuous mathematical model of tumor growth is developed, which considers variables of tumor cells, necrotic tissue, capillaries, and glucose as the crucial nutrient. The model accounts for the intrinsic motility of tumor cells and for the convective motion, arising due to their proliferation, thus allowing considering two types of tumor growth—invasive and compact—as well as their combination. Analytical estimations of tumor growth speed are obtained for compact and invasive tumors. They suggest that antiangiogenic therapy may provide a several times decrease of compact tumor growth speed, but the decrease of growth speed for invasive tumors should be only modest. These estimations are confirmed by numerical simulations, which further allow evaluating the effect of antiangiogenic therapy on tumors with mixed growth type and highlight the non-additive character of the two types of growth.

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.

Gibbs point field model quantifies disorder in microvasculature of U87-glioblastoma

Microvascular proliferation in glioblastoma multiforme is a biological key mechanism to facilitate tumor growth and infiltration and a main target for treatment interventions. The vascular architecture can be obtained by Single Plane Illumination Microscopy (SPIM) to evaluate vascular heterogeneity in tumorous tissue. We make use of the Gibbs point field model to quantify the order of regularity in capillary distributions found in the U87 glioblastoma model in a murine model and to compare tumorous and healthy brain tissue. A single model parameter Γ was assigned that is linked to tissue-specific vascular topology through Monte-Carlo simulations. Distributions of the model parameter Γ differ significantly between glioblastoma tissue with mean 〈Γ_G〉=2.1±0.4, as compared to healthy brain tissue with mean 〈Γ_H〉=4.9±0.4, suggesting that the average Γ-value allows for tissue differentiation. These results may be used for diagnostic magnetic resonance imaging, where it has been shown recently that Γ is linked to tissue-inherent relaxation parameters.

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.

Extension of a multiphase tumour growth model to study nanoparticle delivery to solid tumours

One of the main challenges in increasing the efficacy of conventional chemotherapeutics is the fact that they do not reach cancerous cells at a sufficiently high dosage. In order to remedy this deficiency, nanoparticle-based drugs have evolved as a promising novel approach to more specific tumour targeting. Nevertheless, several biophysical phenomena prevent the sufficient penetration of nanoparticles in order to target the entire tumour. We therefore extend our vascular multiphase tumour growth model, enabling it to investigate the influence of different biophysical factors on the distribution of nanoparticles in the tumour microenvironment. The novel model permits the examination of the interplay between the size of vessel-wall pores, the permeability of the blood-vessel endothelium and the lymphatic drainage on the delivery of particles of different sizes. Solid tumours develop a non-perfused core and increased interstitial pressure. Our model confirms that those two typical features of solid tumours limit nanoparticle delivery. Only in case of small nanoparticles is the transport dominated by diffusion, and particles can reach the entire tumour. The size of the vessel-wall pores and the permeability of the blood-vessel endothelium have a major impact on the amount of delivered nanoparticles. This extended in-silico tumour growth model permits the examination of the characteristics and of the limitations of nanoparticle delivery to solid tumours, which currently complicate the translation of nanoparticle therapy to a clinical stage.

Properties of the Vascular Networks in Malignant Tumors

This work presents an analysis for real and synthetic angiogenic networks using a tomography image that obtains a portrait of a vascular network. After the image conversion into a binary format it is possible to measure various network properties, which includes the average path length, the clustering coefficient, the degree distribution and the fractal dimension. When comparing the observed properties with that produced by the Invasion Percolation algorithm (IPA), we observe that there exist differences between the properties obtained by the real and the synthetic networks produced by the IPA algorithm. Taking into account the former, a new algorithm which models the expansion of an angiogenic network through randomly heuristic rules is proposed. When comparing this new algorithm with the real networks it is observed that now both share some properties. Once creating synthetic networks, we prove the robustness of the network by subjecting the original angiogenic and the synthetic networks to the removal of the most connected nodes, and see to what extent the properties changed. Using this concept of robustness, in a very naive fashion it is possible to launch a hypothetical proposal for a therapeutic treatment based on the robustness of the network.

Notch signaling and taxis mechanisms regulate early stage angiogenesis: A mathematical and computational model

During angiogenesis, new blood vessels sprout and grow from existing ones. This process plays a crucial role in organ development and repair, in wound healing and in numerous pathological processes such as cancer progression or diabetes. Here, we present a mathematical model of early stage angiogenesis that permits exploration of the relative importance of mechanical, chemical and cellular cues. Endothelial cells proliferate and move over an extracellular matrix by following external gradients of Vessel Endothelial Growth Factor, adhesion and stiffness, which are incorporated to a Cellular Potts model with a finite element description of elasticity. The dynamics of Notch signaling involving Delta-4 and Jagged-1 ligands determines tip cell selection and vessel branching. Through their production rates, competing Jagged-Notch and Delta-Notch dynamics determine the influence of lateral inhibition and lateral induction on the selection of cellular phenotypes, branching of blood vessels, anastomosis (fusion of blood vessels) and angiogenesis velocity. Anastomosis may be favored or impeded depending on the mechanical configuration of strain vectors in the ECM near tip cells. Numerical simulations demonstrate that increasing Jagged production results in pathological vasculatures with thinner and more abundant vessels, which can be compensated by augmenting the production of Delta ligands.

Angiogenesis is the process by which new blood vessels grow from existing ones. This process plays a crucial role in organ development, in wound healing and in numerous pathological processes such as cancer growth or in diabetes. Angiogenesis is a complex, multi-step and well regulated process where biochemistry and physics are intertwined. The process entails signaling in vessel cells being driven by both chemical and mechanical mechanisms that result in vascular cell movement, deformation and proliferation. Mathematical models have the ability to bring together these mechanisms in order to explore their relative relevance in vessel growth. Here, we present a mathematical model of early stage angiogenesis that is able to explore the role of biochemical signaling and tissue mechanics. We use this model to unravel the regulating role of Jagged, Notch and Delta dynamics in vascular cells. These membrane proteins have an important part in determining the leading cell in each neo-vascular sprout. Numerical simulations demonstrate that increasing Jagged production results in pathological vasculatures with thinner and more abundant vessels, which can be compensated by augmenting the production of Delta ligands.

Mathematical Modeling to Address Challenges in Pancreatic Cancer

Pancreatic Ductal Adenocarcinoma (PDAC) is regarded as one of the most lethal cancer types for its challenges associated with early diagnosis and resistance to standard chemotherapeutic agents, thereby leading to a poor five-year survival rate. The complexity of the disease calls for a multidisciplinary approach to better manage the disease and improve the status quo in PDAC diagnosis, prognosis, and treatment. To this end, the application of quantitative tools can help improve the understanding of disease mechanisms, develop biomarkers for early diagnosis, and design patient-specific treatment strategies to improve therapeutic outcomes. However, such approaches have only been minimally applied towards the investigation of PDAC, and we review the current status of mathematical modeling works in this field.

Oxygen in the Tumor Microenvironment: Mathematical and Numerical Modeling

There are many reasons to try to achieve a good grasp of the distribution of oxygen in the tumor microenvironment. The lack of oxygen - hypoxia - is a main actor in the evolution of tumors and in their growth and appears to be just as important in tumor invasion and metastasis. Mathematical models of the distribution of oxygen in tumors which are based on reaction-diffusion equations provide partial but qualitatively significant descriptions of the measured oxygen concentrations in the tumor microenvironment, especially when they incorporate important elements of the blood vessel network such as the blood vessel size and spatial distribution and the pulsation of local pressure due to blood circulation. Here, we review our mathematical and numerical approaches to the distribution of oxygen that yield insights both on the role of the distribution of blood vessel density and size and on the fluctuations of blood pressure.

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.

An approach for vascular tumor growth based on a hybrid embedded/homogenized treatment of the vasculature within a multiphase porous medium model

The aim of this work is to develop a novel computational approach to facilitate the modeling of angiogenesis during tumor growth. The preexisting vasculature is modeled as a 1D inclusion and embedded into the 3D tissue through a suitable coupling method, which allows for nonmatching meshes in 1D and 3D domain. The neovasculature, which is formed during angiogenesis, is represented in a homogenized way as a phase in our multiphase porous medium system. This splitting of models is motivated by the highly complex morphology, physiology, and flow patterns in the neovasculature, which are challenging and computationally expensive to resolve with a discrete, 1D angiogenesis and blood flow model. Moreover, it is questionable if a discrete representation generates any useful additional insight. By contrast, our model may be classified as a hybrid vascular multiphase tumor growth model in the sense that a discrete, 1D representation of the preexisting vasculature is coupled with a continuum model describing angiogenesis. It is based on an originally avascular model which has been derived via the thermodynamically constrained averaging theory. The new model enables us to study mass transport from the preexisting vasculature into the neovasculature and tumor tissue. We show by means of several illustrative examples that it is indeed capable of reproducing important aspects of vascular tumor growth phenomenologically.