Accelerated, physics-inspired inference of skeletal muscle microstructure from diffusion-weighted MRI.

Muscle health is a critical component of overall health and quality of life. However, current measures of skeletal muscle health take limited account of microstructural variations within muscle, which play a crucial role in mediating muscle function. To address this, we present a physics-inspired, machine learning-based framework for the non-invasive estimation of microstructural organization in skeletal muscle from diffusion-weighted MRI (dMRI) in an uncertainty-aware manner. To reduce the computational expense associated with direct numerical simulations of dMRI physics, a polynomial meta-model is developed that accurately represents the input/output relationships of a high-fidelity numerical model. This meta-model is used to develop a Gaussian process (GP) model that provides voxel-wise estimates and confidence intervals of microstructure organization in skeletal muscle. Given noise-free data, the GP model accurately estimates microstructural parameters. In the presence of noise, the diameter, intracellular diffusion coefficient, and membrane permeability are accurately estimated with narrow confidence intervals, while volume fraction and extracellular diffusion coefficient are poorly estimated and exhibit wide confidence intervals. A reduced-acquisition GP model, consisting of one-third the diffusion-encoding measurements, is shown to predict parameters with similar accuracy to the original model. The fiber diameter and volume fraction estimated by the reduced GP model is validated via histology, with both parameters accurately estimated, demonstrating the capability of the proposed framework as a promising non-invasive tool for assessing skeletal muscle health and function.

Highly stretchable and customizable microneedle electrode arrays for intramuscular electromyography.

Stretchable three-dimensional (3D) penetrating microelectrode arrays have potential utility in various fields, including neuroscience, tissue engineering, and wearable bioelectronics. These 3D microelectrode arrays can penetrate and conform to dynamically deforming tissues, thereby facilitating targeted sensing and stimulation of interior regions in a minimally invasive manner. However, fabricating custom stretchable 3D microelectrode arrays presents material integration and patterning challenges. In this study, we present the design, fabrication, and applications of stretchable microneedle electrode arrays (SMNEAs) for sensing local intramuscular electromyography signals ex vivo. We use a unique hybrid fabrication scheme based on laser micromachining, microfabrication, and transfer printing to enable scalable fabrication of individually addressable SMNEA with high device stretchability (60 to 90%). The electrode geometries and recording regions, impedance, array layout, and length distribution are highly customizable. We demonstrate the use of SMNEAs as bioelectronic interfaces in recording intramuscular electromyography from various muscle groups in the buccal mass of Aplysia.

Assessment of Mechanically Induced Changes in Helical Fiber Microstructure Using Diffusion Tensor Imaging.

Noninvasive methods to detect microstructural changes in collagen-based fibrous tissues are necessary to differentiate healthy from damaged tissues in vivo but are sparse. Diffusion Tensor Imaging (DTI) is a noninvasive imaging technique used to quantitatively infer tissue microstructure with previous work primarily focused in neuroimaging applications. Yet, it is still unclear how DTI metrics relate to fiber microstructure and function in musculoskeletal tissues such as ligament and tendon, in part because of the high heterogeneity inherent to such tissues. To address this limitation, we assessed the ability of DTI to detect microstructural changes caused by mechanical loading in tissue-mimicking helical fiber constructs of known structure. Using high-resolution optical and micro-computed tomography imaging, we found that static and fatigue loading resulted in decreased sample diameter and a re-alignment of the macro-scale fiber twist angle similar with the direction of loading. However, DTI and micro-computed tomography measurements suggest microstructural differences in the effect of static versus fatigue loading that were not apparent at the bulk level. Specifically, static load resulted in an increase in diffusion anisotropy and a decrease in radial diffusivity suggesting radially uniform fiber compaction. In contrast, fatigue loads resulted in increased diffusivity in all directions and a change in the alignment of the principal diffusion direction away from the constructs' main axis suggesting fiber compaction and microstructural disruptions in fiber architecture. These results provide quantitative evidence of the ability of DTI to detect mechanically induced changes in tissue microstructure that are not apparent at the bulk level, thus confirming its potential as a noninvasive measure of microstructure in helically architected collagen-based tissues, such as ligaments and tendons.

