Scaling rules for diffusive drug delivery in tumor and normal tissues

Identification of established and novel extracellular matrix components in glioblastoma as targets for angiogenesis and prognosis

Glioblastomas (GBM) are aggressive tumors known for their heterogeneity, rapid proliferation, treatment resistance, and extensive vasculature. Angiogenesis, the formation of new vessels, involves endothelial cell (EC) migration and proliferation. Various extracellular matrix (ECM) molecules regulate EC survival, migration, and proliferation. Culturing human brain EC (HBMEC) on GBM-derived ECM revealed a decrease in EC numbers compared to controls. Through in silico analysis, we explored ECM gene expression differences between GBM and brain normal glia cells and the impact of GBM microenvironment on EC ECM transcripts. ECM molecules such as collagen alpha chains (COL4A1, COL4A2, p < 0.0001); laminin alpha (LAMA4), beta (LAMB2), and gamma (LAMC1) chains (p < 0.0005); neurocan (NCAN), brevican (BCAN) and versican (VCAN) (p < 0.0005); hyaluronan synthase (HAS) 2 and metalloprotease (MMP) 2 (p < 0.005); MMP inhibitors (TIMP1-4, p < 0.0005), transforming growth factor beta-1 (TGFB1) and integrin alpha (ITGA3/5) (p < 0.05) and beta (ITGB1, p < 0.0005) chains showed increased expression in GBM. Additionally, GBM-influenced EC exhibited elevated expression of COL5A3, COL6A1, COL22A1 and COL27A1 (p < 0.01); LAMA1, LAMB1 (p < 0.001); fibulins (FBLN1/2, p < 0.01); MMP9, HAS1, ITGA3, TGFB1, and wingless-related integration site 9B (WNT9B) (p < 0.01) compared to normal EC. Some of these molecules: COL5A1/3, COL6A1, COL22/27A1, FBLN1/2, ITGA3/5, ITGB1 and LAMA1/B1 (p < 0.01); NCAN, HAS1, MMP2/9, TIMP1/2 and TGFB1 (p < 0.05) correlated with GBM patient survival. In conclusion, this study identified both established and novel ECM molecules regulating GBM angiogenesis, suggesting NCAN and COL27A1 are new potential prognostic biomarkers for GBM.

Insights into the roles of non-coding RNAs and angiogenesis in glioblastoma: An overview of current research and future perspectives

Glioblastoma (GBM) is a highly aggressive type of primary brain cancer with a poor prognosis, and despite intensive research, survival rates have not significantly improved. Non-coding RNAs (ncRNAs) are emerging as critical regulators of GBM pathogenesis, including angiogenesis, which is essential for tumor growth and invasion. MicroRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) have been identified as regulators of angiogenesis in GBM. miRNAs such as miR-21, miR-10b, and miR-26a promote angiogenesis by targeting anti-angiogenic factors, while lncRNAs such as H19 and MALAT1 inhibit angiogenesis by regulating pro-angiogenic factors. CircRNAs, such as circSMARCA5 and circBACH2, also regulate angiogenesis through various mechanisms. Similarly, signaling pathways such as the vascular endothelial growth factor (VEGF) pathway play critical roles in angiogenesis and have been targeted for GBM therapy. However, resistance to anti-angiogenic therapies is a significant obstacle in clinical practice. Developing novel therapeutic strategies targeting ncRNAs and angiogenesis is a promising approach for GBM. Potential targets include miRNAs, lncRNAs, circRNAs, and downstream signaling pathways that regulate angiogenesis. This review highlights the critical roles of ncRNAs and angiogenesis in GBM pathogenesis and the potential for new therapeutic strategies targeting these pathways to improve the prognosis and quality of life for GBM patients.

Mathematical model of oxygen, nutrient, and drug transport in tuberculosis granulomas

Physiological abnormalities in pulmonary granulomas-pathological hallmarks of tuberculosis (TB)-compromise the transport of oxygen, nutrients, and drugs. In prior studies, we demonstrated mathematically and experimentally that hypoxia and necrosis emerge in the granuloma microenvironment (GME) as a direct result of limited oxygen availability. Building on our initial model of avascular oxygen diffusion, here we explore additional aspects of oxygen transport, including the roles of granuloma vasculature, transcapillary transport, plasma dilution, and interstitial convection, followed by cellular metabolism. Approximate analytical solutions are provided for oxygen and glucose concentration, interstitial fluid velocity, interstitial fluid pressure, and the thickness of the convective zone. These predictions are in agreement with prior experimental results from rabbit TB granulomas and from rat carcinoma models, which share similar transport limitations. Additional drug delivery predictions for anti-TB-agents (rifampicin and clofazimine) strikingly match recent spatially-resolved experimental results from a mouse model of TB. Finally, an approach to improve molecular transport in granulomas by modulating interstitial hydraulic conductivity is tested in silico.

Biomechanical forces and force-triggered drug delivery in tumor neovascularization

Tumor angiogenesis is one of the typical hallmarks of tumor occurrence and development, and tumor neovascularization also exhibits distinct characteristics from normal blood vessels. As the number of cells and matrix inside the tumor increases, the biomechanical force is enhanced, specifically manifested as solid stress, fluid stress, stiffness, and topology. This mechanical microenvironment also provides shelter for tumors and intensifies angiogenesis, providing oxygen and nutritional support for tumor progression. During tumor development, the biomechanical microenvironment also emerges, which in turn feeds back to regulate the tumor progression, including tumor angiogenesis, and biochemical and biomechanical signals can regulate tumor angiogenesis. Blood vessels possess inherent sensitivity to mechanical stimuli, but compared to the extensive research on biochemical signal regulation, the study of the regulation of tumor neovascularization by biomechanical signals remains relatively scarce. Biomechanical forces can affect the phenotypic characteristics and mechanical signaling pathways of tumor blood vessels, directly regulating angiogenesis. Meanwhile, they can indirectly regulate tumor angiogenesis by causing an imbalance in angiogenesis signals and affecting stromal cell function. Understanding the regulatory mechanism of biomechanical forces in tumor angiogenesis is beneficial for better identifying and even taming the mechanical forces involved in angiogenesis, providing new therapeutic targets for tumor vascular normalization. Therefore, we summarized the composition of biomechanical forces and their direct or indirect regulation of tumor neovascularization. In addition, this review discussed the use of biomechanical forces in combination with anti-angiogenic therapies for the treatment of tumors, and biomechanical forces triggered delivery systems.

Sensitivity analysis of a multi-physics model for the vascular microenvironment

The vascular microenvironment is the scale at which microvascular transport, interstitial tissue properties and cell metabolism interact. The vascular microenvironment has been widely studied by means of quantitative approaches, including multi-physics mathematical models as it is a central system for the pathophysiology of many diseases, such as cancer. The microvascular architecture is a key factor for fluid balance and mass transfer in the vascular microenvironment, together with the physical parameters characterizing the vascular wall and the interstitial tissue. The scientific literature of this field has witnessed a long debate about which factor of this multifaceted system is the most relevant. The purpose of this work is to combine the interpretative power of an advanced multi-physics model of the vascular microenvironment with state of the art and robust sensitivity analysis methods, in order to determine the factors that most significantly impact quantities of interest, related in particular to cancer treatment. We are particularly interested in comparing the factors related to the microvascular architecture with the ones affecting the physics of microvascular transport. Ultimately, this work will provide further insight into how the vascular microenvironment affects cancer therapies, such as chemotherapy, radiotherapy or immunotherapy.

The relationship between cancer and biomechanics

