Next Article in Journal
On the Geometric Description of Nonlinear Elasticity via an Energy Approach Using Barycentric Coordinates
Next Article in Special Issue
The Mathematics of Quasi-Diffusion Magnetic Resonance Imaging
Previous Article in Journal
Bounds for the Energy of Graphs
Previous Article in Special Issue
Fractional Order Magnetic Resonance Fingerprinting in the Human Cerebral Cortex
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diffusion in Sephadex Gel Structures: Time Dependency Revealed by Multi-Sequence Acquisition over a Broad Diffusion Time Range

1
Center for Magnetic Resonance Research, University of Illinois at Chicago, Chicago, IL 60612, USA
2
Department of Bioengineering, University of Illinois at Chicago, Chicago, IL 60607, USA
3
Research Resource Center, University of Illinois at Chicago, Chicago, IL 60612, USA
4
Department of Radiology, Northwestern University, Evanston, IL 60611, USA
5
Department of Radiology, University of Illinois College of Medicine, Chicago, IL 60612, USA
6
Department of Neurosurgery, University of Illinois College of Medicine, Chicago, IL 60612, USA
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(14), 1688; https://doi.org/10.3390/math9141688
Submission received: 7 June 2021 / Revised: 8 July 2021 / Accepted: 10 July 2021 / Published: 19 July 2021
(This article belongs to the Special Issue Fractional Calculus in Magnetic Resonance)

Abstract

:
It has been increasingly reported that in biological tissues diffusion-weighted MRI signal attenuation deviates from mono-exponential decay, especially at high b-values. A number of diffusion models have been proposed to characterize this non-Gaussian diffusion behavior. One of these models is the continuous-time random-walk (CTRW) model, which introduces two new parameters: a fractional order time derivative α and a fractional order spatial derivative β. These new parameters have been linked to intravoxel diffusion heterogeneities in time and space, respectively, and are believed to depend on diffusion times. Studies on this time dependency are limited, largely because the diffusion time cannot vary over a board range in a conventional spin-echo echo-planar imaging sequence due to the accompanying T2 decays. In this study, we investigated the time-dependency of the CTRW model in Sephadex gel phantoms across a broad diffusion time range by employing oscillating-gradient spin-echo, pulsed-gradient spin-echo, and pulsed-gradient stimulated echo sequences. We also performed Monte Carlo simulations to help understand our experimental results. It was observed that the diffusion process fell into the Gaussian regime at extremely short diffusion times whereas it exhibited a strong time dependency in the CTRW parameters at longer diffusion times.

1. Introduction

Using water diffusion as a probe, diffusion-weighted MRI (DW-MRI) has become a promising technique to reveal the underlying micrometer-scale structural properties in millimeter-resolution MR images [1,2]. In DW-MRI, two diffusion gradient lobes are employed to dephase and rephase spins, respectively. The displacement of water molecules is quantified over a given time period, known as effective diffusion time (Δeff), which is constrained by the separation and duration of diffusion gradient lobes. For water molecules that diffuse during Δeff, the varying degree of spatial dislocation results in a phase dispersion (Φ) of the magnetization. The probabilistic distribution function (PDF) of the net displacement of diffusing water molecules is related to a probability distribution of Φ, which leads to signal attenuation in an MRI measurement [3,4].
It is widely accepted that the probability distribution of molecular displacement is Gaussian in an isotropic, homogeneous, and unrestricted medium. In that case, the second moment of the distribution, or mean squared displacement (MSD), scales linearly with diffusion time: 〈x2〉 ~ t [5]. In the presence of restricting barriers in complex materials, however, the probability distribution of molecular displacement no longer follows Gaussian distribution. The MSD in the non-Gaussian diffusion case can be characterized as a function of intrinsic diffusion coefficient, restrictive geometry, and diffusion time [6]. One way of characterizing the non-Gaussian diffusion behavior is to employ the continuous time random walk (CTRW) theory, in which the MSD can be expressed by a composite power law as: 〈x2〉 ~ t2α/β, where α and β are the fractional order time and space derivatives, respectively, in fractionalized Fick’s second law [7,8]. This generalized description enables the CTRW model to provide a more realistic description of the complex diffusion pattern in biological tissues [8].
In parallel to the development of non-Gaussian diffusion models, it has been increasingly recognized that diffusion parameters derived from various diffusion models exhibit dependence on diffusion time. Pyatigorskaya et al. [9] and Aggarwal et al. [10] observed substantial time dependency in diffusion kurtosis imaging (DKI) in the mouse brain while noticeable time dependency of intravoxel incoherent motion (IVIM) model was also reported by Wu et al. [11] in a flow phantom and mouse brain. Varying diffusion time enables exploration of the interaction between diffusing water molecules and the surrounding environment at different spatial scales [12], providing a new degree of freedom to estimate parameters that are related to the underlying tissue microstructures [12,13]. Conventional spin-echo-based DWI sequences have limited ability to vary diffusion time. On one hand, a longer diffusion time results in substantial increase in echo time, leading to signal loss due to T2 decay. On the other hand, a shorter diffusion time reduces the b-value, leading to inadequate diffusion-weighting. Therefore, investigation of diffusion time dependency in DWI-MRI over a broad range requires alternative pulse sequences to lengthen or shorten the effective diffusion time. Although previous studies investigated the time dependency of the CTRW parameters in the intermediate to long diffusion time range [8,14,15], the time dependency of the CTRW model at short diffusion time remains unexplored. Furthermore, the consistency of the CTRW parameters under different pulse sequences has not been well studied.
Sephadex gels are structurally heterogeneous swollen polymers with relatively uniform pore size [16], forming an ideal test bed to mimic the diffusion environment in complex biological tissues. The multi-compartment Sephadex beads provide different physical environments, such as restricted diffusion environment within the bead, free liquid space outside the bead, and open pores for water molecules to diffuse through. In this complex multi-compartment environment, water molecules move freely at short diffusion times while they experience hinderance or restriction when interacting with polymer structures at longer diffusion times. Sephadex gel phantoms have been widely used for validating advanced diffusion models, such as the CTRW model’s predecessor–fractional order calculus model [7,14,17].
In this study, we investigated the diffusion time dependency of the CTRW model by employing oscillating-gradient spin-echo (OGSE), pulsed-gradient spin-echo (PGSE), and pulsed-gradient stimulated echo (PGSTE) pulse sequences. Collectively, these sequences enabled investigation of a broad range of diffusion times, spanning from short-, intermediate-, to long-time regime. Two series of Sephadex gels, each with different pore size or bead diameters, were selected as the experimental material to mimic the tissue environment. In addition, Monte Carlo simulations were performed to help understand the experimental results.

