Human whole blood 1 H2 O transverse relaxation with gadolinium-based contrast reagents: Magnetic susceptibility and transmembrane water exchange

Detection of elevated succinate in brain during circulatory arrest in a piglet model: A 3T 1H MR spectroscopy study

PURPOSE

To measure and validate elevated succinate in brain during circulatory arrest in a piglet model of cardiopulmonary bypass.

METHODS

Using data from an archive of 3T 1H MR spectra acquired in previous in-magnet studies, dynamic plots of succinate, spectral simulations and difference spectra were generated for analysis and validation.

RESULTS

Elevation of succinate during circulatory arrest was observed and validated. Fitting bias was evaluated as a function of the line-widths and signal-to-noise ratios of the archived data. Succinate increases were independent of bypass temperature. Succinate elevation was also not observed with antegrade cerebral perfusion.

CONCLUSION

Although spectrally overlapped and at sub-millimolar levels, elevated brain succinate can be reliably measured by dynamic MR spectroscopy at 3T. Noise dependent bias of the stronger overlapping signals did not impact the succinate measurement. Elevated succinate during circulatory arrest and its recovery after reperfusion was observed. This finding is consistent with earlier reports that correlate elevated succinate with ischemic-reperfusion injury.

The Road to a Realistic 3D Model for Estimating R2 and R2* Relaxation Versus Gd-DTPA Concentration in Whole Blood and Brain Tumor Vasculature

Dynamic susceptibility contrast (DSC) MRI is commonly part of brain tumor imaging. For quantitative analysis, measurement of the arterial input function and tissue concentration time curve is required. Usually, a linear relationship between the MR signal changes and contrast agent concentration ([Gd]) is assumed, even though this is a known simplification. The aim of this study was to develop a realistic 3D simulation model as an efficient method to assess the relationship between ΔR2 (*) and [Gd] both in whole blood and brain tissue. We modified an open-source 3D simulation model to study different red blood cell configurations for assessing whole-blood ΔR2 (*) versus [Gd]. The results were validated against previously obtained 2D data and in vitro data. Furthermore, hematocrit levels (30%-50%) and field strengths (1.5-3.0-7.0 T) were varied. Subsequently, realistic tumor vascular networks were derived from intraoperative high framerate Doppler ultrasound data to study the influence of vascular structure and orientation with respect to the main magnetic field (1.5-3.0-7.0 T) for the calculation of ΔR2 (*) versus [Gd] in brain tissue. For whole blood, good agreement of the 3D model was found with in vitro and 2D simulation data when red blood cells were aligned with the blood flow. For brain tissue, minor differences were found between the vascular networks. The effect of vessel direction with respect to B0 was apparent in case of clear directionality of the main vessels. The dependency on field strength agreed with previous reports. In conclusion, we have shown that the relationship between ΔR2 (*) and [Gd] is affected by the organization of red blood cells and orientation of blood vessels with respect to the main magnetic field, as well as the field strength. These findings are important for further optimization of the realistic 3D model that could eventually be used to improve the estimation of hemodynamic parameters from DSC-MRI.

Evolving Characteristics of Gadolinium-Based Contrast Agents for MR Imaging: A Systematic Review of the Importance of Relaxivity

Gadolinium-based contrast agents (GBCAs) are widely and routinely used to enhance the diagnostic performance of magnetic resonance imaging and magnetic resonance angiography examinations. T1 relaxivity (r1) is the measure of their ability to increase signal intensity in tissues and blood on T1-weighted images at a given dose. Pharmaceutical companies have invested in the design and development of GBCAs with higher and higher T1 relaxivity values, and "high relaxivity" is a claim frequently used to promote GBCAs, with no clear definition of what "high relaxivity" means, or general concurrence about its clinical benefit. To understand whether higher relaxivity values translate into a material clinical benefit, well-designed, and properly powered clinical studies are necessary, while mere in vitro measurements may be misleading. This systematic review of relevant peer-reviewed literature provides high-quality clinical evidence showing that a difference in relaxivity of at least 40% between two GBCAs results in superior diagnostic efficacy for the higher-relaxivity agent when this is used at the same equimolar gadolinium dose as the lower-relaxivity agent, or similar imaging performance when used at a lower dose. Either outcome clearly implies a relevant clinical benefit. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 3.

