A Spatial Interpolation Approach to Assign Magnetic Resonance Imaging-Derived Material Properties for Finite Element Models of Adeno-Associated Virus Infusion Into a Recurrent Brain Tumor.

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.

Sequential Dual Alignments Introduce Synergistic Effect on Hexagonal Boron Nitride Platelets for Superior Thermal Performance.

Planarly aligning 2D platelets is challenging due to their additional orientational freedom compared to 1D materials. This study reports a sequential dual-alignment approach, employing an extrusion-printing-induced shear force and rotating-magnetic-field-induced force couple for platelet planarly alignment in a yield-stress support bath. It is hypothesized that the partial alignment induced by a directional shear force facilitates subsequent axial rotation of the platelets for planar alignment under an external force couple, resulting in a synergistic alignment effect. This sequential dual-alignment approach achieves better planar alignment of 2D modified hexagonal boron nitride (mhBN). Specifically, the thermal conductivity of the 40 wt% mhBN/epoxy composite is significantly higher (692%) than that of unaligned composites, surpassing the cumulative effect of individual methods (only 133%) with a 5 times more synergistic effect. For 30, 40, and 50 wt% mhBN composites, the thermal conductivity values (5.9, 9.5, and 13.8 W m-1 K-1 ) show considerable improvement compared to the previously reported highest values (5.3, 6.6, and 8.6 W m-1 K-1 ). Additionally, a 3D mhBN/epoxy heat sink is printed and evaluated to demonstrate the feasibility of device fabrication. The approach enables the planar alignment of electrically or thermally conducting 2D fillers during 3D fabrication.

Estimation of shear stress heterogeneity along capillary segments in angiogenic rat mesenteric microvascular networks.

OBJECTIVE: Fluid shear stress is thought to be a regulator of endothelial cell behavior during angiogenesis. The link, however, requires an understanding of stress values at the capillary level in angiogenic microvascular networks. Critical questions remain. What are the stresses? Do capillaries experience similar stress magnitudes? Can variations explain vessel-specific behavior? The objective of this study was to estimate segment-specific shear stresses in angiogenic networks. METHODS: Images of angiogenic networks characterized by increased vascular density were obtained from rat mesenteric tissues stimulated by compound 48/80-induced mast cell degranulation. Vessels were identified by perfusion of a 40 kDa fixable dextran prior to harvesting and immunolabeling for PECAM. Using a network flow-based segment model with physiologically relevant parameters, stresses were computed per vessel for regions across multiple networks. RESULTS: Stresses ranged from 0.003 to 2328.1 dyne/cm2 and varied dramatically at the capillary level. For all regions, the maximum segmental shear stresses were for capillary segments. Stresses along proximal capillaries branching from arteriole inlets were increased compared to stresses along capillaries in more distal regions. CONCLUSIONS: The results highlight the variability of shear stresses along angiogenic capillaries and motivate new discussions on how endothelial cells may respond in vivo to segment-specific microenvironment during angiogenesis.

Perivascular network segmentations derived from high-field MRI and their implications for perivascular and parenchymal mass transport in the rat brain.

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.

Label-Free Quantification of Microscopic Alignment in Engineered Tissue Scaffolds by Polarized Raman Spectroscopy.

Monitoring of extracellular matrix (ECM) microstructure is essential in studying structure-associated cellular processes, improving cellular function, and for ensuring sufficient mechanical integrity in engineered tissues. This paper describes a novel method to study the microscale alignment of the matrix in engineered tissue scaffolds (ETS) that are usually composed of a variety of biomacromolecules derived by cells. First, a trained loading function was derived from Raman spectra of highly aligned native tissue via principal component analysis (PCA), where prominent changes associated with specific Raman bands (e.g., 1444, 1465, 1605, 1627-1660, and 1665-1689 cm-1) were detected with respect to the polarization angle. These changes were mainly caused by the aligned matrix of many compounds within the tissue relative to the laser polarization, including proteins, lipids, and carbohydrates. Hence this trained function was applied to quantify the alignment within ETS of various matrix components derived by cells. Furthermore, a simple metric called Amplitude Alignment Metric (AAM) was derived to correlate the orientation dependence of polarized Raman spectra of ETS to the degree of matrix alignment. It was found that the AAM was significantly higher in anisotropic ETS than isotropic ones. The PRS method revealed a lower p-value for distinguishing the alignment between these two types of ETS as compared to the microscopic method for detecting fluorescent-labeled protein matrices at a similar microscopic scale. These results indicate that the anisotropy of a complex matrix in engineered tissue can be assessed at the microscopic scale using a PRS-based simple metric, which is superior to the traditional microscopic method. This PRS-based method can serve as a complementary tool for the design and assessment of engineered tissues that mimic the native matrix organizational microstructures.

