Glioblastoma is the most common and deadly primary brain tumor with an average survival of 15 months following diagnosis. Characterized as highly infiltrative with diffuse tumor margins, complete resection and annihilation of tumor cells is impossible following current standard of care therapies. Thus, tumor recurrence is inevitable. Interstitial fluid surrounds all of the cells in the body and has been linked to elevated invasion in glioma, which highlights the importance of this understudied fluid compartment in the brain. The primary objective of this dissertation was to identify specific interstitial fluid transport behaviors associated with elevated invasion surrounding glioma tumors. We first describe our methods to measure interstitial fluid flow in the brain using dynamic contrast enhanced magnetic resonance imaging (DCE-MRI), a clinically used, non-invasive imaging modality. We highlight the versatility of the technique and the possibilities that could arise from widespread adoption into existing perfusion-based imaging protocols. Using this method, we examined transport associated with invasion in a murine GL261 cell line. We found that elevated interstitial fluid velocity magnitudes, decreased diffusion coefficients and regions with accumulating flow were significantly associated with invasion. We tested the validity of our invasive trends by extending our analysis to multiple, clinically-relevant tumor locations in the brain. Interestingly, we found invasion did not follow the same trends across brain regions indicating location-specific structures may drive both interstitial flow and corresponding invasion heterogeneities. Lastly, we aimed to manipulate flow by engaging with the meningeal lymphatics, an established pathway for interstitial fluid drainage. Over-expression of VEGF-C in the tumor microenvironment neither enhanced drainage nor altered invasion in comparison to our control, indicating other tumor-secreted growth factors, such as VEGF-A, may play a larger role in mediating flow and invasion. Taken together, by identifying specific transport factors associated with invasion, we may be better equipped to target and treat infiltrative tumor margins, ultimately extending survival in patients diagnosed with this devastating disease. / Doctor of Philosophy / Glioblastoma is the most common and deadly primary brain tumor with an average survival of 15 months following diagnosis. Characterized as highly infiltrative with diffuse tumor margins, complete resection and annihilation of tumor cells is impossible following current standard of care therapies. Thus, tumor recurrence is inevitable. Interstitial fluid surrounds all of the cells in the body and has been linked to elevated invasion in glioma, which highlights the importance of this understudied fluid compartment in the brain. The primary objective of this dissertation was to identify specific interstitial fluid transport behaviors associated with elevated invasion surrounding glioma tumors. We first describe our methods to measure interstitial fluid flow in the brain using dynamic contrast enhanced magnetic resonance imaging (DCE-MRI), a clinically used, non-invasive imaging modality. We highlight the versatility of the technique and the possibilities that could arise from widespread adoption into existing imaging projects. Using this method, we examined transport associated with cancer cell invasion in a mouse tumor cell line. We found that interstitial fluid speeds were elevated while diffusion was decreased in regions of invasion. Further, regions that had interstitial fluid flow congregation were significantly associated with invasion. We tested the validity of these invasive trends by extending our analysis to multiple, clinically-relevant tumor locations in the brain. Interestingly, we found invasion did not follow the same trends across brain regions, indicating location-specific structures may drive both interstitial flow and invasion differences. Lastly, we aimed to manipulate flow by engaging with the meningeal lymphatics, an established pathway for interstitial fluid drainage. Following administration of a meningeal lymphatic-relevant protein, we saw no changes in flow or invasion in comparison to our untreated control, indicating other tumor-secreted proteins may play a larger role in these responses. Taken together, by identifying specific transport factors associated with invasion, we may be better equipped to target and treat infiltrative tumor margins, ultimately extending survival in patients diagnosed with this devastating disease.