Using Adaptive Imaging Parameters to Improve PEGylated Ultrasmall Iron Oxide Nanoparticles-Enhanced Magnetic Resonance Angiography

The PEGylated ultrasmall iron oxide nanoparticles (PUSIONPs) exhibit longer blood residence time and better biodegradability than conventional gadolinium-based contrast agents (GBCAs), enabling prolonged acquisitions in contrast-enhanced magnetic resonance angiography (CE-MRA) applications. The image quality of CE-MRA is dependent on the contrast agent concentration and the parameters of the pulse sequences. Here, a closed-form mathematical model is demonstrated and validated to automatically optimize the concentration, echo time (TE), repetition time (TR) and flip angle (FA). The pharmacokinetic studies are performed to estimate the dynamic intravascular concentrations within 12 h postinjection, and the adaptive concentration-dependent sequence parameters are determined to achieve optimal signal enhancement during a prolonged measurement window. The presented model is tested on phantom and in vivo rat images acquired from a 3T scanner. Imaging results demonstrate excellent agreement between experimental measurements and theoretical predictions, and the adaptive sequence parameters obtain better signal enhancement than the fixed ones. The low-dose PUSIONPs (0.03 mmol kg-1 and 0.05 mmol kg-1) give a comparable signal intensity to the high-dose one (0.10 mmol kg-1) within 2 h postinjection. The presented mathematical model provides guidance for the optimization of the concentration and sequence parameters in PUSIONPs-enhanced MRA, and has great potential for further clinical translation.

Efficacy of Whole-Blood Model of Gadolinium-Based Contrast Agent Relaxivity in Predicting Vascular MR Signal Intensity In Vivo

BACKGROUND

Previous in vitro studies have described sub-linear longitudinal and heightened transverse H2O relaxivities of gadolinium-based contrast agents (GBCAs) in blood due to their extracellular nature. However, in vivo validation is lacking.

PURPOSE

Validate theory describing blood behavior of R1 and R2* in an animal model.

STUDY TYPE

Prospective, animal.

ANIMAL MODEL

Seven swine (54-65 kg).

FIELD STRENGTH/SEQUENCE

1.5 T; time-resolved 3D spoiled gradient-recalled echo (SPGR) and quantitative Look-Locker and multi-echo fast field echo sequences.

ASSESSMENT

Seven swine were each injected three times with 0.1 mmol/kg intravenous doses of one of three GBCAs: gadoteridol, gadobutrol, and gadobenate dimeglumine. Injections were randomized for rate (1, 2, and 3 mL/s) and order, during which time-resolved aortic 3D SPGR imaging was performed concurrently with aortic blood sampling via an indwelling catheter. Time-varying [GBCA] was measured by mass spectrometry of sampled blood. Predicted signal intensity (SI) was determined from a model incorporating sub-linear R1 and R2* effects (whole-blood model) and simpler models incorporating linear R1, with and without R2* effects. Predicted SIs were compared to measured aortic SI.

STATISTICAL TESTS

Linear correlation (coefficient of determination, R2) and mean errors were compared across the SI prediction models.

RESULTS

There was an excellent correlation between predicted and measured SI across all injections and swine when accounting for the non-linear dependence of R1 and high blood R2* (regression slopes 0.91-1.04, R2 ≥ 0.91). Simplified models (linear R1 with and without R2* effects) showed poorer correlation (slopes 0.67-0.85 and 0.54-0.64 respectively, both R2 ≥ 0.89) and higher averaged mean absolute and mean square errors (128.4 and 177.4 vs. 42.0, respectively, and 5506 and 11,419 vs. 699, respectively).