2. Theory

According to the simple random walk (RW) theory, the one-dimensional Brownian motion of a diffusing particle in a homogeneous and isotropic environment can be described by a second-order partial differential equation,
P ( x , t ) t = D 2 P ( x , t ) | x | 2 ,
where P(x,t) is one-dimensional Brownian motion of a diffusing particle and D is the diffusion coefficient. The solution to Equation (1) yields Gaussian distribution of the displacement where the MSD is proportional to diffusion time t, 〈x2(t)〉 ~ t [8].
In the context of the CTRW theory, where jump length and jump waiting time follow asymptotic power law distributions, the one-dimensional anomalous motion of a diffusing particle in a heterogenous environment can be described with a dual space-time fractional order diffusion equation of the form [7,18,19,20],
D 0 C t α P ( x , t ) = D α , β β P ( x , t ) | x | β ,
where Dα,β is the anomalous diffusion coefficient (in mmβ/sα), D 0 C t α is the αth (0 < α ≤ 1) fractional order time derivative in the Caputo form, given as [21]:
D 0 C t α f ( t ) = 1 Γ ( n α ) 0 t ( t τ ) n α 1 ( d d τ ) n f ( τ ) d τ .
β P ( x , t ) | x | β in Equation (2) is the βth (0 < β ≤ 2) fractional order space derivative in the Riesz form. For n−1 < β < n and the finite interval 0 ≤ xL, the Riesz fractional operator is defined as [22]:
β u ( x , t ) | x | β = 1 2 cos ( π β 2 ) [ D 0   x β + D x   L β ] u ( x , t ) ,
where
D 0   x β u ( x , t ) = 1 Γ ( n β ) n x n 0 x u ( ξ , t ) d ξ ( x ξ ) β + 1 n
D x   L β u ( x , t ) = 1 Γ ( n β ) n x n x L u ( ξ , t ) d ξ ( x ξ ) β + 1 n .
With the representation in Equation (2), the MSD can be represented by 〈x2(t)〉 ~ t(2α/β). When α = 1 and β = 2, this formalism is reduced to the classical Gaussian expression. In comparison, when 2α > β or 2α < β, the anomalous diffusion process is referred to as super-diffusive or sub-diffusive [7,17] respectively, and when 2α = β, the non-Gaussian dynamics is described as quasi-diffusion [23].
For a Stejskal–Tanner diffusion gradient pulse, the solution to Equation (2) can be described as:
S S 0 = P ( q , Δ e f f ) = E α ( D α , β q β Δ e f f α ) .
In Equation (7), S0 is the signal intensity without diffusion weighting and S is the signal intensity at q and Δeff, where q = γGdiffδ and Δeff = Δ − δ/3 in which γ is the gyrometric ratio, Gdiff is the diffusion gradient amplitude, δ is the diffusion gradient pulse width, and Δ is the diffusion gradient separation. Eα is a single-parameter Mittag–Leffler function [8]. For other diffusion gradient waveforms, expressions analogous to Equation (7) can be derived in reference of the methods described in [24,25].

3. Methods

3.1. Sephadex Gel Phantom Preparation

In this study, we used two series of Sephadex gels (GE Healthcare) which were characterized by two numbers (e.g., G25–50), where the first number indicates the macromolecular exclusion limit (in kDaltons, positively correlated with internal pore size) and the second denotes the maximum dry bead diameter (in microns). The first series of gels (G25–50, G50–50, G75–50) had the same dry bead diameter of 50 microns, but with increased macromolecular exclusion limit. The second series of gels (G50–50, G50–80, G50–150) had the same internal pore size of 50 kDaltons, but increased bead diameter. The first series of gels were designed to mimic varying microstructure permeability while the second series to simulate varying microstructural scale.
Sephadex gel phantoms were prepared by gently pouring excess distilled water into the dry power gel in a cylindrical test tube (inner diameter = 13.5 mm) at room temperature; and mixed evenly by using a vortex shaker. The slurry was allowed to settle under the influence of gravity; and the residual water was removed by pipette before sealing.

3.2. Data Acquisition

