Fluid dynamics and leukocyte transit in the lymphatic system.

The lymphatic system plays a vital role in maintaining fluid balance in living tissue and serves as a pathway for the transport of antigen, immune cells, and metastatic cancer cells. In this study, we investigate how the movement of cells through a contracting lymphatic vessel differs from steady flow, using a lattice Boltzmann-based computational model. Our model consists of cells carried by flow in a 2D vessel with regularly spaced, bi-leaflet valves that ensure net downstream flow as the vessel walls contract autonomously in response to calcium and nitric oxide levels regulated by stretch and shear stress levels. The orientation of the vessel with respect to gravity, which may oppose or assist fluid flow, significantly modulates cellular motion due to its effect on the contraction dynamics of the vessel, even when the cells themselves are neutrally buoyant. Additionally, our model shows that cells are carried along with the flow, but when the vessel is actively contracting, they move faster than the average fluid velocity. We also find that the fluid forces cause significant deformation of the compliant cells, especially in the vicinity of the valves. Our study highlights the importance of considering the complex, transient flows near the valves in understanding cellular motion in lymphatic vessels.

Experimental and theoretical model of microvascular network remodeling and blood flow redistribution following minimally invasive microvessel laser ablation.

Overly dense microvascular networks are treated by selective reduction of vascular elements. Inappropriate manipulation of microvessels could result in loss of host tissue function or a worsening of the clinical problem. Here, experimental, and computational models were developed to induce blood flow changes via selective artery and vein laser ablation and study the compensatory collateral flow redistribution and vessel diameter remodeling. The microvasculature was imaged non-invasively by bright-field and multi-photon laser microscopy, and optical coherence tomography pre-ablation and up to 30 days post-ablation. A theoretical model of network remodeling was developed to compute blood flow and intravascular pressure and identify vessels most susceptible to changes in flow direction. The skin microvascular remodeling patterns were consistent among the five specimens studied. Significant remodeling occurred at various time points, beginning as early as days 1-3 and continuing beyond day 20. The remodeling patterns included collateral development, venous and arterial reopening, and both outward and inward remodeling, with variations in the time frames for each mouse. In a representative specimen, immediately post-ablation, the average artery and vein diameters increased by 14% and 23%, respectively. At day 20 post-ablation, the maximum increases in arterial and venous diameters were 2.5× and 3.3×, respectively. By day 30, the average artery diameter remained 11% increased whereas the vein diameters returned to near pre-ablation values. Some arteries regenerated across the ablation sites via endothelial cell migration, while veins either reconnected or rerouted flow around the ablation site, likely depending on local pressure driving forces. In the intact network, the theoretical model predicts that the vessels that act as collaterals after flow disruption are those most sensitive to distant changes in pressure. The model results correlate with the post-ablation microvascular remodeling patterns.

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.

Cancer immunotherapy responses persist after lymph node resection.

Surgical removal of lymph nodes (LNs) to prevent metastatic recurrence, including sentinel lymph node biopsy (SLNB) and completion lymph node dissection (CLND), are performed in routine practice. However, it remains controversial whether removing LNs which are critical for adaptive immune responses impairs immune checkpoint blockade (ICB) efficacy. Here, our retrospective analysis demonstrated that stage III melanoma patients retain robust response to anti-PD1 inhibition after CLND. Using orthotopic murine mammary carcinoma and melanoma models, we show that responses to ICB persist in mice after TDLN resection. Mechanistically, after TDLN resection, antigen can be re-directed to distant LNs, which extends the responsiveness to ICB. Strikingly, by evaluating head and neck cancer patients treated by neoadjuvant durvalumab and irradiation, we show that distant LNs (metastases-free) remain reactive in ICB responders after tumor and disease-related LN resection, hence, persistent anti-cancer immune reactions in distant LNs. Additionally, after TDLN dissection in murine models, ICB delivered to distant LNs generated greater survival benefit, compared to systemic administration. In complete responders, anti-tumor immune memory induced by ICB was systemic rather than confined within lymphoid organs. Based on these findings, we constructed a computational model to predict free antigen trafficking in patients that will undergo LN dissection.