DATA CONCLUSION

Incorporating sub-linear R1 and high first-pass R2* effects in arterial blood models allows accurate SPGR SI prediction in an in vivo animal model, and might be utilized when modeling MR blood SI.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY: Stage 1.

In-cell NMR: Why and how?

NMR spectroscopy has been applied to cells and tissues analysis since its beginnings, as early as 1950. We have attempted to gather here in a didactic fashion the broad diversity of data and ideas that emerged from NMR investigations on living cells. Covering a large proportion of the periodic table, NMR spectroscopy permits scrutiny of a great variety of atomic nuclei in all living organisms non-invasively. It has thus provided quantitative information on cellular atoms and their chemical environment, dynamics, or interactions. We will show that NMR studies have generated valuable knowledge on a vast array of cellular molecules and events, from water, salts, metabolites, cell walls, proteins, nucleic acids, drugs and drug targets, to pH, redox equilibria and chemical reactions. The characterization of such a multitude of objects at the atomic scale has thus shaped our mental representation of cellular life at multiple levels, together with major techniques like mass-spectrometry or microscopies. NMR studies on cells has accompanied the developments of MRI and metabolomics, and various subfields have flourished, coined with appealing names: fluxomics, foodomics, MRI and MRS (i.e. imaging and localized spectroscopy of living tissues, respectively), whole-cell NMR, on-cell ligand-based NMR, systems NMR, cellular structural biology, in-cell NMR… All these have not grown separately, but rather by reinforcing each other like a braided trunk. Hence, we try here to provide an analytical account of a large ensemble of intricately linked approaches, whose integration has been and will be key to their success. We present extensive overviews, firstly on the various types of information provided by NMR in a cellular environment (the "why", oriented towards a broad readership), and secondly on the employed NMR techniques and setups (the "how", where we discuss the past, current and future methods). Each subsection is constructed as a historical anthology, showing how the intrinsic properties of NMR spectroscopy and its developments structured the accessible knowledge on cellular phenomena. Using this systematic approach, we sought i) to make this review accessible to the broadest audience and ii) to highlight some early techniques that may find renewed interest. Finally, we present a brief discussion on what may be potential and desirable developments in the context of integrative studies in biology.

Proton magnetic resonance spectroscopy assessment of neonatal brain metabolism during cardiopulmonary bypass surgery

Here, we report on the development and performance of a robust 3-T single-voxel proton magnetic resonance spectroscopy (¹ H MRS) experimental protocol and data analysis pipeline for quantifying brain metabolism during cardiopulmonary bypass (CPB) surgery in a neonatal porcine model, with the overall goal of elucidating primary mechanisms of brain injury associated with these procedures. The specific aims were to assess which metabolic processes can be reliably interrogated by ¹ H MRS on a 3-T clinical scanner and to provide an initial assessment of brain metabolism during deep hypothermia cardiac arrest (DHCA) surgery and recovery. Fourteen neonatal pigs underwent CPB surgery while placed in a 3-T MRI scanner for 18, 28, and 37°C DHCA studies under hyperglycemic, euglycemic, and hypoglycemic conditions. Total imaging times, including baseline measurements, circulatory arrest (CA), and recovery averaged 3 h/animal, during which 30-40 single-voxel ¹ H MRS spectra (sLASER pulse sequence, TR/TE = 2000/30 ms, 64 or 128 averages) were acquired from a 2.2-cc right midbrain voxel. ¹ H MRS at 3 T was able to reliably quantify (1) anaerobic metabolism via depletion of brain glucose and the associated build-up of lactate during CA, (2) phosphocreatine (PCr) to creatine (Cr) conversion during CA and subsequent recovery upon reperfusion, (3) a robust increase in the glutamine-to-glutamate (Gln/Glu) ratio during the post-CA recovery period, and (4) a broadening of the water peak during CA. In vivo ¹ H MRS at 3 T can reliably quantify subtle metabolic brain changes previously deemed challenging to interrogate, including brain glucose concentrations even under hypoglycemic conditions, ATP usage via the conversion of PCr to Cr, and differential changes in Glu and Gln. Observed metabolic changes during CPB surgery of a neonatal porcine model provide new insights into possible mechanisms for prevention of neuronal injury.