The experiments were performed on an Agilent 9.4 T small animal MRI scanner with a maximum gradient of 1000 mT/m. Prepared Sephadex gels were scanned on the scanner at the room temperature of 22 °C. As illustrated in Figure 1, three different DW pulse sequences were employed to investigate the diffusion time dependency of the CTRW model parameters across a broad range of diffusion times.
(I) Customized cosine-trapezoid OGSE sequence: OGSE sequence enables a short effective diffusion time by periodically varying the polarity of diffusion encoding gradients. Cosine-trapezoid OGSE waveforms start and end with a quarter-period lobe. Δeff and b-value of the cosine-trapezoid OGSE sequence are given by [26]:
Δeff = δ/(3N)
b = γ 2 G d i f f 2 δ 3 6 N 2
where δ is the total waveform duration and N is the number of half oscillation periods (Figure 1a). Three OGSE acquisitions were performed with a constant δ of 30 ms while N was set to 6, 4, and 2, resulting in Δeff values of 1.67, 2.5, and 5 ms, respectively.
(II) PGSE sequence: The effective diffusion time, Δeff, under a Stejskal–Tanner diffusion sensitizing gradient pair in a PGSE sequence is given by [3]:
Δeff = Δ − δ/3
b = γ 2 G d i f f 2 δ 2 Δ e f f
The PGSE experiments were performed with Δ = 11 ms and 35 ms, and δ = 2.5 ms, resulting in Δeff values of 10.17 ms and 34.17 ms, respectively (Figure 1b).
(III) PGSTE sequence: With the same Stejskal–Tanner diffusion sensitizing gradient pair as in PGSE, a PGSTE sequence achieves long diffusion time by taking advantage of the slower T1 relaxation rate during mixing time (TM). The DW-MRI data were acquired with an identical Δ as in the PGSE sequence (35 ms), and three longer Δ values (80 ms, 100 ms, and 150 ms). The corresponding Δeff values were 34.17, 59.17, 99.17, and 149.17 ms, respectively (Figure 1c).
At each diffusion time, DW images with 11 b-values (0, 25, 100, 225, 400, 625, 900, 1225, 1600, 2025, and 2500 s/mm2) were acquired from the Sephadex gel phantoms by varying Gdiff. The other imaging parameters, TR (4000 ms) and TE (75 ms), diffusion gradient direction = R/L, FOV = 36 × 36 mm2, acquisition matrix = 32 × 32, slice thickness = 2 mm and number of repetitions (NEX = 4), were kept the same in all sequences.

3.3. Data Analysis

The DW images acquired by OGSE, PGSE, and PGSTE were first normalized by dividing DW signal, S, at each b-value by S0. The CTRW model in Equation (7) was fit to the DW images voxel-by-voxel by using an iterative non-linear Levenberg–Marquardt algorithm in MATLAB. To improve the fitting stability and alleviate the degeneracy, Dα,β was first estimated by a mono-exponential model at lower b-values, followed by a simultaneous estimation of other parameters with appropriate constraints (0 < α ≤ 1 and 0 < β ≤ 2) by using all b-values [27]. Measurements were made from each quantitative parameter map (Dα,β, α, and β) by computing the mean value over a ~16 mm2 region-of-interest (ROI) within each vial of the Sephadex gel. Representative DW images and ROIs are shown in Figure 2.

3.4. Monte Carlo Simulations

Monte Carlo simulations of the time-dependent MR signals were performed with random walkers implemented by using the Monte Carlo Diffusion Simulator of Camino Diffusion MRI Toolkit (UCL) [28,29]. In our simulations, we modeled Sephadex beads with permeable pores as the perpendicular sections of parallelly packed, non-overlapping cylinders with permeable membrane. The permeability was defined as a fixed probability of a random walker stepping through the membrane. Different Sephadex bead sizes and macromolecular exclusion limits were simulated by adjusting the circle radius and permeability as detailed below.
Random walkers were randomly seeded inside and between the hexagonally packed 3D parallel cylinders with permeable membranes; and updated the positions following the rules described in a previous study [28]. Phase change of each random walker was calculated under the custom-specific diffusion gradients. The synthetic DWI signals were then generated by summing the contributions from all random walkers at the echo time. To investigate the effect of varying microstructure scale and membrane permeability on the anomalous diffusion signal behavior and its time dependency, two sets of simulation data were generated to simulate the two Sephadex gel series in the MRI experiments. In the first simulation dataset, substrates with different permeability (p = 0.1%, 0.2%, and 0.4%) were chosen with fixed cylinder radius (r = 8 µm). In the second simulation dataset, substrates with different cylinder radii (r = 6 µm, 7 µm, and 8 µm) were selected with fixed membrane permeability (p = 0.2%). Synthetic DWI signals were generated with diffusion gradients that were designed to be perpendicular to the cylinder long axes. For each substrate, three different synthetic OGSE signals were simulated with the oscillating diffusion gradients at three Δeff values (3.33, 5, and 10 ms). Eleven different synthetic PGSE/PGSTE signals were simulated with Stejskal–Tanner diffusion gradients at eleven Δeff (25, 30, 35, 40, 45, 50, 60, 70, 80, 90, and 100 ms).
All simulations were performed on an 8-core Intel i7–2600 CPU with 100,000 random walkers and 20,000 time-steps with intrinsic diffusivity of 2 × 10−3 mm2/s, intracellular volume ratio of 0.5 and 7 b-values (0, 200, 500, 1000, 1500, 3000, and 6000 s/mm2). The normalized simulated signal intensities over all b-values were fit to Equation (7) with the same fitting procedure as for the experimental data.

4. Results

