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.

Dynamic heterogeneity in COVID-19: Insights from a mathematical model.

Critical illness, such as severe COVID-19, is heterogenous in presentation and treatment response. However, it remains possible that clinical course may be influenced by dynamic and/or random events such that similar patients subject to similar injuries may yet follow different trajectories. We deployed a mechanistic mathematical model of COVID-19 to determine the range of possible clinical courses after SARS-CoV-2 infection, which may follow from specific changes in viral properties, immune properties, treatment modality and random external factors such as initial viral load. We find that treatment efficacy and baseline patient or viral features are not the sole determinant of outcome. We found patients with enhanced innate or adaptive immune responses can experience poor viral control, resolution of infection or non-infectious inflammatory injury depending on treatment efficacy and initial viral load. Hypoxemia may result from poor viral control or ongoing inflammation despite effective viral control. Adaptive immune responses may be inhibited by very early effective therapy, resulting in viral load rebound after cessation of therapy. Our model suggests individual disease course may be influenced by the interaction between external and patient-intrinsic factors. These data have implications for the reproducibility of clinical trial cohorts and timing of optimal treatment.

A Multiresolution Approach with Method-Informed Statistical Analysis for Quantifying Lymphatic Pumping Dynamics.

Despite significant strides in lymphatic system imaging, the timely diagnosis of lymphatic disorders remains elusive. One main cause for this is the absence of standardized, quantitative methods for real-time analysis of lymphatic contractility. Here, we address this unmet need by combining near-infrared lymphangiography imaging with an innovative analytical workflow. We combined data acquisition, signal processing, and statistical analysis to integrate traditional peak and-valley with advanced wavelet time-frequency analyses. Decision theory was used to evaluate the primary drivers of attributable variance in lymphangiography measurements to generate a strategy for optimizing the number of repeat measurements needed per subject to increase measurement reliability. This approach not only offers detailed insights into lymphatic pumping behaviors across species, sex and age, but also significantly boosts the reliability of these measurements by incorporating multiple regions of interest and evaluating the lymphatic system under various gravitational loads. By addressing the critical need for improved imaging and quantification methods, our study offers a new standard approach for the imaging and analysis of lymphatic function that can improve our understanding, diagnosis, and treatment of lymphatic diseases. The results highlight the importance of comprehensive data acquisition strategies to fully capture the dynamic behavior of the lymphatic system.

Mechanics of Lymphatic Pumping and Lymphatic Function.

The lymphatic system plays a crucial role in maintaining tissue fluid balance, immune surveillance, and the transport of lipids and macromolecules. Lymph is absorbed by initial lymphatics and then driven through lymph nodes and to the blood circulation by the contraction of collecting lymphatic vessels. Intraluminal valves in collecting lymphatic vessels ensure the unidirectional flow of lymph centrally. The lymphatic muscle cells that invest in collecting lymphatic vessels impart energy to propel lymph against hydrostatic pressure gradients and gravity. A variety of mechanical and biochemical stimuli modulate the contractile activity of lymphatic vessels. This review focuses on the recent advances in our understanding of the mechanisms involved in regulating and collecting lymphatic vessel pumping in normal tissues and the association between lymphatic pumping, infection, inflammatory disease states, and lymphedema.

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.

In silico clinical studies for optimal COVID-19 vaccination schedules in patients with cancer.

This study introduces a tailored COVID-19 model for patients with cancer, incorporating viral variants and immune-response dynamics. The model aims to optimize vaccination strategies, contributing to personalized healthcare for vulnerable groups.

Reprogramming the Intrahepatic Cholangiocarcinoma Immune Microenvironment by Chemotherapy and CTLA-4 Blockade Enhances Anti-PD-1 Therapy.

Intrahepatic cholangiocarcinoma (ICC) has limited therapeutic options and a dismal prognosis. Adding blockade of the anti-programmed cell death protein (PD)-1 pathway to gemcitabine/cisplatin chemotherapy has recently shown efficacy in biliary tract cancers but with low response rates. Here, we studied the effects of anti-cytotoxic T lymphocyte antigen (CTLA)-4 when combined with anti-PD-1 and gemcitabine/cisplatin in orthotopic murine models of ICC. This combination therapy led to substantial survival benefits and reduction of morbidity in two aggressive ICC models that were resistant to immunotherapy alone. Gemcitabine/cisplatin treatment increased tumor-infiltrating lymphocytes and normalized the ICC vessels and, when combined with dual CTLA-4/PD-1 blockade, increased the number of activated CD8+Cxcr3+IFNγ+ T cells. CD8+ T cells were necessary for the therapeutic benefit because the efficacy was compromised when CD8+ T cells were depleted. Expression of Cxcr3 on CD8+ T cells is necessary and sufficient because CD8+ T cells from Cxcr3+/+ but not Cxcr3-/- mice rescued efficacy in T cell‒deficient mice. Finally, rational scheduling of anti-CTLA-4 "priming" with chemotherapy followed by anti-PD-1 therapy achieved equivalent efficacy with reduced overall drug exposure. These data suggest that this combination approach should be clinically tested to overcome resistance to current therapies in ICC patients.

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.

