Evaluation and comparison of diffusion MR methods for measuring apparent transcytolemmal water exchange rate constant

Comparison of water exchange measurements between filter-exchange imaging and diffusion time-dependent kurtosis imaging in the human brain

PURPOSE

Filter-exchange imaging (FEXI) and diffusion time (t)-dependent kurtosis imaging (DKI(t)) are two diffusion-based methods that have been proposed for in vivo measurements of water exchange rates. Few studies have directly compared these methods. We aimed to investigate whether FEXI and DKI(t) yield comparable water exchange measurements in the human brain in vivo.

METHODS

Eight healthy volunteers underwent multiple-direction FEXI and DKI(t) acquisitions on a 3T scanner. We performed region of interest (ROI) analysis to determine correlations between FEXI-derived apparent exchange rate (AXR) and DKI(t)-derived reciprocal of exchange time (1 / τ ex $$ 1/{\tau}_{ex} $$ ).

RESULTS

In both white matter (WM) and gray matter (GM), DKI(t) revealed substantial diffusion-time dependence of diffusivity and kurtosis. However, at t ≥ 100 ms, the diffusivity showed weak time dependence. In WM, this time dependence may be due to water exchange between myelin water and "free" water with different T1 values, although other factors, such as remaining restrictive effects from microstructural barriers, cannot be excluded. We found a significant correlation between DKI(t)-derived1 / τ ex $$ 1/{\tau}_{ex} $$ and FEXI-derived AXR in the axial direction within WM. No such correlation was present in GM, although both values showed similar ranges.

CONCLUSION

These results suggest that FEXI and DKI(t) could be sensitive to the same water exchange process only when the diffusion time in DKI(t) is sufficiently long, and only in WM. In both GM and WM, the restrictive effect of microstructure is non-negligible, especially at short diffusion times (<100 ms).

MATI: A GPU-accelerated toolbox for microstructural diffusion MRI simulation and data fitting with a graphical user interface

PURPOSE

To introduce MATI (Microstructural Analysis Toolbox for Imaging), a versatile MATLAB-based toolbox that combines both simulation and data fitting capabilities for microstructural dMRI research.

METHODS

MATI provides a user-friendly, graphical user interface that enables researchers, including those without much programming experience, to perform advanced simulations and data analyses for microstructural MRI research. For simulation, MATI supports arbitrary microstructural tissues and pulse sequences. For data fitting, MATI supports a range of fitting methods, including traditional non-linear least squares, Bayesian approaches, machine learning, and dictionary matching methods, allowing users to tailor analyses based on specific research needs.

RESULTS

Optimized with vectorized matrix operations and high-performance numerical libraries, MATI achieves high computational efficiency, enabling rapid simulations and data fitting on CPU and GPU hardware. While designed for microstructural dMRI, MATI's generalized framework can be extended to other imaging methods, making it a flexible and scalable tool for quantitative MRI research.

CONCLUSION

MATI offers a significant step toward translating advanced microstructural MRI techniques into clinical applications.

tDKI-Net: A Joint q-t Space Learning Network for Diffusion-Time-Dependent Kurtosis Imaging

Time-dependent diffusion magnetic resonance imaging (TDDMRI) is useful for the non-invasive characterization of tissue microstructure. These models require densely sampled q-t space data for microstructural fitting, leading to very time-consuming acquisition protocols. To overcome this problem, we present a joint q-t space model-tDKI-Net to estimate diffusion-time dependent kurtosis and the transmembrane exchange, using downsampled q-t space data. The tDKI-Net is composed of several q-Encoders and a t-Encoder, designed based on the extragradient mechanism, each integrated with their respective mapping networks. In the tDKI-Net, two types of encoders along with their mapping networks are employed sequentially to generate kurtosis at individual diffusion times and to fit the transmembrane exchange time () using the time-dependent kurtosis according to the Kärger's model. Meanwhile, we proposed a three-stage training strategy, including physics-informed self-supervised pretraining, DKI warm-up, and joint training, to match the network structure. Our results demonstrated that the proposed tDKI-Net could effectively accelerate tDKI scans, resulting in lower estimation error compared with other methods. Our proposed three-stage training strategy demonstrated superior results than those training from scratch, e.g., the normalized root mean square error (NRMSE) of decreased by up to 1.4%. We also investigated the training data size effects and found that although we used one-subject training, the network achieved lower NRMSEs for , and (2.50%, 3.04%, 10.86%) than previous work that used three-subject training (3.8%, 9.5%, 12.1%). tDKI-Net can considerably reduce the scan time by 10.5-fold by joint downsampling the q-t space data without compromising the estimation accuracy.

