Time-dependent diffusion coefficient as a probe of the permeability of the pore wall

The Diffusion Exchange Ratio (DEXR): A minimal sampling of diffusion exchange spectroscopy to probe exchange, restriction, and time-dependence

Water exchange is increasingly recognized as an important biological process that can affect the study of biological tissue using diffusion MR. Methods to measure exchange, however, remain immature as opposed to those used to characterize restriction, with no consensus on the optimal pulse sequence (s) or signal model (s). In general, the trend has been towards data-intensive fitting of highly parameterized models. We take the opposite approach and show that a judicious sub-sample of diffusion exchange spectroscopy (DEXSY) data can be used to robustly quantify exchange, as well as restriction, in a data-efficient manner. This sampling produces a ratio of two points per mixing time: (i) one point with equal diffusion weighting in both encoding periods, which gives maximal exchange contrast, and (ii) one point with the same total diffusion weighting in just the first encoding period, for normalization. We call this quotient the Diffusion EXchange Ratio (DEXR). Furthermore, we show that it can be used to probe time-dependent diffusion by estimating the velocity autocorrelation function (VACF) over intermediate to long times (∼2-500ms). We provide a comprehensive theoretical framework for the design of DEXR experiments in the case of static or constant gradients. Data from Monte Carlo simulations and experiments acquired in fixed and viable ex vivo neonatal mouse spinal cord using a permanent magnet system are presented to test and validate this approach. In viable spinal cord, we report the following apparent parameters from just 6 data points: τk=17±4ms, fNG=0.72±0.01, Reff=1.05±0.01μm, and κeff=0.19±0.04μm/ms, which correspond to the exchange time, restricted or non-Gaussian signal fraction, an effective spherical radius, and permeability, respectively. For the VACF, we report a long-time, power-law scaling with ≈t-2.4, which is approximately consistent with disordered domains in 3-D. Overall, the DEXR method is shown to be highly efficient, capable of providing valuable quantitative diffusion metrics using minimal MR data.

Limitations and generalizations of the first order kinetics reaction expression for modeling diffusion-driven exchange: Implications on NMR exchange measurements

The study and modeling of water exchange in complex media using different applications of diffusion and relaxation magnetic resonance (MR) have been of interest in recent years. Most models attempt to describe this process using a first order kinetics expression, which is appropriate to describe chemical exchange; however, it may not be suitable to describe diffusion-driven exchange since it has no direct relationship to diffusion dynamics of water molecules. In this paper, these limitations are addressed through a more general exchange expression that does consider such important properties. This exchange fraction expression features a multi-exponential recovery at short times and a mono-exponential decay at long times, both of which are not captured by the first order kinetics expression. Furthermore, simplified exchange expressions containing partial information of the analyzed system's diffusion and relaxation processes and geometry are proposed, which can potentially be employed in already established estimation protocols. Finally, exchange fractions estimated from simulated MR data and derived here were compared, showing qualitative similarities but quantitative differences, suggesting that the features of the derived exchange fraction in this paper can be partially recovered by employing an existing estimation framework.

A single-shot measurement of time-dependent diffusion over sub-millisecond timescales using static field gradient NMR

Time-dependent diffusion behavior is probed over sub-millisecond timescales in a single shot using a nuclear magnetic resonance static gradient time-incremented echo train acquisition (SG-TIETA) framework. The method extends the Carr-Purcell-Meiboom-Gill cycle under a static field gradient by discretely incrementing the π-pulse spacings to simultaneously avoid off-resonance effects and probe a range of timescales (50-500 µs). Pulse spacings are optimized based on a derived ruleset. The remaining effects of pulse inaccuracy are examined and found to be consistent across pure liquids of different diffusivities: water, decane, and octanol-1. A pulse accuracy correction is developed. Instantaneous diffusivity, Dinst(t), curves (i.e., half of the time derivative of the mean-squared displacement in the gradient direction) are recovered from pulse accuracy-corrected SG-TIETA decays using a model-free log-linear least squares inversion method validated by Monte Carlo simulations. A signal-averaged 1-min experiment is described. A flat Dinst(t) is measured on pure dodecamethylcyclohexasiloxane, whereas decreasing Dinst(t) is measured on yeast suspensions, consistent with the expected short-time Dinst(t) behavior for confining microstructural barriers on the order of micrometers.

From diffusion in compartmentalized media to non-Gaussian random walks

In this work we establish a link between two different phenomena that were studied in a large and growing number of biological, composite and soft media: the diffusion in compartmentalized environment and the non-Gaussian diffusion that exhibits linear or power-law growth of the mean square displacement joined by the exponential shape of the positional probability density. We explore a microscopic model that gives rise to transient confinement, similar to the one observed for hop-diffusion on top of a cellular membrane. The compartmentalization of the media is achieved by introducing randomly placed, identical barriers. Using this model of a heterogeneous medium we derive a general class of random walks with simple jump rules that are dictated by the geometry of the compartments. Exponential decay of positional probability density is observed and we also quantify the significant decrease of the long time diffusion constant. Our results suggest that the observed exponential decay is a general feature of the transient regime in compartmentalized media.

Slow polymer diffusion on brush-patterned surfaces in aqueous solution

A model system for the investigation of diffusional transport in compartmentalized nanosystems is described. Arrays of "corrals" enclosed within poly[oligo(ethylene glycol)methyl ether methacrylate] (POEGMA) "walls" were fabricated using double-exposure interferometric lithography to deprotect aminosilane films protected by a nitrophenyl group. In exposed regions, removal of the nitrophenyl group enabled attachment of an initiator for the atom-transfer radical polymerization of end-grafted POEGMA (brushes). Diffusion coefficients for poly(ethylene glycol) in these corrals were obtained by fluorescence correlation spectroscopy. Two modes of surface diffusion were observed: one which is similar to diffusion on the unpatterned surface and a very slow mode of surface diffusion that becomes increasingly important as confinement increases. Diffusion within the POEGMA brushes does not significantly contribute to the results.

Tissue microstructure features derived from anomalous diffusion measurements in magnetic resonance imaging

OBJECTIVES

Tissue microstructure features, namely axon radius and volume fraction, provide important information on the function of white matter pathways. These parameters vary on the scale much smaller than imaging voxels (microscale) yet influence the magnetic resonance imaging diffusion signal at the image voxel scale (macroscale) in an anomalous manner. Researchers have already mapped anomalous diffusion parameters from magnetic resonance imaging data, but macroscopic variations have not been related to microscale influences. With the aid of a tissue model, we aimed to connect anomalous diffusion parameters to axon radius and volume fraction using diffusion-weighted magnetic resonance imaging measurements.

EXPERIMENTAL DESIGN

An ex vivo human brain experiment was performed to directly validate axon radius and volume fraction measurements in the human brain. These findings were validated using electron microscopy. Additionally, we performed an in vivo study on nine healthy participants to map axon radius and volume fraction along different regions of the corpus callosum projecting into various cortical areas identified using tractography.

PRINCIPAL OBSERVATIONS

We found a clear relationship between anomalous diffusion parameters and axon radius and volume fraction. We were also able to map accurately the trend in axon radius along the corpus callosum, and in vivo findings resembled the low-high-low-high behaviour in axon radius demonstrated previously.

CONCLUSIONS

Axon radius and volume fraction measurements can potentially be used in brain connectivity studies and to understand the implications of white matter structure in brain diseases and disorders. Hum Brain Mapp 38:1068-1081, 2017. © 2016 Wiley Periodicals, Inc.

Spatially dependent diffusion coefficient as a model for pH sensitive microgel particles in microchannels

Solute transport and intermixing in microfluidic devices is strongly dependent on diffusional processes. Brownian Dynamics simulations of pressure-driven flow of model microgel particles in microchannels have been carried out to explore these processes and the factors that influence them. The effects of a pH-field that induces a spatial dependence of particle size and consequently the self-diffusion coefficient and system thermodynamic state were focused on. Simulations were carried out in 1D to represent some of the cross flow dependencies, and in 2D and 3D to include the effects of flow and particle concentration, with typical stripe-like diffusion coefficient spatial variations. In 1D, the mean square displacement and particle displacement probability distribution function agreed well with an analytically solvable model consisting of infinitely repulsive walls and a discontinuous pH-profile in the middle of the channel. Skew category Brownian motion and non-Gaussian dynamics were observed, which follows from correlations of step lengths in the system, and can be considered to be an example of so-called "diffusing diffusivity." In Poiseuille flow simulations, the particles accumulated in regions of larger diffusivity and the largest particle concentration throughput was found when this region was in the middle of the channel. The trends in the calculated cross-channel diffusional behavior were found to be very similar in 2D and 3D.

Multiple-echo diffusion tensor acquisition technique (MEDITATE) on a 3T clinical scanner

This article describes the concepts and implementation of an MRI method, the multiple-echo diffusion tensor acquisition technique (MEDITATE), which is capable of acquiring apparent diffusion tensor maps in two scans on a 3T clinical scanner. In each MEDITATE scan, a set of RF pulses generates multiple echoes, the amplitudes of which are diffusion weighted in both magnitude and direction by a pattern of diffusion gradients. As a result, two scans acquired with different diffusion weighting strengths suffice for accurate estimation of diffusion tensor imaging (DTI) parameters. The MEDITATE variation presented here expands previous MEDITATE approaches to adapt to the clinical scanner platform, such as exploiting longitudinal magnetization storage to reduce T2 weighting. Fully segmented multi-shot Cartesian encoding is used for image encoding. MEDITATE was tested on isotropic (agar gel), anisotropic diffusion phantoms (asparagus) and in vivo skeletal muscle in healthy volunteers with cardiac gating. Comparisons of accuracy were performed with standard twice-refocused spin echo (TRSE) DTI in each case and good quantitative agreement was found between diffusion eigenvalues, mean diffusivity and fractional anisotropy derived from TRSE DTI and from the MEDITATE sequence. Orientation patterns were correctly reproduced in both isotropic and anisotropic phantoms, and approximately for in vivo imaging. This illustrates that the MEDITATE method of compressed diffusion encoding is feasible on the clinical scanner platform. With future development and employment of appropriate view-sharing image encoding, this technique may be used in clinical applications requiring time-sensitive acquisition of DTI parameters such as dynamical DTI in muscle.

Accounting for adsorption and desorption in lattice Boltzmann simulations

We report a Lattice-Boltzmann scheme that accounts for adsorption and desorption in the calculation of mesoscale dynamical properties of tracers in media of arbitrary complexity. Lattice Boltzmann simulations made it possible to solve numerically the coupled Navier-Stokes equations of fluid dynamics and Nernst-Planck equations of electrokinetics in complex, heterogeneous media. With the moment propagation scheme, it became possible to extract the effective diffusion and dispersion coefficients of tracers, or solutes, of any charge, e.g., in porous media. Nevertheless, the dynamical properties of tracers depend on the tracer-surface affinity, which is not purely electrostatic and also includes a species-specific contribution. In order to capture this important feature, we introduce specific adsorption and desorption processes in a lattice Boltzmann scheme through a modified moment propagation algorithm, in which tracers may adsorb and desorb from surfaces through kinetic reaction rates. The method is validated on exact results for pure diffusion and diffusion-advection in Poiseuille flows in a simple geometry. We finally illustrate the importance of taking such processes into account in the time-dependent diffusion coefficient in a more complex porous medium.

Random walk with barriers

Restrictions to molecular motion by barriers (membranes) are ubiquitous in porous media, composite materials and biological tissues. A major challenge is to characterize the microstructure of a material or an organism nondestructively using a bulk transport measurement. Here we demonstrate how the long-range structural correlations introduced by permeable membranes give rise to distinct features of transport. We consider Brownian motion restricted by randomly placed and oriented membranes (d - 1 dimensional planes in d dimensions) and focus on the disorder-averaged diffusion propagator using a scattering approach. The renormalization group solution reveals a scaling behavior of the diffusion coefficient for large times, with a characteristically slow inverse square root time dependence for any d. Its origin lies in the strong structural fluctuations introduced by the spatially extended random restrictions, representing a novel universality class of the structural disorder. Our results agree well with Monte Carlo simulations in two dimensions. They can be used to identify permeable barriers as restrictions to transport, and to quantify their permeability and surface area.

Pulsed-gradient spin-echo monitoring of restricted diffusion in multilayered structures

A general mathematical basis is developed for computation of the pulsed-gradient spin-echo signal attenuated due to restricted diffusion in multilayered structures (e.g., multiple slabs, cylindrical or spherical shells). Individual layers are characterized by (different) diffusion coefficients and relaxation times, while boundaries between adjacent layers are characterized by (different) permeabilities. Arbitrary temporal profile of the applied magnetic field can be incorporated. The signal is represented in a compact matrix form and the explicit analytical formulas for the elements of the underlying matrices are derived. The implemented algorithm is faster and much more accurate than classical techniques such as Monte Carlo simulations or numerical resolutions of the Bloch-Torrey equation. The algorithm can be applied for studying restricted diffusion in biological systems which exhibit a multilayered structure such as composite tissues, axons and living cells.

Quantitative permeability imaging of plant tissues

A method for mapping tissue permeability based on time-dependent diffusion measurements is presented. A pulsed field gradient sequence to measure the diffusion encoding time dependence of the diffusion coefficients based on the detection of stimulated spin echoes to enable long diffusion times is combined with a turbo spin echo sequence for fast NMR imaging (MRI). A fitting function is suggested to describe the time dependence of the apparent diffusion constant in porous (bio-)materials, even if the time range of the apparent diffusion coefficient is limited due to relaxation of the magnetization. The method is demonstrated by characterizing anisotropic cell dimensions and permeability on a subpixel level of different tissues of a carrot (Daucus carota) taproot in the radial and axial directions.

Spectral characterization of diffusion with chemical shift resolution: highly concentrated water-in-oil emulsion