Experimental and Computational Models of Transport of Galectin-3 Through Glycosylated Matrix.

Altered extracellular matrix (ECM) production is a hallmark of many fibroproliferative diseases, including certain cancers. The high incidence of glycan-rich components within altered ECM makes the use of glycan-binding proteins such as Galectin-3 (G3) a promising therapeutic strategy. The complexity of ECM as a rich 3D network of proteins with varied glycosylation states makes it challenging to determine the retention of glycan-binding proteins in altered ECM environments. Computational models capable of predicting the transport of glycan-binding proteins in altered ECM can benefit the design and testing of such proteins and associated novel therapeutic strategies. However, such computational models require many kinetic parameters that cannot be estimated from traditional 2D pharmacokinetic assays. To validate transport properties of G3 in 3D ECM constructs, we developed a species transport model that includes diffusion and matrix-binding components to predict retention of G3 fusion proteins in glycan-rich ECM. By iteratively comparing our computational model to experimental results, we are able to determine a reasonable range of parameters for a robust computational model of G3 transport. We anticipate this overall approach to building a data-driven model is translatable to other ECM-targeting therapeutic strategies.

Estimation of shear stress values along endothelial tip cells past the lumen of capillary sprouts.

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.

Label-free quantification of soft tissue alignment by polarized Raman spectroscopy.

The organization of proteins is an important determinant of functionality in soft tissues. However, such organization is difficult to monitor over time in soft tissue with complex compositions. Here, we establish a method to determine the alignment of proteins in soft tissues of varying composition by polarized Raman spectroscopy (PRS). Unlike most conventional microscopy methods, PRS leverages non-destructive, label-free sample preparation. PRS data from highly aligned muscle layers were utilized to derive a weighting function for aligned proteins via principal component analysis (PCA). This trained weighting function was used as a master loading function to calculate a principal component score (PC1 Score) as a function of polarized angle for tendon, dermis, hypodermis, and fabricated collagen gels. Since the PC1 Score calculated at arbitrary angles was insufficient to determine level of alignment, we developed an Amplitude Alignment Metric by fitting a sine function to PC1 Score with respect to polarized angle. We found that our PRS-based Amplitude Alignment Metric can be used as an indicator of level of protein alignment in soft tissues in a non-destructive manner with label-free preparation and has similar discriminatory capacity among isotropic and anisotropic samples compared to microscopy-based image processing method. This PRS method does not require a priori knowledge of sample orientation nor composition and appears insensitive to changes in protein composition among different tissues. The Amplitude Alignment Metric introduced here could enable convenient and adaptable evaluation of protein alignment in soft tissues of varying protein and cell composition. STATEMENT OF SIGNIFICANCE: Polarized Raman spectroscopy (PRS) has been used to characterize the of organization of soft tissues. However, most of the reported applications of PRS have been on collagen-rich tissues and reliant on intensities of collagen-related vibrations. This work describes a PRS method via a multivariate analysis to characterize alignment in soft tissues composed of varying proteins. Of note, the highly aligned muscle layer of mouse skin was used to train a master function then applied to other soft tissue samples, and the degree of anisotropy in the PRS response was evaluated to obtain the level of alignment in tissues. We have demonstrated that this method supports convenient and adaptable evaluation of protein alignment in soft tissues of varying protein and cell composition.

A computational model of glioma reveals opposing, stiffness-sensitive effects of leaky vasculature and tumor growth on tissue mechanical stress and porosity.

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.

Imaging acute effects of bevacizumab on tumor vascular kinetics in a preclinical orthotopic model of U251 glioma.

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.

Raman Spectroscopy Methods to Characterize the Mechanical Response of Soft Biomaterials.

Raman spectroscopy has been used extensively to characterize the influence of mechanical deformation on microstructure changes in biomaterials. While traditional piezo-spectroscopy has been successful in assessing internal stresses of hard biomaterials by tracking prominent peak shifts, peak shifts due to applied loads are near or below the resolution limit of the spectrometer for soft biomaterials with moduli in the kilo- to mega-Pascal range. In this Review, in addition to peak shifts, other spectral features (e.g., polarized intensity and intensity ratio) that provide quantitative assessments of microstructural orientation and secondary structure in soft biomaterials and their strain dependence are discussed. We provide specific examples for each method and classify sensitive Raman characteristic bands common across natural (e.g., soft tissue) and synthetic (e.g., polymeric scaffolds) soft biomaterials upon mechanical deformation. This Review can provide guidance for researchers aiming to analyze micromechanics of soft tissues and engineered tissue constructs by Raman spectroscopy.