Tuned exchange imaging: Can the filter exchange imaging pulse sequence be adapted for applications with thin slices and restricted diffusion?

Filter exchange imaging (FEXI) is a double diffusion-encoding (DDE) sequence that is specifically sensitive to exchange between sites with different apparent diffusivities. FEXI uses a diffusion-encoding filtering block followed by a detection block at varying mixing times to map the exchange rate. Long mixing times enhance the sensitivity to exchange, but they pose challenges for imaging applications that require a stimulated echo sequence with crusher gradients. Thin imaging slices require strong crushers, which can introduce significant diffusion weighting and bias exchange rate estimates. Here, we treat the crushers as an additional encoding block and consider FEXI as a triple diffusion-encoding sequence. This allows the bias to be corrected in the case of multi-Gaussian diffusion, but not easily in the presence of restricted diffusion. Our approach addresses challenges in the presence of restricted diffusion and relies on the ability to independently gauge sensitivities to exchange and restricted diffusion for arbitrary gradient waveforms. It follows two principles: (i) the effects of crushers are included in the forward model using signal cumulant expansion; and (ii) timing parameters of diffusion gradients in filter and detection blocks are adjusted to maintain the same level of restriction encoding regardless of the mixing time. This results in the tuned exchange imaging (TEXI) protocol. The accuracy of exchange mapping with TEXI was assessed through Monte Carlo simulations in spheres of identical sizes and gamma-distributed sizes, and in parallel hexagonally packed cylinders. The simulations demonstrate that TEXI provides consistent exchange rates regardless of slice thickness and restriction size, even with strong crushers. However, the accuracy depends on b-values, mixing times, and restriction geometry. The constraints and limitations of TEXI are discussed, including suggestions for protocol adaptations. Further studies are needed to optimize the precision of TEXI and assess the approach experimentally in realistic, heterogeneous substrates.

Filter-exchange spectroscopy is sensitive to gradual cell membrane degradation

Transmembrane water permeability changes occur after initialization of necrosis and are a mechanism for early detection of cell death. Filter-exchange spectroscopy (FEXSY) is sensitive to transmembrane water permeability and enables its quantification by magnetic resonance via the apparent exchange rate (AXR). In this study, we investigate AXR changes during necrotic cell death. FEXSY measurements of yeast cells in different necrotic stages were performed and compared with established fluorescence cell death markers and pulsed gradient spin echo measurements. Furthermore, the influence of T2 relaxation on AXR was examined in a two-compartment system. The AXR of yeast cells increased slightly after incubation with 20% isopropanol, whereas it peaked sharply after incubation with 25% isopropanol. At this point, almost all the yeast cells were vital but showed compromised membranes. After incubation with 30% isopropanol, AXR measurements showed high variability, at a point corresponding to a majority of the yeast cells being in late-stage necrosis with disrupted cell membranes. Simulations revealed that, for FEXSY measurements in a two-compartment system, a long filter echo time (TEf), compared with the T2 of the slow-diffusing compartment, filters out a fraction of the slow-diffusing compartment signal and leads to overestimation of apparent diffusion coefficient (ADC) and underestimation of AXR. Our results demonstrate that AXR is sensitive to gradual permeabilization of the cell membrane of living cells in different permeabilization stages without exogenous contrast agents. AXR measurements were sensitive to permeability changes induced by relatively low concentrations of isopropanol, at levels for which no measurable effect was detectable by ADC measurements. TEf may act as a signal filter that affects the estimated AXR value of a system consisting of a variety of local diffusivities and a range of T2 that includes T2 values shorter or comparable with the TEf.

Restriction-induced time-dependent transcytolemmal water exchange: Revisiting the Kӓrger exchange model

The Kӓrger model and its derivatives have been widely used to incorporate transcytolemmal water exchange rate, an essential characteristic of living cells, into analyses of diffusion MRI (dMRI) signals from tissues. The Kӓrger model consists of two homogeneous exchanging components coupled by an exchange rate constant and assumes measurements are made with sufficiently long diffusion time and slow water exchange. Despite successful applications, it remains unclear whether these assumptions are generally valid for practical dMRI sequences and biological tissues. In particular, barrier-induced restrictions to diffusion produce inhomogeneous magnetization distributions in relatively large-sized compartments such as cancer cells, violating the above assumptions. The effects of this inhomogeneity are usually overlooked. We performed computer simulations to quantify how restriction effects, which in images produce edge enhancements at compartment boundaries, influence different variants of the Kӓrger-model. The results show that the edge enhancement effect will produce larger, time-dependent estimates of exchange rates in e.g., tumors with relatively large cell sizes (>10 μm), resulting in overestimations of water exchange as previously reported. Moreover, stronger diffusion gradients, longer diffusion gradient durations, and larger cell sizes, all cause more pronounced edge enhancement effects. This helps us to better understand the feasibility of the Kärger model in estimating water exchange in different tissue types and provides useful guidance on signal acquisition methods that may mitigate the edge enhancement effect. This work also indicates the need to correct the overestimated transcytolemmal water exchange rates obtained assuming the Kärger-model.

