IHMT: Is it misnamed? A simple theoretical description of "inhomogeneous" MT.
Alan P Manning1, Kimberley L Chang2, Alex MacKay1,3, and Carl A Michal1

1Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada, 2Department of Neurology, University of British Columbia, Vancouver, BC, Canada, 3UBC MRI Research Centre, Department of Radiology, University of British Columbia, Vancouver, BC, Canada

Synopsis

Inhomogeneous MT (IHMT) shows promise for myelin-selectivity. Images acquired with soft prepulses at positive and negative offsets simultaneously show a reduced intensity compared to images with a single positive or negative offset prepulse. The leading hypothesis is that this works due to inhomogeneous broadening of the lipid proton line. Our results contradict this. We show that IHMT can be explained by a simple spin-1 model of a coupled methylene pair, and that it occurs in homogeneously-broadened systems (hair and wood). We propose the relevant timescales for IHMT are the dipolar coupling correlation time and the prepulse nutation period.

Introduction

Inhomogeneous Magnetization Transfer (IHMT) appears to be myelin-selective[1,2,3]. In IHMT, images are acquired with four different soft prepulses: no prepulse (S0), prepulse at offset +Δ (S+) or -Δ (S-), and prepulses at both ±Δ simultaneously (Sboth). The IHMT Ratio, typically measured from aqueous proton intensity, is

$$\text{IHMTR}=\frac{S_+ + S_- - 2S_{\text{both}}}{2S_0} \le 1$$

IHMTR≠0 in systems with lipid lamellar structures (e.g. myelin, hair conditioner)[1,2], but ≈0 in many other systems (e.g. gelatin, agarose)[4].

Lipid hydrocarbon chain spin systems are thermally-averaged, resulting in spectral behaviour dominated by intra-methylene dipolar couplings. These methylenes resemble an ensemble of spin-1 systems[5,6], and others have suggested this underlies the IHMT mechanism[1,7]. Recent work has also highlighted IHMT's sensitivity to T1D, the dipolar order relaxation time[2]. Nonetheless, the accepted view remains that IHMT's lipid-selectivity results from lipid chains having inhomogeneous lineshapes[1,2,3,7].

Here, we present experimental results on homogeneously-broadened systems with IHMTR≠0, contradicting the inhomogeneous explanation. The following simple spin-1 model shows how this is possible.

We consider the non-aqueous protons only, since their behaviour under the prepulses determines if IHMT will occur. The Hamiltonian of a dipolar-coupled methylene pair is $$${\cal H} = -\omega_0 I_z -\omega_D I_z^2$$$ (ω0=Larmor frequency, ωD=dipolar interaction strength, ω0»ωD). The equilibrium density matrix is

$$\rho = M_0 I_z = M_0\, \text{diag}(1,0,-1).$$

Then, assuming the prepulses are calibrated for inversion, S+ exchanges ρ11/ρ22, S- exchanges ρ22/ρ33, and Sboth exchanges ρ11/ρ33 (using 2× the power to invert both transitions). Following the prepulse, $$$\langle I_z \rangle$$$ is the magnetization, and $$$\langle I_z^2 \rangle$$$ the dipolar order (nonzero only in S+/S-).

We form the "non-aqueous IHMTR":

$$\begin{align}\text{non-aqueous IHMTR} &=\frac{(\langle I_z \rangle_+ -\langle I_z \rangle_0)+(\langle I_z \rangle_- -\langle I_z \rangle_0) - (\langle I_z \rangle_{\text{both}} - \langle I_z \rangle_0)}{2\langle I_z \rangle_0} \\&= \frac{ (-\frac{1}{2}) + (-\frac{1}{2}) - (-2) }{2} =\frac{1}{2}.\end{align}$$

We don't multiply $$$(\langle I_z \rangle_0 - \langle I_z \rangle_{\text{both}})$$$ by 2 because we assumed 2× the power was used to invert both transitions.

Following a hard observe pulse of flip angle α, the spectral line amplitudes A± at ω0±ωD are

$$A_\pm = [\frac{1}{2}(\rho_{11}-\rho_{33})]\sin\alpha\pm [\frac{1}{2}(\rho_{11}+\rho_{33}) -\rho_{22}] \cos\alpha \sin\alpha. ~~~~~\text{(1)}$$

Evidently, irradiating one transition affects the other's amplitude, and both amplitudes are equal when α=π/2. The lines cannot be independently saturated, as would be the case in an inhomogeneously-broadened system, yet it has a nonzero IHMTR. This equation's second term produces spectral asymmetry from dipolar order, and can only be observed when α≠nπ/2.

Methods

10% (w/w) Prolipid-161 (PL161, Ashland, DE, USA) in 99% D2O was used for spin-1 model comparison, T1D, and IHMT measurements. PL161 is a suitable myelin model for IHMT[4,8]. IHMT and T1D was also measured in air-dried Douglas Fir wood (sapwood), human hair, and lamb Achilles tendon. Wood and hair have homogeneously-broadened non-aqueous proton spectra, primarily from cellulose and α-keratin crystallites, respectively[9,10,11]. Both have exchangeable aqueous protons. The NMR experiments utilized a home-built spectrometer (200 MHz / 4.7 T). Figure 1 shows the IHMT and Adiabatic Demagnetization/Remagnetization in the Rotating Frame (ADRF/ARRF) sequences used.

Results and Discussion

The PL161 and fir non-aqueous proton spectra following IHMT prepulses (Figure 2) showed no evidence of hole-burning. However, with α=33o, the characteristic spectral asymmetry from dipolar order in S+/S- is obvious.

As a function of α, the ω>0/ω<0 spectral integrals of PL161 were fit to (1) (Figure 3). Inhomogeneously-broadened spectra would only have a sin(α) dependence, but we clearly see more complex behaviour consistent with dipolar order.

Our samples' non-aqueous and aqueous IHMTRs vs. Δ and their T1Ds are given in Figures 4 and 5, respectively. Fir and hair have IHMTR≠0, despite being homogeneously-broadened systems. Also, although PL161's T1D is 20-40× larger than the fir and hair's, the corresponding IHMTRs are only 2-3× larger. Consequently, T1D is likely not the only relevant non-aqueous timescale.

Many tissues and phantoms in which IHMTR≈0 are highly dynamic systems, so dipolar couplings are transient. For IHMT to occur, we believe the coupling lifetime must be long enough for dipolar order to develop: τD$$$\gtrsim$$$1/(4ν1) (τD = dipolar coupling correlation time, ν1 = prepulse nutation frequency). ν1 must be small enough to only excite one side of the spectrum, setting a lower limit for τD. In a system with transient dipolar couplings, T1D is constrained to be no greater than τD.

Conclusions

We have given evidence that IHMT doesn't rely upon inhomogeneously-broadened non-aqueous proton lines. Our nonzero IHMTR results in homogeneously-broadened systems are explained with a simple spin-1 model. Furthermore, we have proposed that nonzero IHMTRs depend upon the dipolar coupling correlation time (≈T1D in motionally-dynamic systems), which must be long enough for appreciable dipolar order development during the prepulse.

Acknowledgements

APM acknowledges an NSERC postgraduate award. AM and CAM acknowledge NSERC Discovery grants.

References

[1] Varma, G., Duhamel, G., de Bazelaire, C. & Alsop, D. C. Magnetization transfer from inhomogeneously broadened lines: A potential marker for myelin. Magn. Reson. Med. 73, 614–622 (2015)

[2] Varma, G. et al. Interpretation of magnetization transfer from inhomogeneously broadened lines (ihMT) in tissues as a dipolar order effect within motion restricted molecules. J. Magn. Reson. 260, 67–76 (2015)

[3] Girard, O. M. et al. Magnetization transfer from inhomogeneously broadened lines (ihMT): Experimental optimization of saturation parameters for human brain imaging at 1.5 Tesla. Magn. Reson. Med. 73, 2111–2121 (2015)

[4] Swanson, S. D., Alsop, D. C. & Malyarenko, D. I. The physical basis of inhomogeneous magnetization transfer. In Proceedings of the 22nd Annual Meeting of ISMRM, 0391 (Milano, Italy, 2014)

[5] Higgs, T. & Mackay, A. Determination of the complete order parameter tensor for a lipid methylene group from 1H- and 2H-NMR spin labels. Chem. Phys. Lipids 20, 105–114 (1977)

[6] Bloom, M., Burnell, E. E., Roeder, S. B. W. & Valic, M. I. Nuclear magnetic resonance line shapes in lyotropic liquid crystals and related systems. J. Chem. Phys. 66, 3012 (1977)

[7] Swanson, S. D., Malyarenko, D. I. & Fabiilli, M. L. Towards a quantitative theory for inhomogeneous magnetization transfer. In Proceedings of the 23th Annual Meeting of ISMRM, 0994 (Toronto, Canada, 2015)

[8] Zhao, F., Nielsen, J.-F., Swanson, S. D., Fessler, J. A. & Noll, D. C. Simultaneous fat saturation and magnetization transfer contrast imaging with steady-state incoherent sequences. Magn. Reson. Med. 74, 739–746 (2015)

[9] Peemoeller, H. et al. NMR detection of liquid-like wood polymer component in dryaspen wood. Polymer 54, 1524–1529 (2013)

[10] Wolfram, L. J. Human hair: A unique physicochemical composite. Journal of the American Academy of Dermatology 48, S106–S114 (2003)

[11] Melian, C. et al. Morphology and side-chain dynamics in hydrated hard α-keratin fibres by 1H solid-state NMR. Chem. Phys. Lett. 480, 300–304 (2009)

Figures

Pulse sequences. A/B: Single soft Gaussian/rectangular prepulses at +Δ or -Δ, (giving S+ or S-). C: Via sine modulation, the prepulse is applied simultaneously at ±Δ (giving Sboth). D: ADRF/ARRF to measure T1D. (B1=pulse amplitude, τ=pulse time, δ=delay time, α=hard observe pulse flip angle.)

Non-aqueous PL161 and Douglas Fir proton spectra following IHMT prepulses. S+ and S- spectra show dipolar order (manifesting as spectral asymmetry), but not S0 or Sboth. No evidence of isolated "hole burning" characteristic of inhomogeneous broadening is seen. Parameters: τ=20 ms, |Δ|=10 kHz, B1=1 kHz, α=33o.

PL161 spectrum vs. observe pulse flip angle α, compared to the spin-1 model. S+ used sequence A (Δ=+2.5 kHz, B1 = 2.5 kHz, α=0.5 ms). I< / I> are non-aqueous intensities of (S+ - S0), integrated between (-80,-3)/(3,80) kHz. An exp(-Cα) multiplicative term (C=parameter) was included to account for B1 inhomogeneities.

Aqueous and non-aqueous IHMTRs vs. prepulse offset. Samples known to have homogeneously broadened spectra (wood and hair) still have IHMTR≠0. Sequences A and B were used (τ=500 ms, B1=415 Hz, δ=0.5 ms, α=33o). The non-aqueous and aqueous intensities were integrated between ((-80,3) and (3,80) kHz) and (-1,1) kHz, respectively.

T1D values measured with the sequence C. A stretched exponential fit was used to account for a distribution of T1Ds from multiple orientations. Even with δ=0.1 ms, there was no detectable signal from lamb tendon. Ramp times were τ=1 ms.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0305