Dissecting the Impact of the Gut Microbiome on Cancer Immunotherapy.

The gut microbiome has emerged as a key regulator of response to cancer immunotherapy. However, there is a gap in our understanding of the underlying mechanisms by which the microbiome influences immunotherapy. To this end, we developed a mathematical model based on i) gut microbiome data derived from preclinical studies on melanomas after fecal microbiota transplant, ii) mechanistic modeling of antitumor immune response, and iii) robust association analysis of murine and human microbiome profiles with model-predicted immune profiles. Using our model, we could distill the complexity of these murine and human studies on microbiome modulation in terms of just two model parameters: the activation and killing rate constants of immune cells. We further investigated associations between specific bacterial taxonomies and antitumor immunity and immunotherapy efficacy. This model can guide the design of studies to refine and validate mechanistic links between the microbiome and immune system.

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.

Lymphatic muscle cells are unique cells that undergo aging induced changes.

Lymphatic muscle cells (LMCs) within the wall of collecting lymphatic vessels exhibit tonic and autonomous phasic contractions, which drive active lymph transport to maintain tissue-fluid homeostasis and support immune surveillance. Damage to LMCs disrupts lymphatic function and is related to various diseases. Despite their importance, knowledge of the transcriptional signatures in LMCs and how they relate to lymphatic function in normal and disease contexts is largely missing. We have generated a comprehensive transcriptional single-cell atlas-including LMCs-of collecting lymphatic vessels in mouse dermis at various ages. We identified genes that distinguish LMCs from other types of muscle cells, characterized the phenotypical and transcriptomic changes in LMCs in aged vessels, and uncovered a pro-inflammatory microenvironment that suppresses the contractile apparatus in advanced-aged LMCs. Our findings provide a valuable resource to accelerate future research for the identification of potential drug targets on LMCs to preserve lymphatic vessel function as well as supporting studies to identify genetic causes of primary lymphedema currently with unknown molecular explanation.

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.

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.

A Self-Assembly Method for Creating Vascularized Tumor Explants Using Biomaterials for 3D Culture.

Validation of potential therapeutic targets in cancer requires functional live assays that recapitulate the biology, anatomy, and physiology of human tumors. We present a methodology for maintaining mouse and patient tumor samples ex vivo for in vitro drug-screening as well as for the guidance of patient-specific chemotherapies. The harvested tumor biopsy, excised from mice or patients, is integrated into a support tissue that includes extended stroma and vasculature. The methodology is more representative than tissue culture assays, faster than patient-derived xenograft models, easy to implement, amenable to high-throughput assays and does not carry the ethical issues or expense associated with animal studies. Our physiologically relevant model can be successfully used for high-throughput drug screening.

Endothelial mechanobiology in atherosclerosis.

Cardiovascular disease (CVD) is a serious health challenge, causing more deaths worldwide than cancer. The vascular endothelium, which forms the inner lining of blood vessels, plays a central role in maintaining vascular integrity and homeostasis and is in direct contact with the blood flow. Research over the past century has shown that mechanical perturbations of the vascular wall contribute to the formation and progression of atherosclerosis. While the straight part of the artery is exposed to sustained laminar flow and physiological high shear stress, flow near branch points or in curved vessels can exhibit 'disturbed' flow. Clinical studies as well as carefully controlled in vitro analyses have confirmed that these regions of disturbed flow, which can include low shear stress, recirculation, oscillation, or lateral flow, are preferential sites of atherosclerotic lesion formation. Because of their critical role in blood flow homeostasis, vascular endothelial cells (ECs) have mechanosensory mechanisms that allow them to react rapidly to changes in mechanical forces, and to execute context-specific adaptive responses to modulate EC functions. This review summarizes the current understanding of endothelial mechanobiology, which can guide the identification of new therapeutic targets to slow or reverse the progression of atherosclerosis.

Reprogramming Intrahepatic Cholangiocarcinoma Immune Microenvironment by Chemotherapy and CTLA-4 Blockade Enhances Anti-PD1 Therapy.

Intrahepatic cholangiocarcinoma (ICC) has limited therapeutic options and a dismal prognosis. Anti-PD-L1 immunotherapy combined with gemcitabine/cisplatin chemotherapy has recently shown efficacy in biliary tract cancers, but responses are seen only in a minority of patients. Here, we studied the roles of anti-PD1 and anti-CTLA-4 immune checkpoint blockade (ICB) therapies when combined with gemcitabine/cisplatin and the mechanisms of treatment benefit in orthotopic murine ICC models. We evaluated the effects of the combined treatments on ICC vasculature and immune microenvironment using flow cytometry analysis, immunofluorescence, imaging mass cytometry, RNA-sequencing, qPCR, and in vivo T-cell depletion and CD8+ T-cell transfer using orthotopic ICC models and transgenic mice. Combining gemcitabine/cisplatin with anti-PD1 and anti-CTLA-4 antibodies led to substantial survival benefits and reduction of morbidity in two aggressive ICC models, which were ICB-resistant. Gemcitabine/cisplatin treatment increased the frequency of tumor-infiltrating lymphocytes and normalized the ICC vessels, and when combined with dual CTLA-4/PD1 blockade, increased the number of activated CD8+Cxcr3+IFN-γ+ T-cells. Depletion of CD8+ but not CD4+ T-cells compromised efficacy. Conversely, CD8+ T-cell transfer from Cxcr3-/- versus Cxcr3+/+ mice into Rag1-/- immunodeficient mice restored the anti-tumor effect of gemcitabine/cisplatin/ICB combination therapy. Finally, rational scheduling of the ICBs (anti-CTLA-4 "priming") with chemotherapy and anti-PD1 therapy achieved equivalent efficacy with continuous dosing while reducing overall drug exposure. In summary, gemcitabine/cisplatin chemotherapy normalizes vessel structure, increases activated T-cell infiltration, and enhances anti-PD1/CTLA-4 immunotherapy efficacy in aggressive murine ICC. This combination approach should be clinically tested to overcome resistance to current therapies in ICC patients.

Inhibition of CXCR4 Enhances the Efficacy of Radiotherapy in Metastatic Prostate Cancer Models.

Radiotherapy (RT) is a standard treatment for patients with advanced prostate cancer (PCa). Previous preclinical studies showed that SDF1α/CXCR4 axis could mediate PCa metastasis (most often to the bones) and cancer resistance to RT. We found high levels of expression for both SDF1α and its receptor CXCR4 in primary and metastatic PCa tissue samples. In vitro analyses using PCa cells revealed an important role of CXCR4 in cell invasion but not radiotolerance. Pharmacologic inhibition of CXCR4 using AMD3100 showed no efficacy in orthotopic primary and bone metastatic PCa models. However, when combined with RT, AMD3100 potentiated the effect of local single-dose RT (12 Gy) in both models. Moreover, CXCR4 inhibition also reduced lymph node metastasis from primary PCa. Notably, CXCR4 inhibition promoted the normalization of bone metastatic PCa vasculature and reduced tissue hypoxia. In conclusion, the SDF1α/CXCR4 axis is a potential therapeutic target in metastatic PCa patients treated with RT.

Mechanistic model for booster doses effectiveness in healthy, cancer, and immunosuppressed patients infected with SARS-CoV-2.

SARS-CoV-2 vaccines are effective at limiting disease severity, but effectiveness is lower among patients with cancer or immunosuppression. Effectiveness wanes with time and varies by vaccine type. Moreover, previously prescribed vaccines were based on the ancestral SARS-CoV-2 spike-protein that emerging variants may evade. Here, we describe a mechanistic mathematical model for vaccination-induced immunity. We validate it with available clinical data and use it to simulate the effectiveness of vaccines against viral variants with lower antigenicity, increased virulence, or enhanced cell binding for various vaccine platforms. The analysis includes the omicron variant as well as hypothetical future variants with even greater immune evasion of vaccine-induced antibodies and addresses the potential benefits of the new bivalent vaccines. We further account for concurrent cancer or underlying immunosuppression. The model confirms enhanced immunogenicity following booster vaccination in immunosuppressed patients but predicts ongoing booster requirements for these individuals to maintain protection. We further studied the impact of variants on immunosuppressed individuals as a function of the interval between multiple booster doses. Our model suggests possible strategies for future vaccinations and suggests tailored strategies for high-risk groups.

Experimental and Theoretical Model of Single Vessel Minimally Invasive Micro-Laser Ablation: Inducing Microvascular Network Remodeling and Blood Flow Redistribution Without Compromising Host Tissue Function.

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.5x and 3.3x, 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 match the post-ablation microvascular remodeling patterns.

Wnt inhibition alleviates resistance to immune checkpoint blockade in glioblastoma.

Wnt signaling plays a critical role in the progression and treatment outcome of glioblastoma (GBM). Here, we identified WNT7b as a heretofore unknown mechanism of resistance to immune checkpoint inhibition (αPD1) in GBM patients and murine models. Acquired resistance to αPD1 was found to be associated with the upregulation of Wnt7b and ß-catenin protein levels in GBM in patients and in a clinically relevant, stem-rich GBM model. Combining the porcupine inhibitor WNT974 with αPD1 prolonged the survival of GBM-bearing mice. However, this combination had a dichotomous response, with a subset of tumors showing refractoriness. WNT974 and αPD1 expanded a subset of DC3-like dendritic cells (DCs) and decreased the granulocytic myeloid-derived suppressor cells (gMDSCs) in the tumor microenvironment (TME). By contrast, monocytic MDSCs (mMDSCs) increased, while T-cell infiltration remained unchanged, suggesting potential TME-mediated resistance. Our preclinical findings warrant the testing of Wnt7b/ß-catenin combined with αPD1 in GBM patients with elevated Wnt7b/ß-catenin signaling.