Comprehensive characterization of tumor therapeutic response with simultaneous mapping cell size, density, and transcytolemmal water exchange

Early assessment of tumor therapeutic response is an important topic in precision medicine to optimize personalized treatment regimens and reduce unnecessary toxicity, cost, and delay. Although diffusion MRI (dMRI) has shown potential to address this need, its predictive accuracy is limited, likely due to its unspecific sensitivity to overall pathological changes. In this work, we propose a new quantitative dMRI-based method dubbed EXCHANGE (MRI of water Exchange, Confined and Hindered diffusion under Arbitrary Gradient waveform Encodings) for simultaneous mapping of cell size, cell density, and transcytolemmal water exchange. Such rich microstructural information comprehensively evaluates tumor pathologies at the cellular level. Validations using numerical simulations and in vitro cell experiments confirmed that the EXCHANGE method can accurately estimate mean cell size, density, and water exchange rate constants. The results from in vivo animal experiments show the potential of EXCHANGE for monitoring tumor treatment response. Finally, the EXCHANGE method was implemented in breast cancer patients with neoadjuvant chemotherapy, demonstrating its feasibility in assessing tumor therapeutic response in clinics. In summary, a new, quantitative dMRI-based EXCHANGE method was proposed to comprehensively characterize tumor microstructural properties at the cellular level, suggesting a unique means to monitor tumor treatment response in clinical practice.

A Magnetic Resonance Imaging-Chemical Exchange Saturation Transfer (MRI-CEST) Method for the Detection of Water Cycling across Cellular Membranes

Water cycling across the membrane transporters is considered a hallmark of cellular metabolism and it could be of high diagnostic relevance in the characterization of tumors and other diseases. The method relies on the response of intracellular proton exchanging molecules to the presence of extracellular Gd-based contrast agents (GBCAs). Paramagnetic GBCAs enhances the relaxation rate of water molecules in the extracellular compartment and, through membrane exchange, the relaxation enhancement is transferred to intracellular molecules. The effect is detected at the MRI-CEST (Magnetic Resonance Imaging - Chemical Exchange Saturation Transfer) signal of intracellular proton exchanging molecules. The magnitude of the change in the CEST response reports on water cycling across the membrane. The method has been tested on Red Blood Cells and on orthotopic murine models of breast cancer with different degree of malignancy (4T1, TS/A and 168FARN). The distribution of voxels reporting on membrane permeability fits well with the cells' aggressiveness and acts as an early reporter to monitor therapeutic treatments.

Measuring water exchange on a preclinical MRI system using filter exchange and diffusion time dependent kurtosis imaging

PURPOSE

Filter exchange imaging (FEXI) and diffusion time (t)-dependent diffusion kurtosis imaging (DKI(t)) are both sensitive to water exchange between tissue compartments. The restrictive effects of tissue microstructure, however, introduce bias to the exchange rate obtained by these two methods, as their interpretation conventionally rely on the Kärger model of barrier limited exchange between Gaussian compartments. Here, we investigated whether FEXI and DKI(t) can provide comparable exchange rates in ex vivo mouse brains.

THEORY AND METHODS

FEXI and DKI(t) data were acquired from ex vivo mouse brains on a preclinical MRI system. Phase cycling and negative slice prewinder gradients were used to minimize the interferences from imaging gradients.

RESULTS

In the corpus callosum, apparent exchange rate (AXR) from FEXI correlated with the exchange rate (the inverse of exchange time, 1/τex ) from DKI(t) along the radial direction. In comparison, discrepancies between FEXI and DKI(t) were found in the cortex due to low filter efficiency and confounding effects from tissue microstructure.

CONCLUSION

The results suggest that FEXI and DKI(t) are sensitive to the same exchange processes in white matter when separated from restrictive effects of microstructure. The complex microstructure in gray matter, with potential exchange among multiple compartments and confounding effects of microstructure, still pose a challenge for FEXI and DKI(t).

Metabolic activity diffusion imaging (MADI): I. Metabolic, cytometric modeling and simulations