The CTRW parameter, Dα,β, obtained from the first Sephadex gel series was plotted as a function of Δeff in Figure 3. Two trends were observed. First, a downward trend reaching a plateau was seen for all gels, suggesting increased hinderance at longer diffusion times. Second, the gels with larger pore sizes (G50–50 and G75–50 in Figure 3b,c) exhibited higher Dα,β values at all diffusion times.
In Figure 4 and Figure 5, the fractional order time and space derivatives, α and β, were plotted against Δeff, respectively. α and β exhibited a similar trend to each other. For all the gels, as Δeff decreased to 0, α and β values approached to 1 and 2, respectively, indicating that the diffusion signal behavior approaches to the Gaussian regime in the limit of short diffusion times. The gels with larger pore sizes exhibited higher α and β values (G50–50 and G75–50 in Figure 4a,b and Figure 5a,b), suggesting less deviation from Gaussian diffusion dynamics.
In Figure 6, Dα,β is plotted as a function of Δeff for the second Sephadex gel series, G50–50, G50–80, and G50–150. Similar to the first Sephadex gel series, Dα,β followed a downward trend reaching a plateau for all the Sephadex gels. G50–50, G50–80, and G50–150 in Figure 6a–c show similar Dα,β values at short diffusion times. However, at longer diffusion times, the gels with larger bead sizes (G50–80 and G50–150 in Figure 6b,c) exhibited higher Dα,β values, similar to what is shown in Figure 3.
In Figure 7 and Figure 8, α and β are plotted against Δeff for gels in the second Sephadex series (G50–50, G50–80, and G50–150). Similar to the first Sephadex gel series, α and β showed a decreasing trend against Δeff in all gels. As the dry bead size increased, higher α and β values were observed in general. 2α/β < 1 was observed in all Sephadex gels, indicating the diffusion dynamics fell into sub-diffusion regime.
Plots in Figure 9 and Figure 10 show the time-dependent changes observed in the CTRW parameters in the Monte Carlo simulations. In the first simulation dataset (fixed r and varying p of 0.1%, 0.2%, and 0.4%), the CTRW parameters, Dα,β (Figure 9a), α (Figure 9b), and β (Figure 9c), yielded higher values in the data with higher permeability. In the second simulation dataset (fixed p and varying r of 6, 7, and 8 µm), the simulated data with a larger cylinder radius yielded higher Dα,β (Figure 10a), α (Figure 10b), and β (Figure 10c). The simulation results exhibited a good agreement with experimental results. In both simulations, Dα,β, α, and β exhibited a monotonically decreasing trend.

5. Discussion

By employing OGSE, PGSE, and PGSTE sequences to span a broad range of diffusion times, we investigated the diffusion time dependency of the CTRW parameters in Sephadex gel phantoms and correlated our results with Monte Carlo simulations. We observed a monotonic decrease in the CTRW parameters, Dα,β, α, and β, as the diffusion time increased. Our Monte Carlo simulations exhibited a similar trend with the experimental results. To the best of our knowledge, this is the first study which investigates the time dependency behavior of the CTRW model over a wide range of diffusion times using a multi-sequence acquisition scheme.
In the classical mono-exponential model, where water molecules diffuse freely without hinderance and restriction, the MR signal attenuation function is concisely characterized by a commonly used exponential function, exp(-bD). In the CTRW model, α and β quantitatively describe the deviation of diffusion dynamics from the mono-exponential decay [30]. At short diffusion times where MSD is much smaller than the obstacle scales, the water molecules can diffuse freely in all directions, leading to a process that follows Gaussian dynamics. As expected, at this short diffusion time extreme, α and β values approached to 1 and 2, respectively, while Dα,β converged to the diffusion coefficient of pure water, D0. These outcomes were clearly observed in our experimental data and confirmed in Monte Carlo simulations. In contrast, the long diffusion time provides water molecules a greater opportunity to explore the heterogeneity of the surrounding environment, resulting in reduced α and β. The increased hinderance and restriction experienced by water molecules at long diffusion times yielded reduced Dα,β, α, and β values, which allows us to infer information on microstructures and micro-environment.
The lower Dα,β values observed in the gels with smaller pore sizes in the first Sephadex gel series (Figure 3a,b) is consistent with the general belief that diffusion coefficient is lower in materials with increased micro-structural barriers [31]. Sephadex gels with larger dry bead sizes exhibited similar Dα,β at low diffusion times, but higher Dα,β at long diffusion times (Figure 6b,c), suggesting that the influence of microstructure scale on diffusion dynamics is more evident at longer diffusion times.
The CTRW model incorporates temporal and spatial diffusion heterogeneities through fractional order time and space parameters, α and β, respectively. The α parameter describes the likelihood of water molecules to be “trapped” or “released” in complex materials, which reflects temporal heterogeneity of diffusion process. The β parameter is mathematically equivalent to the heterogeneity parameter in the stretched-exponential model [32], which has been linked with the non-Gaussian distribution of diffusion displacement and shown to be related to intravoxel spatial heterogeneity [18]. In our experimental and simulation results, the lower permeability and smaller microstructural scale of the structural barriers led to a higher likelihood for water molecules to interact with the surrounding structures. This explains the lower α and β values observed in Sephadex gels with lower macromolecular exclusion limit (Figure 4a and Figure 5a) and smaller bead diameter (Figure 7a and Figure 8a). The experimental observations were further reinforced by the simulation results with lower permeability (Figure 9a) and smaller cylinder radius (Figure 10a).
In this study, some discontinuities were observed in the parameter values when different pulse sequences were employed. For example, Dα,β values obtained from the experiments performed with the PGSTE were higher than those observed by the PGSE at Δeff = 35 ms (Figure 3 and Figure 6). Also, α and β values estimated from the data acquired with the PGSE at Δ = 11 ms were higher than those with the OGSE at Δeff = 5 ms in the first Sephadex series, as shown in Figure 4 and Figure 5. The root cause of the discontinuities is unknown and requires further investigation. Nevertheless, the overall monotonic trend across a broad diffusion time range is consistent in all CTRW parameters. Two Sephadex gel series with varying pore size or bead diameter exhibited the same monotonic trend, which was consistent with the trends revealed by the Monte Carlo simulations.
The accuracy of our simulation depends upon the Monte Carlo Diffusion Simulator of Camino Diffusion MRI Toolkit. The assumption is that this simulator is capable of simulating diverse diffusion processes across a broad range of environments (i.e., substrates), from simple to exceedingly complex. In our simulations, the varying radius and permeability were employed to mimic the varying dry bead diameter and macromolecular exclusion limit of the Sephadex gels, respectively. Although we did not attempt to explicitly evaluate the accuracy and precision of the simulations in this study, a previous study [28] has illustrated the accuracy of a similar simulation approach. Furthermore, the accuracy was likely enhanced by the large number of random walkers (100,000) and time steps (20,000) employed in our study.
Varying diffusion time provides a viable approach to investigating the length scale of tissue microstructures using water diffusion as a probe [10]. This has a direct impact on investigating a range of clinical problems. For example, Lemberskiy et al. [33] utilized time dependent MD and FA for prostate cancer grading while Iima et al. [34] distinguished malignant from benign head and neck tumors by using time dependent apparent diffusion coefficient. Several other studies investigated the time dependency of alternative diffusion models [12,35]. For example, Iima et al. [12] observed significant time dependence of IVIM and DKI parameters in breast cancer and hepatocellular carcinoma xenograft; and Zhou et al. [35] reported time-dependence of the parameters in a fractional order calculus model in the human brain. Although the present study did not focus on a specific clinical problem, the results from the phantom study provided a useful guide to future investigations involving pathological tissue specimens, animal models, and clinical patients.
Our study has several limitations. First, the longest Δeff in our experiments was limited to 149.17 ms. This was largely due to the inadequate signal-to-noise ratio in the PGSTE acquisition. Additionally, a moderate TE of 75 ms was chosen to match the TE in the OGSE acquisition, thereby mitigating the potential issue with the TE-dependence in diffusion characterization, which can be particularly evident in a multi-compartmental environment [36,37,38]. If the PGSTE sequence is employed alone without the need to match parameters in other sequences, then studies on diffusion time dependency at long diffusion time regime can be conducted with a shorter TE. Second, our simulations did not cover the low diffusion time regime (e.g., <3.33 ms for OGSE and <25 ms for PGSE/PGSTE). This was because the limited total step size of 20,000 could lead to unstable results at shorter diffusion times. Optimized algorithms and more powerful computational platforms may help overcome this limitation. Third, although significant time dependency of CTRW model parameters was observed in this study, the analytical expressions [39] of this time dependency in a multi-compartmental environment requires further investigation. Finally, Sephadex gel phantoms provide a simple diffusion environment with spherical beads and permeable pores. Although they helped provide valuable insights into understanding of the complex diffusion processes, their limitations in adequately mimicking actual biological tissue structures must be recognized.

