Monte Carlo Simulations of Diffusion Weighted MRI in Myocardium: Validation and Sensitivity Analysis

Characterization of anthracycline-induced cardiotoxicity by diffusion tensor magnetic resonance imaging

Anthracyclines are highly potent anti-cancer drugs, but their clinical use is limited by severe cardiotoxic side effects. The impact of anthracycline-induced cardiotoxicity (AIC) on left ventricular (LV) microarchitecture and diffusion properties remains unknown. This study sought to characterize AIC by cardiovascular magnetic resonance diffusion tensor imaging (DTI). Mice were treated with Doxorubicin (DOX; n = 16) for induction of AIC or saline as corresponding control (n = 15). Cardiac function was assessed via echocardiography at the end of the study period. Whole hearts (n = 8 per group) were scanned ex vivo by high-resolution DTI at 7 T. Results were correlated with histopathology and mass spectrometry imaging. Mice with AIC demonstrated systolic dysfunction (LVEF 52 ± 3% vs. 43 ± 6%, P < 0.001), impaired global longitudinal strain (-19.6 ± 2.0% vs. -16.6 ± 3.0%, P < 0.01), and cardiac atrophy (LV mass index [mg/mm], 4.3 ± 0.1 vs. 3.6 ± 0.2, P < 0.01). Regional sheetlet angles were significantly lower in AIC, whereas helix angle and relative helicity remained unchanged. In AIC, fractional anisotropy was increased (0.12 ± 0.01 vs. 0.14 ± 0.02, P < 0.05). DOX-treated mice displayed higher planar and less spherical anisotropy (CPlanar 0.07 ± 0.01 vs. 0.09 ± 0.01, P < 0.01; CSpherical 0.89 ± 0.01 vs. 0.87 ± 0.02, P < 0.05). CPlanar and CSpherical yielded good discriminatory power to distinguish between mice with and without AIC (c-index 0.91 and 0.84, respectively, P for both < 0.05). AIC is associated with regional changes in sheetlet angle but no major abnormalities of global LV microarchitecture. The geometric shape of the diffusion tensor is altered in AIC. DTI may provide a new tool for myocardial characterization in patients with AIC, which warrants future clinical studies to evaluate its diagnostic utility.

Particle-based MR modeling with diffusion, microstructure, and enzymatic reactions

PURPOSE

To develop a novel particle-based in silico MR model and demonstrate applications of this model to signal mechanisms which are affected by the spatial organization of particles, including metabolic reaction kinetics, microstructural effects on diffusion, and radiofrequency (RF) refocusing effects in gradient-echo sequences.

METHODS

The model was developed by integrating a forward solution of the Bloch equations with a Brownian dynamics simulator. Simulation configurations were then designed to model MR signal dynamics of interest, with a primary focus on hyperpolarized 13C MRI methods. Phantom scans and spectrophotometric assays were conducted to validate model results in vitro.

RESULTS

The model accurately reproduced the reaction kinetics of enzyme-mediated conversion of pyruvate to lactate. When varying proportions of restrictive structure were added to the reaction volume, nonlinear changes in the reaction rate measured in vitro were replicated in silico. Modeling of RF refocusing effects characterized the degree of diffusion-weighted contribution from preserved residual magnetization in nonspoiled gradient-echo sequences.

CONCLUSIONS

These results show accurate reproduction of a range of MR signal mechanisms, establishing the model's capability to investigate the multifactorial signal dynamics such as those underlying hyperpolarized 13C MRI data.

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.

CMRsim-A python package for cardiovascular MR simulations incorporating complex motion and flow

PURPOSE

To present an open-source MR simulation framework that facilitates the incorporation of complex motion and flow for studying cardiovascular MR (CMR) acquisition and reconstruction.

METHODS

CMRsim is a Python package that allows simulation of CMR images using dynamic digital phantoms with complex motion as input. Two simulation paradigms are available, namely, numerical and analytical solutions to the Bloch equations, using a common motion representation. Competitive simulation speeds are achieved using TensorFlow for GPU acceleration. To demonstrate the capability of the package, one introductory and two advanced CMR simulation experiments are presented. The latter showcase phase-contrast imaging of turbulent flow downstream of a stenotic section and cardiac diffusion tensor imaging on a contracting left ventricle. Additionally, extensive documentation and example resources are provided.

RESULTS

The Bloch simulation with turbulent flow using approximately 1.5 million particles and a sequence duration of 710 ms for each of the seven different velocity encodings took a total of 29 min on a NVIDIA Titan RTX GPU. The results show characteristic phase contrast and magnitude modulation present in real data. The analytical simulation of cardiac diffusion tensor imaging with bulk-motion phase sensitivity took approximately 10 s per diffusion-weighted image, including preparation and loading steps. The results exhibit the expected alteration of diffusion metrics due to strain.

CONCLUSION

CMRsim is the first simulation framework that allows one to feasibly incorporate complex motion, including turbulent flow, to systematically study advanced CMR acquisition and reconstruction approaches. The open-source package features modularity and transparency, facilitating maintainability and extensibility in support of reproducible research.

Monte Carlo simulation of water diffusion through cardiac tissue models

Monte Carlo diffusion simulations are commonly used to establish a reliable ground truth of tissue microstructure, including for the validation of diffusion-weighted MRI. However, selecting simulation parameters is challenging and affects validity and reproducibility. We conducted experiments to investigate critical conditions in Monte Carlo simulations, such as tissue representation complexity, simulated molecules, update duration, and compartment size. Results show significant changes in microstructure characteristics when parameters are altered, emphasizing the importance of careful control for a reliable ground truth.