Evidence mounts that the steady-state cellular water efflux (unidirectional) first-order rate constant (k_io [s^-1 ]) magnitude reflects the ongoing, cellular metabolic rate of the cytolemmal Na⁺ , K⁺ -ATPase (NKA), ^c MR_NKA (pmol [ATP consumed by NKA]/s/cell), perhaps biology's most vital enzyme. Optimal ¹ H₂ O MR k_io determinations require paramagnetic contrast agents (CAs) in model systems. However, results suggest that the homeostatic metabolic k_io biomarker magnitude in vivo is often too large to be reached with allowable or possible CA living tissue distributions. Thus, we seek a noninvasive (CA-free) method to determine k_io in vivo. Because membrane water permeability has long been considered important in tissue water diffusion, we turn to the well-known diffusion-weighted MRI (DWI) modality. To analyze the diffusion tensor magnitude, we use a parsimoniously primitive model featuring Monte Carlo simulations of water diffusion in virtual ensembles comprising water-filled and -immersed randomly sized/shaped contracted Voronoi cells. We find this requires two additional, cytometric properties: the mean cell volume (V [pL]) and the cell number density (ρ [cells/μL]), important biomarkers in their own right. We call this approach metabolic activity diffusion imaging (MADI). We simulate water molecule displacements and transverse MR signal decays covering the entirety of b-space from pure water (ρ = V = 0; k_io undefined; diffusion coefficient, D₀ ) to zero diffusion. The MADI model confirms that, in compartmented spaces with semipermeable boundaries, diffusion cannot be described as Gaussian: the nanoscopic D (D_n ) is diffusion time-dependent, a manifestation of the "diffusion dispersion". When the "well-mixed" (steady-state) condition is reached, diffusion becomes limited, mainly by the probabilities of (1) encountering (ρ, V), and (2) permeating (k_io ) cytoplasmic membranes, and less so by D_n magnitudes. Importantly, for spaces with large area/volume (A/V; claustrophobia) ratios, this can happen in less than a millisecond. The model matches literature experimental data well, with implications for DWI interpretations.

What does FEXI measure?

Filter-exchange imaging (FEXI) has already been utilized in several biomedical studies for evaluating the permeability of cell membranes. The method relies on suppressing the extracellular signal using strong diffusion weighting (the mobility filter causing a reduction in the overall diffusivity) and monitoring the subsequent diffusivity recovery. Using Monte Carlo simulations, we demonstrate that FEXI is sensitive not uniquely to the transcytolemmal exchange but also to the geometry of involved compartments: complex geometry offers locations where spins remain unaffected by the mobility filter; moving to other locations afterwards, such spins contribute to the diffusivity recovery without actually permeating any membrane. This exchange mechanism is a warning for those who aim to use FEXI in complex media such as brain gray matter and opens wide scope for investigation towards crystallizing the genuine membrane permeation and characterizing the compartment geometry.

Improving MR cell size imaging by inclusion of transcytolemmal water exchange

The goal of the current study is to include transcytolemmal water exchange in MR cell size imaging using the IMPULSED model for more accurate characterization of tissue cellular properties (e.g., apparent volume fraction of intracellular space v in ) and quantification of indicators of transcytolemmal water exchange. We propose a heuristic model that incorporates transcytolemmal water exchange into a multicompartment diffusion-based method (IMPULSED) that was developed previously to extract microstructural parameters (e.g., mean cell size d and apparent volume fraction of intracellular space v in ) assuming no water exchange. For t diff ≤ 5 ms, the water exchange can be ignored, and the signal model is the same as the IMPULSED model. For t diff ≥ 30 ms, we incorporated the modified Kärger model that includes both restricted diffusion and exchange between compartments. Using simulations and previously published in vitro cell data, we evaluated the accuracy and precision of model-derived parameters and determined how they are dependent on SNR and imaging parameters. The joint model provides more accurate d values for cell sizes ranging from 10 to 12 microns when water exchange is fast (e.g., intracellular water pre-exchange lifetime τ in ≤ 100 ms) than IMPULSED, and reduces the bias of IMPULSED-derived estimates of v in , especially when water exchange is relatively slow (e.g., τ in &gt; 200 ms). Indicators of transcytolemmal water exchange derived from the proposed joint model are linearly correlated with ground truth τ in values and can detect changes in cell membrane permeability induced by saponin treatment in murine erythroleukemia cancer cells. Our results suggest this joint model not only improves the accuracy of IMPULSED-derived microstructural parameters, but also provides indicators of water exchange that are usually ignored in diffusion models of tissues.

Disentangling the effects of restriction and exchange with diffusion exchange spectroscopy