Dependency of R2 and R2 * relaxation on Gd-DTPA concentration in arterial blood: Influence of hematocrit and magnetic field strength

Dynamic susceptibility contrast (DSC) MRI is clinically used to measure brain perfusion by monitoring the dynamic passage of a bolus of contrast agent through the brain. For quantitative analysis of the DSC images, the arterial input function is required. It is known that the original assumption of a linear relation between the R₂^(*) relaxation and the arterial contrast agent concentration is invalid, although the exact relation is as of yet unknown. Studying this relation in vitro is time-consuming, because of the widespread variations in field strengths, MRI sequences, contrast agents, and physiological conditions. This study aims to simulate the R₂^(*) versus contrast concentration relation under varying physiological and technical conditions using an adapted version of an open-source simulation tool. The approach was validated with previously acquired data in human whole blood at 1.5 T by means of a gradient-echo sequence (proof-of-concept). Subsequently, the impact of hematocrit, field strength, and oxygen saturation on this relation was studied for both gradient-echo and spin-echo sequences. The results show that for both gradient-echo and spin-echo sequences, the relaxivity increases with hematocrit and field strength, while the hematocrit dependency was nonlinear for both types of MRI sequences. By contrast, oxygen saturation has only a minor effect. In conclusion, the simulation setup has proven to be an efficient method to rapidly calibrate and estimate the relation between R₂^(*) and gadolinium concentration in whole blood. This knowledge will be useful in future clinical work to more accurately retrieve quantitative information on brain perfusion.

Relationship Between Simulated Gadolinium-Based Contrast Agent Injection Profile and Achievable Resolution Metrics in Contrast-Enhanced Magnetic Resonance Angiography

BACKGROUND

Contrast bolus variation during contrast-enhanced magnetic resonance angiography (CE-MRA) acquisition may lead to vessel blurring.

PURPOSE

To combine knowledge of how contrast signal intensity (SI) evolves for different injection strategies with anatomically familiar parametric computer models to measure and visually assess the effects of a wide range of variables on modeled CE-MRA, and in doing so develop contrast rate injection guidelines.

STUDY TYPE

Computer modeling.

PHANTOM

Digital three-dimensional phantom consisting of orthogonal "aorta," 7 mm diameter "renal arteries" (with 57% and 86% diameter stenoses), and 7 mm diameter "superior mesenteric artery" (with 57% diameter stenosis).

FIELD STRENGTH/SEQUENCE

One millimeter in-plane resolution arterial CE-MRA imaging at 3 T.

ASSESSMENT

"Background" (time invariant) and "vascular" (time varying) components of the phantom were each Fourier transformed into the spatial frequency domain, the latter modulated by the SI evolution of a contrast bolus of varying "plateau" lengths and "tail" heights. Data are presented as surface plots of stenosis measurement error and blurring vs. a reference-standard injection.

STATISTICAL TESTS

Descriptive.

RESULTS

Shorter plateau lengths and lower tail heights resulted in increased measured stenosis error and blurring vs. the reference standard. Under a 44-second acquisition, full width half maximum stenosis error of the 86% stenosis with 25% plateau length and 25% tail height is 24% as compared to that from the reference standard. As plateau length and tail height approach 100%, stenosis error and blurring approach a floor defined by the MR acquisition's limitations.

DATA CONCLUSION

We propose that to achieve minimal degradation with CE-MRA, one can create a contrast bolus with either 60% plateau and 50% tail height or 80% plateau with any tail. These considerations may well prove to be of practical importance, possibly via manipulating the tail by means of multiphasic contrast injections.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 1.

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.

Quantitative 3D myocardial perfusion with an efficient arterial input function

PURPOSE