Friction modulation in limbless, three-dimensional gaits and heterogeneous terrains.

Motivated by a possible convergence of terrestrial limbless locomotion strategies ultimately determined by interfacial effects, we show how both 3D gait alterations and locomotory adaptations to heterogeneous terrains can be understood through the lens of local friction modulation. Via an effective-friction modeling approach, compounded by 3D simulations, the emergence and disappearance of a range of locomotory behaviors observed in nature is systematically explained in relation to inhabited environments. Our approach also simplifies the treatment of terrain heterogeneity, whereby even solid obstacles may be seen as high friction regions, which we confirm against experiments of snakes 'diffracting' while traversing rows of posts, similar to optical waves. We further this optic analogy by illustrating snake refraction, reflection and lens focusing. We use these insights to engineer surface friction patterns and demonstrate passive snake navigation in complex topographies. Overall, our study outlines a unified view that connects active and passive 3D mechanics with heterogeneous interfacial effects to explain a broad set of biological observations, and potentially inspire engineering design.

Sensitivity analysis of effective transverse shear viscoelastic and diffusional properties of myelinated white matter.

Motivated by the need to interpret the results from a combined use of in vivo brain Magnetic Resonance Elastography (MRE) and Diffusion Tensor Imaging (DTI), we developed a computational framework to study the sensitivity of single-frequency MRE and DTI metrics to white matter microstructure and cell-level mechanical and diffusional properties. White matter was modeled as a triphasic unidirectional composite, consisting of parallel cylindrical inclusions (axons) surrounded by sheaths (myelin), and embedded in a matrix (glial cells plus extracellular matrix). Only 2D mechanics and diffusion in the transverse plane (perpendicular to the axon direction) was considered, and homogenized (effective) properties were derived for a periodic domain containing a single axon. The numerical solutions of the MRE problem were performed with ABAQUS and by employing a sophisticated boundary-conforming grid generation scheme. Based on the linear viscoelastic response to harmonic shear excitation and steady-state diffusion in the transverse plane, a systematic sensitivity analysis of MRE metrics (effective transverse shear storage and loss moduli) and DTI metric (effective radial diffusivity) was performed for a wide range of microstructural and intrinsic (phase-based) physical properties. The microstructural properties considered were fiber volume fraction, and the myelin sheath/axon diameter ratio. The MRE and DTI metrics are very sensitive to the fiber volume fraction, and the intrinsic viscoelastic moduli of the glial phase. The MRE metrics are nonlinear functions of the fiber volume fraction, but the effective diffusion coefficient varies linearly with it. Finally, the transverse metrics of both MRE and DTI are insensitive to the axon diameter in steady state. Our results are consistent with the limited anisotropic MRE and co-registered DTI measurements, mainly in the corpus callosum, available in the literature. We conclude that isotropic MRE and DTI constitutive models are good approximations for myelinated white matter in the transverse plane. The unidirectional composite model presented here is used for the first time to model harmonic shear stress under MRE-relevant frequency on the cell level. This model can be extended to 3D in order to inform the solution of the inverse problem in MRE, establish the biological basis of MRE metrics, and integrate MRE/DTI with other modalities towards increasing the specificity of neuroimaging.

Lattice Boltzmann method for simulation of diffusion magnetic resonance imaging physics in multiphase tissue models.