Three-dimensional micro-structurally informed in silico myocardium-Towards virtual imaging trials in cardiac diffusion weighted MRI

In silico tissue models (viz. numerical phantoms) provide a mechanism for evaluating quantitative models of magnetic resonance imaging. This includes the validation and sensitivity analysis of imaging biomarkers and tissue microstructure parameters. This study proposes a novel method to generate a realistic numerical phantom of myocardial microstructure. The proposed method extends previous studies by accounting for the variability of the cardiomyocyte shape, water exchange between the cardiomyocytes (intercalated discs), disorder class of myocardial microstructure, and four sheetlet orientations. In the first stage of the method, cardiomyocytes and sheetlets are generated by considering the shape variability and intercalated discs in cardiomyocyte-cardiomyocyte connections. Sheetlets are then aggregated and oriented in the directions of interest. The morphometric study demonstrates no significant difference (p>0.01) between the distribution of volume, length, and primary and secondary axes of the numerical and real (literature) cardiomyocyte data. Moreover, structural correlation analysis validates that the in-silico tissue is in the same class of disorderliness as the real tissue. Additionally, the absolute angle differences between the simulated helical angle (HA) and input HA (reference value) of the cardiomyocytes (4.3°±3.1°) demonstrate a good agreement with the absolute angle difference between the measured HA using experimental cardiac diffusion tensor imaging (cDTI) and histology (reference value) reported by (Holmes et al., 2000) (3.7°±6.4°) and (Scollan et al. 1998) (4.9°±14.6°). Furthermore, the angular distance between eigenvectors and sheetlet angles of the input and simulated cDTI is much smaller than those between measured angles using structural tensor imaging (as a gold standard) and experimental cDTI. Combined with the qualitative results, these results confirm that the proposed method can generate richer numerical phantoms for the myocardium than previous studies.

Anderson-fabry disease: role of traditional and new cardiac MRI techniques

Anderson-Fabry (FD) disease is a rare X-linked disorder caused by different mutations in the Galactosidase α (GLA) gene, which leads to α-galactosidase A enzyme deficiency and the storage of glycosphingolipids in different kinds of organs, included the heart. This results in myocardial inflammation and left ventricular hypertrophy (LVH) and fibrosis. Echocardiography and cardiac magnetic resonance (C-MRI), in particular with new techniques, such as mapping analysis, late gadolinium enhancement (LGE) assessment and strain imaging, are important tools that allow a correct diagnosis, discriminating FD from other hypertrophic heart conditions. C-MRI is able to detect tissue alterations in the early stages of the disease, when an appropriate treatment could be more effective, and it has a fundamental role in monitoring therapy.

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.

Apparent diffusion coefficient measured by diffusion MRI of moving and deforming domains

The modeling of the diffusion MRI signal from moving and deforming organs such as the heart is challenging due to significant motion and deformation of the imaged medium during the signal acquisition. Recently, a mathematical formulation of the Bloch-Torrey equation, describing the complex transverse magnetization due to diffusion-encoding magnetic field gradients, was developed to account for the motion and deformation. In that work, the motivation was to cancel the effect of the motion and deformation in the MRI image and the space scale of interest spans multiple voxels. In the present work, we adapt the mathematical equation to study the diffusion MRI signal at the much smaller scale of biological cells. We start with the Bloch-Torrey equation defined on a cell that is moving and deforming and linearize the equation around the magnitude of the diffusion-encoding gradient. The result is a second order signal model in which the linear term gives the imaginary part of the diffusion MRI signal and the quadratic term gives the apparent diffusion coefficient (ADC) attributable to the biological cell. We numerically validate this model for a variety of motions and deformations.

Novel insights into in-vivo diffusion tensor cardiovascular magnetic resonance using computational modeling and a histology-based virtual microstructure

Purpose

To develop histology‐informed simulations of diffusion tensor cardiovascular magnetic resonance (DT‐CMR) for typical in‐vivo pulse sequences and determine their sensitivity to changes in extra‐cellular space (ECS) and other microstructural parameters.

Methods

We synthesised the DT‐CMR signal from Monte Carlo random walk simulations. The virtual tissue was based on porcine histology. The cells were thickened and then shrunk to modify ECS. We also created idealised geometries using cuboids in regular arrangement, matching the extra‐cellular volume fraction (ECV) of 16–40%. The simulated voxel size was 2.8 × 2.8 × 8.0 mm³ for pulse sequences covering short and long diffusion times: Stejskal–Tanner pulsed‐gradient spin echo, second‐order motion‐compensated spin echo, and stimulated echo acquisition mode (STEAM), with clinically available gradient strengths.

Results

The primary diffusion tensor eigenvalue increases linearly with ECV at a similar rate for all simulated geometries. Mean diffusivity (MD) varies linearly, too, but is higher for the substrates with more uniformly distributed ECS. Fractional anisotropy (FA) for the histology‐based geometry is higher than the idealised geometry with low sensitivity to ECV, except for the long mixing time of the STEAM sequence. Varying the intra‐cellular diffusivity (D_IC) results in large changes of MD and FA. Varying extra‐cellular diffusivity or using stronger gradients has minor effects on FA. Uncertainties of the primary eigenvector orientation are reduced using STEAM.

Conclusions

We found that the distribution of ECS has a measurable impact on DT‐CMR parameters. The observed sensitivity of MD and FA to ECV and D_IC has potentially interesting applications for interpreting in‐vivo DT‐CMR parameters.