The purpose of this study was to further develop and combine several innovative sequence designs to achieve quantitative 3D myocardial perfusion. These developments include an optimized 3D stack-of-stars readout (150 ms per beat), efficient acquisition of a 2D arterial input function, tailored saturation pulse design, and potential whole heart coverage during quantitative stress perfusion.

THEORY AND METHODS

All studies were performed free-breathing on a Prisma 3T MRI scanner. Phantom validation was used to verify sequence accuracy. A total of 21 subjects (3 patients with known disease) were scanned, 12 with a rest only protocol and 9 with both stress (regadenoson) and rest protocols. First pass quantitative perfusion was performed with gadoteridol (0.075 mmol/kg).

RESULTS

Implementation and quantitative perfusion results are shown for healthy subjects and subjects with known coronary disease. Average rest perfusion for the 15 included healthy subjects was 0.79 ± 0.19 mL/g/min, the average stress perfusion for 6 healthy subject studies was 2.44 ± 0.61 mL/g/min, and the average global myocardial perfusion reserve ratio for 6 healthy subjects was 3.10 ± 0.24. Perfusion deficits for 3 patients with ischemia are shown. Average resting heart rate was 59 ± 7 bpm and the average stress heart rate was 81 ± 10 bpm.

CONCLUSION

This work demonstrates that a quantitative 3D myocardial perfusion sequence with the acquisition of a 2D arterial input function is feasible at high stress heart rates such as during stress. T1 values and gadolinium concentrations of the sequence match the reference standard well in a phantom, and myocardial rest and stress perfusion and myocardial perfusion reserve values are consistent with those published in literature.

Quantitative theory for the transverse relaxation time of blood water

An integrative model is proposed to describe the dependence of the transverse relaxation rate of blood water protons (R_2blood = 1/T_2blood ) on hematocrit fraction and oxygenation fraction (Y). This unified model takes into account (a) the diamagnetic effects of albumin, hemoglobin and the cell membrane; (b) the paramagnetic effect of hemoglobin; (c) the effect of compartmental exchange between plasma and erythrocytes under both fast and slow exchange conditions that vary depending on field strength and compartmental relaxation rates and (d) the effect of diffusion through field gradients near the erythrocyte membrane. To validate the model, whole-blood and lysed-blood R₂ data acquired previously using Carr-Purcell-Meiboom-Gill measurements as a function of inter-echo spacing τ_cp at magnetic fields of 3.0, 7.0, 9.4 and 11.7 T were fitted to determine the lifetimes (field-independent physiological constants) for water diffusion and exchange, as well as several physical constants, some of which are field-independent (magnetic susceptibilities) and some are field-dependent (relaxation rates for water protons in solutions of albumin and oxygenated and deoxygenated hemoglobin, ie, blood plasma and erythrocytes, respectively). This combined exchange-diffusion model allowed excellent fitting of the curve of the τ_cp -dependent relaxation rate dispersion at all four fields using a single average erythrocyte water lifetime, τ_ery = 9.1 ± 1.4 ms, and an averaged diffusional correlation time, τ_D = 3.15 ± 0.43 ms. Using this model and the determined physiological time constants and relaxation parameters, blood T₂ values published by multiple groups based on measurements at magnetic field strengths of 1.5 T and higher could be predicted correctly within error. Establishment of this theory is a fundamental step for quantitative modeling of the BOLD effect underlying functional MRI.

In vitro T2 relaxivities of the Gd-based contrast agents (GBCAs) in human blood at 1.5 and 3 T

BACKGROUND

The availability of data in the medical literature for the T₂ relaxivities of the Gd-based contrast agents (GBCAs) is limited. A comprehensive comparison between the agents available commercially (other than in Europe) is lacking, with no data available that most closely reflect the clinic, which is in human whole blood at body temperature.

PURPOSE

To complement the existing literature by determining T₂ relaxivity data for eight GBCAs in vitro.

MATERIAL AND METHODS

The relaxivities of eight GBCAs diluted in human whole blood at 1.5 and 3 T were determined at 37 ± 0.5 °C. Gd was in the range of 0-4 mM. Multi-echo sequences with variable echo times were acquired using a phantom containing a dilution series with each agent, and SigmaPlot 12.0 was used to calculate the R₂ relaxation rate and finally r₂. Statistical comparisons between agents and field strengths were conducted.

RESULTS

The relationship between R₂ vs. Gd was observed to be linear at 1.5 and 3 T, with a mild increase in r₂ from 1.5 to 3 T for all GBCAs. T₂ relaxivity data were compared with prior results. The GBCAs are closely clustered into two groups, with higher r₂ noted for the two lipophilic (those with partial hepatobiliary excretion) compounds.

CONCLUSION

The r₂ values at 1.5 and 3 T, determined for the eight GBCAs still clinically available (other than in Europe), provide a definitive baseline for future evaluations, including theoretical calculations of signal intensity and their clinical impact on T2-weighted scans.

Evaluation of four injection profiles for uniform contrast-enhanced signal intensity profiles in MR angiography

BACKGROUND

Gadolinium concentration variation during acquisition of contrast-enhanced MR angiography (CE-MRA) may lead to artifacts.

PURPOSE

To compare signal intensity (SI) profiles of four different contrast agent injection strategies during CE-MRA with the goal of minimizing SI variation during acquisition.

STUDY TYPE

Prospective.

SUBJECTS

Forty subjects randomized to receive one of four injection profiles of gadobenate dimeglumine (0.1 mmol/kg), either undiluted (0.5 M) or diluted to 40 ml total volume. Tested profiles: 1) nondiluted single-phase ("standard" NS; 1.6 ml/s), 2) diluted single-phase (DS; 1.6 ml/s), 3) diluted biphasic (DB; 9 ml @ 3.3 ml/s, 29 ml @ 1.4 ml/s), 4) patient-tailored protocol using linear prediction (DT).

FIELD STRENGTH/SEQUENCE

Time-resolved SI measured at 3T with spoiled gradient echo sequences having analogous parameters to those of CE-MRA.

ASSESSMENT

Plateau arrival time, rise time, duration, peak and tail SI, plateau quality (sum of squared residuals; SSR), average SI for each injection type derived were used.

STATISTICAL TEST

Two-tailed t-test.

RESULTS

Peak SI, arrival, and rise times were not significantly different between groups, excepting peak SI DB slightly > DS (P = 0.042). Duration of NS vs. the diluted groups was significantly shorter (all P < 0.0001), and DS duration was significantly shorter than that of DT and DB (NS 11.4 ± 3.5 vs. DS 22.9 ± 4.3, DB 25.4 ± 2.3, DT 28.3 ± 4.1 sec). Quality (SSR) of the 20-second plateau was significantly better for DS, DB, DT as compared with NS (all P < 0.001).

DATA CONCLUSION

Three different strategies to power-inject diluted gadobenate dimeglumine targeting a 20-second plateau produced SI profiles with longer duration, more consistent plateau, and no significant loss in peak SI. Such injection profiles may provide more uniform SI during CE-MRA, potentially reducing blurring artifacts.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:1808-1816.

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.

NMR shutter-speed elucidates apparent population inversion of 1 H2 O signals due to active transmembrane water cycling

Purpose

The desire to quantitatively discriminate the extra‐ and intracellular tissue ¹H₂O MR signals has gone hand‐in‐hand with the continual, historic increase in MRI instrument magnetic field strength [B₀]. However, recent studies have indicated extremely valuable, novel metabolic information can be readily accessible at ultra–low B₀. The two signals can be distinguished, and the homeostatic activity of the cell membrane sodium/potassium pump (Na⁺,K⁺,ATPase) detected. The mechanism allowing ¹H₂O MRI to do this is the newly discovered active transmembrane water cycling (AWC) phenomenon, which we found using paramagnetic extracellular contrast agents at clinical B₀ values. AWC is important because Na⁺,K⁺,ATPase can be considered biology’s most vital enzyme, and its in vivo steady‐state activity has not before been measurable, let alone amenable to mapping with high spatial resolution. Recent reports indicate AWC correlates with neuronal firing rate, with malignant tumor metastatic potential, and inversely with cellular reducing equivalent fraction. We wish to systematize the ways AWC can be precisely measured.

Methods

We present a theoretical longitudinal relaxation analysis of considerable scope: it spans the low‐ and high–field situations.

Results

We show the NMR shutter‐speed organizing principle is pivotal in understanding how trans–membrane steady–state water exchange kinetics are manifest throughout the range. Our findings illuminate an aspect, apparent population inversion, which is crucial in understanding ultra‐low field results.

Conclusions

Without an appreciation of apparent population inversion, significant misinterpretations of future data are likely. These could have unfortunate diagnostic consequences.

Relaxivity of Ferumoxytol at 1.5 T and 3.0 T

OBJECTIVES

The aim of this study was to determine the relaxation properties of ferumoxytol, an off-label alternative to gadolinium-based contrast agents, under physiological conditions at 1.5 T and 3.0 T.

MATERIALS AND METHODS

Ferumoxytol was diluted in gradually increasing concentrations (0.26-4.2 mM) in saline, human plasma, and human whole blood. Magnetic resonance relaxometry was performed at 37°C at 1.5 T and 3.0 T. Longitudinal and transverse relaxation rate constants (R1, R2, R2*) were measured as a function of ferumoxytol concentration, and relaxivities (r1, r2, r2*) were calculated.

RESULTS

A linear dependence of R1, R2, and R2* on ferumoxytol concentration was found in saline and plasma with lower R1 values at 3.0 T and similar R2 and R2* values at 1.5 T and 3.0 T (1.5 T: r1saline = 19.9 ± 2.3 smM; r1plasma = 19.0 ± 1.7 smM; r2saline = 60.8 ± 3.8 smM; r2plasma = 64.9 ± 1.8 smM; r2*saline = 60.4 ± 4.7 smM; r2*plasma = 64.4 ± 2.5 smM; 3.0 T: r1saline = 10.0 ± 0.3 smM; r1plasma = 9.5 ± 0.2 smM; r2saline = 62.3 ± 3.7 smM; r2plasma = 65.2 ± 1.8 smM; r2*saline = 57.0 ± 4.7 smM; r2*plasma = 55.7 ± 4.4 smM). The dependence of relaxation rates on concentration in blood was nonlinear. Formulas from second-order polynomial fittings of the relaxation rates were calculated to characterize the relationship between R1blood and R2 blood with ferumoxytol.

CONCLUSIONS

Ferumoxytol demonstrates strong longitudinal and transverse relaxivities. Awareness of the nonlinear relaxation behavior of ferumoxytol in blood is important for ferumoxytol-enhanced magnetic resonance imaging applications and for protocol optimization.

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.

Transverse signal decay under the weak field approximation: Theory and validation

PURPOSE

To derive an expression for the transverse signal time course from systems in the motional narrowing regime, such as water diffusing in blood. This was validated in silico and experimentally with ex vivo blood samples.

METHODS

A closed-form solution (CFS) for transverse signal decay under any train of refocusing pulses was derived using the weak field approximation. The CFS was validated via simulations of water molecules diffusing in the presence of spherical perturbers, with a range of sizes and under various pulse sequences. The CFS was compared with more conventional fits assuming monoexponential decay, including chemical exchange, using ex vivo blood Carr-Purcell-Meiboom-Gill data.

RESULTS

From simulations, the CFS was shown to be valid in the motional narrowing regime and partially into the intermediate dephasing regime, with increased accuracy with increasing Carr-Purcell-Meiboom-Gill refocusing rate. In theoretical calculations of the CFS, fitting for the transverse relaxation rate (R₂ ) gave excellent agreement with the weak field approximation expression for R₂ for Carr-Purcell-Meiboom-Gill sequences, but diverged for free induction decay. These same results were confirmed in the ex vivo analysis.

