Materials: Carboxylic acid-coated iron oxide nanoparticles were purchased from Sigma Aldrich and MK Nano. Specifically, the 4.9 nm, 8.5 nm and 15.7 nm particles were acquired from Sigma Aldrich (the labels report diameters of 5 nm, 10 nm and 20 nm, respectively). The 12.9 nm particles were acquired from MK Nano (listed as 15 nm particles). Agarose (BioReagent for molecular biology, low electroendosmosis) was also purchased from Sigma Aldrich.
MRI Sample Preparation: Samples were prepared by diluting contrast agents in an agarose gel (prepared with mass fraction of 1 % agarose in H2O) based on the concentrations listed on the manufacturers label. Samples were prepared with concentrations of 0.03 mmol/L, 0.06 mmol/L, 0.12 mmol/L, 0.25 mmol/L, 0.5 mmol/L and 1 mmol/L of metal (Fe for the SPIONs and Gd for Gd-BOPTA).
Transmission Electron Microscopy (TEM): TEM was performed using a Tecnai T12 Spirit BT microscope with a LaB6 filament. Samples were prepared by diluting nanoparticles in a 1:1 volume ratio mixture of isopropyl alcohol and H2O. The diluted samples were added dropwise to the top of Formvar/carbon-coated copper grids purchased from Ted Pella. Physical core sizes of particles were analyzed using TEM images and ImageJ46.
3 Tesla MRI: An Agilent preclinical scanner was used for 3 T MRI. T1 measurements were made using an inversion recovery sequence with a 256 x 256 matrix and a 128 mm x 128 mm field of view. Images were acquired from 6 axial slices with a thickness of 2 mm and a gap spacing of 4 mm. The inversion times used for the sequence were 50 ms, 75 ms, 100 ms, 125 ms, 250 ms, 500 ms, 1000 ms, 1500 ms, 2000 ms, 3000 ms and 6000 ms. The repetition time was 10000 ms and the echo time was 13.92 ms. T2 measurements were performed using a spin echo sequence with the same resolution, field of view and slice parameters as the T1 measurements. Echo times used for T2 acquisition were 14 ms, 28 ms, 56 ms, 112 ms, and 224 ms. The repetition time was 10000 ms. Measurements were made using a temperature-controlled phantom. Temperature control was achieved using a closed flow loop of a perfluorocarbon coolant and monitored by a fiber optic sensor. The temperature was set to 21.5° C in order to match the lab temperature of the 64 mT MRI measurements.
64 mT MRI: A Hyperfine Swoop scanner with hardware version 1.8 and software version rc8.3.1, was used to acquire T1 and T2 measurements at 64 mT. All scans used an 8-channel receive, 1-channel transmit head coil. The. T1 measurements were made using a research version of the Hyperfine proprietary T1-weighted inversion recovery 3D fast spin echo (FSE) sequence with a 22 cm x 18 cm x 18 cm field of view, an in-plane resolution of 1.6 mm x 1.6 mm, and a slice thickness of 5 mm. The inversion times used for the sequence were 100 ms, 200 ms, 300 ms, 400 ms, 500 ms, 600 ms, 700 ms, 800 ms, 900 ms, 1100 ms, 1300 ms, 1500 ms, 1800 ms, 2100 ms, 2500 ms. The repetition time was 3000 ms and the echo time was 5.96 ms. T2 measurements were performed using a research version of the Hyperfine proprietary T2-weighted 3D FSE sequence with a 22 cm x 18 cm x 18 cm field of view, an in-plane resolution of 1.5 mm x 1.5 mm, and a slice thickness of 5 mm. Echo times used for T2 acquisition were 37 ms, 111 ms, 184 ms, 259 ms, 333 ms, 407 ms, 480 ms, 554 ms, 628 ms, 702 ms. The repetition time was 3000 ms. Measurements were performed at an ambient lab temperature of 21.5 °C.
Concentration Measurements: The Fe concentration of each SPION-based contrast agent was measured using inductively coupled plasma optical emission spectroscopy (ICP-OES). ICP-OES was performed using a Perkin Elmer Optima 8300 ICP-OES optical system with a segmented-array charge-coupled device detector. Briefly, samples were digested using nitric acid, then further diluted to generate samples for measurement with ICP-OES. A detailed description of sample preparation is given in the Supplemental Information.
T1/T2 Mapping and Relaxivity Analysis: Samples were prepared with nominal concentrations (diluted according to concentration on manufacturer’s label) of 0 mmol/L, 0.03 mmol/L, 0.06 mmol/L, 0.12 mmol/L, 0.25 mmol/L, 0.5 mmol/L and 1 mmol/L. For measurements at 3 T, the 0 mmol/L sample was not measured because the 3 T sample holder accommodates one fewer sample than 64 mT. The concentration of the SPIONs were measured using inductively coupled plasma optical emission spectroscopy and the nominal concentrations were rescaled to reflect the actual concentration of the samples (see Supplemental Information for more details). Samples were prepared by embedding contrast agents in an agarose solution (prepared with mass fraction of 1 % agarose in H2O), which was allowed to set into a semi-solid gel at room temperature. An agarose medium was chosen for imaging because it provided long term colloidal stability, so that the same samples could be scanned at 64 mT and 3 T. Diluting the samples in other solutions did not confer the same degree of colloidal stability (see Supplemental Information). T1 was calculated for each voxel using LMFIT in Python for the inversion-recovery model,
with T1 the target value to fit, inversion time TI, repetition time TR, scale factor for imperfect inversion d, the nominal signal intensity for a voxel S0, and measured signal intensity Si. At 64 mT, T2 maps were acquired using a research version of the Hyperfine proprietary T2-weighted 3D FSE sequence. The T2 map is calculated in the Hyperfine protocol using SciPy optimize curve_fit in Python for the model,
with T2 the target value for the fit, echo time TE, nominal signal intensity without relaxation S0, and measured signal intensity Si. At 3 T, the T2 maps were found by using LMFIT for each voxel using Equation 3.
After T1 and T2 were measured for all concentrations and field strengths, the data was fit using Eqn. 1 to calculate the relaxivities. Unfortunately, the experimental protocol at 64 mT could not accurately measure T1 less than the first inversion time of 100 ms. At high concentrations (0.25 L mmol-1 s-1 – 1 L mmol-1 s-1), some of the SPION-based agents with larger physical core sizes had T1 values lower than 100 ms. So, we instituted a selection rule to exclude cases where the concentration of contrast agents was too high for accurate determination of measurement times. The rule states that if the T1 value calculated using a fit of Eqn. 2 is less than the first inversion time (100 ms for 64 mT), then the data was excluded from the relaxivity fit (see Supplemental Information for more details). The data corresponding to these concentrations was then excluded from fitting the relaxivities at other conditions (r2 at 64 mT, r1 at 3 T and r2 at 3 T) to maintain consistency across all field strengths and relaxation types. Out of 42 measurements (6 concentrations for 7 different agents), this rule eliminated 11 data sets from relaxivity analysis and fitting. The data sets that were eliminated corresponded to samples with large nanoparticle sizes at high concentrations, since these tended to relax signal very rapidly. A table has been included in the Supplemental Information to indicate the samples that were excluded on this basis.
SQUID Magnetometry: The magnetization of samples was measured using a Quantum Design SQUID MPMS 3 magnetometer. Samples were immobilized in a water-soluble hydrogel and cured with UV-light to prevent magnetic-field induced chaining and aggregation of particles during measurement.
NMRD: Water relaxation profiles were acquired with a Stelar Spinmaster FFC2000-1 T relaxometer by measuring the water proton relaxation rates as a function of the applied magnetic field (0.01 MHz – 40 MHz proton Larmor frequency). The field-cycling technique is used to detect the magnetization decay/recovery curves by acquiring the free induction decay signals from the sample after exposure to a given magnetic field for 16 different intervals of time47. By changing the magnetic field, longitudinal relaxation rates can be determined at all frequencies permitted by the instrument. For magnetic fields lower than 0.35T, a 1 T pre-polarization field is applied before the measurement to generate sufficiently large signals for NMRD acquisition. The relaxation measurements, obtained from the fit of the magnetization decay/recovery curves against a mono-exponential function, were affected by an error below ±1 %. Samples were prepared with a nominal metal concentration of 1 mmol/L (using concentration specified by manufacturer) in an agarose gel (prepared with mass fraction of 1 % agarose in H2O). The relaxivity profiles were obtained by normalization of the measured relaxation data to the metal concentration (measured by ICP-OES) after subtracting the diamagnetic relaxation rate contribution. The measurements were performed at 21.5 °C for comparison with the MRI data and 37 °C to evaluate relaxivity at physiological temperature. Additional measurements were performed at 15 °C and 25 °C and are included in the Supplemental Information.