An MRI-based switched gradient impulse response characterization method with uniform eigenmode excitation.

Switching gradients generate eddy currents and mechanical vibrations of the gradient assembly causing errors in the gradient time integrals. This results in image distortions in k-space and inaccuracies in q-space imaging. The purpose of this work is to develop an MRI based unbiased measurement of the switched gradient impulse response function (sGIRF). A new gradient pattern, called the Tukey windowed Shifted Sine-Integral (Tw-SSI) pulse, is introduced to excite the gradient eigenmodes uniformly over a user-defined bandwidth. A 3D MRI-based method with Hadamard encoding was developed to map the spatiotemporal magnetic field generated after the excitation pulse to obtain the sGIRF for all the three gradient axes simultaneously. Compared to an energy-equivalent traditional trapezoidal pulse, the Tw-SSI pulse is able to excite the weak bandlimited cross-terms of the sGIRF by uniformly distributing the energy across eigenmodes. The developed field mapping method is sensitive enough to capture both the direct and cross-terms in the sGIRF. The various mechanical resonant modes of the gradient coils are also revealed, which were found to last longer than eddy currents in the shielded gradient coil studied. Tunable Tw-SSI pulse offers the flexibility to perform unbiased sGIRF measurements over a bandwidth of interest. Rapid MRI field mapping can be easily implemented in any MRI system. The method may be used to perform gradient pre-emphasis, to evaluate new gradient coil designs, and to characterize higher order shims.

Cochlear basal turn patency in unrecognized perilymph gushers.

OBJECTIVES: Perilymph gusher (PLG), an uncommon complication of otologic surgery, has been attributed to communication between the cochlea and the internal auditory canal (IAC). Subtle osseous defects may be missed on routine review of computed tomography (CT). This study aimed to quantify cochlear basal turn patency not seen on axial CT in patients with PLG and compare those against patients without intraoperative PLG. METHODS: Ears that underwent cochlear implantation or stapedotomy with preoperative helical CT that was interpreted as "normal" at a tertiary referral center. An otologist and a radiologist independently and in a blinded fashion measured the dimensions of cochlear basal turn patency on CT images in oblique plane and parasagittal planes along the interface of the cochlea and IAC fundus. RESULTS: Sixty-one ears were reviewed, including 3 with surgically confirmed PLGs and 12 with apparent dehiscence without a PLG. Mean defect width with PLG was 0.83 mm (range 0.75-0.9 mm) and without PLG was 0.43 mm (range 0.3-0.65 mm, p = 0.011). A greater proportion of PLGs occurred in ears with defects (3 of 15) than in ears without (0 of 46, p = 0.013). Using a cutoff of 0.75 mm, a greater proportion of PLGs occurred with defect width >0.75 mm (3 of 3) than in defects <0.75 mm (0 of 12, p = 0.022). CONCLUSIONS: CT dehiscence between the IAC and cochlear basal turn, particularly with a width > 0.75 mm, should be considered a risk for PLG with stapedotomy or cochlear implantation.

MRI of Whole Rat Brain Perivascular Network Reveals Role for Ventricles in Brain Waste Clearance.

Investigating the mechanisms by which metabolic wastes are cleared from nervous tissue is important for understanding natural function and the pathophysiology of several neurological disorders including Alzheimer's disease. Recent evidence suggests clearance may be the function of annular spaces around cerebral blood vessels, called perivascular spaces (PVS), through which cerebrospinal fluid (CSF) is transported from the subarachnoid space into brain parenchyma to exchange with interstitial fluid (also known as the glymphatic system). In this work, an MRI-based methodology was developed to reconstruct the PVS network in whole rat brain to better elucidate both PVS uptake and clearance pathways. MR visible tracer (Gd-albumin) was infused in vivo into the CSF-filled lateral ventricle followed by ex vivo high-resolution MR imaging at 17.6 T with an image voxel volume two orders of magnitude smaller than previously reported. Imaged tracer distribution patterns were reconstructed to obtain a more complete brain PVS network. Several PVS connections were repeatedly highlighted across different animals, and new PVS connections between ventricles and different parts of the brain parenchyma were revealed suggesting a possible role for the ventricles as a source or sink for solutes in the brain. In the future, this methodology may be applied to understand changes in the PVS network with disease.

Longitudinal evaluation of tumor microenvironment in rat focal brainstem glioma using diffusion and perfusion MRI.

