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.

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.

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.

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.

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.

Voxelized Model of Brain Infusion That Accounts for Small Feature Fissures: Comparison With Magnetic Resonance Tracer Studies.

Convection enhanced delivery (CED) is a promising novel technology to treat neural diseases, as it can transport macromolecular therapeutic agents greater distances through tissue by direct infusion. To minimize off-target delivery, our group has developed 3D computational transport models to predict infusion flow fields and tracer distributions based on magnetic resonance (MR) diffusion tensor imaging data sets. To improve the accuracy of our voxelized models, generalized anisotropy (GA), a scalar measure of a higher order diffusion tensor obtained from high angular resolution diffusion imaging (HARDI) was used to improve tissue segmentation within complex tissue regions of the hippocampus by capturing small feature fissures. Simulations were conducted to reveal the effect of these fissures and cerebrospinal fluid (CSF) boundaries on CED tracer diversion and mistargeting. Sensitivity analysis was also conducted to determine the effect of dorsal and ventral hippocampal infusion sites and tissue transport properties on drug delivery. Predicted CED tissue concentrations from this model are then compared with experimentally measured MR concentration profiles. This allowed for more quantitative comparison between model predictions and MR measurement. Simulations were able to capture infusate diversion into fissures and other CSF spaces which is a major source of CED mistargeting. Such knowledge is important for proper surgical planning.

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.

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.

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.

Refractive index measurement of acute rat brain tissue slices using optical coherence tomography.

An optical coherence tomography (OCT) system employing a microelectromechanical system (MEMS) mirror was used to measure the refractive index (RI) of anatomically different regions in acute brain tissue slices, in which viability was maintained. RI was measured in white-matter and grey-matter regions, including the cerebral cortex, putamen, hippocampus, thalamus and corpus callosum. The RI in the corpus callosum was found to be ~4% higher than the RIs in other regions. Changes in RI with tissue deformation were also measured in the cerebral cortex and corpus callosum under uniform compression (20-80% strain). For 80% strain, measured RIs increased nonlinearly by up to 70% and 90% in the cerebral cortex and corpus callosum respectively. Knowledge of RI in heterogeneous tissues can be used to correct distorted optical images caused by RI variations between different regions. Also deformation-dependent changes in RI can be applied to OCT elastography or to mechanical tests based on optical imaging such as indentation tests.

Role of convection and diffusion on DCE-MRI parameters in low leakiness KHT sarcomas.

The solid tumor is an abnormal environment that is resistant to systemically delivered drugs. Increased plasma leakiness and extracellular matrix density along with poor lymphatic function can result in interstitial flow that attenuates the effectiveness of therapeutics. This study expands upon a previously presented magnetic resonance (MR) imaging-based porous media model by investigating low permeability tumors, where interstitial flow may have increased effect on systemically delivered solutes. The solute transport of the porous media model is compared to that of experiment and the two-compartment model. Small non-necrotic tumors (n=3) were MR-imaged, serially, for 90 min after a bolus injection of Gd-based contrast agent (CA). These data provided for the calculation of experimental CA concentration over 90 min, while only early time points (15 min) were used to create vascular permeability, K(trans), maps for the porous media model. A K(trans) scale factor (range=1.3-2.5) in the porous media model was found to account for the reduction of permeability (measured by two-compartment model) due to interstitial flow. The optimized porous media simulations showed: 1) better dynamic CA behavior agreement with the experimental data than the two-compartment model (>33% reduction of RMS error); 2) similar spatial CA distribution trends across tumor with increased uptake at the tumor boundary.

Influence of needle insertion speed on backflow for convection-enhanced delivery.