6. Conclusions

We have investigated time dependency of the CTRW model parameters in Sephadex gel phantoms across a broad range of diffusion times by using a set of pulse sequences comprising OGSE, PGSE, and PGSTE. We have experimentally observed monotonic decreases in Dα,β, α, and β as the diffusion time increased. These experimental results were reinforced by the Monte Carlo simulations. The present study provides valuable insights into probing microstructures by characterizing the time dependency of the CTRW model parameters, paving the way for future investigations on biological systems.

Author Contributions

Conceptualization, M.M.K. and X.J.Z.; methodology, G.D., R.L.M., X.J.Z. and M.M.K.; software, G.D. and Z.Z.; validation, G.D., X.J.Z. and M.M.K.; formal analysis, G.D. and M.M.K.; investigation, G.D., M.M.K., K.S. and Q.L.; resources, W.L.; data curation, G.D. and W.L.; writing—original draft preparation, G.D.; writing—review and editing, M.M.K. and X.J.Z.; visualization, G.D.; supervision, M.M.K. and X.J.Z.; project administration, M.M.K.; funding acquisition, X.J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Institutes of Health (5R01EB026716–01). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available on request.

Acknowledgments

We thank Yanyu Qian for assistance with phantom preparation.

Conflicts of Interest

X.J.Z. was a consultant to General Electric Healthcare. The other authors declare no conflict of interest.