BACKGROUND: Brainstem gliomas are aggressive and difficult to treat. Growth of these tumors may be characterized with MRI methods. PURPOSE: To visualize longitudinal changes in tumor volume, vascular leakiness, and tissue microstructure in an animal model of brainstem glioma. STUDY TYPE: Prospective animal model. ANIMAL MODEL: Male Sprague-Dawley rats (n = 9) were imaged with 9L gliosarcoma cells infused into the pontine reticular formation of the brainstem. The MRI tumor microenvironment was studied at 3 and 10 days postimplantation of tumor cells. FIELD STRENGTH/SEQUENCE: Diffusion tensor imaging (DTI) and dynamic contrast-enhanced (DCE)-MRI were performed at 4.7T using spin-echo multislice echo planar imaging and gradient echo multislice imaging, respectively. ASSESSMENT: Tumor leakiness was assessed by the forward volumetric transfer constant, Ktrans , estimated from DCE-MRI data. Tumor structure was evaluated with fractional anisotropy (FA) obtained from DTI. Tumor volumes, delineated by a T1 map, T2 -weighted image, FA, and DCE signal enhancement were compared. STATISTICAL TESTS: Changes in the assessed parameters within and across the groups (ie, rats 3 and 10 days post tumor cell implantation) were evaluated with Wilcoxon rank-sum tests. RESULTS: Day 3 tumors were visible mainly on contrast-enhanced images, while day 10 tumors were visible in both contrast-enhanced and diffusion-weighted images. Mean Ktrans at day 10 was 41% lower than at day 3 (P = 0.23). In day 10 tumors, FA was regionally lower in the tumor compared to normal tissue (P = 0.0004), and tumor volume, segmented based on FA map, was significantly smaller (P ≤ 0.05) than that obtained from other contrasts. DATA CONCLUSION: Contrast-enhanced MRI was found to be more sensitive in detecting early-stage tumor boundaries than other contrasts. Areas of the tumor outlined by DCE-MRI and DTI were significantly different. Over the observed period of tumor growth, average vessel leakiness decreased with tumor progression. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2019;49:1322-1332.

Biphasic analysis of rat brain slices under creep indentation shows nonlinear tension-compression behavior.

Biphasic theory can provide a mechanistic description of deformation and transport phenomena in soft tissues, and has been used to model surgery and drug delivery in the brain for decades. Knowledge of corresponding mechanical properties of the brain is needed to accurately predict tissue deformation and flow transport in these applications. Previously in our group, creep indentation tests were conducted for multiple anatomical regions in acute rat brain tissue slices. In the current study, a biphasic finite element model of creep indentation was developed with which to compare these data. Considering the soft tissue structure of brain, the solid matrix was assumed to be composed of a neo-Hookean ground matrix reinforced by continuously distributed fibers that exhibits tension-compression nonlinearity during deformation. By fixing Poisson's ratio of the ground matrix, Young's modulus, fiber modulus and hydraulic permeability were estimated. Hydraulic permeability was found to be nearly independent of the properties of the solid matrix. Estimated modulus (40 Pa to 1.1 kPa for the ground matrix, 3.2-18.2 kPa for fibers) and hydraulic permeability (1.2-5.5×10-13m4/N s) fell within an acceptable range compared with those in previous studies. Instantaneous indentation depth was dominated by tension provided by fibers, while the tissue response at equilibrium was sensitive to Poisson's ratio. Results of sensitivity analysis also point to the necessity of considering tension-compression nonlinearity in the solid phase when the biphasic material undergoes large creep deformation.

Phase contrast MRI of creeping flows using stimulated echo.

Creeping flows govern many important physiological phenomena such as elevated interstitial fluid flows in tumors, glymphatic flows in the brain, among other applications. However, few methods exist to measure such slow flows non-invasively in optically opaque biological tissues in vivo. Phase-contrast MRI is a velocimetry technique routinely used in the clinic to measure fast flows in biological tissues, such as blood and cerebrospinal fluid (CSF), in the order of cm/s. Use of this technique to encode slower flows is hampered by diffusion weighting and phase error introduced by gradient hardware imperfections. In this study, a new PC-MRI technique is developed using stimulated echo preparation to overcome these challenges. Flows as slow as 1 µm/s are measured and validated using controlled water flow through a pipe at 4.7 T. The error in measured flow rate obtained by integrating the measured velocity over the cross-sectional area of the pipe is less than 10%. The developed method was also able to capture slow natural convection flows appearing in liquids placed inside a horizontal bore magnet. Monitoring the 4D velocity vector field revealed that the natural convection flows decay exponentially with time. This method could be applied in future to study creeping flows, e.g. in tissue.

Pulsatile flow drivers in brain parenchyma and perivascular spaces: a resistance network model study.