Transport Barriers Influence the Activation of Anti-Tumor Immunity: A Systems Biology Analysis.

Effective anti-cancer immune responses require activation of one or more naïve T cells. If the correct naïve T cell encounters its cognate antigen presented by an antigen presenting cell, then the T cell can activate and proliferate. Here, mathematical modeling is used to explore the possibility that immune activation in lymph nodes is a rate-limiting step in anti-cancer immunity and can affect response rates to immune checkpoint therapy. The model provides a mechanistic framework for optimizing cancer immunotherapy and developing testable solutions to unleash anti-tumor immune responses for more patients with cancer. The results show that antigen production rate and trafficking of naïve T cells into the lymph nodes are key parameters and that treatments designed to enhance tumor antigen production can improve immune checkpoint therapies. The model underscores the potential of radiation therapy in augmenting tumor immunogenicity and neoantigen production for improved ICB therapy, while emphasizing the need for careful consideration in cases where antigen levels are already sufficient to avoid compromising the immune response.

Regeneration of collecting lymphatic vessels following injury.

Secondary lymphedema is a debilitating condition driven by impaired regeneration of lymphatic vasculature following lymphatic injury, surgical removal of lymph nodes in cancer patients or infection. However, the extent to which collecting lymphatic vessels regenerate following injury remains unclear. Here, we employed a novel mouse model of lymphatic injury in combination with state-of-the-art lymphatic imaging to demonstrate that the implantation of an optimized fibrin gel following lymphatic vessel injury leads to the growth and reconnection of the injured lymphatic vessel network, resulting in the restoration of lymph flow to the draining node. Intriguingly, we found that fibrin implantation elevates the tissue levels of CCL5, a potent macrophage-recruiting chemokine. Notably, CCL5-KO mice displayed a reduced ability to reconnect injured vessels following fibrin gel implantation. These novel findings shed light on the mechanisms underlying lymphatic regeneration and suggest that enhancing CCL5 signaling may be a promising therapeutic strategy for enhancing lymphatic regeneration.

Normalizing tumor microenvironment with nanomedicine and metronomic therapy to improve immunotherapy.

Nanomedicine offered hope for improving the treatment of cancer but the survival benefits of the clinically approved nanomedicines are modest in many cases when compared to conventional chemotherapy. Metronomic therapy, defined as the frequent, low dose administration of chemotherapeutics - is being tested in clinical trials as an alternative to the conventional maximum tolerated dose (MTD) chemotherapy schedule. Although metronomic chemotherapy has not been clinically approved yet, it has shown better survival than MTD in many preclinical studies. When beneficial, metronomic therapy seems to be associated with normalization of the tumor microenvironment including improvements in tumor perfusion, tissue oxygenation and drug delivery as well as activation of the immune system. Recent preclinical studies suggest that nanomedicines can cause similar changes in the tumor microenvironment. Here, by employing a mathematical framework, we show that both approaches can serve as normalization strategies to enhance treatment. Furthermore, employing murine breast and fibrosarcoma tumor models as well as ultrasound shear wave elastography and contrast-enhanced ultrasound, we provide evidence that the approved nanomedicine Doxil can induce normalization in a dose-dependent manner by improving tumor perfusion as a result of tissue softening. Finally, we show that pretreatment with a normalizing dose of Doxil can improve the efficacy of immune checkpoint inhibition.

Combining microenvironment normalization strategies to improve cancer immunotherapy.