Diffusion exchange spectroscopy (DEXSY) is a multidimensional NMR technique that can reveal how water molecules exchange between compartments within heterogeneous media, such as biological tissue. Data from DEXSY experiments is typically processed using numerical inverse Laplace transforms (ILTs) to produce a diffusion-diffusion spectrum. A tacit assumption of this ILT approach is that the signal behavior is Gaussian - i.e., the spin echo intensity decays exponentially with the degree of diffusion weighting. The assumptions that underlie Gaussian signal behavior may be violated, however, depending on the gradient strength applied and the sample under study. We argue that non-Gaussian signal behavior due to restrictions is to be expected in the study of biological tissue using diffusion NMR. Further, we argue that this signal behavior can produce confounding features in the diffusion-diffusion spectra obtained from numerical ILTs of DEXSY data - entangling the effects of restriction and exchange. Specifically, restricted signal behavior can result in broadening of peaks and in the appearance of illusory exchanging compartments with distributed diffusivities, which pearl into multiple peaks if not highly regularized. We demonstrate these effects on simulated data. That said, we suggest the use of features in the signal acquisition domain that can be used to rapidly probe exchange without employing an ILT. We also propose a means to characterize the non-Gaussian signal behavior due to restrictions within a sample using DEXSY measurements with a near zero mixing time or storage interval. We propose a combined acquisition scheme to independently characterize restriction and exchange with various DEXSY measurements, which we term Restriction and Exchange from Equally-weighted Double and Single Diffusion Encodings (REEDS-DE). We test this method on ex vivo neonatal mouse spinal cord - a sample consisting primarily of gray matter - using a low-field, static gradient NMR system. In sum, we highlight critical shortcomings of prevailing DEXSY analysis methods that conflate the effects of restriction and exchange, and suggest a viable experimental approach to disentangle them.

Measurement of cellular-interstitial water exchange time in tumors based on diffusion-time-dependent diffusional kurtosis imaging

PURPOSE

To assess the feasibility of using diffusion-time-dependent diffusional kurtosis imaging (tDKI) to measure cellular-interstitial water exchange time (τex ) in tumors, both in animals and in humans.

METHODS

Preclinical tDKI studies at 7 T were performed with the GL261 glioma model and the 4T1 mammary tumor model injected into the mouse brain. Clinical studies were performed at 3 T with women who had biopsy-proven invasive ductal carcinoma. tDKI measurement was conducted using a diffusion-weighted STEAM pulse sequence with multiple diffusion times (20-800 ms) at a fixed echo time, while keeping the b-values the same (0-3000 s/mm2 ) by adjusting the diffusion gradient strength. The tDKI data at each diffusion time t were used for a weighted linear least-squares fit method to estimate the diffusion-time-dependent diffusivity, D(t), and diffusional kurtosis, K(t).

RESULTS

Both preclinical and clinical studies showed that, when diffusion time t ≥ 200 ms, D(t) did not have a noticeable change while K(t) decreased monotonically with increasing diffusion time in tumors and t ≥ 100 ms for the cortical ribbon of the mouse brain. The estimated τex averaged median and interquartile range (IQR) of GL261 and 4T1 tumors were 93 (IQR = 89) ms and 68 (78) ms, respectively. For the cortical ribbon, the estimated τex averaged median and IQR were 41 (34) ms for C57BL/6 and 30 (17) ms for BALB/c. For invasive ductal carcinoma, the estimated τex median and IQR of the two breast cancers were 70 (94) and 106 (92) ms.

CONCLUSION

The results of this proof-of-concept study substantiate the feasibility of using tDKI to measure cellular-interstitial water exchange time without using an exogenous contrast agent.

Microstructure Modeling of High b-Value Diffusion-Weighted Images in Glioblastoma

Apparent diffusion coefficient has limits to differentiate solid tumor from normal tissue or edema in glioblastoma (GBM). This study investigated a microstructure model (MSM) in GBM using a clinically available diffusion imaging technique. The MSM was modified to integrate with bi-polar diffusion gradient waveforms, and applied to 30 patients with newly diagnosed GBM. Diffusion-weighted (DW) images acquired on a 3 T scanner with b-values from 0 to 2500 s/mm² were fitted in volumes of interest (VOIs) of solid tumor to obtain the apparent restriction size of intracellular water (ARS), the fractional volume of intracellular water (V_in), and extracellular (D_ex) water diffusivity. The parameters in solid tumor were compared with those of other tissue types by Students’ t test. For comparison, DW images were fitted by conventional mono-exponential and bi-exponential models. ARS, D_ex, and V_in from the MSM in tumor VOIs were significantly greater than those in WM, GM, and edema (P values of .01–.001). ARS values in solid tumors (from 21.6 to 34.5 um) had absolutely no overlap with those in all other tissue types (from 0.9 to 3.5 um). V_in values showed a descending order from solid tumor (from 0.32 to 0.52) to WM, GM, and edema (from 0.05 to 0.25), consisting with the descending cellularity in these tissue types. The parameters from mono-exponential and bi-exponential models could not significantly differentiate solid tumor from all other tissue types, particularly from edema. Further development and histopathological validation of the MSM will warrant its role in clinical management of GBM.