The onset, development, diagnosis, and treatment of cancer involve intricate interactions among various factors, spanning the realms of mechanics, physics, chemistry, and biology. Within our bodies, cells are subject to a variety of forces such as gravity, magnetism, tension, compression, shear stress, and biological static force/hydrostatic pressure. These forces are perceived by mechanoreceptors as mechanical signals, which are then transmitted to cells through a process known as mechanical transduction. During tumor development, invasion and metastasis, there are significant biomechanical influences on various aspects such as tumor angiogenesis, interactions between tumor cells and the extracellular matrix (ECM), interactions between tumor cells and other cells, and interactions between tumor cells and the circulatory system and vasculature. The tumor microenvironment comprises a complex interplay of cells, ECM and vasculature, with the ECM, comprising collagen, fibronectins, integrins, laminins and matrix metalloproteinases, acting as a critical mediator of mechanical properties and a key component within the mechanical signaling pathway. The vasculature exerts appropriate shear forces on tumor cells, enabling their escape from immune surveillance, facilitating their dissemination in the bloodstream, dictating the trajectory of circulating tumor cells (CTCs) and playing a pivotal role in regulating adhesion to the vessel wall. Tumor biomechanics plays a critical role in tumor progression and metastasis, as alterations in biomechanical properties throughout the malignant transformation process trigger a cascade of changes in cellular behavior and the tumor microenvironment, ultimately culminating in the malignant biological behavior of the tumor.

Mechanistically Weighted Metric to Predict In Vivo Antibody-Receptor Occupancy: An Analytical Approach

In situ clinical measurement of receptor occupancy (RO) is challenging, particularly for solid tumors, necessitating the use of mathematical models that predict tumor receptor occupancy to guide dose decisions. A potency metric, average free tissue target to initial target ratio (AFTIR), was previously described based on a mechanistic compartmental model and is informative for near-saturating dose regimens. However, the metric fails at clinically relevant subsaturating antibody doses, as compartmental models cannot capture the spatial heterogeneity of distribution faced by some antibodies in solid tumors. Here we employ a partial differential equation (PDE) Krogh cylinder model to simulate spatiotemporal receptor occupancy and derive an analytical solution, a mechanistically weighted global AFTIR, that can better predict receptor occupancy regardless of dosing regimen. In addition to the four key parameters previously identified, a fifth key parameter, the absolute receptor density (targets/cell), is incorporated into the mechanistic AFTIR metric. Receptor density can influence equilibrium intratumoral drug concentration relative to whether the dose is saturating or not, thereby influencing the tumor penetration depth of the antibody. We derive mechanistic RO predictions based on distinct patterns of antibody tumor penetration, presented as a global AFTIR metric guided by a Thiele Modulus and a local saturation potential (drug equivalent of binding potential for positron emissions tomography imaging) and validate the results using rigorous global and local sensitivity analysis. This generalized AFTIR serves as a more accurate analytical metric to aid clinical dose decisions and rational design of antibody-based therapeutics without the need for extensive PDE simulations. SIGNIFICANCE STATEMENT: Determining antibody-receptor occupancy (RO) is critical for dosing decisions in pharmaceutical development, but direct clinical measurement of RO is often challenging and invasive, particularly for solid tumors. Significant efforts have been made to develop mathematical models and simplified analytical metrics of RO, but these often require complex computer simulations. Here we present a mathematically rigorous but simplified analytical model to accurately predict RO across a range of affinities, doses, drug, and tumor properties.

Positive Chemotaxis of CREKA-Modified Ceria@Polydopamine Biomimetic Nanoswimmers for Enhanced Penetration and Chemo-photothermal Tumor Therapy

Tumor interstitial pressure represents the greatest barrier against drug diffusion into the depth of the tumor. Biometric nanomotors highlight the possibility of enhanced deep penetration and improve cellular uptake. However, control of their directionality remains difficult to achieve. Herein, we report cysteine-arginine-glutamic acid-lysine-alanine (CREKA)-modified ceria@polydopamine nanobowls as tumor microenvironment-fueled nanoscale motors for positive chemotaxis into the tumor depth or toward tumor cells. Upon laser irradiation, this nanoswimmer rapidly depletes the tumor microenvironment-specific hydrogen peroxide (H2O2) in the nanobowl, contributing to a self-generated gradient and subsequently propulsion (9.5 μm/s at 46 °C). Moreover, the asymmetrical modification of CREKA on nanobowls could automatically reconfigure the motion direction toward tumor depth or tumor cells in response to receptor-ligand interaction, leading to a deep penetration (70 μm in multicellular spheroids) and enhanced antitumor effects over conventional nanomedicine-induced chemo-photothermal therapy (tumor growth inhibition rate: 84.2% versus 56.9%). Thus, controlling the direction of nanomotors holds considerable potential for improved antitumor responses, especially in solid tumors with high tumor interstitial pressure.

Investigating the Mechanism of Intravascular Bubble Formation in Designed Arrays of Vascularized Systems on a Chip

Vascular gas embolism is a rare medical condition, resulting from the existence of air or gas in the venous or arterial system. Gas embolism is associated with a wide range of circulatory, cardiovascular, and neurological complications that can lead to sudden and unexplained death. Despite the recent increase in related studies, gas embolism remains under-reported with a poor understanding of its genesis and pathophysiology. In this work, intravascular bubble formation is investigated in an array of biomimetic microscale systems, where the endogenous generation of gas bubbles is induced by variations in the surrounding pressure. Microfluidic devices, based on polydimethylsiloxane, are designed and fabricated as vascularized systems on a chip with one main channel at two different diameters (30 µm, and 40 µm), surrounded by a pressure chamber (200 µm) on each side, at a separation of 50 µm. Two blood-equivalent solutions, at 20% and 46% hematocrit concentrations were prepared from a glycerin and xanthan gum mixture to mimic the physicochemical characteristics of the blood. As the volume of injected air increased, the events related to gas embolism were occurring at shorter timespans with more significant characteristics, i.e., length and number of bubbles. Additionally, correlations were established between the input parameters, i.e., the vascular diameter and equivalent hematocrit concentration, and the output parameters, i.e., the bubble size, velocity, frequency, and nucleation sites.Clinical Relevance- The reported results constitute a reproducible observation and quantification of intravascular bubble formation induced by global pressure variations, where the emergence of bubbles exhibits different patterns depending on biological characteristics related to gender and medical history.

Computational simulations of tumor growth and treatment response: Benefits of high-frequency, low-dose drug regimens and concurrent vascular normalization

Implementation of effective cancer treatment strategies requires consideration of how the spatiotemporal heterogeneities within the tumor microenvironment (TME) influence tumor progression and treatment response. Here, we developed a multi-scale three-dimensional mathematical model of the TME to simulate tumor growth and angiogenesis and then employed the model to evaluate an array of single and combination therapy approaches. Treatments included maximum tolerated dose or metronomic (i.e., frequent low doses) scheduling of anti-cancer drugs combined with anti-angiogenic therapy. The results show that metronomic therapy normalizes the tumor vasculature to improve drug delivery, modulates cancer metabolism, decreases interstitial fluid pressure and decreases cancer cell invasion. Further, we find that combining an anti-cancer drug with anti-angiogenic treatment enhances tumor killing and reduces drug accumulation in normal tissues. We also show that combined anti-angiogenic and anti-cancer drugs can decrease cancer invasiveness and normalize the cancer metabolic microenvironment leading to reduced hypoxia and hypoglycemia. Our model simulations suggest that vessel normalization combined with metronomic cytotoxic therapy has beneficial effects by enhancing tumor killing and limiting normal tissue toxicity.

Investigation of the evolution of tumor-induced microvascular network under the inhibitory effect of anti-angiogenic factor, angiostatin: A mathematical study