BACKGROUND: In animal models, dissolved compounds in the subarachnoid space and parenchyma have been found to preferentially transport through the cortex perivascular spaces (PVS) but the transport phenomena involved are unclear. METHODS: In this study two hydraulic network models were used to predict fluid motion produced by blood vessel pulsations and estimate the contribution made to solute transport in PVS and parenchyma. The effect of varying pulse amplitude and timing, PVS dimensions, and tissue hydraulic conductivity on fluid motion was investigated. RESULTS: Periodic vessel pulses resulted in oscillatory fluid motion in PVS and parenchyma but no net flow over time. For baseline parameters, PVS and parenchyma peak fluid velocity was on the order of 10 µm/s and 1 nm/s, with corresponding Peclet numbers below 103 and 10-1 respectively. Peak fluid velocity in the PVS and parenchyma tended to increase with increasing pulse amplitude and vessel size, and exhibited asymptotic relationships with hydraulic conductivity. CONCLUSIONS: Solute transport in parenchyma was predicted to be diffusion dominated, with a negligible contribution from convection. In the PVS, dispersion due to oscillating flow likely plays a significant role in PVS rapid transport observed in previous in vivo experiments. This dispersive effect could be more significant than convective solute transport from net flow that may exist in PVS and should be studied further.

Controlled single bubble cavitation collapse results in jet-induced injury in brain tissue.

Multiscale damage due to cavitation is considered as a potential mechanism of traumatic brain injury (TBI) associated with explosion. In this study, we employed a TBI relevant hippocampal ex vivo slice model to induce bubble cavitation. Placement of single reproducible seed bubbles allowed control of size, number, and tissue location to visualize and measure deformation parameters. Maximum strain value was measured at 45 µs after bubble collapse, presented with a distinct contour and coincided temporally and spatially with the liquid jet. Composite injury maps combined this maximum strain value with maximum measured bubble size and location along with histological injury patterns. This facilitated the correlation of bubble location and subsequent jet direction to the corresponding regions of high strain which overlapped with regions of observed injury. A dynamic threshold strain range for tearing of cerebral cortex was estimated to be between 0.5 and 0.6. For a seed bubble placed underneath the hippocampus, cavitation induced damage was observed in hippocampus (local), proximal cerebral cortex (marginal) and the midbrain/forebrain (remote) upon histological evaluation. Within this test model, zone of cavitation injury was greater than the maximum radius of the bubble. Separation of apposed structures, tissue tearing, and disruption of cellular layers defined early injury patterns that were not detected in the blast-exposed half of the brain slice. Ultrastructural pathology of the neurons exposed to cavitation was characterized by disintegration of plasma membrane along with loss of cellular content. The developed test system provided a controlled experimental platform to study cavitation induced high strain deformations on brain tissue slice. The goal of the future studies will be to lower underpressure magnitude and cavitation bubble size for more sensitive evaluation of injury.

Simulated blast overpressure induces specific astrocyte injury in an ex vivo brain slice model.

Exposure to explosive blasts can produce functional debilitation in the absence of brain pathology detectable at the scale of current diagnostic imaging. Transient (ms) overpressure components of the primary blast wave are considered to be potentially damaging to the brain. Astrocytes participate in neuronal metabolic maintenance, blood-brain barrier, regulation of homeostatic environment, and tissue remodeling. Damage to astrocytes via direct physical forces has the potential to disrupt local and global functioning of neuronal tissue. Using an ex vivo brain slice model, we tested the hypothesis that viable astrocytes within the slice could be injured simply by transit of a single blast wave consisting of overpressure alone. A polymer split Hopkinson pressure bar (PSHPB) system was adapted to impart a single positive pressure transient with a comparable magnitude to those that might be present inside the head. A custom built test chamber housing the brain tissue slice incorporated revised design elements to reduce fluid space and promote transit of a uniform planar waveform. Confocal microscopy, stereology, and morphometry of glial fibrillary acidic protein (GFAP) immunoreactivity revealed that two distinct astrocyte injury profiles were identified across a 4 hr post-test survival interval: (a) presumed conventional astrogliosis characterized by enhanced GFAP immunofluorescence intensity without significant change in tissue area fraction and (b) a process comparable to clasmatodendrosis, an autophagic degradation of distal processes that has not been previously associated with blast induced neurotrauma. Analysis of astrocyte branching revealed early, sustained, and progressive differences distinct from the effects of slice incubation absent overpressure testing. Astrocyte vulnerability to overpressure transients indicates a potential for significant involvement in brain blast pathology and emergent dysfunction. The testing platform can isolate overpressure injury phenomena to provide novel insight on physical and biological mechanisms.