Advances in immunotherapy have revolutionized the treatment of multiple cancers. Unfortunately, tumors usually have impaired blood perfusion, which limits the delivery of therapeutics and cytotoxic immune cells to tumors and also results in hypoxia-a hallmark of the abnormal tumor microenvironment (TME)-that causes immunosuppression. We proposed that normalization of TME using antiangiogenic drugs and/or mechanotherapeutics can overcome these challenges. Recently, immunotherapy with checkpoint blockers was shown to effectively induce vascular normalization in some types of cancer. Although these therapeutic approaches have been used in combination in preclinical and clinical studies, their combined effects on TME are not fully understood. To identify strategies for improved immunotherapy, we have developed a mathematical framework that incorporates complex interactions among various types of cancer cells, immune cells, stroma, angiogenic molecules, and the vasculature. Model predictions were compared with the data from five previously reported experimental studies. We found that low doses of antiangiogenic treatment improve immunotherapy when the two treatments are administered sequentially, but that high doses are less efficacious because of excessive vessel pruning and hypoxia. Stroma normalization can further increase the efficacy of immunotherapy, and the benefit is additive when combined with vascular normalization. We conclude that vessel functionality dictates the efficacy of immunotherapy, and thus increased tumor perfusion should be investigated as a predictive biomarker of response to immunotherapy.

Solid stress in brain tumours causes neuronal loss and neurological dysfunction and can be reversed by lithium.

The compression of brain tissue by a tumour mass is believed to be a major cause of the clinical symptoms seen in patients with brain cancer. However, the biological consequences of these physical stresses on brain tissue are unknown. Here, via imaging studies in patients and by using mouse models of human brain tumours, we show that a subgroup of primary and metastatic brain tumours, classified as nodular on the basis of their growth pattern, exert solid stress on the surrounding brain tissue, causing a decrease in local vascular perfusion as well as neuronal death and impaired function. We demonstrate a causal link between solid stress and neurological dysfunction by applying and removing cerebral compression, which respectively mimic the mechanics of tumour growth and of surgical resection. We also show that, in mice, treatment with lithium reduces solid-stress-induced neuronal death and improves motor coordination. Our findings indicate that brain-tumour-generated solid stress impairs neurological function in patients, and that lithium as a therapeutic intervention could counter these effects.

Experimental and computational analyses reveal dynamics of tumor vessel cooption and optimal treatment strategies.

Cooption of the host vasculature is a strategy that some cancers use to sustain tumor progression without-or before-angiogenesis or in response to antiangiogenic therapy. Facilitated by certain growth factors, cooption can mediate tumor infiltration and confer resistance to antiangiogenic drugs. Unfortunately, this mode of tumor progression is difficult to target because the underlying mechanisms are not fully understood. Here, we analyzed the dynamics of vessel cooption during tumor progression and in response to antiangiogenic treatment in gliomas and brain metastases. We followed tumor evolution during escape from antiangiogenic treatment as cancer cells coopted, and apparently mechanically compressed, host vessels. To gain deeper understanding, we developed a mathematical model, which incorporated compression of coopted vessels, resulting in hypoxia and formation of new vessels by angiogenesis. Even if antiangiogenic therapy can block such secondary angiogenesis, the tumor can sustain itself by coopting existing vessels. Hence, tumor progression can only be stopped by combination therapies that judiciously block both angiogenesis and cooption. Furthermore, the model suggests that sequential blockade is likely to be more beneficial than simultaneous blockade.

Lymph node effective vascular permeability and chemotherapy uptake.

OBJECTIVE: Lymph node metastases are a poor prognostic factor. Additionally, responses of lymph node metastasis to therapy can be different from the primary tumor. Investigating the physiologic lymph node blood vasculature might give insight into the ability of systemic drugs to penetrate the lymph node, and thus into the differential effect of therapy between lymph node metastasis and primary tumors. Here, we measured effective vascular permeability of lymph node blood vessels and attempted to increase chemotherapy penetration by increasing effective vascular permeability. METHODS: We developed a novel three-dimensional method to measure effective vascular permeability in murine lymph nodes in vivo. VEGF-A was systemically administered to increase effective vascular permeability. Validated high-performance liquid chromatography protocols were used to measure chemotherapeutic drug concentrations in untreated and VEGF-A-treated lymph nodes, liver, spleen, brain, and blood. RESULTS: VEGF-A-treated lymph node blood vessel effective vascular permeability (mean 3.83 × 10-7  cm/s) was significantly higher than untreated lymph nodes (mean 9.87 × 10-8  cm/s). No difference was found in lymph node drug accumulation in untreated versus VEGF-A-treated mice. CONCLUSIONS: Lymph node effective vascular permeability can be increased (~fourfold) by VEGF-A. However, no significant increase in chemotherapy uptake was measured by pretreatment with VEGF-A.