Anti-angiogenesis as a treatment strategy for normalizing the microvascular network of tumors is of great interest among researchers, especially in combination with chemotherapy or radiotherapy. According to the vital role that angiogenesis plays in tumor growth and in exposing the tumor to therapeutic agents, this work develops a mathematical framework to study the influence of angiostatin, a plasminogen fragment that shows the anti-angiogenic function, in the evolutionary behavior of tumor-induced angiogenesis. Angiostatin-induced microvascular network reformation is investigated in a two-dimensional space by considering two parent vessels around a circular tumor by a modified discrete angiogenesis model in different tumor sizes. The effects of imposing modifications on the existing model, i.e., the matrix-degrading enzyme effect, proliferation and death of endothelial cells, matrix density function, and a more realistic chemotactic function, are investigated in this study. Results show a decrease in microvascular density in response to the angiostatin. A functional relationship exists between angiostatin's ability to normalize the capillary network and tumor size or progression stage, such that capillary density decreases by 55%, 41%, 24%, and 13% in tumors with a non-dimensional radius of 0.4, 0.3, 0.2, and 0.1, respectively, after angiostatin administration.

Bioactive Glasses as Carriers of Cancer-Targeted Drugs: Challenges and Opportunities in Bone Cancer Treatment

The treatment of bone cancer involves tumor resection followed by bone reconstruction of the defect caused by the tumor using biomaterials. Additionally, post-surgery protocols cover chemotherapy, radiotherapy, or drug administration, which are employed as adjuvant treatments to prevent tumor recurrence. In this work, we reviewed new strategies for bone cancer treatment based on bioactive glasses as carriers of cancer-targeted and other drugs that are intended for bone regeneration in conjunction with adjuvant treatments. Drugs used in combination with bioactive glasses can be classified into cancer-target, osteoclast-target, and new therapies (such as gene delivery and bioinorganic). Microparticulated, nanoparticulated, or mesoporous bioactive glasses have been used as drug-delivery systems. Additionally, surface modification through functionalization or the production of composites based on polymers and hydrogels has been employed to improve drug-release kinetics. Overall, although different drugs and drug delivery systems have been developed, there is still room for new studies involving kinase inhibitors or antibody-conjugated drugs, as these drugs have been poorly explored in combination with bioactive glasses.

Combination of Carriers with Complementary Intratumoral Microdistributions of Delivered -Particles May Realize the Promise for 225Ac in Large, Solid Tumors

α-particle radiotherapy has already been shown to be impervious to most resistance mechanisms. However, in established (i.e., large, vascularized) soft-tissue lesions, the diffusion-limited penetration depths of radiolabeled antibodies or nanocarriers (≤50-80 μm) combined with the short range of α-particles (4-5 cell diameters) may result in only partial tumor irradiation, potentially limiting treatment efficacy. To address this challenge, we combined carriers with complementary intratumoral microdistributions of the delivered α-particles. We used the α-particle generator ²²⁵Ac, and we combined a tumor-responsive liposome (which, on tumor uptake, releases into the interstitium a highly diffusing form of its radioactive payload [²²⁵Ac-DOTA], potentially penetrating the deeper parts of tumors where antibodies do not reach) with a separately administered, less-penetrating radiolabeled antibody (irradiating the tumor perivascular regions where liposome contents clear too quickly). Methods: In a murine model with orthotopic human epidermal growth factor receptor 2-positive BT474 breast cancer xenografts, the biodistributions of each carrier were evaluated, and the control of tumor growth was monitored after administration of the same total radioactivity of ²²⁵Ac delivered by the ²²⁵Ac-DOTA-encapsulating liposomes, by the ²²⁵Ac-DOTA-SCN--labeled trastuzumab, and by both carriers at equally split radioactivities. Results: Tumor growth was significantly more inhibited when the same total injected radioactivity was divided between the 2 separate carriers than when delivered by either of the carriers alone. The combined carriers enabled more uniform intratumoral microdistributions of α-particles, at a tumor dose that was lower than the dose delivered by the antibody alone. Conclusion: This strategy demonstrates that more uniform microdistributions of the delivered α-particles within established solid tumors improve efficacy even at lower tumor doses. Augmentation of antibody-targeted α-particle therapies with tumor-responsive liposomes may address partial tumor irradiation, improving therapeutic effects.

Design of an Integrated Microvascularized Human Skin-on-a-Chip Tissue Equivalent Model

Tissue-engineered skin constructs have been under development since the 1980s as a replacement for human skin tissues and animal models for therapeutics and cosmetic testing. These have evolved from simple single-cell assays to increasingly complex models with integrated dermal equivalents and multiple cell types including a dermis, epidermis, and vasculature. The development of micro-engineered platforms and biomaterials has enabled scientists to better recreate and capture the tissue microenvironment in vitro, including the vascularization of tissue models and their integration into microfluidic chips. However, to date, microvascularized human skin equivalents in a microfluidic context have not been reported. Here, we present the design of a novel skin-on-a-chip model integrating human-derived primary and immortalized cells in a full-thickness skin equivalent. The model is housed in a microfluidic device, in which a microvasculature was previously established. We characterize the impact of our chip design on the quality of the microvascular networks formed and evidence that this enables the formation of more homogenous networks. We developed a methodology to harvest tissues from embedded chips, after 14 days of culture, and characterize the impact of culture conditions and vascularization (including with pericyte co-cultures) on the stratification of the epidermis in the resulting skin equivalents. Our results indicate that vascularization enhances stratification and differentiation (thickness, architecture, and expression of terminal differentiation markers such as involucrin and transglutaminase 1), allowing the formation of more mature skin equivalents in microfluidic chips. The skin-on-a-chip tissue equivalents developed, because of their realistic microvasculature, may find applications for testing efficacy and safety of therapeutics delivered systemically, in a human context.

Multiphoton Phosphorescence Quenching Microscopy Reveals Kinetics of Tumor Oxygenation during Antiangiogenesis and Angiotensin Signaling Inhibition

PURPOSE

The abnormal function of tumor blood vessels causes tissue hypoxia, promoting disease progression and treatment resistance. Although tumor microenvironment normalization strategies can alleviate hypoxia globally, how local oxygen levels change is not known because of the inability to longitudinally assess vascular and interstitial oxygen in tumors with sufficient resolution. Understanding the spatial and temporal heterogeneity should help improve the outcome of various normalization strategies.

EXPERIMENTAL DESIGN

We developed a multiphoton phosphorescence quenching microscopy system using a low-molecular-weight palladium porphyrin probe to measure perfused vessels, oxygen tension, and their spatial correlations in vivo in mouse skin, bone marrow, and four different tumor models. Further, we measured the temporal and spatial changes in oxygen and vessel perfusion in tumors in response to an anti-VEGFR2 antibody (DC101) and an angiotensin-receptor blocker (losartan).

RESULTS

We found that vessel function was highly dependent on tumor type. Although some tumors had vessels with greater oxygen-carrying ability than those of normal skin, most tumors had inefficient vessels. Further, intervessel heterogeneity in tumors is associated with heterogeneous response to DC101 and losartan. Using both vascular and stromal normalizing agents, we show that spatial heterogeneity in oxygen levels persists, even with reductions in mean extravascular hypoxia.

CONCLUSIONS

High-resolution spatial and temporal responses of tumor vessels to two agents known to improve vascular perfusion globally reveal spatially heterogeneous changes in vessel structure and function. These dynamic vascular changes should be considered in optimizing the dose and schedule of vascular and stromal normalizing strategies to improve the therapeutic outcome.

Modeling the impact of spatial oxygen heterogeneity on radiolytic oxygen depletion during FLASH radiotherapy