Fluid flow back along the outer surface of a needle (backflow) can be a significant problem during the direct infusion of drugs into brain tissues for procedures such as convection-enhanced delivery (CED). This study evaluates the effects of needle insertion speed (0.2 and 1.8 mm/s) as well as needle diameter and flow rate on the extent of backflow and local damage to surrounding tissues. Infusion experiments were conducted on a transparent tissue phantom, 0.6% (w/v) agarose hydrogel, to visualize backflow. Needle insertion experiments were also performed to evaluate local damage at the needle tip and to back out the prestress in the surrounding media for speed conditions where localized damage was not excessive. Prestress values were then used in an analytical model of backflow. At the higher insertion speed (1.8 mm/s), local insertion damage was found to be reduced and backflow was decreased. The compressive prestress at the needle-tissue interface was estimated to be approximately constant (0.812 kPa), and backflow distances were similar regardless of needle gauge (22, 26, and 32 gauge). The analytical model underestimated backflow distances at low infusion flow rates and overestimated backflow at higher flow rates. At the lower insertion speed (0.2 mm/s), significant backflow was measured. This corresponded to an observed accumulation of material at the needle tip which produced a gap between the needle and the surrounding media. Local tissue damage was also evaluated in excised rat brain tissues, and insertion tests show similar rate-dependent accumulation of tissue at the needle tip at the lower insertion speed. These results indicate that local tissue damage and backflow may be avoided by using an appropriate insertion speed.

Evaluation of a voxelized model based on DCE-MRI for tracer transport in tumor.

Recent advances in the treatment of cancer involving therapeutic agents have shown promising results. However, treatment efficacy can be limited due to inadequate and uneven uptake in solid tumors, thereby making the prediction of drug transport important for developing effective therapeutic strategies. In this study, a patient-specific computational porous media model (voxelized model) was developed for predicting the interstitial flow field and distribution of a systemically delivered magnetic resonance (MR) visible tracer in a tumor. The benefits of a voxel approach include less labor and less computational time (approximately an order of magnitude reduction compared to the traditional computational fluid dynamics (CFD) approach developed earlier by our group). The model results were compared with that obtained from a previous approach based on unstructured meshes along with MR-measured tracer concentration data within tumors, using statistical analysis and qualitative representations. The statistical analysis indicated the similarity between the structured and unstructured models' results with a low root mean square error (RMS) and a high correlation coefficient. The voxelized model captured features of the flow field and tracer distribution such as high interstitial fluid pressure inside the tumor and the heterogeneous distribution of the tracer. Predictions of tracer distribution by the voxelized approach also resulted in low RMS error when compared with MR-measured data over a 1 h time course. The similarity in the voxelized model results with experiment and the nonvoxelized model predictions were maintained across three different tumors. Overall, the voxelized model serves as a reliable and swift alternative to approaches using unstructured meshes in predicting extracellular transport within tumors.

Polymer-coated cannulas for the reduction of backflow during intraparenchymal infusions.

Infusate backflow or leak-back along the cannula track can occur during intraparenchymal infusions resulting in non-specific targeting of therapeutic agents. The occurrence of backflow depends on several variables including cannula radius, infusate flow rate, and tip location. In this study, polymer coatings that swell in situ were developed and tested with in vitro hydrogel experiments for backflow reduction. Coatings were applied to the external cannula surface in a dual layer arrangement with a poly(vinyl alcohol) outer layer atop an inner poly(ethylene oxide) and alginate layer. Once these coated cannulas were inserted and allotted an 8-10 min waiting period for hydration, backflow during infusions of 4.0 µl of a macromolecular tracer (Evans Blue labeled albumin) was reduced significantly under flow rates of 0.3-0.6 µl/min, allowing for more effective distribution within targeted regions. Polymer coating thicknesses before and after hydrations were 0.035 and 0.370 mm, respectively. Also, backflow data was fit to a model to estimate the effective local compressive stress caused by the hydrated polymers. After withdrawal of the cannula from the insertion site, the hydrated polymer coatings remained within the cavity left in the hydrogel tissue phantom and formed a seal at the infusion site that prevented further backflow during needle withdrawal. Ex vivo infusions in excised porcine brain tissues also showed significant backflow reduction while also demonstrating the ability to leave a polymer seal in the tissue cavity after cannula removal. Thus, application of these polymers as needle or cannula coatings offers a potentially simple method to improve targeting for local drug delivery.

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.

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.

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.

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.

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.