Apparent exchange rate imaging: On its applicability and the connection to the real exchange rate

PURPOSE

Water exchange between the intracellular and extracellular space can be measured using apparent exchange rate (AXR) imaging. The aim of this study was to investigate the relationship between the measured AXR and the geometry of diffusion restrictions, membrane permeability, and the real exchange rate, as well as to explore the applicability of AXR for typical human measurement settings.

METHODS

The AXR measurements and the underlying exchange rates were simulated using the Monte Carlo method with different geometries, size distributions, packing densities, and a broad range of membrane permeabilities. Furthermore, the influence of SNR and sequence parameters was analyzed.

RESULTS

The estimated AXR values correspond to the simulated values and show the expected proportionality to membrane permeability, except for fast exchange (ie, AXR > 20 - 30 s - 1 ) and small packing densities. Moreover, it was found that the duration of the filter gradient must be shorter than 2 · AX R - 1 . In cell size and permeability distributions, AXR depends on the average surface-to-volume ratio, permeability, and the packing density. Finally, AXR can be reliably determined in the presence of orientation dispersion in axon-like structures with sufficient gradient sampling (ie, 30 gradient directions).

CONCLUSION

Currently used experimental settings for in vivo human measurements are well suited for determining AXR, with the exception of single-voxel analysis, due to limited SNR. The detection of changes in membrane permeability in diseased tissue is nonetheless challenging because of the AXR dependence on further factors, such as packing density and geometry, which cannot be disentangled without further knowledge of the underlying cell structure.

Real-time measurement of diffusion exchange rate in biological tissue

Diffusion exchange spectroscopy (DEXSY) provides a means to isolate the signal attenuation associated with exchange from other sources of signal loss. With the total diffusion weighting b₁+b₂=b_s held constant, DEXSY signals acquired with b₁=0 or b₂=0 have no exchange weighting, while a DEXSY signal acquired with b₁=b₂ has maximal exchange weighting. The exchange rate can be estimated by fitting a diffusion exchange model to signals acquired with variable mixing times. Conventionally, acquired signals are normalized by a signal with b₁=0 and b₂=0 to remove the decay due to spin-lattice relaxation. Instead, division by a signal with equal b_s but b₁=0 or b₂=0 reduces spin-lattice relaxation weighting of the apparent exchange rate (AXR). Furthermore, apparent diffusion-weighted R₁ relaxation rates can be estimated from non-exchange-weighted DEXSY signals. Estimated R₁ values are utilized to remove signal decay due to spin-lattice relaxation from exchange-weighted signals, permitting a more precise estimate of AXR with less data. Data reduction methods are proposed and tested with regards to statistical accuracy and precision of AXR estimates on simulated and experimental data. Simulations show that the methods are capable of accurately measuring the ground-truth exchange rate. The methods remain accurate even when the assumption that DEXSY signals attenuate with b is violated, as occurs for restricted diffusion. Experimental data was collected from fixed neonatal mouse spinal cord samples at 25 and 7°C using the strong static magnetic field gradient produced by a single-sided permanent magnet (i.e., an NMR MOUSE). The most rapid method for exchange measurements requires only five data points (an 80 s experiment as implemented) and achieves a similar level of accuracy and precision to the baseline method using 44 data points. This represents a significant improvement in acquisition speed, overcoming a barrier which has limited the use of DEXSY on living specimen.

In silico simulations of erythrocyte aquaporins with quantitative in vitro validation

Modelling water and membrane lipids is an essential element in the computational research of biophysical/biochemical processes such as water transport across the cell membrane. In this study, we examined the accuracies of two popular water models, TIP3P and TIP4P, in the molecular dynamics simulations of erythrocyte aquaporins (AQP1 and AQP3). We modelled the erythrocyte membrane as an asymmetric lipid bilayer with appropriate lipid compositions of its inner and outer leaflet, in comparison with a symmetric lipid bilayer of a single lipid type. We computed the AQP1/3 permeabilities with the transition state theory with full correction for recrossing events. We also conducted cell swelling assays for water transport across the erythrocyte membrane. The experimental results agree with the TIP3P water–erythrocyte membrane model, in confirmation of the expected accuracy of the erythrocyte membrane model, the TIP3P water model, and the CHARMM parameters for water–protein interactions.