Role of vascular normalization in benefit from metronomic chemotherapy.

Metronomic dosing of chemotherapy-defined as frequent administration at lower doses-has been shown to be more efficacious than maximum tolerated dose treatment in preclinical studies, and is currently being tested in the clinic. Although multiple mechanisms of benefit from metronomic chemotherapy have been proposed, how these mechanisms are related to one another and which one is dominant for a given tumor-drug combination is not known. To this end, we have developed a mathematical model that incorporates various proposed mechanisms, and report here that improved function of tumor vessels is a key determinant of benefit from metronomic chemotherapy. In our analysis, we used multiple dosage schedules and incorporated interactions among cancer cells, stem-like cancer cells, immune cells, and the tumor vasculature. We found that metronomic chemotherapy induces functional normalization of tumor blood vessels, resulting in improved tumor perfusion. Improved perfusion alleviates hypoxia, which reprograms the immunosuppressive tumor microenvironment toward immunostimulation and improves drug delivery and therapeutic outcomes. Indeed, in our model, improved vessel function enhanced the delivery of oxygen and drugs, increased the number of effector immune cells, and decreased the number of regulatory T cells, which in turn killed a larger number of cancer cells, including cancer stem-like cells. Vessel function was further improved owing to decompression of intratumoral vessels as a result of increased killing of cancer cells, setting up a positive feedback loop. Our model enables evaluation of the relative importance of these mechanisms, and suggests guidelines for the optimal use of metronomic therapy.

Measuring Vascular Permeability In Vivo.

Over the past decades, in vivo vascular permeability measurements have provided significant insight into vascular functions in physiological and pathophysiological conditions such as the response to pro- and anti-angiogenic signaling, abnormality of tumor vasculature and its normalization, and delivery and efficacy of therapeutic agents. Different approaches for vascular permeability measurements have been established. Here, we describe and discuss a conventional 2D imaging method to measure vascular permeability, which was originally documented by Gerlowski and Jain in 1986 (Microvasc Res 31:288-305, 1986) and further developed by Yuan et al. in the early 1990s (Microvasc Res 45:269-289, 1993; Cancer Res 54:352-3356, 1994), and our recently developed 3D imaging method, which advances the approach originally described by Brown et al. in 2001 (Nat Med 7:864-868, 2001).

Mathematical Model of Oxygen Transport in Tuberculosis Granulomas.

Pulmonary granulomas--the hallmark of Mycobacterium tuberculosis (MTB) infection--are dense cellular lesions that often feature regions of hypoxia and necrosis, partially due to limited transport of oxygen. Low oxygen in granulomas can impair the host immune response, while MTB are able to adapt and persist in hypoxic environments. Here, we used a physiologically based mathematical model of oxygen diffusion and consumption to calculate oxygen profiles within the granuloma, assuming Michaelis-Menten kinetics. An approximate analytical solution--using a priori and newly estimated parameters from experimental data in a rabbit model of tuberculosis--was able to predict the size of hypoxic and necrotic regions in agreement with experimental results from the animal model. Such quantitative understanding of transport limitations can inform future tuberculosis therapeutic strategies that may include adjunct host-directed therapies that facilitate oxygen and drug delivery for more effective treatment.

Solid stress and elastic energy as measures of tumour mechanopathology.

Solid stress and tissue stiffness affect tumour growth, invasion, metastasis and treatment. Unlike stiffness, which can be precisely mapped in tumours, the measurement of solid stresses is challenging. Here, we show that two-dimensional spatial mappings of solid stress and the resulting elastic energy in excised or in situ tumours with arbitrary shapes and wide size ranges can be obtained via three distinct and quantitative techniques that rely on the measurement of tissue displacement after disruption of the confining structures. Application of these methods in models of primary tumours and metastasis revealed that: (i) solid stress depends on both cancer cells and their microenvironment; (ii) solid stress increases with tumour size; and (iii) mechanical confinement by the surrounding tissue significantly contributes to intratumoural solid stress. Further study of the genesis and consequences of solid stress, facilitated by the engineering principles presented here, may lead to significant discoveries and new therapies.