Purpose.It has been postulated that the delivery of radiotherapy at ultra-high dose rates ('FLASH') reduces normal tissue toxicities by depleting them of oxygen. The fraction of normal tissue and cancer cells surviving radiotherapy depends on dose and oxygen levels in an exponential manner and even a very small fraction of tissue at low oxygen levels can determine radiotherapy response. To quantify the differential impact of FLASH radiotherapy on normal and tumour tissues, the spatial heterogeneity of oxygenation in tissue should thus be accounted for.Methods.The effect of FLASH on radiation-induced normal and tumour tissue cell killing was studied by simulating oxygen diffusion, metabolism, and radiolytic oxygen depletion (ROD) over domains with simulated capillary architectures. To study the impact of heterogeneity, two architectural models were used: (1) randomly distributed capillaries and (2) capillaries forming a regular square lattice array. The resulting oxygen partial pressure distribution histograms were used to simulate normal and tumour tissue cell survival using the linear quadratic model of cell survival, modified to incorporate oxygen-enhancement ratio effects. The ratio ('dose modifying factors') of conventional low-dose-rate dose and FLASH dose at iso-cell survival was computed and compared with empirical iso-toxicity dose ratios.Results.Tumour cell survival was found to be increased by FLASH as compared to conventional radiotherapy, with a 0-1 order of magnitude increase for expected levels of tumour hypoxia, depending on the relative magnitudes of ROD and tissue oxygen metabolism. Interestingly, for the random capillary model, the impact of FLASH on well-oxygenated (normal) tissues was found to be much greater, with an estimated increase in cell survival by up to 10 orders of magnitude, even though reductions in mean tissue partial pressure were modest, less than ∼7 mmHg for the parameter values studied. The dose modifying factor for normal tissues was found to lie in the range 1.2-1.7 for a representative value of normal tissue oxygen metabolic rate, consistent with preclinical iso-toxicity results.Conclusions.The presence of very small nearly hypoxic regions in otherwise well-perfused normal tissues with high mean oxygen levels resulted in a greater proportional sparing of normal tissue than tumour cells during FLASH irradiation, possibly explaining empirical normal tissue sparing and iso-tumour control results.

Investigation of cancer response to chemotherapy: a hybrid multi-scale mathematical and computational model of the tumor microenvironment

Tumor microenvironment (TME) is a multi-scale biological environment that can control tumor dynamics with many biomechanical and biochemical factors. Investigating the physiology of TME with a heterogeneous structure and abnormal functions not only can achieve a deeper understanding of tumor behavior but also can help develop more efficient anti-cancer strategies. In this work, we develop a hybrid multi-scale mathematical model of TME to simulate the progression of a three-dimensional tumor and elucidate its response to different chemotherapy approaches. The chemotherapy approaches include multiple low dose (MLD) of anti-cancer drug, maximum tolerated dose (MTD) of anti-cancer drug, combination therapy of MLD and anti-angiogenic drug, and combination therapy of MTD and anti-angiogenic drug. The results show that combining anti-angiogenic agent with anti-cancer drug increases the performance of cancer treatment and decreases side effects for normal tissue. Indeed, the vascular normalization caused by anti-angiogenic therapy improves anti-cancer drug delivery for both MLD and MTD approaches. The results show that anti-cancer drug administered in a lower dose than the maximum tolerated dose repetitively over a long period treats cancer with a considerable performance and fewer side effects. We also show that tumor morphology and distribution of cancer cell phenotypes can be considered as the characteristics to distinguish different chemotherapy approaches. This robust model can be applied to predict the behavior of any type of cancer and quantify cancer response to different chemotherapy approaches. The computational results of cancer response to chemotherapy are in good agreement with experimental measurements.

Functions and clinical significance of mechanical tumor microenvironment: cancer cell sensing, mechanobiology and metastasis

Dynamic and heterogeneous interaction between tumor cells and the surrounding microenvironment fuels the occurrence, progression, invasion, and metastasis of solid tumors. In this process, the tumor microenvironment (TME) fractures cellular and matrix architecture normality through biochemical and mechanical means, abetting tumorigenesis and treatment resistance. Tumor cells sense and respond to the strength, direction, and duration of mechanical cues in the TME by various mechanotransduction pathways. However, far less understood is the comprehensive perspective of the functions and mechanisms of mechanotransduction. Due to the great therapeutic difficulties brought by the mechanical changes in the TME, emerging studies have focused on targeting the adverse mechanical factors in the TME to attenuate disease rather than conventionally targeting tumor cells themselves, which has been proven to be a potential therapeutic approach. In this review, we discussed the origins and roles of mechanical factors in the TME, cell sensing, mechano‐biological coupling and signal transduction, in vitro construction of the tumor mechanical microenvironment, applications and clinical significance in the TME.

Longitudinal in-vivo quantification of tumour microvascular heterogeneity by optical coherence angiography in pre-clinical radiation therapy

Stereotactic body radiotherapy (SBRT) is an emerging cancer treatment due to its logistical and potential therapeutic benefits as compared to conventional radiotherapy. However, its mechanism of action is yet to be fully understood, likely involving the ablation of tumour microvasculature by higher doses per fraction used in SBRT. In this study, we hypothesized that longitudinal imaging and quantification of the vascular architecture may elucidate the relationship between the microvasculature and tumour response kinetics. Pancreatic human tumour xenografts were thus irradiated with single doses of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10$$\end{document}10, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$20$$\end{document}20 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$30$$\end{document}30 Gy to simulate the first fraction of a SBRT protocol. Tumour microvascular changes were monitored with optical coherence angiography for up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8$$\end{document}8 weeks following irradiation. The temporal kinetics of two microvascular architectural metrics were studied as a function of time and dose: the diffusion-limited fraction, representing poorly vascularized tissue \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$>150$$\end{document}>150 μm from the nearest detected vessel, and the vascular distribution convexity index, a measure of vessel aggregation at short distances. These biological metrics allowed for dose dependent temporal evaluation of tissue (re)vascularization and vessel aggregation after radiotherapy, showing promise for determining the SBRT dose–response relationship.

Brain tumor classification of virtual NMR voxels based on realistic blood vessel-induced spin dephasing using support vector machines

Remodeling of tissue microvasculature commonly promotes neoplastic growth; however, there is no imaging modality in oncology yet that noninvasively quantifies microvascular changes in clinical routine. Although blood capillaries cannot be resolved in typical magnetic resonance imaging (MRI) measurements, their geometry and distribution influence the integral nuclear magnetic resonance (NMR) signal from each macroscopic MRI voxel. We have numerically simulated the expected transverse relaxation in NMR voxels with different dimensions based on the realistic microvasculature in healthy and tumor-bearing mouse brains (U87 and GL261 glioblastoma). The 3D capillary structure in entire, undissected brains was acquired using light sheet fluorescence microscopy to produce large datasets of the highly resolved cerebrovasculature. Using this data, we trained support vector machines to classify virtual NMR voxels with different dimensions based on the simulated spin dephasing accountable to field inhomogeneities caused by the underlying vasculature. In prediction tests with previously blinded virtual voxels from healthy brain tissue and GL261 tumors, stable classification accuracies above 95% were reached. Our results indicate that high classification accuracies can be stably attained with achievable training set sizes and that larger MRI voxels facilitated increasingly successful classifications, even with small training datasets. We were able to prove that, theoretically, the transverse relaxation process can be harnessed to learn endogenous contrasts for single voxel tissue type classifications on tailored MRI acquisitions. If translatable to experimental MRI, this may augment diagnostic imaging in oncology with automated voxel-by-voxel signal interpretation to detect vascular pathologies.