References

  1. Le Bihan, D.; Poupon, C.; Amadon, A.; Lethimonnier, F. Artifacts and pitfalls in diffusion MRI. J. Magn. Reson. Imaging 2006, 24, 478–488. [Google Scholar] [CrossRef]
  2. Novikov, D.S.; Kiselev, V.G.; Jespersen, S.N. On modeling. Magn. Reson. Med. 2018, 79, 3172–3193. [Google Scholar] [CrossRef]
  3. Stejskal, E.O.; Tanner, J.E. Spin diffusion measurements: Spin echoes in the presence of a time-dependent field gradient. J. Chem. Phys. 1965, 42, 288–292. [Google Scholar] [CrossRef] [Green Version]
  4. Bernstein, M.A.; King, K.F.; Zhou, X.J. Handbook of MRI Pulse Sequences; Elsevier: Burlington, MA, USA, 2004. [Google Scholar]
  5. Mazo, R.M. Brownian Motion: Fluctuations, Dynamics, and Applications; Oxford University Press: New York, NY, USA, 2009. [Google Scholar]
  6. Sen, P.N. Time-dependent diffusion coefficient as a probe of geometry. Concepts Magn. Reson. Part A Bridg. Educ. Res. 2004, 23, 1–21. [Google Scholar] [CrossRef]
  7. Magin, R.L.; Abdullah, O.; Baleanu, D.; Zhou, X.J. Anomalous diffusion expressed through fractional order differential operators in the Bloch-Torrey equation. J. Magn. Reson. 2008, 190, 255–270. [Google Scholar] [CrossRef] [PubMed]
  8. Ingo, C.; Magin, R.L.; Colon-Perez, L.; Triplett, W.; Mareci, T.H. On random walks and entropy in diffusion-weighted magnetic resonance imaging studies of neural tissue. Magn. Reson. Med. 2014, 71, 617–627. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Pyatigorskaya, N.; Le Bihan, D.; Reynaud, O.; Ciobanu, L. Relationship between the diffusion time and the diffusion mri signal observed at 17.2 tesla in the healthy rat brain cortex. Magn. Reson. Med. 2014, 72, 492–500. [Google Scholar] [CrossRef]
  10. Aggarwal, M.; Smith, M.D.; Calabresi, P.A. Diffusion-time dependence of diffusional kurtosis in the mouse brain. Magn. Reson. Med. 2020, 84, 1564–1578. [Google Scholar] [CrossRef]
  11. Wu, D.; Zhang, J. Evidence of the diffusion time dependence of intravoxel incoherent motion in the brain. Magn. Reson. Med. 2019, 82, 2225–2235. [Google Scholar] [CrossRef] [PubMed]
  12. Iima, M.; Nobashi, T.; Imai, H.; Koyasu, S.; Saga, T.; Nakamoto, Y.; Kataoka, M.; Yamamoto, A.; Matsuda, T.; Togashi, K.; et al. Effects of diffusion time on non-gaussian diffusion and Intravoxel Incoherent Motion (IVIM) MRI parameters in breast cancer and hepatocellular carcinoma xenograft models. Acta Radiol. Open 2018, 7, 2058460117751565. [Google Scholar] [CrossRef]
  13. Reynaud, O. Time-dependent diffusion mri in cancer: Tissue modeling and applications. Front. Phys. 2017, 5, 1–16. [Google Scholar] [CrossRef] [Green Version]
  14. Magin, R.L.; Akpa, B.S.; Neuberger, T.; Webb, A.G. Fractional order analysis of sephadex gel structures: Nmr measurements reflecting anomalous diffusion. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 4581–4587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Liang, Y.; Chen, W.; Akpa, B.S.; Neuberger, T.; Webb, A.G.; Magin, R.L. Using spectral and cumulative spectral entropy to classify anomalous diffusion in sephadex™ gels. Comput. Math. Appl. 2017, 73, 765–774. [Google Scholar] [CrossRef]
  16. Li, H.W.; Ai, K.; Wu, Y. Fluorescence visual gel-separation of dansylated bsa-protected gold-nanoclusters. Chem. Commun. 2011, 47, 9852–9854. [Google Scholar] [CrossRef]
  17. Zhou, X.J.; Gao, Q.; Abdullah, O.; Magin, R.L. Studies of anomalous diffusion in the human brain using fractional order calculus. Magn. Reson. Med. 2010, 63, 562–569. [Google Scholar] [CrossRef] [PubMed]
  18. Karaman, M.M.; Sui, Y.; Wang, H.; Magin, R.L.; Li, Y.; Zhou, X.J. Differentiating low- and high-grade pediatric brain tumors using a continuous-time random-walk diffusion model at high b-values. Magn. Reson. Med. 2016, 76, 1149–1157. [Google Scholar] [CrossRef] [Green Version]
  19. Saichev, A.I.; Zaslavsky, G.M. Fractional kinetic equations: Solutions and applications. Chaos 1997, 7, 753–764. [Google Scholar] [CrossRef] [Green Version]
  20. Gorenflo, R.; Iskenderov; Luchko, Y. Mapping between solutions of fractional diffusion-wave equations. Fract. Calc. Appl. Anal. 2000, 3, 75–86. [Google Scholar]
  21. Magin, R.; Feng, X.; Baleanu, D. Solving the fractional order bloch equation. Concepts Magn. Reson. Part A Bridg. Educ. Res. 2009, 34, 16–23. [Google Scholar] [CrossRef]
  22. Ray, S.S.; Sahoo, S. Analytical approximate solutions of Riesz fractional diffusion equation and riesz fractional advection-dispersion equation involving nonlocal space fractional derivatives. Math. Methods Appl. Sci. 2015, 38, 2840–2849. [Google Scholar] [CrossRef]
  23. Barrick, T.R.; Spilling, C.A.; Ingo, C.; Madigan, J.; Isaacs, J.D.; Rich, P.; Jones, T.L.; Magin, R.L.; Hall, M.G.; Howe, F.A. Quasi-Diffusion Magnetic Resonance Imaging (QDI): A fast, high b-value diffusion imaging technique. Neuroimage 2020, 211, 116606. [Google Scholar] [CrossRef] [PubMed]
  24. Gao, Q.; Srinivasan, G.; Magin, R.L.; Joe Zhou, X. Anomalous diffusion measured by a twice-refocused spin echo pulse sequence: Analysis using fractional order calculus. J. Magn. Reson. Imaging 2011, 33, 1177–1183. [Google Scholar] [CrossRef] [Green Version]
  25. Karaman, M.M.; Zhou, X.J. A Fractional motion diffusion model for a twice-refocused spin-echo pulse sequence. NMR Biomed. 2018, 31, e3960. [Google Scholar] [CrossRef]
  26. Baron, C.A.; Beaulieu, C. Oscillating Gradient Spin-Echo (OGSE) diffusion tensor imaging of the human brain. Magn. Reson. Med. 2014, 72, 726–736. [Google Scholar] [CrossRef] [PubMed]
  27. Zhong, Z.; Merkitch, D.; Karaman, M.M.; Zhang, J.; Sui, Y.; Goldman, J.G.; Zhou, X.J. High-spatial-resolution diffusion mri in Parkinson disease: Lateral asymmetry of the substantia nigra. Radiology 2019, 291, 149–157. [Google Scholar] [CrossRef] [PubMed]
  28. Hall, M.G.; Alexander, D.C. Convergence and parameter choice for monte-carlo simulations of diffusion mri. IEEE Trans. Med. Imaging 2009, 28, 1354–1364. [Google Scholar] [CrossRef]
  29. Cook, P.A.; Bai, Y.; Hall, M.G.; Nedjati-Gilani, S.; Seunarine, K.K.; Alexander, D.C. Camino: Diffusion-MRI Reconstruction and Processing. In Proceedings of the 14th Scientific Meeting of the International Society for Magnetic Resonance in Medicine, Seattle, WA, USA, 6–12 May 2006; p. 2759. [Google Scholar]
  30. Magin, R.L.; Karani, H.; Wang, S.; Liang, Y. Fractional order complexity model of the diffusion signal decay in MRI. Mathematics 2019, 7, 348. [Google Scholar] [CrossRef] [Green Version]
  31. Novikov, D.S.; Fieremans, E.; Jensen, J.H.; Helpern, J.A. Random walks with barriers. Nat. Phys. 2011, 7, 508–514. [Google Scholar] [CrossRef] [Green Version]
  32. Bennett, K.M.; Schmainda, K.M.; Bennett, R.; Rowe, D.B.; Lu, H.; Hyde, J.S. Characterization of continuously distributed cortical water diffusion rates with a stretched-exponential model. Magn. Reson. Med. 2003, 50, 727–734. [Google Scholar] [CrossRef]
  33. Lemberskiy, G.; Rosenkrantz, A.B.; Veraart, J.; Taneja, S.S.; Novikov, D.S.; Fieremans, E. Time-dependent diffusion in prostate cancer. Investig. Radiol. 2017, 52, 405–411. [Google Scholar] [CrossRef]
  34. Iima, M.; Yamamoto, A.; Kataoka, M.; Yamada, Y.; Omori, K.; Feiweier, T.; Togashi, K. Time-dependent diffusion mri to distinguish malignant from benign head and neck tumors. J. Magn. Reson. Imaging 2019, 50, 88–95. [Google Scholar] [CrossRef] [PubMed]
  35. Zhou, X.J.; Gao, Q.; Srinivasan, G.; Magin, R.L. Dependence of Fractional Order Diffusion Model Parameters on Diffusion Time. In Proceedings of the 17th Annual Meeting of the International Society for Magnetic Resonance in Medicine, Honolulu, HI, USA, 18–24 April 2009; p. 1369. [Google Scholar]
  36. Veraart, J.; Novikov, D.S.; Fieremans, E. TE Dependent Diffusion Imaging (TEdDI) distinguishes between compartmental T2 relaxation times. Neuroimage 2018, 182, 360–369. [Google Scholar] [CrossRef] [PubMed]
  37. Lin, M.; He, H.; Tong, Q.; Ding, Q.; Yan, X.; Feiweier, T.; Zhong, J. Effect of myelin water exchange on DTI-derived parameters in diffusion mri: Elucidation of te dependence. Magn. Reson. Med. 2018, 79, 1650–1660. [Google Scholar] [CrossRef] [PubMed]
  38. Chatterjee, A.; Bourne, R.M.; Wang, S.; Devaraj, A.; Gallan, A.J.; Antic, T.; Karczmar, G.S.; Oto, A. Diagnosis of prostate cancer with noninvasive estimation of prostate tissue composition by using hybrid multidimensional MR imaging: A feasibility study. Radiology 2018, 287, 863–873. [Google Scholar] [CrossRef] [PubMed]
  39. Stanisz, G.J.; Szafer, A.; Wright, G.A.; Henkelman, R.M. An Analytical Model of Restricted Diffusion in Bovine Optic Nerve. Magn. Reson. Med. 1997, 37, 103–111. [Google Scholar] [CrossRef]
