Simulating the Impact of Tumor Mechanical Forces on Glymphatic Networks in the Brain Parenchyma

Abstract

Background: The brain glymphatic system is currently being explored in the context of many neurological disorders and diseases, including traumatic brain injury, Alzheimer's disease, and ischemic stroke. However, little is known about the impact of brain tumors on glymphatic function. Mechanical forces generated during tumor development and growth may be responsible for compromised glymphatic transport pathways, reducing waste clearance and cerebrospinal fluid (CSF) transport in the brain parenchyma. One such force is solid stress, i.e., growth-induced forces from cell hyperproliferation and excess matrix deposition. Because there are no prior studies assessing the impact of tumor-derived solid stress on glymphatic system structure and performance in the brain parenchyma, this study serves to fill an important gap in the field. Methods: We adapted a previously developed Electrical Analog Model using MATLAB Simulink for glymphatic transport coupled with Finite Element Analysis for tumor mechanical stresses and strains in COMSOL. This allowed simulation of the impact of tumor mechanical force generation on fluid transport within brain parenchymal glymphatic units, which include paravascular spaces, astrocytic networks, interstitial spaces, and capillary basement membranes. We conducted a parametric analysis to compare the contributions of tumor size, tumor proximity, and ratio of glymphatic subunits to the stress and strain experienced by the glymphatic unit and corresponding reduction in flow rate of CSF. Results:. Mechanical stresses intensify with proximity to the tumor and increasing tumor size, highlighting the vulnerability of nearby glymphatic units to tumor-derived forces. Our stress and strain profiles reveal compressive deformation of these surrounding glymphatics and demonstrate that varying the relative contributions of astrocytes vs. interstitial spaces impact the resulting glymphatic structure significantly under tumor mechanical forces. Increased tumor size and proximity caused increased stress and strain across all glymphatic subunits, as does decreased astrocyte composition. Indeed, our model reveals an inverse correlation between extent of astrocyte contribution to the composition of the glymphatic unit and the resulting mechanical stress. This increased mechanical strain across the glymphatic unit decreases the venous efflux rate of CSF, dependent on the degree of strain and the specific glymphatic subunit of interest. For example, a 20% mechanical strain on capillary basement membranes does not significantly decrease venous efflux (2% decrease in flow rates), while the same magnitude of strain on astrocyte networks and interstitial spaces decreases efflux flow rates by 7% and 22%, respectively. Conclusion: Our simulations reveal that solid stress from brain tumors directly reduces glymphatic fluid transport, independently from biochemical effects from cancer cells. Understanding these pathophysiological implications is crucial for developing targeted interventions aimed at restoring effective waste clearance mechanisms in the brain. This study opens potential avenues for future experimental research in brain tumor-related glymphatic dysfunction.