Model studies of the role of oxygen in the FLASH effect

Current radiotherapy facilities are standardized to deliver dose rates around 0.1-0.4 Gy/s in 2 Gy daily fractions, designed to deliver total accumulated doses to reach the tolerance limit of normal tissues undergoing irradiation. FLASH radiotherapy (FLASH-RT), on the other hand, relies on facilities capable of delivering ultrahigh dose rates in large doses in a single microsecond pulse, or in a few pulses given over a very short time sequence. For example, most studies to date have implemented 4-6 MeV electrons with intra-pulse dose rates in the range 10⁶ -10⁷  Gy/s. The proposed dependence of the FLASH effect on oxygen tension has stimulated several theoretical models based on three different hypotheses: (i) Radiation-induced transient oxygen depletion; (ii) cell-specific differences in the ability to detoxify and/or recover from injury caused by reactive oxygen species; (iii) self-annihilation of radicals by bimolecular recombination. This article focuses on the observations supporting or refuting these models in the frame of the chemical-biological bases of the impact of oxygen on the radiation response of cell free, in vitro and in vivo model systems.

Bridging the macro to micro resolution gap with angiographic optical coherence tomography and dynamic contrast enhanced MRI

Dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) is emerging as a valuable tool for non-invasive volumetric monitoring of the tumor vascular status and its therapeutic response. However, clinical utility of DCE-MRI is challenged by uncertainty in its ability to quantify the tumor microvasculature (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu \mathrm{m}$$\end{document}μm scale) given its relatively poor spatial resolution (mm scale at best). To address this challenge, we directly compared DCE-MRI parameter maps with co-registered micron-scale-resolution speckle variance optical coherence tomography (svOCT) microvascular images in a window chamber tumor mouse model. Both semi and fully quantitative (Toft’s model) DCE-MRI metrics were tested for correlation with microvascular svOCT biomarkers. svOCT’s derived vascular volume fraction (VVF) and the mean distance to nearest vessel (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm{DNV} }$$\end{document}DNV¯) metrics were correlated with DCE-MRI vascular biomarkers such as time to peak contrast enhancement (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r=-0.81$$\end{document}r=-0.81 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.83$$\end{document}0.83 respectively, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P<0.0001$$\end{document}P<0.0001 for both), the area under the gadolinium-time concentration curve (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r=0.50$$\end{document}r=0.50 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-0.48$$\end{document}-0.48 respectively, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P<0.0001$$\end{document}P<0.0001 for both) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${k}_{trans}$$\end{document}ktrans (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r=0.64$$\end{document}r=0.64 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-0.61$$\end{document}-0.61 respectively, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P<0.0001$$\end{document}P<0.0001 for both). Several other correlated micro–macro vascular metric pairs were also noted. The microvascular insights afforded by svOCT may help improve the clinical utility of DCE-MRI for tissue functional status assessment and therapeutic response monitoring applications.

ROS-Responsive Dexamethasone Micelles Normalize the Tumor Microenvironment to Enhance Hypericin in Cancer Photodynamic Therapy

Efficacy of photodynamic therapy (PDT) for cancer is limited due to the abnormality of tumor microenvironment (TME), such as dysfunctional tumor vascular system leading to restrict the drug distribution in...

Modified microvessel density based on perfusion distance: a preferable NSCLC prognostic factor

Background

Despite the vital role of blood perfusion in tumor progression, the prognostic value of typical blood perfusion markers, such as microvessel density (MVD) or microvessel area (MVA), in patients with non-small cell lung cancer (NSCLC) is still unclear. This study established a modified MVD (mMVD) measurement based on perfusion distance and determined its prognostic value in patients with NSCLC.

Methods

A total of 100 patients with NSCLC were enrolled in this retrospective study. The intratumor microvessels of NSCLC patients were visualized using immunohistochemical staining for CD31. The blood perfusion distance was evaluated as the distance from each vessel to its nearest cancer cell (D_mvcc), and the cutoff value for prognosis was determined. Apart from the total MVD (tMVD), microvessels near cancer cells within the cutoff-D_mvcc were counted as mMVD. Predictive values for mortality and recurrence were evaluated and compared.

Results

The D_mvcc ranged from 1.6 to 269.8 µm (median, 13.1 µm). The mMVD (range: 2-70; median 23) was counted from tMVD according to the cutoff-D_mvcc (~20 µm). Compared with tMVD, a larger fraction of mMVD (80% vs. 2.9%) played a significant role in overall survival, with an improved area under the receiver operating characteristic (ROC) curve (AUC) (0.74 vs. 0.56). A high mMVD was an independent positive indicator of overall survival (OS) and progression-free survival (PFS). In contrast, tMVD was only related to PFS at the optimal cutoff.

Conclusions

Perfusion-distance-based mMVD is a promising prognostic factor for NSCLC patients with superior sensitivity, specificity, and clinical applicability compared to tMVD. This study provides novel insights into the prognostic role of tumor vessel perfusion in patients with NSCLC.

Engineered Bacteria Enhance Immunotherapy and Targeted Therapy through Stromal Remodeling of Tumors

Desmoplastic solid tumors are characterized by the rapid build-up of extracellular matrix (ECM) macromolecules, such as hyaluronic acid (HA). The resulting physiological barrier prevents the infiltration of immune cells and also impedes the delivery of anticancer agents. The development of a hypervesiculating Escherichia coli Nissle (ΔECHy) based tumor targeting bacterial system capable of distributing a fusion peptide, cytolysin A (ClyA)-hyaluronidase (Hy) via outer membrane vesicles (OMVs) is reported. The capability of targeting hypoxic tumors, manufacturing recombinant proteins in situ and the added advantage of an on-site OMV based distribution system makes the engineered bacterial vector a unique candidate for peptide delivery. The HA degrading potential of Hy for stromal modulation is combined with the cytolytic activity of ClyA followed by testing it within syngeneic cancer models. ΔECHy is combined with immune checkpoint antibodies and tyrosine kinase inhibitors (TKIs) to demonstrate that remodeling the tumor stroma results in the improvement of immunotherapy outcomes and enhancing the efficacy of biological signaling inhibitors. The biocompatibility of ΔECHy is also investigated to show that the engineered bacteria are effectively cleared, elicit minimal inflammatory and immune responses, and therefore could be a reliable candidate as a live biotherapeutic.

How translational modeling in oncology needs to get the mechanism just right

Translational model‐based approaches have played a role in increasing success in the development of novel anticancer treatments. However, despite this, significant translational uncertainty remains from animal models to patients. Optimization of dose and scheduling (regimen) of drugs to maximize the therapeutic utility (maximize efficacy while avoiding limiting toxicities) is still predominately driven by clinical investigations. Here, we argue that utilizing pragmatic mechanism‐based translational modeling of nonclinical data can further inform this optimization. Consequently, a prototype model is demonstrated that addresses the required fundamental mechanisms.

Roles of Non-coding RNAs and Angiogenesis in Glioblastoma

One of the significant hallmarks of cancer is angiogenesis. It has a crucial function in tumor development and metastasis. Thus, angiogenesis has become one of the most exciting targets for drug development in cancer treatment. Here we discuss the regulatory effects on angiogenesis in glioblastoma (GBM) of non-coding RNAs (ncRNAs), including long ncRNA (lncRNA), microRNA (miRNA), and circular RNA (circRNA). These ncRNAs may function in trans or cis forms and modify gene transcription by various mechanisms, including epigenetics. NcRNAs may also serve as crucial regulators of angiogenesis-inducing molecules. These molecules include, metalloproteinases, cytokines, several growth factors (platelet-derived growth factor, vascular endothelial growth factor, fibroblast growth factor, hypoxia-inducible factor-1, and epidermal growth factor), phosphoinositide 3-kinase, mitogen-activated protein kinase, and transforming growth factor signaling pathways.