CONCLUSION

Transverse signal decay in the motional narrowing regime can be accurately described analytically. This theory has applications in areas such as tissue iron imaging, relaxometry of blood, and contrast agent imaging. Magn Reson Med 80:341-350, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Fast, Na+ /K+ pump driven, steady-state transcytolemmal water exchange in neuronal tissue: A study of rat brain cortical cultures

PURPOSE

Water homeostasis and transport play important roles in brain function (e.g., ion homeostasis, neuronal excitability, cell volume regulation, etc.). However, specific mechanisms of water transport across cell membranes in neuronal tissue have not been completely elaborated.

METHODS

The kinetics of transcytolemmal water exchange were measured in neuronal tissue using simultaneous, real-time fluorescence and nuclear magnetic resonance (NMR) measurements of perfused, active brain organotypic cortical cultures. Perfusion with a paramagnetic MRI contrast agent, gadoteridol, allows NMR determination of the unidirectional rate constant for steady-state cellular water efflux (k_io ), and the mole fraction of intracellular water ( pi), related to the average cell volume (V). Changes in intracellular calcium concentration [Cai2+] were used as a proxy for neuronal activity and were monitored by fluorescence imaging.

RESULTS

The k_io value, averaged over all cultures (N = 99) at baseline, was 2.02 (±1.72) s^-1 , indicating that on average, the equivalent of the entire intracellular water volume turns over twice each second. To probe possible molecular pathways, the specific Na⁺ -K⁺ -ATPase (NKA) inhibitor, ouabain (1 mM), was transiently introduced into the perfusate. This caused significant transient changes (N = 8): [Cai2+] rose ∼250%, V rose ∼89%, and k_io fell ∼45%, with a metabolically active k_io contribution probably eliminated by ouabain saturation.

CONCLUSIONS

These results suggest that transcytolemmal water exchange in neuronal tissue involves mechanisms affected by NKA activity as well as passive pathways. The active pathway may account for half of the basal homeostatic water flux. Magn Reson Med 79:3207-3217, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

T1-refBlochi: high resolution 3D post-contrast T1 myocardial mapping based on a single 3D late gadolinium enhancement volume, Bloch equations, and a reference T1

BACKGROUND

High resolution 3D T1 mapping is important for assessment of diffuse myocardial fibrosis in left atrium or other thin-walled structures. In this work, we investigated a fast single-TI 3D high resolution T1 mapping method that directly transforms a 3D late gadolinium enhancement (LGE) volume to a 3D T1 map.

METHODS

The proposed method, T1-refBlochi, is based on Bloch equation modeling of the LGE signal, a single-point calibration, and assumptions that proton density and T2* are relatively uniform in the heart. Several sources of error of this method were analyzed mathematically and with simulations. Imaging was performed in phantoms, eight swine and five patients, comparing T1-refBlochi to a standard spin-echo T1 mapping, 3D multi-TI T1 mapping, and 2D ShMOLLI, respectively.

RESULTS

The method has a good accuracy and adequate precision, even considering various sources of error. In phantoms, over a range of protocols, heart-rates and T1 s, the bias ±1SD was -3 ms ± 9 ms. The porcine studies showed excellent agreement between T1-refBlochi and the multi-TI method (bias ±1SD = -6 ± 22 ms). The proton density and T2* weightings yielded ratios for scar/blood of 0.94 ± 0.01 and for myocardium/blood of 1.03 ± 0.02 in the eight swine, confirming that sufficient uniformity of proton density and T2* weightings exists among heterogeneous tissues of the heart. In the patients, the mean T1 bias ±1SD in myocardium and blood between T1-refBlochi and ShMOLLI was -9 ms ± 21 ms.

CONCLUSION

T1-refBlochi provides a fast single-TI high resolution 3D T1 map of the heart with good accuracy and adequate precision.