Quantitatively predictive study of aquaporins in model erythrocyte membrane validated with cellular experiments.

Constant gradient FEXSY: A time-efficient method for measuring exchange

Filter-Exchange NMR Spectroscopy (FEXSY) is a method for measurement of apparent transmembranal water exchange rates. The experiment is comprised of two co-linear sequential pulsed-field gradient (PFG) blocks, separated by a mixing period in which exchange takes place. The first block remains constant and serves as a diffusion-based filter that removes signal coming from fast-diffusing water. The mixing time and the gradient area (q-value) of the second block are varied on repeated iterations to produce a 2D data set that is analyzed using a bi-compartmental model which assumes that intra- and extra-cellular water are slow and fast diffusing, respectively. Here we suggest a variant of the FEXSY method in which measurements for different mixing times are taken at a constant gradient. This Constant Gradient FEXSY (CG-FEXSY) allows for the determination of the exchange rate by using a smaller 1D data set, so that the same information can be gathered during a considerably shorter scan time. Furthermore, in the limit of high diffusion weighting, such that the extra-cellular water signal is removed while the intra-cellular signal is retained, CG-FEXSY also allows for determination of the intra-cellular mean residence time (MRT). The theoretical results are validated on a living yeast cells sample and on a fixed porcine optic nerve, where the values obtained from the two methods are shown to be in agreement.

Higher-order diffusion MRI characterization of mesorectal lymph nodes in rectal cancer

PURPOSE

Mesorectal lymph node staging plays an important role in treatment decision making. Here, we explore the benefit of higher-order diffusion MRI models accounting for non-Gaussian diffusion effects to classify mesorectal lymph nodes both 1) ex vivo at ultrahigh field correlated with histology and 2) in vivo in a clinical scanner upon patient staging.

METHODS

The preclinical investigation included 54 mesorectal lymph nodes, which were scanned at 16.4 T with an extensive diffusion MRI acquisition. Eight diffusion models were compared in terms of goodness of fit, lymph node classification ability, and histology correlation. In the clinical part of this study, 10 rectal cancer patients were scanned with diffusion MRI at 1.5 T, and 72 lymph nodes were analyzed with Apparent Diffusion Coefficient (ADC), Intravoxel Incoherent Motion (IVIM), Kurtosis, and IVIM-Kurtosis.

RESULTS

Compartment models including restricted and anisotropic diffusion improved the preclinical data fit, as well as the lymph node classification, compared to standard ADC. The comparison with histology revealed only moderate correlations, and the highest values were observed between diffusion anisotropy metrics and cell area fraction. In the clinical study, the diffusivity from IVIM-Kurtosis was the only metric showing significant differences between benign (0.80 ± 0.30 μm² /ms) and malignant (1.02 ± 0.41 μm² /ms, P = .03) nodes. IVIM-Kurtosis also yielded the largest area under the receiver operating characteristic curve (0.73) and significantly improved the node differentiation when added to the standard visual analysis by experts based on T₂ -weighted imaging.

CONCLUSION

Higher-order diffusion MRI models perform better than standard ADC and may be of added value for mesorectal lymph node classification in rectal cancer patients.

Advances in Diffusion and Perfusion MRI for Quantitative Cancer Imaging

PURPOSE OF REVIEW

This article is to review recent technical developments and their clinical applications in cancer imaging quantitative measurement of cellular and vascular properties of the tumors.

RECENT FINDINGS

Rapid development of fast Magnetic Resonance Imaging (MRI) technologies over last decade brought new opportunities in quantitative MRI methods to measure both cellular and vascular properties of tumors simultaneously.

SUMMARY

Diffusion MRI (dMRI) and dynamic contrast enhanced (DCE)-MRI have become widely used to assess the tissue structural and vascular properties, respectively. However, the ultimate potential of these advanced imaging modalities has not been fully exploited. The dependency of dMRI on the diffusion weighting gradient strength and diffusion time can be utilized to measure tumor perfusion, cellular structure, and cellular membrane permeability. Similarly, DCE-MRI can be used to measure vascular and cellular membrane permeability along with cellular compartment volume fractions. To facilitate the understanding of these potentially important methods for quantitative cancer imaging, we discuss the basic concepts and recent developments, as well as future directions for further development.

INDIANA: An in-cell diffusion method to characterize the size, abundance and permeability of cells