Targeting tumor vascularization: promising strategies for vascular normalization

Tumor recurrence after the clinical cure of tumor often results from the presence of an abnormal microenvironment, including an aberrant vasculature. The tumor microenvironment is rich in pro-angiogenic factors but lacks pro-maturation factors. Pro-angiogenic conditions in the tumor microenvironment, such as hypoxia, are double-edged swords, promoting both the repair of normal tissues and the development of an abnormal blood vessel network. The coexistence of perfusion and hypoxic zones and uneven blood vessel distribution in tumor tissues profoundly influence tumor deterioration, recurrence, and metastasis. Traditional anti-angiogenic therapies have shown limited efficacy, and promote drug resistance, and even metastasis. In contrast, vascular normalization therapy induces a more physiological-like state, leading to better outcomes and fewer side effects. Vascular normalization entails modifying the tumor vascular system to improve tumor oxygenation and substance transport, thereby contributing to improving the efficacy of radiotherapy, chemotherapy, and immunotherapy. This review mainly focuses on the process of tumor vascularization; potential therapeutic targets, including cells, metabolism, signaling pathways, and angiogenesis-related genes; and possible strategies to normalize blood vessels through regulating tumor vessel generation, the development of tumor vessels, and blood vessel fusion and pruning.

Chemo-mechanistic multi-scale model of a three-dimensional tumor microenvironment to quantify the chemotherapy response of cancer

Exploring efficient chemotherapy would benefit from a deeper understanding of the tumor microenvironment (TME) and its role in tumor progression. As in vivo experimental methods are unable to isolate or control individual factors of the TME, and in vitro models often cannot include all the contributing factors, some questions are best addressed with mathematical models of systems biology. In this study, we establish a multi-scale mathematical model of the TME to simulate three-dimensional tumor growth and angiogenesis and then implement the model for an array of chemotherapy approaches to elucidate the effect of TME conditions and drug scheduling on controlling tumor progression. The hyperglycemic condition as the most common disorder for cancer patients is considered to evaluate its impact on cancer response to chemotherapy. We show that combining antiangiogenic and anticancer drugs improves the outcome of treatment and can decrease accumulation of the drug in normal tissue and enhance drug delivery to the tumor. Our results demonstrate that although both concurrent and neoadjuvant combination therapies can increase intratumoral drug exposure and therapeutic accuracy, neoadjuvant therapy surpasses this, especially against hyperglycemia. Our model provides mechanistic explanations for clinical observations of tumor progression and response to treatment and establishes a computational framework for exploring better treatment strategies.

Remodeling tumor microenvironment with nanomedicines

The tumor microenvironment (TME) has been recognized as a major contributor to cancer malignancy and therapeutic resistance. Thus, strategies directed to re-engineer the TME are emerging as promising approaches for improving the efficacy of antitumor therapies by enhancing tumor perfusion and drug delivery, as well as alleviating the immunosuppressive TME. In this regard, nanomedicine has shown great potential for developing effective treatments capable of re-modeling the TME by controlling drug action in a spatiotemporal manner and allowing long-lasting modulatory effects on the TME. Herein, we review recent progress on TME re-engineering by using nanomedicine, particularly focusing on formulations controlling TME characteristics through targeted interaction with cellular components of the TME. Importantly, the TME should be re-engineering to a quiescent phenotype rather than be destroyed. Finally, immediate challenges and future perspectives of TME-re-engineering nanomedicines are discussed, anticipating further innovation in this growing field. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease.

Topological data analysis distinguishes parameter regimes in the Anderson-Chaplain model of angiogenesis

Angiogenesis is the process by which blood vessels form from pre-existing vessels. It plays a key role in many biological processes, including embryonic development and wound healing, and contributes to many diseases including cancer and rheumatoid arthritis. The structure of the resulting vessel networks determines their ability to deliver nutrients and remove waste products from biological tissues. Here we simulate the Anderson-Chaplain model of angiogenesis at different parameter values and quantify the vessel architectures of the resulting synthetic data. Specifically, we propose a topological data analysis (TDA) pipeline for systematic analysis of the model. TDA is a vibrant and relatively new field of computational mathematics for studying the shape of data. We compute topological and standard descriptors of model simulations generated by different parameter values. We show that TDA of model simulation data stratifies parameter space into regions with similar vessel morphology. The methodologies proposed here are widely applicable to other synthetic and experimental data including wound healing, development, and plant biology.

Agent-Based Learning Model for the Obesity Paradox in RCC

A recent study on the immunotherapy treatment of renal cell carcinoma reveals better outcomes in obese patients compared to lean subjects. This enigmatic contradiction has been explained, in the context of the debated obesity paradox, as the effect produced by the cell-cell interaction network on the tumor microenvironment during the immune response. To better understand this hypothesis, we provide a computational framework for the in silico study of the tumor behavior. The starting model of the tumor, based on the cell-cell interaction network, has been described as a multiagent system, whose simulation generates the hypothesized effects on the tumor microenvironment. The medical needs in the immunotherapy design meet the capabilities of a multiagent simulator to reproduce the dynamics of the cell-cell interaction network, meaning a reaction to environmental changes introduced through the experimental data.

Phenotypic Heterogeneity in Tumor Progression, and Its Possible Role in the Onset of Cancer

Heterogeneity among isogenic cells/individuals has been known for at least 150 years. Even Mendel, working on pea plants, realized that not all tall plants were identical. However, Mendel was more interested in the discontinuous variation between genetically distinct individuals. The concept of environment dictating distinct phenotypes among isogenic individuals has since been shown to impact the evolution of populations in numerous examples at different scales of life. In this review, we discuss how phenotypic heterogeneity and its evolutionary implications exist at all levels of life, from viruses to mammals. In particular, we discuss how a particular disease condition (cancer) is impacted by heterogeneity among isogenic cells, and propose a potential role that phenotypic heterogeneity might play toward the onset of the disease.

Enhancer-bound Nrf2 licenses HIF-1α transcription under hypoxia to promote cisplatin resistance in hepatocellular carcinoma cells

Tumor microenvironment is hypoxic, which can cause resistance to chemotherapy, but the detailed mechanisms remain elusive. Here we find that mild hypoxia (5% O₂) further increases cisplatin resistance in the already resistant HepG2/DDP but not the sensitive HepG2 cells. We find that Nrf2 is responsible for cisplatin resistance under hypoxia, as Nrf2 knockdown sensitizes HepG2/DDP cells while Nrf2 hyper-activation (though KEAP1 knockdown) increases resistance of HepG2 cells to cisplatin. Nrf2 binds to an enhancer element in the upstream of HIF-1α gene independently of hypoxia, promoting HIF-1α mRNA synthesis under hypoxic condition. As a result, Nrf2-dependent transcription counteracts HIF-1α degradation under mild hypoxia condition, leading to preferential cisplatin-resistance in HepG2/DDP cells. Our data suggest that Nrf2 regulation of HIF-1α could be an important mechanism for chemotherapy resistance in vivo.

Bicomponent polymeric micelles for pH-controlled delivery of doxorubicin