We report an implementation of the lattice Boltzmann method (LBM) to integrate the Bloch-Torrey equation, which describes the evolution of the transverse magnetization vector and the fate of the signal of diffusion magnetic resonance imaging (dMRI). Motivated by the need to interpret dMRI experiments in biological tissues, and to offset the small time-step limitation of classical LBM, a hybrid LBM scheme is introduced and implemented to solve the Bloch-Torrey equation. A membrane boundary condition is presented which is able to accurately represent the effects of thin curvilinear membranes typically found in biological tissues. As implemented, the hybrid LBM scheme accommodates piece-wise uniform transport, dMRI parameters, periodic and mirroring outer boundary conditions, and finite membrane permeabilities on non-boundary-conforming inner boundaries. By comparing with analytical solutions of limiting cases, we demonstrate that the hybrid LBM scheme is more accurate than the classical LBM scheme. The proposed explicit LBM scheme maintains second-order spatial accuracy, stability, and first-order temporal accuracy for a wide range of parameters. The parallel implementation of the hybrid LBM code in a multi-CPU computer system, as well as on GPUs, is straightforward and efficient. Along with offering certain advantages over finite element or Monte Carlo schemes, the proposed hybrid LBM constitutes a flexible scheme that can by easily adapted to model more complex interfacial conditions and physics in heterogeneous multiphase tissue models and to accommodate sophisticated dMRI sequences.

Global sensitivity analysis of skeletal muscle dMRI metrics: Effects of microstructural and pulse parameters.

PURPOSE: Estimating microstructural parameters of skeletal muscle from diffusion MRI (dMRI) signal requires understanding the relative importance of both microstructural and dMRI sequence parameters on the signal. This study seeks to determine the sensitivity of dMRI signal to variations in microstructural and dMRI sequence parameters, as well as assess the effect of noise on sensitivity. METHODS: Using a cylindrical myocyte model of skeletal muscle, numerical solutions of the Bloch-Torrey equation were used to calculate global sensitivity indices of dMRI metrics (FA, RD, MD, λ1 , λ2 , λ3 ) for wide ranges of microstructural and dMRI sequence parameters. The microstructural parameters were: myocyte diameter, volume fraction, membrane permeability, intra- and extracellular diffusion coefficients, and intra- and extracellular T2 times. Two separate pulse sequences were examined, a PGSE and a generalized diffusion-weighted sequence that accommodates a larger range of diffusion times. The effect of noise and signal averaging on the sensitivity of the dMRI metrics was examined by adding synthetic noise to the simulated signal. RESULTS: Among the examined parameters, the intracellular diffusion coefficient has the strongest effect, and myocyte diameter is more influential than permeability for FA and RD. The sensitivity indices do not vary significantly between the two pulse sequences. Also, noise strongly affects the sensitivity of the dMRI signal to microstructural variations. CONCLUSIONS: With the identification of key microstructural features that affect dMRI measurements, the reported sensitivity results can help interpret dMRI measurements of skeletal muscle in terms of the underlying microstructure and further develop parsimonious dMRI models of skeletal muscle.

Comparison of two-compartment exchange and continuum models of dMRI in skeletal muscle.

Clinical diffusion MRI (dMRI) is sensitive to micrometer scale spin displacements, but the image resolution is ∼mm, so the biophysical interpretation of the signal relies on establishing appropriate subvoxel tissue models. A class of two-compartment exchange models originally proposed by Kärger have been used successfully in neural tissue dMRI. Their use to interpret the signal in skeletal muscle dMRI is challenging because myocyte diameters are comparable to the root-mean-square of spin displacement and their membrane permeability is high. A continuum tissue model consisting of the Bloch-Torrey equation integrated by a hybrid lattice Boltzmann scheme is used for comparison. The validity domain of a classical two-compartment tissue model is probed by comparing it with the prediction of the continuum model for a 2D unidirectional composite continuum model of myocytes embedded in a uniform matrix. This domain is described in terms of two dimensionless parameters inspired by mass transfer phenomena, the Fourier (F) and Biot (B) numbers. The two-compartment model is valid when [Formula: see text] and [Formula: see text], or when [Formula: see text] and [Formula: see text]. The model becomes less appropriate for muscle dMRI as the cell diameter and volume fraction increase, with the primary source of error associated with modeling diffusion in the extracellular matrix.