The human brain is a complex organ that is responsible for regulating all the physiological processes in the body, ranging from memory to movement. As humans age, the brain goes through a variety of changes including a reduction in glymphatic waste clearance and increase in glial reactivity. Two neurological conditions that affect individuals over the age of 65 include glioblastoma (GBM) and Alzheimer's disease (AD). Interestingly, patients with GBM do not present with AD and vice versa. Both conditions are characterized by a disruption in interstitial fluid flow (IFF) and an increase in neuroinflammation. Throughout the following dissertation, we examined the role of IFF in AD and GBM progression using a three-sided approach (in vivo, in vitro, and in silico). Increased IFF underlies glioma invasion into the surrounding tumor microenvironment (TME) in GBM. We used a 3D hydrogel model of the GBM TME to examine potential pathways by which astrocytes and microglia contribute to glioma invasion. A reduction in IFF contributes to accumulation of the toxic protein amyloid beta (Aβ) in AD. We sought to create a novel, patient-inspired model of the AD hippocampus for examination of the relationship between IFF and Aβ clearance. Human AD and unaffected control hippocampal brain samples were stained for markers of neurons, astrocytes, microglia and Aβ. The percentage of each cell population in the CA1 region of the hippocampus was calculated. We also analyzed the amount and characteristics of the Aβ aggregates present in this hippocampal region. Pearson correlation analysis was completed to assess the relationships between the various cell populations, Aβ load, and patient descriptors. The cell ratios gleaned from the patient samples were incorporated into a novel, 3D hydrogel model of the AD hippocampus. This model features a hydrogel mixture like the native brain extracellular matrix (ECM) and allows for the application of IFF and Aβ. To our knowledge, we are the first group to create a patient-specific triculture model of the AD hippocampus, which is the main site of Aβ aggregation in the AD brain. We used this model to examine the relationship between IFF-mediated Aβ clearance and glial reactivity. The last aim of this dissertation was to create a computational model for examining Aβ binding within the ECM and the effects of IFF on Aβ clearance. In vitro experiments were conducted to generate 3D renderings of glial cells and to determine relevant parameters for our model. Throughout this work, we discuss the relationship between disruption in IFF and glial reactivity in the context of GBM and AD. / Doctor of Philosophy / The human brain is a complex organ that is responsible for regulating all the physiological processes in the body, ranging from memory to movement. As humans age, the brain goes through a variety of changes including a reduction in brain waste removal and an increase in inflammation. Two neurological conditions that affect individuals over the age of 65 include glioblastoma (GBM) and Alzheimer's disease (AD). Interestingly, patients with GBM do not present with AD and vice versa. Both conditions are characterized by a disruption in brain interstitial fluid flow (IFF) and an increase in neuroinflammation. Throughout the following dissertation, we examined the role of IFF in AD and GBM progression using a three-sided approach including analysis of mouse and human tissues, engineered cell models, and computational methods. Specific interactions between brain cell types and their relationships with glioma invasion were examined using a 3D cell model that mimics the brain. Through the work presented here, we also sought to create a novel cell model of the hippocampus region located in the AD brain. We quantified the various cell types in the hippocampus of AD patient samples and incorporated this information into our hydrogel model. The resulting model features three brain cell types (astrocytes, microglia, and neurons) that are added at patient relevant ratios, a matrix that mimics the native brain scaffold, and allows for the application of IFF. In the AD brain there is a reduction in brain waste removal that leads to accumulation of the toxic protein amyloid beta (Aβ). We were successfully able to incorporate this protein within our model so we could assess the relationship between IFF and Aβ removal from the brain. We further studied this relationship using a new computational model of Aβ accumulation in the brain. Throughout this work, we discuss the connection between disrupted IFF and neuroinflammation in the context of GBM and AD.

<p>Subcutaneous administration has become a common approach for drug delivery of biotherapeutics, such as monoclonal antibodies, which is achieved mainly by absorption through the lymphatic system. This dissertation focuses on the computational modeling of the fluid flow and solute transport in the skin tissue and the quantitative investigation of lymphatic uptake. First, the various mechanisms governing drug transport and lymphatic uptake of biotherapeutics through subcutaneous injection are investigated quantitatively through high-fidelity numerical simulations, including lymphatic drainage, blood perfusion, binding, and metabolism. The tissue is modeled as a homogeneous porous medium using both a single-layered domain and a multi-layered domain, which includes the epidermis, dermis, hypodermis (subcutaneous tissue), and muscle layers. A systematic parameter study is conducted to understand the roles of different properties of the tissue in terms of permeability, porosity, and vascular permeability. The role of binding and metabolism on drug absorption is studied by varying the binding parameters for different macromolecules after coupling the transport equation with a pharmacokinetic equation. The interstitial pressure plays an essential role in regulating the absorption of unbound drug proteins during the injection, while the binding and metabolism of drug molecules reduce the total free drugs. </p> <p>  </p> <p>The lymphatic vessel network is essential to achieve the functions of the lymphatic system. Thus, the drug transport and lymphatic uptake through a three-dimensional hybrid discrete-continuum vessel network in the skin tissue are investigated through high-fidelity numerical simulations. The explicit heterogeneous vessel network is embedded into the continuum model to investigate the role of explicit heterogeneous vessel network in drug transport and absorption. The solute transport across the vessel wall is investigated under various transport conditions. The diffusion of the drug solutes through the explicit vessel wall affects the drug absorption after the injection, while the convection under large interstitial pressure dominates the drug absorption during the injection. The effect of diffusion cannot be captured by the previously developed continuum model. Furthermore, the effects of injection volume and depth on the lymphatic uptake are investigated in a multi-layered domain. The injection volume significantly affects lymphatic uptake through the heterogeneous vessel network, while the injection depth has little influence. At last, the binding and metabolism of drug molecules are studied to bridge the simulation to the experimentally measured drug clearance. </p> <p><br></p> <p>Convective transport of drug solutes in biological tissues is regulated by the interstitial fluid pressure, which plays a crucial role in drug absorption into the lymphatic system through the subcutaneous (SC) injection.  An approximate continuum poroelasticity model is developed to simulate the pressure evolution in the soft porous tissue during an SC injection. This poroelastic model mimics the deformation of the tissue by introducing the time variation of the interstitial fluid pressure. The advantage of this method lies in its computational time efficiency and simplicity, and it can accurately model the relaxation of pressure. The interstitial fluid pressure obtained using the proposed model is validated against both the analytical and the numerical solution of the poroelastic tissue model. The decreasing elasticity elongates the relaxation time of pressure, and the sensitivity of pressure relaxation to elasticity decreases with the hydraulic permeability, while the increasing porosity and permeability due to deformation alleviate the high pressure. </p> <p><br></p> <p>At last, an improved Kedem-Katchalsky model is developed to study solute transport across the lymphatic vessel network, including convection and diffusion in the multi-layered poroelastic tissue with a hybrid discrete-continuum vessel network embedded inside. The effect of different drug solutes with different Stokes radii and different structures of the lymphatic vessel network, such as fractal trees and Voronoi structure, on the lymphatic uptake is investigated. The drug solute with a small size has a larger partition coefficient and diffusivity across the openings of the lymphatic vessel wall, which favors drug absorption. The Voronoi structure is found to be more efficient in lymphatic uptake. </p> <p><br></p> <p>The systematic and quantitative investigation of subcutaneous absorption based on high-fidelity numerical simulations can provide guidance on the optimization of drug delivery systems and is valuable for the translation of bioavailability from the pre-clinical species to humans. We provide a novel approach to studying the diffusion and convection of drug molecules into the lymphatic system by developing the hybrid discrete-continuum vessel network. The study of the solute transport across the discrete lymphatic vessel walls further improves our understanding of lymphatic uptake. The novel and time-efficient computational model for solute transport across the lymphatic vasculature connects the microscopic properties of the lymphatic vessel membrane to macroscopic drug absorption. The comprehensive hybrid vessel network model developed here can be further used to improve our understanding of the diseases caused by the disturbed lymphatic system, such as lymphedema, and provide insights into the treatment of diseases caused by the malfunction of lymphatics.</p>