We present a modulated gradient spin-echo method, which uses a train of sinusoidally shaped gradient pulses separated by 180 degrees radio-frequency (RF) pulses. The RF pulses efficiently refocus chemical shifts and de-phasing due to susceptibility differences, resulting in undistorted, high-resolution diffusion weighted spectra. This allows for the simultaneous spectral characterization of the diffusion of several molecular species with different chemical shifts. The technique is robust against susceptibility artifacts, field inhomogeneity and imperfections in the gradient generating equipment. The feasibility of the technique is demonstrated by measuring the diffusion of water, oil, and water-soluble salt in a highly concentrated water-in-oil emulsion. The diffusion of water and salt reveal precise information about the droplet size distribution below the mum-range. Common droplet size distribution explains both the data for water with finite long-range diffusion and the data for salt with negligible long-range diffusion. The results of water diffusion show that the technique is efficient in deconvolving the effects of molecular exchange between droplets and restricted diffusion within droplets. The effects of water exchange suggest that droplets of different sizes are uniformly distributed within the sample.

Spatial fluorescence cross-correlation spectroscopy by means of a spatial light modulator

Spatial fluorescence cross-correlation spectroscopy is a rarely investigated version of fluorescence correlation spectroscopy, in which the fluorescence signals from different observation volumes are cross-correlated. In the reported experiments, two observation volumes, typically shifted by a few microm, are produced, with a spatial light modulator and two adjustable pinholes. We illustrated the feasibility and potentiality of this technique by: i) measuring molecular flows, in the range 0.2-1.5 microm/ms, of solutions seeded with fluorescent nanobeads or rhodamine molecules (simulating active transport phenomenons); ii) investigating the permeability of the phospholipidic membrane of giant unilamellar vesicles versus hydrophilic or hydrophobic molecules (in that case the laser spots were set on both sides of the membrane). Theoretical descriptions are proposed together with a discussion about fluorescence-correlation-spectroscopy-based, alternative methods.

Probing restrictive diffusion dynamics at short time scales

Diffusion imaging gradients serve to spectrally filter the temporally evolving diffusion tensor. In this framework, the design of diffusion sensitizing gradients is reduced to the problem of adequately sampling q-space in the spectral domain. The practical limitations imposed by the requirement for delta-function type diffusion-sensitizing gradients to adequately sample q-space, can be relaxed if these impulse gradients are replaced with chirped oscillatory gradients. It is well known that in many systems of interest, dispersion of velocity will itself produce a peak in the velocity correlation function near w=0, while restricted diffusion will manifest itself in the dispersion spectrum at higher frequencies. In this paper, chirped diffusion-sensitizing gradients are proposed and analytically shown to yield an efficient sampling of q-space in a manner that asymptotically approaches that using delta-function diffusion-sensitizing gradient. The challenge is the consequent reduction in diffusion sensitivity as one probes higher frequency dynamics. This problem is addressed by restricting the gradient power to a spectral bandwidth corresponding to the diffusion spectral range of the underlying restrictive geometry. Simultaneous imaging of diffusion and flow at microscopic resolution and at temporally resolvable diffusion time scales thus becomes possible in vivo. Simulations and experiments validate the proposed approach.

Displacement power spectrum measurement by CPMG in constant gradient

The modulation of spin phase produced by Carr-Purcell-Meiboom-Gill (CPMG) sequence in combination with constant magnetic field gradient is appropriate to probe the displacement power spectrum (DPS). The spin-echo attenuation is directly proportional to the DPS value at the applied modulation frequency. Relaxation and selective excitation effects can be factored out while probing the DPS. The modulation frequency is adjusted by varying the pulse separation time while the gradient strength and the time of acquisition are kept constant. In designing the experiment gradient strength limitations, imposed by off-resonance effects, as well as limitations arising from using Gaussian phase approximation must be considered. An effective experimental strategy is presented, supported by experimental results for free and restricted diffusion.

Acquisition of pure-phase diffusion spectra using oscillating-gradient spin echo

Oscillating-gradient spin echo (OGSE) diffusion experiments have long been used to measure the short-time apparent diffusion coefficient, D(app)(t), in the presence of restricted diffusion, as well as the spectrum of the slow-motion velocity autocorrelation function. In this work, we focus on two previously unexplored aspects of OGSE experiments: convection compensation and acquisition of pure-phase diffusion spectra in the presence of homonuclear scalar couplings. We demonstrate that convection compensation afforded by single-echo OGSE compares well with that in double-echo convection-compensated PGSE experiments. We also show that, in the presence of homonuclear scalar couplings, setting the OGSE echo time to 1/2J enables acquisition of pure-phase diffusion spectra and yields more reliable D estimates than mixed-phase PGSE or OGSE spectra. Pure-phase OGSE acquisition is also compatible with measurements of the apparent diffusion coefficient at an arbitrary diffusion time. These features of OGSE can be valuable in diffusion measurements of scalar-coupled small-molecule probes in cellular and other heterogeneous systems.