Figure 1. Pulse sequences employed in this study. (a) Cosine-trapezoid OGSE where N is the number of half oscillation period and δ is the total waveform duration (N = 4 in the sequence diagram). (b) PGSE where δ is the diffusion lobe duration and Δ is the diffusion lobe separation. (c) PGSTE where δ and Δ are defined similarly as in (b).
Figure 1. Pulse sequences employed in this study. (a) Cosine-trapezoid OGSE where N is the number of half oscillation period and δ is the total waveform duration (N = 4 in the sequence diagram). (b) PGSE where δ is the diffusion lobe duration and Δ is the diffusion lobe separation. (c) PGSTE where δ and Δ are defined similarly as in (b).
Mathematics 09 01688 g001
Figure 2. DW images acquired by using the PGSE sequence (Δ = 35 ms and b = 0). (a) The first Sephadex series: G25–50 (bottom right), G50–50 (bottom left), and G75–50 (top). (b) The second Sephadex series: G50–50 (bottom left), G50–80 (bottom right), and G50–150 (top). The rhombus-shaped ROIs indicate the regions used to calculate the mean parameter values.
Figure 2. DW images acquired by using the PGSE sequence (Δ = 35 ms and b = 0). (a) The first Sephadex series: G25–50 (bottom right), G50–50 (bottom left), and G75–50 (top). (b) The second Sephadex series: G50–50 (bottom left), G50–80 (bottom right), and G50–150 (top). The rhombus-shaped ROIs indicate the regions used to calculate the mean parameter values.
Mathematics 09 01688 g002
Figure 3. Plots of Dα,β versus effective diffusion time, Δeff, for gels G25–50 (a) and G50–50 (b), and G75–50 (c) with increased macromolecular exclusion limit. The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Figure 3. Plots of Dα,β versus effective diffusion time, Δeff, for gels G25–50 (a) and G50–50 (b), and G75–50 (c) with increased macromolecular exclusion limit. The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Mathematics 09 01688 g003
Figure 4. Plots of temporal fractional order (α) versus effective diffusion time, Δeff, for gels G25–50 (a) and G50–50 (b), and G75–50 (c). The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Figure 4. Plots of temporal fractional order (α) versus effective diffusion time, Δeff, for gels G25–50 (a) and G50–50 (b), and G75–50 (c). The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Mathematics 09 01688 g004
Figure 5. Plots of spatial fractional order (β) versus effective diffusion time, Δeff, for gels G25–50 (a) and G50–50 (b), and G75–50 (c). The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Figure 5. Plots of spatial fractional order (β) versus effective diffusion time, Δeff, for gels G25–50 (a) and G50–50 (b), and G75–50 (c). The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Mathematics 09 01688 g005
Figure 6. Plot of Dα,β versus effective diffusion time, Δeff, for gels G50–50 (a), G50–80 (b), and G50–150 (c). The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Figure 6. Plot of Dα,β versus effective diffusion time, Δeff, for gels G50–50 (a), G50–80 (b), and G50–150 (c). The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Mathematics 09 01688 g006
Figure 7. Plot of temporal fractional order (α) versus effective diffusion time, Δeff, for gels G50–50 (a), G50–80 (b), and G50–150 (c). The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Figure 7. Plot of temporal fractional order (α) versus effective diffusion time, Δeff, for gels G50–50 (a), G50–80 (b), and G50–150 (c). The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Mathematics 09 01688 g007
Figure 8. Plot of spatial fractional order (β) versus effective diffusion time, Δeff, for gels G50–50 (a), G50–80 (b), and G50–150 (c). The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Figure 8. Plot of spatial fractional order (β) versus effective diffusion time, Δeff, for gels G50–50 (a), G50–80 (b), and G50–150 (c). The data acquired by using the OGSE, PGSE, and PGSTE pulse sequences are marked in black, red, and blue, respectively.
Mathematics 09 01688 g008
Figure 9. Plots of Dα,β (a), α (b), and β (c) versus Δeff obtained from the Monte Carlo simulations with fixed r = 8 µm and varying p of 0.1% (red), 0.2% (green) and 0.4% (blue). The rhombi and circles represent the simulation results with oscillating diffusion gradient (OGSE) and Stejskal–Tanner diffusion gradient (PGSE/PGSTE), respectively.
Figure 9. Plots of Dα,β (a), α (b), and β (c) versus Δeff obtained from the Monte Carlo simulations with fixed r = 8 µm and varying p of 0.1% (red), 0.2% (green) and 0.4% (blue). The rhombi and circles represent the simulation results with oscillating diffusion gradient (OGSE) and Stejskal–Tanner diffusion gradient (PGSE/PGSTE), respectively.
Mathematics 09 01688 g009
Figure 10. Plots of Dα,β (a), α (b), and β (c) versus Δeff from Monte Carlo simulation with fixed p = 0.2% and varying r of 6 µm (red), 7 µm (green) and 8 µm (blue). The rhombi and circles represent simulation results with oscillating diffusion gradient (OGSE) and Stejskal–Tanner diffusion gradient (PGSE/PGSTE), respectively.
Figure 10. Plots of Dα,β (a), α (b), and β (c) versus Δeff from Monte Carlo simulation with fixed p = 0.2% and varying r of 6 µm (red), 7 µm (green) and 8 µm (blue). The rhombi and circles represent simulation results with oscillating diffusion gradient (OGSE) and Stejskal–Tanner diffusion gradient (PGSE/PGSTE), respectively.
Mathematics 09 01688 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dan, G.; Li, W.; Zhong, Z.; Sun, K.; Luo, Q.; Magin, R.L.; Zhou, X.J.; Karaman, M.M. Diffusion in Sephadex Gel Structures: Time Dependency Revealed by Multi-Sequence Acquisition over a Broad Diffusion Time Range. Mathematics 2021, 9, 1688. https://doi.org/10.3390/math9141688

AMA Style

Dan G, Li W, Zhong Z, Sun K, Luo Q, Magin RL, Zhou XJ, Karaman MM. Diffusion in Sephadex Gel Structures: Time Dependency Revealed by Multi-Sequence Acquisition over a Broad Diffusion Time Range. Mathematics. 2021; 9(14):1688. https://doi.org/10.3390/math9141688

Chicago/Turabian Style

Dan, Guangyu, Weiguo Li, Zheng Zhong, Kaibao Sun, Qingfei Luo, Richard L. Magin, Xiaohong Joe Zhou, and M. Muge Karaman. 2021. "Diffusion in Sephadex Gel Structures: Time Dependency Revealed by Multi-Sequence Acquisition over a Broad Diffusion Time Range" Mathematics 9, no. 14: 1688. https://doi.org/10.3390/math9141688

APA Style

Dan, G., Li, W., Zhong, Z., Sun, K., Luo, Q., Magin, R. L., Zhou, X. J., & Karaman, M. M. (2021). Diffusion in Sephadex Gel Structures: Time Dependency Revealed by Multi-Sequence Acquisition over a Broad Diffusion Time Range. Mathematics, 9(14), 1688. https://doi.org/10.3390/math9141688

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop