RESUMO
Mechanical stress and fluid flow influence glioma cell phenotype in vitro, but measuring these quantities in vivo continues to be challenging. The purpose of this study was to predict these quantities in vivo, thus providing insight into glioma physiology and potential mechanical biomarkers that may improve glioma detection, diagnosis, and treatment. Image-based finite element models of human U251N orthotopic glioma in athymic rats were developed to predict structural stress and interstitial flow in and around each animal's tumor. In addition to accounting for structural stress caused by tumor growth, our approach has the advantage of capturing fluid pressure-induced structural stress, which was informed by in vivo interstitial fluid pressure (IFP) measurements. Because gliomas and the brain are soft, elevated IFP contributed substantially to tumor structural stress, even inverting this stress from compressive to tensile in the most compliant cases. The combination of tumor growth and elevated IFP resulted in a concentration of structural stress near the tumor boundary where it has the greatest potential to influence cell proliferation and invasion. MRI-derived anatomical geometries and tissue property distributions resulted in heterogeneous interstitial fluid flow with local maxima near cerebrospinal fluid spaces, which may promote tumor invasion and hinder drug delivery. In addition, predicted structural stress and interstitial flow varied markedly between irradiated and radiation-naïve animals. Our modeling suggests that relative to tumors in stiffer tissues, gliomas experience unusual mechanical conditions with potentially important biological (e.g., proliferation and invasion) and clinical consequences (e.g., drug delivery and treatment monitoring).
RESUMO
Adeno-associated virus (AAV) is a clinically useful gene delivery vehicle for treating neurological diseases. To deliver AAV to focal targets, direct infusion into brain tissue by convection-enhanced delivery (CED) is often needed due to AAV's limited penetration across the blood-brain-barrier and its low diffusivity in tissue. In this study, computational models that predict the spatial distribution of AAV in brain tissue during CED were developed to guide future placement of infusion catheters in recurrent brain tumors following primary tumor resection. The brain was modeled as a porous medium, and material property fields that account for magnetic resonance imaging (MRI)-derived anatomical regions were interpolated and directly assigned to an unstructured finite element mesh. By eliminating the need to mesh complex surfaces between fluid regions and tissue, mesh preparation was expedited, increasing the model's clinical feasibility. The infusion model predicted preferential fluid diversion into open fluid regions such as the ventricles and subarachnoid space (SAS). Additionally, a sensitivity analysis of AAV delivery demonstrated that improved AAV distribution in the tumor was achieved at higher tumor hydraulic conductivity or lower tumor porosity. Depending on the tumor infusion site, the AAV distribution covered 3.67-70.25% of the tumor volume (using a 10% AAV concentration threshold), demonstrating the model's potential to inform the selection of infusion sites for maximal tumor coverage.
Assuntos
Neoplasias Encefálicas , Dependovirus , Análise de Elementos Finitos , Imageamento por Ressonância Magnética , Neoplasias Encefálicas/diagnóstico por imagem , Neoplasias Encefálicas/terapia , Imageamento por Ressonância Magnética/métodos , Humanos , Modelos Biológicos , Porosidade , Recidiva Local de Neoplasia/diagnóstico por imagemRESUMO
A custom segmentation workflow was applied to ex vivo high-field MR images of rat brains acquired following in vivo intraventricular contrast agent infusion to generate maps of the perivascular spaces (PVS). The resulting perivascular network segmentations enabled analysis of perivascular connections to the ventricles, parenchymal solute clearance, and dispersive solute transport within PVS. Numerous perivascular connections between the brain surface and the ventricles suggest the ventricles integrate into a PVS-mediated clearance system and raise the possibility of cerebrospinal fluid (CSF) return from the subarachnoid space to the ventricles via PVS. Assuming rapid solute exchange between the PVS and CSF spaces primarily by advection, the extensive perivascular network decreased the mean clearance distance from parenchyma to the nearest CSF compartment resulting in an over 21-fold reduction in the estimated diffusive clearance time scale, irrespective of solute diffusivity. This corresponds to an estimated diffusive clearance time scale under 10 min for amyloid-beta which suggests that the widespread distribution of PVS may render diffusion an effective parenchymal clearance mechanism. Additional analysis of oscillatory solute dispersion within PVS indicates that advection rather than dispersion is likely the primary transport mechanism for dissolved compounds greater than 66 kDa in the long (> 2 mm) perivascular segments identified here, although dispersion may be significant for smaller compounds in shorter perivascular segments.
Assuntos
Encéfalo , Imageamento por Ressonância Magnética , Ratos , Animais , Encéfalo/diagnóstico por imagem , Encéfalo/irrigação sanguínea , Imageamento por Ressonância Magnética/métodos , Espaço Subaracnóideo , Meios de Contraste , DifusãoRESUMO
Shear stress is recognized as a regulator of angiogenesis. However, the shear stress experienced by the endothelial cells of capillary sprouts remains unknown. The objective of this study was to estimate shear stress due to local interstitial flow along endothelial tip cells at the end of the capillary sprout lumen. Computational fluid dynamics were used to model flow within a blind-ended vessel, transendothelial flow across the vessel wall, and flow within the surrounding perivascular/interstitial space. Shear stress along the wall of the tip cells was calculated while varying sprout length, perivascular space channel width, and vessel wall hydraulic conductivity. Increasing sprout length, increasing wall hydraulic conductivity, and decreasing perivascular space width increased shear stress magnitude. Wall shear stress magnitude within the lumen ranged from 0.015 to 0.55 dyne/cm2 at the sprout entrance and linearly decreased to near zero at the base of the tip cells. Tip cell wall shear stress magnitude due to interstitial flow ranged from 0.009 to 4.65 dyne/cm2. In 3 out of 8 cases, shear stress magnitude was above 1 dyne/cm2 and considered physiologically relevant. The results provide a framework for discussing the role of local mechanical cues in regulating endothelial cell dynamics involved in angiogenesis. Mainly, interstitial flows may generate physiologically relevant shear stresses on tip cells in certain scenarios. This source of tip cell shear stress has not been previously considered or modeled.
Assuntos
Capilares , Células Endoteliais , Capilares/fisiologia , Hidrodinâmica , Estresse Mecânico , VeiasRESUMO
A biphasic computational model of a growing, vascularized glioma within brain tissue was developed to account for unique features of gliomas, including soft surrounding brain tissue, their low stiffness relative to brain tissue, and a lack of draining lymphatics. This model is the first to couple nonlinear tissue deformation with porosity and tissue hydraulic conductivity to study the mechanical interaction of leaky vasculature and solid growth in an embedded glioma. The present model showed that leaky vasculature and elevated interstitial fluid pressure produce tensile stress within the tumor in opposition to the compressive stress produced by tumor growth. This tensile effect was more pronounced in softer tissue and resulted in a compressive stress concentration at the tumor rim that increased when tumor was softer than host. Aside from generating solid stress, fluid pressure-driven tissue deformation decreased the effective stiffness of the tumor while growth increased it, potentially leading to elevated stiffness in the tumor rim. A novel prediction of reduced porosity at the tumor rim was corroborated by direct comparison with estimates from our in vivo imaging studies. Antiangiogenic and radiation therapy were simulated by varying vascular leakiness and tissue hydraulic conductivity. These led to greater solid compression and interstitial pressure in the tumor, respectively, the former of which may promote tumor infiltration of the host. Our findings suggest that vascular leakiness has an important influence on in vivo solid stress, stiffness, and porosity fields in gliomas given their unique mechanical microenvironment.
Assuntos
Neoplasias Encefálicas/fisiopatologia , Líquido Extracelular/fisiologia , Glioma/fisiopatologia , Microambiente Tumoral , Animais , Encéfalo , Força Compressiva , Simulação por Computador , Humanos , Linfa/fisiologia , Modelos Biológicos , Modelos Teóricos , Porosidade , Pressão , Estresse Mecânico , Resistência à TraçãoRESUMO
The effect of a human vascular endothelial growth factor antibody on the vasculature of human tumor grown in rat brain was studied. Using dynamic contrast-enhanced magnetic resonance imaging, the effects of intravenous bevacizumab (Avastin; 10 mg/kg) were examined before and at postadministration times of 1, 2, 4, 8, 12 and 24 h (N = 26; 4-5 per time point) in a rat model of orthotopic, U251 glioblastoma (GBM). The commonly estimated vascular parameters for an MR contrast agent were: (i) plasma distribution volume (vp ), (ii) forward volumetric transfer constant (Ktrans ) and (iii) reverse transfer constant (kep ). In addition, extracellular distribution volume (VD ) was estimated in the tumor (VD-tumor ), tumor edge (VD-edge ) and the mostly normal tumor periphery (VD-peri ), along with tumor blood flow (TBF), peri-tumoral hydraulic conductivity (K) and interstitial flow (Flux) and tumor interstitial fluid pressure (TIFP). Studied as % changes from baseline, the 2-h post-treatment time point began showing significant decreases in vp , VD-tumor, VD-edge and VD-peri , as well as K, with these changes persisting at 4 and 8 h in vp , K, VD-tumor, -edge and -peri (t-tests; p < 0.05-0.01). Decreases in Ktrans were observed at the 2- and 4-h time points (p < 0.05), while interstitial volume fraction (ve ; = Ktrans /kep ) showed a significant decrease only at the 2-h time point (p < 0.05). Sustained decreases in Flux were observed from 2 to 24 h (p < 0.01) while TBF and TIFP showed delayed responses, increases in the former at 12 and 24 h and a decrease in the latter only at 12 h. These imaging biomarkers of tumor vascular kinetics describe the short-term temporal changes in physical spaces and fluid flows in a model of GBM after Avastin administration.