Angiopoietin-4 increases permeability of blood vessels and promotes lymphatic dilation.

The angiopoietin (Ang) ligands are potential therapeutic targets for lymphatic related diseases, which include lymphedema and cancer. Ang-1 and Ang-2 functions are established, but those of Ang-4 are poorly understood. We used intravital fluorescence microscopy to characterize Ang-4 actions on T241 murine fibrosarcoma-associated vessels in mice. The diameters of lymphatic vessels draining Ang-4- or VEGF-C (positive control)-expressing tumors increased to 123 and 135 µm, respectively, and parental, mock-transduced (negative controls) and tumors expressing Ang-1 or Ang-2 remained at baseline (∼60 µm). Ang-4 decreased human dermal lymphatic endothelial cell (LEC) monolayer permeability by 27% while increasing human dermal blood endothelial cell (BEC) monolayer permeability by 200%. In vivo, Ang-4 stimulated a 4.5-fold increase in tumor-associated blood vessel permeability compared with control when measured using intravital quantitative multiphoton microscopy. Ang-4 activated receptor signaling in both LECs and BECs, evidenced by tyrosine kinase with Ig and endothelial growth factor homology domains-2 (TIE2) receptor, protein kinase B, and Erk1,2 phosphorylation detectable by immunoblotting. These data suggest that Ang-4 actions are mediated through cell-type-specific networks and that lymphatic vessel dilation occurs secondarily to increased vascular leakage. Ang-4 also promoted survival of LECs. Thus, blocking Ang-4 may prune the draining lymphatic vasculature and decrease interstitial fluid pressure (IFP) by reducing vascular permeability.

Towards Optimal Design of Cancer Nanomedicines: Multi-stage Nanoparticles for the Treatment of Solid Tumors.

Conventional drug delivery systems for solid tumors are composed of a nano-carrier that releases its therapeutic load. These two-stage nanoparticles utilize the enhanced permeability and retention (EPR) effect to enable preferential delivery to tumor tissue. However, the size-dependency of the EPR, the limited penetration of nanoparticles into the tumor as well as the rapid binding of the particles or the released cytotoxic agents to cancer cells and stromal components inhibit the uniform distribution of the drug and the efficacy of the treatment. Here, we employ mathematical modeling to study the effect of particle size, drug release rate and binding affinity on the distribution and efficacy of nanoparticles to derive optimal design rules. Furthermore, we introduce a new multi-stage delivery system. The system consists of a 20-nm primary nanoparticle, which releases 5-nm secondary particles, which in turn release the chemotherapeutic drug. We found that tuning the drug release kinetics and binding affinities leads to improved delivery of the drug. Our results also indicate that multi-stage nanoparticles are superior over two-stage nano-carriers provided they have a faster drug release rate and for high binding affinity drugs. Furthermore, our results suggest that smaller nanoparticles achieve better treatment outcome.

Scaling rules for diffusive drug delivery in tumor and normal tissues.

Delivery of blood-borne molecules and nanoparticles from the vasculature to cells in the tissue differs dramatically between tumor and normal tissues due to differences in their vascular architectures. Here we show that two simple measures of vascular geometry--δ(max) and λ--readily obtained from vascular images, capture these differences and link vascular structure to delivery in both tissue types. The longest time needed to bring materials to their destination scales with the square of δ(max), the maximum distance in the tissue from the nearest blood vessel, whereas λ, a measure of the shape of the spaces between vessels, determines the rate of delivery for shorter times. Our results are useful for evaluating how new therapeutic agents that inhibit or stimulate vascular growth alter the functional efficiency of the vasculature and more broadly for analysis of diffusion in irregularly shaped domains.