NMR and MRI diffusion experiments contain information describing the shape, size, abundance, and membrane permeability of cells although extracting this information can be challenging. Here we present the INDIANA (IN-cell DIffusion ANAlysis) method to simultaneously and non-invasively measure cell abundance, effective radius, permeability and intrinsic relaxation rates and diffusion coefficients within the inter- and intra-cellular populations. The method couples an experimental dataset comprising stimulated-echo diffusion measurements, varying both the gradient strength and the diffusion delay, together with software to fit a model based on the Kärger equations to robustly extract the relevant parameters. A detailed error analysis is presented by comparing the results from fitting simulated data from Monte Carlo simulations, establishing its effectiveness. We note that for parameters typical of mammalian cells the approach is particularly effective, and the shape of the underlying cells does not unduly affect the results. Finally, we demonstrate the performance of the experiment on systems of suspended yeast and mammalian cells. The extracted parameters describing cell abundance, size, permeability and relaxation are independently validated.

Single-channel permeability and glycerol affinity of human aquaglyceroporin AQP3

For its fundamental relevance, transport of water and glycerol across the erythrocyte membrane has long been investigated before and after the discovery of aquaporins (AQPs), the membrane proteins responsible for water and glycerol transport. AQP1 is abundantly expressed in the human erythrocyte for maintaining its hydrohomeostasis where AQP3 is also expressed (at a level ~30-folds lower than AQP1) facilitating glycerol transport. This research is focused on two of the remaining questions: How permeable is AQP3 to water? What is the glycerol-AQP3 affinity under near-physiological conditions? Through atomistic modelling and large-scale simulations, we found that AQP3 is two to three times more permeable to water than AQP1 and that the glycerol-AQP3 affinity is approximately 500/M. Using these computed values along with the data from the latest literature on AQP1 and on erythrocyte proteomics, we estimated the water and glycerol transport rates across the membrane of an entire erythrocyte. We used these rates to predict the time courses of erythrocyte swelling-shrinking in response to inward and outward osmotic gradients. Experimentally, we monitored the time course of human erythrocytes when subject to an osmotic or glycerol gradient with light scattering in a stopped-flow spectrometer. We observed close agreement between the experimentally measured and the computationally predicted time courses of erythrocytes, which corroborated our computational conclusions on the AQP3 water-permeability and the glycerol-AQP3 affinity.

Rapid detection of the presence of diffusion exchange

Diffusion exchange spectroscopy (DEXSY) provides a detailed picture of how fluids in different microenvironments communicate with one another but requires a large amount of data. For DEXSY MRI, a simple measure of apparent exchanging fractions may suffice to characterize and differentiate materials and tissues. Reparameterizing signal intensity from a PGSE-storage-PGSE experiment as a function of the sum, b_s=b₁+b₂, and difference b_d=b₂-b₁ of the diffusion encodings separates diffusion weighting from exchange weighting. Exchange leads to upward curvature along a slice of constant b_s. Exchanging fractions can be measured rapidly by a finite difference approximation of the curvature using four data points. The method is generalized for non-steady-state and multi-site exchange. We apply the method to image exchanging fractions and calculate exchange rates of water diffusing across the bulk water interface of a glass capillary array.

Using 1H2O MR to measure and map sodium pump activity in vivo

The cell plasma membrane Na⁺,K⁺-ATPase [NKA] is one of biology's most [if not the most] significant enzymes. By actively transporting Na⁺ out [and K⁺ in], it maintains the vital trans-membrane ion concentration gradients and the membrane potential. The forward NKA reaction is shown in the Graphical Abstract [which is elaborated in the text]. Crucially, NKA does not operate in isolation. There are other transporters that conduct K⁺ back out of [II, Graphical Abstract] and Na⁺ back into [III, Graphical Abstract] the cell. Thus, NKA must function continually. Principal routes for ATP replenishment include mitochondrial oxidative phosphorylation, glycolysis, and creatine kinase [CrK] activity. However, it has never been possible to measure, let alone map, this integrated, cellular homeostatic NKA activity in vivo. Active trans-membrane water cycling [AWC] promises a way to do this with ¹H₂O MR. Inthe Graphical Abstract, the AWC system is characterized by active contributions totheunidirectional rate constants for steady-state water efflux and influx, respectively, k_io(a) and k_oi(a). The discovery, validation, and initial exploration of active water cycling are reviewed here. Promising applications in cancer, cardiological, and neurological MRI are covered. This initial work employed paramagnetic Gd(III)chelate contrast agents [CAs]. However, the significant problems associated with in vivo CA use are also reviewed. A new analysis of water diffusion-weighted MRI [DWI] is presented. Preliminary results suggest a non-invasive way to measure the cell number density [ρ (cells/μL)], the mean cell volume [V (pL)], and the cellular NKA metabolic rate [^cMR_NKA(fmol(ATP)/s/cell)] with high spatial resolution. These crucial cell biology properties have not before been accessible invivo. Furthermore, initial findings indicate their absolute values can be determined.