Stimuli-responsive drug delivery systems (DDSs) are expected to realize site-specific drug release and kill cancer cells selectively. In this study, a pH-responsive micelle was designed utilizing the pH-sensitivity of borate bonds formed between dopamine and boronic acid. First, methyl (polyethylene glycol)-block-polycaprolactone (mPEG-PCL) was conjugated with 4-cyano-4-(thiobenzoylthio)pentanoic acid (CTP) to obtain a macroinitiator. Two different segments poly(dopamine methacrylamide) (PDMA) and poly(vinylphenylboronic acid) (PVBA) were then grafted to the end of mPEG-PCL. Two triblock copolymers, mPEG-PCL-PDMA and mPEG-PCL-PVBA, were then obtained by reversible addition-fragmentation transfer (RAFT) polymerization. These copolymers and their mixture self-assembled in aqueous solution to form micelles that were able to load hydrophobic anticancer drug doxorubicin (DOX). These two-component micelles were found to be pH-sensitive, in contrast to the one-component micelles. Furthermore, MTT studies showed that the micelles were almost nontoxic. The DOX-loaded micelles showed cytotoxicity equivalent to that of DOX at high concentration. In vivo antitumor experiments showed that this pH-sensitive polymeric micellar system had an enhanced therapeutic effect on tumors. These two-component boronate-based pH micelles are universally applicable to the delivery of anticancer drugs, showing great potential for cancer therapy.

Renal clearable nanocarriers: Overcoming the physiological barriers for precise drug delivery and clearance

Physiological barriers encountered in the clinical translation of cancer nanomedicines inspire the community to more deeply understand nano-bio interactions in not only tumor microenvironment but also entire body and develop new nanocarriers to tackle these barriers. Renal clearable nanocarriers are one kind of these newly emerged drug delivery systems (DDSs), which enable drugs to rapidly penetrate into the tumor cores with no need of long blood retention and escape macrophage uptake in the meantime they can also enhance body elimination of non-targeted anticancer drugs. As a result, they can improve therapeutic efficacies and reduce side effects of anticancer drugs. Not limited to anticancer drugs, diagnostic agents can also be achieved with these renal clearable DDSs, which might also be applied to improve the precision in the gene editing and immunotherapy in the future.

Efficient vasculature investment in tissues can be determined without global information

Cells are the fundamental building blocks of organs and tissues. Information and mass flow through cellular contacts in these structures is vital for the orchestration of organ function. Constraints imposed by packing and cell immobility limit intercellular communication, particularly as organs and organisms scale up to greater sizes. In order to transcend transport limitations, delivery systems including vascular and respiratory systems evolved to facilitate the movement of matter and information. The construction of these delivery systems has an associated cost, as vascular elements do not perform the metabolic functions of the organs they are part of. This study investigates a fundamental trade-off in vascularization in multicellular tissues: the reduction of path lengths for communication versus the cost associated with producing vasculature. Biologically realistic generative models, using multicellular templates of different dimensionalities, revealed a limited advantage to the vascularization of two-dimensional tissues. Strikingly, scale-free improvements in transport efficiency can be achieved even in the absence of global knowledge of tissue organization. A point of diminishing returns in the investment of additional vascular tissue to the increased reduction of path length in 2.5- and three-dimensional tissues was identified. Applying this theory to experimentally determined biological tissue structures, we show the possibility of a co-dependency between the method used to limit path length and the organization of cells it acts upon. These results provide insight as to why tissues are or are not vascularized in nature, the robustness of developmental generative mechanisms and the extent to which vasculature is advantageous in the support of organ function.

Tumor Development and Angiogenesis in Adult Brain Tumor: Glioblastoma

Angiogenesis is the growth of new capillaries from the preexisting blood vessels. Glioblastoma (GBM) tumors are highly vascularized tumors, and glioma growth depends on the formation of new blood vessels. Angiogenesis is a complex process involving proliferation, migration, and differentiation of vascular endothelial cells (ECs) under the stimulation of specific signals. It is controlled by the balance between its promoting and inhibiting factors. Various angiogenic factors and genes have been identified that stimulate glioma angiogenesis. Therefore, attention has been directed to anti-angiogenesis therapy in which glioma proliferation is inhibited by inhibiting the formation of new tumor vessels using angiogenesis inhibitory factors and drugs. Here, in this review, we highlight and summarize the various molecular mediators that regulate GBM angiogenesis with focus on recent clinical research on the potential of exploiting angiogenic pathways as a strategy in the treatment of GBM patients.

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.

Multiscale modeling of solid stress and tumor cell invasion in response to dynamic mechanical microenvironment

Mathematical models can provide a quantitatively sophisticated description of tumor cell (TC) behaviors under mechanical microenvironment and help us better understand the role of specific biophysical factors based on their influences on the TC behaviors. To this end, we propose an off-lattice cell-based multiscale mathematical model to describe the dynamic growth-induced solid stress during tumor progression and investigate the influence of the mechanical microenvironment on TC invasion. At the cellular level, cell-cell and cell-matrix interactive forces depend on the mechanical properties of the cells and the cancer-associated fibroblasts in the stroma, respectively. The constitutive relationship between the interactive forces and cell migrations obeys the Hooke's law and damping effects. At the tissue level, the integrated growth-induced forces caused by proliferating cells within the simulation region are balanced by the external forces applied by the surrounding host tissues. Then, the cell movements are calculated according to the Newton's second law of motion, and the morphology of TC invasion is updated. The simulation results reveal the continuous changes of the macroscopic mechanical forces due to the interactions among the structural components and the microscopic environmental factors. Moreover, the simulation results demonstrate the adverse effect of the stiffness of tumor tissue on tumor growth and invasion. A decrease in the stiffness of tumor and matrix can promote TCs to proliferate at a much faster rate and invade into the surrounding healthy tissue more easily, whereas an increase in the stiffness can lead to an aggressive morphology of tumor invasion. We envision that the proposed model can be served as a quantitative theoretical platform to study the underlying biophysical role of the mechanical microenvironmental factors during tumor invasion and metastasis.

Optical-resolution photoacoustic microscopy for monitoring vascular normalization during anti-angiogenic therapy

Monitoring the changes in tumor vascularity is important for anti-angiogenic therapy assessment with therapeutic implications. However, monitoring vascularity is quite challenging due to the lack of appropriate imaging techniques. Here, we describe a non-invasive imaging technique using optical-resolution photoacoustic microscopy (OR-PAM) to track vascular changes in prostate cancer treated with an anti-angiogenic agent, DC101, on a mouse ear xenograft model. Approximately 1-3 days after the initial therapy, OR-PAM imaging detected tumor vascular changes such as reduced vessel tortuosity, decreased vessel diameter and homogenized intratumoral vessel distribution. These observations indicated vessel normalization, which was pathologically validated as increased fractional pericyte coverage, functional perfusion and drug delivery of the vessels. After four DC101 interventions, OR-PAM imaging eventually revealed intratumoral vessel regression. Therefore, OR-PAM imaging of the vasculature offers a promising method to study anti-angiogenic drug mechanisms of action in vivo and holds potential in monitoring and guiding anti-angiogenic therapy.

Glioblastoma multiforme restructures the topological connectivity of cerebrovascular networks

Glioblastoma multiforme alters healthy tissue vasculature by inducing angiogenesis and vascular remodeling. To fully comprehend the structural and functional properties of the resulting vascular network, it needs to be studied collectively by considering both geometric and topological properties. Utilizing Single Plane Illumination Microscopy (SPIM), the detailed capillary structure in entire healthy and tumor-bearing mouse brains could be resolved in three dimensions. At the scale of the smallest capillaries, the entire vascular systems of bulk U87- and GL261-glioblastoma xenografts, their respective cores, and healthy brain hemispheres were modeled as complex networks and quantified with fundamental topological measures. All individual vessel segments were further quantified geometrically and modular clusters were uncovered and characterized as meta-networks, facilitating an analysis of large-scale connectivity. An inclusive comparison of large tissue sections revealed that geometric properties of individual vessels were altered in glioblastoma in a relatively subtle way, with high intra- and inter-tumor heterogeneity, compared to the impact on the vessel connectivity. A network topology analysis revealed a clear decomposition of large modular structures and hierarchical network organization, while preserving most fundamental topological classifications, in both tumor models with distinct growth patterns. These results augment our understanding of cerebrovascular networks and offer a topological assessment of glioma-induced vascular remodeling. The findings may help understand the emergence of hypoxia and necrosis, and prove valuable for therapeutic interventions such as radiation or antiangiogenic therapy.

Mechanical Forces in Tumor Angiogenesis

A defining hallmark of cancer and cancer development is upregulated angiogenesis. The vasculature formed in tumors is structurally abnormal, not organized in the conventional hierarchical arrangement, and more permeable than normal vasculature. These features contribute to leaky, tortuous, and dilated blood vessels, which act to create heterogeneous blood flow, compression of vessels, and elevated interstitial fluid pressure. As such, abnormalities in the tumor vasculature not only affect the delivery of nutrients and oxygen to the tumor, but also contribute to creating an abnormal tumor microenvironment that further promotes tumorigenesis. The role of chemical signaling events in mediating tumor angiogenesis has been well researched; however, the relative contribution of physical cues and mechanical regulation of tumor angiogenesis is less understood. Growing research indicates that the physical microenvironment plays a significant role in tumor progression and promoting abnormal tumor vasculature. Here, we review how mechanical cues found in the tumor microenvironment promote aberrant tumor angiogenesis. Specifically, we discuss the influence of matrix stiffness and mechanical stresses in tumor tissue on tumor vasculature, as well as the mechanosensory pathways utilized by endothelial cells to respond to the physical cues found in the tumor microenvironment. We also discuss the impact of the resulting aberrant tumor vasculature on tumor progression and therapeutic treatment.

Dexamethasone Increases Cisplatin-Loaded Nanocarrier Delivery and Efficacy in Metastatic Breast Cancer by Normalizing the Tumor Microenvironment

Dexamethasone is a glucocorticoid steroid with anti-inflammatory properties used to treat many diseases, including cancer, in which it helps manage various side effects of chemo-, radio-, and immunotherapies. Here, we investigate the tumor microenvironment (TME)-normalizing effects of dexamethasone in metastatic murine breast cancer (BC). Dexamethasone normalizes vessels and the extracellular matrix, thereby reducing interstitial fluid pressure, tissue stiffness, and solid stress. In turn, the penetration of 13 and 32 nm dextrans, which represent nanocarriers (NCs), is increased. A mechanistic model of fluid and macromolecule transport in tumors predicts that dexamethasone increases NC penetration by increasing interstitial hydraulic conductivity without significantly reducing the effective pore diameter of the vessel wall. Also, dexamethasone increases the tumor accumulation and efficacy of ∼30 nm polymeric micelles containing cisplatin (CDDP/m) against murine models of primary BC and spontaneous BC lung metastasis, which also feature a TME with abnormal mechanical properties. These results suggest that pretreatment with dexamethasone before NC administration could increase efficacy against primary tumors and metastases.

Tuning the In Vivo Transport of Anticancer Drugs Using Renal-Clearable Gold Nanoparticles

Precise control of in vivo transport of anticancer drugs in normal and cancerous tissues with engineered nanoparticles is key to the future success of cancer nanomedicines in clinics. This requires a fundamental understanding of how engineered nanoparticles impact the targeting-clearance and permeation-retention paradoxes in the anticancer-drug delivery. Herein, we systematically investigated how renal-clearable gold nanoparticles (AuNPs) affect the permeation, distribution, and retention of the anticancer drug doxorubicin in both cancerous and normal tissues. Renal-clearable AuNPs retain the advantages of the free drug, including rapid tumor targeting and high tumor vascular permeability. The renal-clearable AuNPs also accelerated body clearance of off-target drug via renal elimination. These results clearly indicate that diverse in vivo transport behaviors of engineered nanoparticles can be used to reconcile long-standing paradoxes in the anticancer drug delivery.

Tumor Ensemble-Based Modeling and Visualization of Emergent Angiogenic Heterogeneity in Breast Cancer

There is a critical need for new tools to investigate the spatio-temporal heterogeneity and phenotypic alterations that arise in the tumor microenvironment. However, computational investigations of emergent inter- and intra-tumor angiogenic heterogeneity necessitate 3D microvascular data from 'whole-tumors' as well as "ensembles" of tumors. Until recently, technical limitations such as 3D imaging capabilities, computational power and cost precluded the incorporation of whole-tumor microvascular data in computational models. Here, we describe a novel computational approach based on multimodality, 3D whole-tumor imaging data acquired from eight orthotopic breast tumor xenografts (i.e. a tumor 'ensemble'). We assessed the heterogeneous angiogenic landscape from the microvascular to tumor ensemble scale in terms of vascular morphology, emergent hemodynamics and intravascular oxygenation. We demonstrate how the abnormal organization and hemodynamics of the tumor microvasculature give rise to unique microvascular niches within the tumor and contribute to inter- and intra-tumor heterogeneity. These tumor ensemble-based simulations together with unique data visualization approaches establish the foundation of a novel 'cancer atlas' for investigators to develop their own in silico systems biology applications. We expect this hybrid image-based modeling framework to be adaptable for the study of other tissues (e.g. brain, heart) and other vasculature-dependent diseases (e.g. stroke, myocardial infarction).

Epigallocatechin-3-gallate-induced vascular normalization in A549-cell xenograft-bearing nude mice: therapeutic efficacy in combination with chemotherapy

Purpose

Large-scale studies have revealed that appropriate antiangiogenic treatment enables the recovery of the normal structure and function of solid tumor vessels. Epigallocatechin-3-gallate (EGCG), a natural extract of green tea, has multiple effects on angiogenesis. However, normalization of blood vessels due to natural ingredients has not yet been reported. Therefore, we examined the microvasculature, microenvironment, and efficacy of EGCG combined with chemotherapy in a xenograft model.

Methods

We treated A549 cell (human lung adenocarcinoma cell line) xenograft-bearing nude mice with EGCG in vivo. CD31, αSMA, and collagen IV were labeled and detected using quantum-dot double-labeled immunofluorescence to measure microvessel density, microvessel pericyte-coverage index, and collagen IV expression. Vessel-perfusion function was determined by lectin injection, permeability by Evans blue extravasation, interstitial fluid pressure using the wick-in-needle technique, and hypoxia levels using a polarographic electrode and immunohistochemical pimonidazole labeling. Cisplatin concentration in tumor tissue was detected using graphite-furnace atomic absorption spectrophotometry. Xenograft mice were randomized into five groups: treated with saline, cisplatin, EGCG, EGCG + cisplatin on day 1, or EGCG + cisplatin during the vascular normalization window. Tumor-growth delay and tumor-suppression rate were measured to evaluate tumor growth.

Results

EGCG treatment in vivo caused temporary changes, including transient depression of microvessel density, microvessel pericyte-coverage index, and collagen IV expression, transient elevation of vessel perfusion and permeability, and decreased interstitial fluid pressure and hypoxia. During vascular normalization, pretreatment with EGCG increased cisplatin concentration in tumor tissue compared with treatment with cisplatin only. Tumor-growth delay after treatment in the five groups during the vascular normalization window was 6.3±1.51, 7.5±1.57, 8.3±1.79, 12.1±1.35, and 15.4±1.99 days, indicating synergistic EGCG–cisplatin effects, especially during the vascular normalization window (P<0.01).

Conclusion

EGCG-induced vascular normalization in human lung adenocarcinoma may be a novel modality for enhancing chemotherapy effects.