Next Article in Journal
Research on the Migration and Adsorption Mechanism Applied to Microplastics in Porous Media: A Review
Next Article in Special Issue
Application of Nanomaterials in Biomedical Imaging and Cancer Therapy II
Previous Article in Journal
The Effect of Surface Oxygen Coverage on the Oxygen Evolution Reaction over a CoFeNiCr High-Entropy Alloy
Previous Article in Special Issue
Magnet-Guided Temozolomide and Ferucarbotran Loaded Nanoparticles to Enhance Therapeutic Efficacy in Glioma Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Magnetic Particle Imaging-Guided Thermal Simulations for Magnetic Particle Hyperthermia

by
Hayden Carlton
1,
Nageshwar Arepally
2,
Sean Healy
1,
Anirudh Sharma
1,
Sarah Ptashnik
3,
Maureen Schickel
3,
Matt Newgren
4,
Patrick Goodwill
4,
Anilchandra Attaluri
2 and
Robert Ivkov
1,5,6,7,*
1
Department of Radiation Oncology and Molecular Radiation Sciences, The Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
2
Department of Mechanical Engineering, School of Science, Engineering, and Technology, The Pennsylvania State University—Harrisburg, Middletown, PA 17057, USA
3
Materialise NV, 3001 Leuven, Belgium
4
Magnetic Insight Inc., Alameda, CA 94502, USA
5
Department of Oncology, Sydney Kimmel Comprehensive Cancer Center, School of Medicine, Johns Hopkins University, Baltimore, MD 21205, USA
6
Department of Mechanical Engineering, Whiting School of Engineering, Johns Hopkins University, Baltimore, MD 21218, USA
7
Department of Materials Science and Engineering, Whiting School of Engineering, Johns Hopkins University, Baltimore, MD 21218, USA
*
Author to whom correspondence should be addressed.
Nanomaterials 2024, 14(12), 1059; https://doi.org/10.3390/nano14121059
Submission received: 20 April 2024 / Revised: 13 June 2024 / Accepted: 17 June 2024 / Published: 20 June 2024

Abstract

:
Magnetic particle hyperthermia (MPH) enables the direct heating of solid tumors with alternating magnetic fields (AMFs). One challenge with MPH is the unknown particle distribution in tissue after injection. Magnetic particle imaging (MPI) can measure the nanoparticle content and distribution in tissue after delivery. The objective of this study was to develop a clinically translatable protocol that incorporates MPI data into finite element calculations for simulating tissue temperatures during MPH. To verify the protocol, we conducted MPH experiments in tumor-bearing mouse cadavers. Five 8–10-week-old female BALB/c mice bearing subcutaneous 4T1 tumors were anesthetized and received intratumor injections of Synomag®-S90 nanoparticles. Immediately following injection, the mice were euthanized and imaged, and the tumors were heated with an AMF. We used the Mimics Innovation Suite to create a 3D mesh of the tumor from micro-computerized tomography data and spatial index MPI to generate a scaled heating function for the heat transfer calculations. The processed imaging data were incorporated into a finite element solver, COMSOL Multiphysics®. The upper and lower bounds of the simulated tumor temperatures for all five cadavers demonstrated agreement with the experimental temperature measurements, thus verifying the protocol. These results demonstrate the utility of MPI to guide predictive thermal calculations for MPH treatment planning.

Graphical Abstract

1. Introduction

Magnetic particle hyperthermia (MPH) is an interstitial thermal therapy approved for treating recurrent glioblastoma with external beam radiation therapy [1]. First proposed in the mid-20th century, clinical MPH involves the intratumor delivery of magnetic nanoparticles (MNPs), which are heated via magnetic hysteresis loss by exposing the region to an alternating magnetic field (AMF) [2]. The magnetocaloric effect has garnered interest for applications in hyperthermia; however, it is prominent in doped manganite magnetic materials and not magnetite/maghemite [3]. For hyperthermia, treatment effectiveness requires the control of thermal energy to heat the tumor to a mild temperature (41–45 °C) for a prescribed time while minimizing the temperature rise in normal tissue. The duration of exposure at an elevated temperature (or “time-at-temperature”) defines the thermal dose [4], which is often expressed by the isoeffect dose metric cumulative equivalent minutes of exposure referenced against 43 °C (CEM43), close to the thermal breakpoint temperature (42.5 °C) of human cells [5,6]. Depending on the thermal dose (typically >15 min), hyperthermia can sensitize cells to radiation [7,8,9] or chemotherapy [10,11], increase blood perfusion to reduce tumor hypoxia [12,13,14], stimulate anti-tumor immune signaling [15,16,17], or it can be directly cytotoxic [18].
While MPH offers substantial advantages for intervention, technical challenges continue to inhibit wider clinical acceptance. Among these, an inability to accurately measure the MNP concentration and its distribution in tissues presents barriers to developing reliable MPH clinical workflows that provide adequate quality assurance measures to compare the delivered thermal doses with the initial prescriptive treatment plan [19]. As sources of heat, the nanoparticle content and distribution determine the thermal dose for MPH; therefore, knowing these in both tumor and surrounding tissues becomes essential when addressing the clinical requirements for reliability and quality in patient care. Previous efforts to quantify the nanoparticles’ distribution have relied on the analyses of computerized tomography (CT) scans [1,20], but there are some limitations. Foremost, the CT signal from the MNPs must be differentiated from the tissue signal; this can prove difficult for intratumor concentrations of less than ~10 g/L [1,19,21].
Magnetic particle imaging (MPI) [22] is a tracer imaging modality that measures time-varying responses of the magnetization vector of a sample and offers advantages for MPH over anatomical imaging. Early efforts demonstrated that the MPI signal correlated with the intratumor nanoparticle content, which can be used to predict a temperature rise when a region is exposed to AMFs [23]. These initial studies were extended to demonstrate how MPI can be used to monitor changes in MNP distribution after MPH [24]. Further advances integrated MPI with an AMF heater and used the MPI gradient field to restrict heating to a particular area within the AMF coil [25,26]. One group developed a computational model for predicting the spatio-thermal resolution of combined MPI/MPH systems [27]. A fully automated prototype AMF heating platform that enables spatially confined heating in a user-selected region of interest has also been described [28]. Buchholz et al. integrated MPI/MPH functionality into a platform that included MPI-based thermometry to monitor safety [29]; the same group also discussed proposed hyperthermia platforms that can integrate with commercial MPI scanners [30].
The technology exists to implement MPI-guided MPH; however, we must be able to incorporate MPI data into predictive thermal simulations to accurately estimate the intratumor temperature. Finite element (FEM) and finite difference methods (FDM) applications are available for other hyperthermia modalities. In those cases, the software enables the user to estimate the thermal dose within specified volumes to create thermal dose contours of the tumor and margins [31,32]. Treatment planning applications for hyperthermia include HyperPlan [33,34,35], SEMCAD X [36,37], and Plan2Heat [38,39]. In some cases, commercial finite element analysis (FEA) applications have been validated for hyperthermia [40]. We have previously validated COMSOL Multiphysics® for preclinical MPH [41]. Lacking in those previous efforts was a direct knowledge of the MNP distribution in tissue or a tissue-mimicking phantom [42]. MNP distribution in tumors and surrounding tissue is subject to considerable individual variability, irrespective of the method used to deliver MNPs [9]. This fact is understudied and often underappreciated, yet it can be a singular source of uncertainty and unreliability in predictive computations of tissue temperature with MPH. It thus becomes essential to develop robust methods to ascertain the MNP tissue content and distribution in each tumor after delivery. Previous work demonstrated integrating MPI data into heat transfer simulations; however, approaches and experimental verification were limited [43,44]. Here, we describe the results of an effort to develop and verify a clinically translatable thermal simulation workflow that uses 3D MPI data as input for finite element computations to predict tumor temperature during a simulated MPH (Figure 1). To process the MPI and micro CT data, we used an imaging-data processing suite that has clinical utility. Verification was achieved by comparing temperatures measured from tumors in mouse cadavers that were heated by AMF-activated MNPs with predicted temperatures of simulated MPHs, where we used experimental temperature data as the initial and boundary conditions for the simulations.

2. Materials and Methods

2.1. Postmortem Animal Studies

2.1.1. Mouse Models

All animal study procedures were conducted according to the protocol approved by the Johns Hopkins University Animal Care and Use Committee (JHU ACUC). The project/protocol identification number is not publicly available due to the JHU institutional policy and can be provided on request from the corresponding author. Five 8–10-week-old female BALB/c mice (Jackson Laboratory, Bar Harbor, ME, USA) were used. All mice were fed a normal diet and water ad libitum, maintained at a 12 h light/12 h dark cycle, and monitored daily for signs of distress or pain.

2.1.2. Cell Line and Tumor Implantation

A vial of the 4T1 murine mammary carcinoma cell line [ER/PR/HER2 negative] was purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA) and maintained according to the supplier’s recommendations. Roswell Park Memorial Institute (RPMI) 1640 media with 10% heat-inactivated fetal bovine serum (FBS) was used to grow the cells. For each mouse, approximately 50,000 cells were suspended in 100 µL of PBS and injected subcutaneously into the right thigh. Once the tumors were palpable, we measured them daily using calipers until they reached a volume between 100–500 mm3.

2.1.3. MNPs

Hydroxyethyl starch-coated Synomag®-S90 iron oxide (mixture of magnetite and maghemite) nanoparticles (Lot #: 14422105–01; micromod Partikeltechnologie GmbH, Rostock, Germany) suspended in water and 50 mg of Fe/mL were used as received [45]. Photon correlation spectroscopy data from the manufacturer indicated a z-average diameter of 109.8 nm with a polydispersity index of 0.092. The particle concentration was verified using a Ferene-s assay [46]. In brief, we first digested the nanoparticles in an acetate buffer with ascorbic acid for at least 20 h. The iron concentration was determined by comparison with reference standards using UV/vis spectrophotometry. The heating rate of the MNPs was estimated using our previously published transient pulse analysis and was reported as specific loss power (SLP) [47]. A sample of nanoparticles at a concentration of 1 mg of Fe/mL H2O was heated at a peak AMF amplitude of 15 mT and at 50% duty (60 s ON/60 s OFF). Each pulse was analyzed and fitted to a non-adiabatic lumped mass heat transfer model, from which the SLP was calculated.

2.1.4. Intratumor MNP Injections

For the MNP dose (i.e., MNP mass) calculation, we used a heat transfer approach to model the MNPs as a spherical uniform heat source [48,49,50] to identify the minimum concentration (c) needed to increase the temperature, ΔT, of a tumor having radius R (Equation (1)):
c = 3 T k S L P × R 2 .
Here, k is the isotropic coefficient of the thermal conductivity of the tumor. Using the measured SLP (496 ± 28 W/g) for the nanoparticle lot used (Figure S4) and assuming idealized conditions, we estimated the minimum MNP concentration required to raise the temperature of a tumor having radius R = 0.5 cm and k = 0.6 W/(m°C) from 37 °C to 43 °C (Equation (2)):
c = 3 × 6 × 0.6 W m 496 W g ×   0.005 m 2 = 0.870 m g F e m L T u m o r .
Then, we adjusted the injection volume and concentration to achieve a target intratumor concentration of ~2 mg Fe/mL tumor to compensate for particle loss and heterogeneous intratumor distribution. Tumor size, for dose determination, was measured with calipers, and the MNP injected volume was adjusted to achieve 2 mg of Fe/mL tumor concentration (Table 1).
Mice were anesthetized via inhalation of 1–2% isoflurane mixed with O2 delivered through a nose cone. The nose cone remained in place for the duration of the MNP injection. MNPs were injected into each tumor at a rate of 2.5 µL/min using a syringe pump (Pump 11 Elite, Harvard Apparatus, Holliston, MA, USA), followed immediately by euthanasia. We decided to perform the experiments with mouse cadavers in order to remove the effects of temperature-dependent perfusion and tissue damage accumulation to better isolate the temperature change induced by the MNPs.

2.1.5. MPI Scanner and Imaging

We used a Momentum® MPI scanner (Magnetic Insight, Inc., Alameda, CA, USA) to generate 3D images of the MNPs within each tumor. The scanner was equipped with a custom mouse holder to limit lateral movement. We used three ~1 µL aliquots of Synomag®-S90 MNPs in microcentrifuge tubes secured to the sample holder within the scanning region as fiducial markers. Each isotropic scan was measured using the ‘Standard mode’ scanner configuration and an excitation field amplitude of 5 mT, with a gradient field of 5.7 T/m and a drive frequency of 45 kHz. Each 3D image consisted of 21 radial scans.

2.1.6. The microCT

We used an IVIS SpectrumCT In Vivo Imaging System (Perkin-Elmer, Shelton, CT, USA) for the anatomical imaging, which was co-registered with the MPI scans. The scans (50 kV at 1 mA) consisted of 720 projections, each with an exposure time of 20 ms. We performed microCT immediately after MPI using the MPI sample holder. A 3D-printed adapter enabled the integration of the MPI holder into the IVIS scanner. We adjusted the scan area to include the mouse cadaver as well as the fiducials within the area of interest.

2.1.7. AMF Heating

All tissue heating was performed on the HYPER device previously described [28]. The selected AMF amplitude for MNP heating by the solenoid was 15 mT ± 10% (12 kA/m) at a set frequency of 341.25 kHz, which is within the Hergt–Dutz biological limit of 5 × 109 A/(m-s) [51]. We verified the amplitude with a 1D magnetic field probe (AMF LifeSystems, LLC, Auburn Hills, MI, USA). The temperature in the sample area encompassed by the solenoid was maintained at 37 ± 0.1 °C by a closed-loop water circulating system.
Before heating the tumors, mouse cadavers were individually placed into plastic bags, which were sealed and immersed into a circulating water bath (Polyscience, Niles, IL, USA) set at 37 °C for 15 min. Mouse cadavers were then removed from the water bath and plastic bags and then placed into a custom 3D-printed sample holder, also maintained at 37.0 ± 0.1 °C by circulated heated water. After securing the cadaver to the sample holder, a single probe was inserted into the tumor to the approximate geometric center. Temperatures were recorded at 1 s intervals. An additional probe was inserted into the rectum to monitor the core body temperature at similar intervals. Each tumor underwent two separate, consecutive heating trials. The first trial consisted of continuous heating at 15 mT for 30 min, or until the tumor exceeded 51 °C. After the trial, the tumor was allowed to cool until 37 °C. When the tumor temperature returned to 37 °C, it was heated at 15 mT at a 67% duty cycle (60 s ON/30 s OFF) for 20 cycles. We performed pulsed heating in addition to continuous heating to mimic the power modulation that may occur during MPH therapy.

2.2. Image Analysis and Computational Modeling

2.2.1. Co-Registration

We used the Mimics Innovation Suite (MIS) Research v.25 (Materialise NV, Leuven, Belgium) with the FEA module for all image analyses, including segmentation, co-registration, conversion to 3D parts, and meshing. While we used the Research version of the software, MIS Medical has FDA 501k medical clearance. The suite consists of two main software: (1) Mimics v26.0®, which we used to create masks and 3D parts from the data, and (2) 3-matic® v18.0 for smoothing and re-meshing. For the MPI-guided simulations, MPI and anatomical imaging data were co-registered to ensure alignment of the shared scale and coordinate system. Mathematically, co-registration transforms each voxel, or 3D pixel, at point x (x, y, z) to point y (x′, y′, z′), as shown in Equation (3):
a 11 a 12 a 13 a 14 a 21 a 22 a 23 a 24 a 31 a 32 a 33 a 34 0 0 0 1 A × x y z 1 x = x y z 1 y .
Within matrix A, the top left 3 × 3 values scale or rotate the image, while the top right 3 × 1 values translate the image. Without this co-registration, the positions of the MNPs (heat generators) within the tumor would be unknown. Practically, we performed the co-registration using the Landmark Registration tool in Mimics®, where the locations of the fiducial markers in 3D space were matched in both the MPI and microCT image stacks (Figure 2). The application automatically generated a transformation matrix (matrix A), which we applied to the MPI data to spatially match the MNP locations within the tumor. The co-registered images for all 5 mice can be found in Figure S1 in the Supplementary Materials.

2.2.2. Imaging Data Calibration

To correlate the MPI signal within a single voxel to MNP thermal output, we prepared serial dilutions of the MNP suspension in 20 µL aliquots in plastic centrifuge tubes having concentrations 0, 0.25, 0.5, 2.5, 5.0, 25, and 50 mg of Fe/mL H2O. The tubes were placed onto the MPI sample holder to ensure separate and distinguishable samples. The calibration samples were then measured with the MPI using the same scanning parameters used to image the tumors. We then used Mimics® to extract the maximum voxel intensity from each calibration sample. The MPI signal intensity of each voxel was mapped to an arbitrary grayscale by the Mimics® software; we will refer to the units simply as a grayscale value (GV). Afterward, we plotted the maximum GV against the estimated volumetric thermal output from each sample (SLP × concentration) and used linear least squares regression to fit the data points (Figure 3). We assumed the SLP did not vary with the concentration or aggregation effects. A representative MPI scan and segmentations of the calibration samples can be found in Figure S2 in the Supplementary Materials. The tabulated calibration curve can also be found in Table S1.

2.2.3. Mesh Generation

We created masks for both the tumor and MNP distribution using the built-in functions within Mimics®. The tumor boundary was segmented, and the created mask was then converted to a .STL file (3D part). We then imported the tumor-part file into 3-matic® (Figure 4a), where the geometry surface was smoothed and re-meshed; a volume mesh was calculated (Figure 4b), and the final mesh was exported as a COMSOL® mesh-part file (Figure 4c), using the built-in functions within 3-matic®. Similarly, the MNP distribution was segmented (Figure 4d), except that the distribution was converted directly to a voxel mesh in Mimics from the image data rather than being converted to a 3D part before meshing (Figure 4e). The MPI mask was created using a threshold segmentation, where the upper bound was the highest GV, and the lower bound was selected manually for each tumor, such that the majority of the signal within the tumor bounds was encompassed. Then, under material assignment, we used the calibration curve to convert the GV of each voxel to units of volumetric thermal output (W/m3), using the samples with the highest and lowest concentrations as references. A material property table was created, which contained the spatial coordinates of each voxel and the corresponding volumetric thermal energy generation (Figure 4f). A more detailed experimental procedure is described in the Supplementary Materials.

2.2.4. FEA Software and Mathematical Models

We used the COMSOL Multiphysics® 6.2 (COMSOL, Burlington, MA, USA) advanced numerical methods software for all heat transfer calculations. The geometry of the computational model is illustrated schematically in Figure 5. The tumor mesh was modeled as a subcutaneous mass, where approximately 1/3 of the tumor volume was embedded within a cuboid geometry representing muscle. The computational phantom tumor and muscle resided within a cylindrical AMF coil, having dimensions similar to that of the experimental coil on the HYPER device.
We used two physics modules available in COMSOL, Heat Transfer in Solids and Magnetic Fields, as well as the Electromagnetic Heating multiphysics module to simulate eddy current heating. The material properties used for the simulation are summarized in Table 2. The Heat Transfer in Solids module implemented the Fourier conduction equation (Equation (4)),
ρ c p 𝜕 T 𝜕 t = k 2 T + Q M N P + Q e d d y
where ρ is the density, cp is the specific heat at constant pressure, k is the coefficient of thermal conductivity, T is the temperature, t is the time, and QMNP and Qeddy are the volumetric thermal power outputs from the hysteresis heating of the MNPs and Joule heating from induced eddy currents, respectively. For QMNP, we imported the material properties table created in Mimics® from the MPI data as a 3D interpolation function, which was used as a spatially variant heat source.
Qeddy was calculated by coupling the frequency domain of Maxwell’s equations (Equations (4)–(7)) with the conduction equation (Equation (4)),
× H = J ,
B = × A ,
J = σ E + j ω D ,
E = j ω A ,
where J is the current density, H is the magnetic field amplitude, B is the magnetic flux density, A is the magnetic vector potential, σ is the electrical conductivity, E is the electric field amplitude, ω is the angular frequency, and D is the electric flux density. As an initial condition, A was assumed zero in all dimensions, and Ampere’s law was applied to the phantom muscle and tumor geometries. The AMF coil was modeled as free space with a magnetic field of amplitude 12 kA/m, applied along the z-axis (along the coil length). Using the equation for Qeddy (Equation (9)), COMSOL solved the Joule heating contribution as follows:
Q e d d y = 0.5 R e J · E + 0.5 R e i ω B · H .
On the top surface of the tumor and muscle, a convective boundary condition (Equation (10)) was implemented using Newton’s law of cooling, with a constant ambient temperature (T) of 37 °C and a heat transfer coefficient (h) of 20 W/(m2·K).
q = h T T ,
where q is the convective heat flux.
A uniform temperature boundary condition was used on the remaining five rectangular faces of the phantom muscle, where the temperature value was interpolated from the experimental cadaver rectal temperature at a given time point. The rectal temperature data can be seen in Figure S3 of the Supplementary Materials. Additionally, we used experimental data for our initial conditions, where the initial temperature of the phantom tumor(s) and muscle geometries corresponded with the respective initial experimental temperatures measured, i.e., at t = 0. A default physics-based “fine” mesh was selected in COMSOL. We performed a frequency-transient study at 341.25 kHz for each phantom tumor to mirror the experimental heating. Specifically, simulated heating was conducted as one continuous heating trial and one pulsed heating trial (60 s ON/30 s OFF). Predicted temperatures from tumor phantoms were reported as a volumetric maximum, minimum, and average of the entire tumor geometry.

2.3. Uncertainty Quantitation

We used the Type A uncertainty evaluation, where valid statistical methods were used to treat our data, as defined by the National Institutes of Standards and Technology [53]. The MPI calibration was performed in triplicate, where the values used to convert the MPI GV to a volumetric thermal output consisted of the average calculated from the 3 replicate measurements. All statistical testing, linear regression, and analyses were performed using Prism 6 software (GraphPad, Boston, MA, USA).

3. Results

All simulated tumor temperature values showed general agreement with the measured intratumor temperatures (Figure 6). For tumors 1, 3, and 4, the calculated average tumor phantom temperatures were ±3 °C of the experimental temperatures measured at the tumor center. Deviations between the simulated and experimental temperatures were observed with tumors 2 and 5, with the latter showing the largest discrepancy. In this case, most of the injected MNPs were concentrated on the left tumor boundary, while the remainder of the tumor (notably the tumor center and location of the probe) was virtually MNP-free. However, the experimental temperature plot for tumor 5 matches the minimum simulated temperature reasonably well, supporting the accuracy of the simulation. We observed significant variances of intratumor MNP and temperature distributions among the tumors after single-point injection, consistent with the previous observations from other tumor models (Figure 6) [9].

4. Discussion

4.1. MNP Distribution and Thermal Probe Placement

The results demonstrate the challenges to achieving a uniform distribution of MNPs within any tumor (constant concentration of MNPs throughout), even with percutaneous delivery, thus highlighting the importance of accounting for the MNP content and distribution in tissues prior to activating them with an AMF. One strategy used to mitigate the effects of unknown and heterogeneous MNP distributions was by manual or automated temperature control mechanisms [42]; however, their success relies on the appropriate placement of the thermal probe within the tumor. Consider tumor 5 (Figure 6). Based solely on the experimental data, we observed temperatures in the range of 40–50 °C for the duration of testing, which is ideal for MPH. Yet, based on the MPI data, the largest concentration of MNPs was along the tumor periphery, away from the position of the probe at the tumor center, resulting in the predicted temperatures being >twofold what was measured experimentally. The MPI data imply that such conditions in a live subject could have produced significant ablative tissue damage at the tumor boundary, assuming a single thermal probe placed in the tumor center as a reference. During treatment, such a disparate probe placement from the concentration of MNPs provides little useful information for accurate thermal dosimetry in the tumor and, worse, fails to anticipate the damage to surrounding normal tissues. MPI, with its enhanced sensitivity to changes in MNP concentration, provides the necessary critical information to locate such potential “hot spot” regions and avoid damaging healthy tissue when compared to works that use anatomical imaging alone.

4.2. Implications for Treatment Planning

Our aim in this work was to develop and verify the capability of MPI to guide thermal simulations using clinically relevant software. To achieve this, we conducted the computational analysis after experimental heating trials, where experimental temperature data was used to provide both boundary and initial conditions for simulations. If simulations are to be used for MPH treatment planning, initial and boundary conditions would be required as initial inputs prior to treatments commencing. Our results here show that reasonable approximations, e.g., subject temperature, could suffice, provided the simulations are used to estimate average, maximum, and minimum temperatures in the target volume (tumor). In the context of treatment planning, these values could be used to generate thermal dose contours to guide AMF power management to ensure treatments achieve an outcome within a prescribed thermal dose range. Certainly, the accuracy of correspondence between planned and achieved treatments would improve if temperature probe placement is performed with imaging guidance. MPI-guided simulations thus can provide the basis of MPH quality control to the researcher or clinician who is able to follow a pre-determined, in silico treatment plan. Additional investigation is needed to ascertain the limits of MPI-guided simulations for treatment planning.

4.3. Limitations

For combined MPI/MPH treatment, the MNPs are both therapeutic and tracer imaging agents; therefore, the selected nanoparticles must generate adequate heat with acceptable MPI tracer quality. For this study, we use commercially available magnetic nanoflowers, Synomag®, that exhibit suitable properties for both MPI and MPH [54]; however, the effect of the MNP dose on MPI image quality bears consideration. The MPI signal varies linearly with the MNP spatial density, i.e., concentration. A high-intensity signal generally produces quality images, but if the concentration is too high, the signal saturates the MPI detector, resulting in image aberrations. Corrective actions can be implemented, such as lowering the excitation field amplitude and increasing the gradient field magnitude, but relatively high MNP concentrations (>1 mg Fe/mL) associated with MPH therapies will likely produce image artifacts when compared to tumor models injected with MNPs at tracer concentrations (~µg Fe/mL).
MPH requires accurate active thermometry to ensure that the prescribed thermal dose is achieved. Further, in vivo thermometry with MPH typically uses invasive probes (usually one) implanted directly within the tumor to measure the intratumor temperature. From a heat transfer perspective, single-point thermometry introduces severe limitations, especially for larger tumors, as seen with tumor 5 (Figure 6). To draw useful conclusions from a single temperature probe, one must assume the tumor is a lumped mass with a uniform spatial temperature distribution, which is only true for sufficiently small tumors that are uniformly loaded with MNPs. Future development of multi-point or volumetric thermometry would enhance the accuracy of MPI-guided simulations of MPH and reduce the associated uncertainties.

5. Conclusions

We developed and verified an MPI-guided thermal simulation workflow for use in MPH preclinical experiments. We conclude that our methodology successfully integrated MPI data into FEA software for thermal simulations. Verification of the proposed methodology was achieved by comparing the simulation results with corresponding experimentally heated tumor-bearing mouse cadavers that had been previously injected with MNPs. The general agreement between the experimental and simulated results raises confidence that this approach can overcome many of the technological limitations encountered with clinical MPH. Achieving further in vivo preclinical validation of this methodology will establish its potential for use in human patients.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano14121059/s1.

Author Contributions

Conceptualization: H.C., N.A., S.H. and R.I.; Resources: H.C., S.P., M.S., M.N., P.G. and R.I.; Funding acquisition: R.I.; Methodology: H.C., N.A., S.H., A.S., S.P., M.S., M.N., P.G. and A.A.; Investigation: H.C.; Data curation: H.C.; Formal analysis: H.C.; Validation: H.C., N.A., S.H., A.A. and R.I.; Visualization: H.C.; Project administration: A.A. and R.I.; Supervision: R.I.; Writing—original draft: H.C. and R.I.; Writing—review and editing: H.C., N.A., S.H., A.S., S.P., M.S., M.N., P.G., A.A. and R.I. All authors have read and agreed to the published version of the manuscript.

Funding

Funding for the experimental portion of this project was provided by the National Cancer Institute of the National Institutes of Health under Award Numbers 1R01 CA257557, 1R01 CA247290, 1S10 OD026740, 5P30 CA006973, 1R44CA285064 and R37 CA229417. The content is solely the responsibility of the authors and does not necessarily represent the official views of Johns Hopkins University, Pennsylvania State University, or the National Institutes of Health.

Institutional Review Board Statement

No human subjects were used. The animal study procedures were conducted according to the Johns Hopkins Institutional Animal Care and Use Committee (IACUC) approved protocol. The protocol/project identification number is not publicly available due to the JHU institutional policy and can be provided on request from the corresponding author.

Data Availability Statement

Data are publicly available on the Johns Hopkins Data Repository.

Conflicts of Interest

M.N. and P.G. are employees of Magnetic Insight, Inc., a company that manufactures and sells MPI scanners and technology used in this study. S.P. and M.S. are employees of Materialise, the company that manufactures and sells the imaging data processing software used in this study. R.I. is an inventor listed on several nanoparticle patents. All patents are assigned to either The Johns Hopkins University or Aduro Biosciences, Inc. All other authors report no other conflicts of interest.

References

  1. Maier-Hauff, K.; Ulrich, F.; Nestler, D.; Niehoff, H.; Wust, P.; Thiesen, B.; Orawa, H.; Budach, V.; Jordan, A. Efficacy and safety of intratumoral thermotherapy using magnetic iron-oxide nanoparticles combined with external beam radiotherapy on patients with recurrent glioblastoma multiforme. J. Neuro-Oncol. 2011, 103, 317–324. [Google Scholar] [CrossRef]
  2. Gilchrist, R.K.; Medal, R.; Shorey, W.D.; Hanselman, R.C.; Parrott, J.C.; Taylor, C.B. Selective inductive heating of lymph nodes. Ann. Surg. 1957, 146, 596–606. [Google Scholar] [CrossRef] [PubMed]
  3. Rmili, N.; Riahi, K.; M’Nassri, R.; Ouertani, B.; Cheikhrouhou-Koubaa, W.; Hlil, E.K. Magnetocaloric and induction heating characteristics of La0.71Sr0.29Mn0.95Fe0.05O3 nanoparticles. J. Sol-Gel Sci. Technol. 2024. ahead of print. [Google Scholar] [CrossRef]
  4. Dewhirst, M.W.; Viglianti, B.L.; Lora-Michiels, M.; Hanson, M.; Hoopes, P.J. Basic principles of thermal dosimetry and thermal thresholds for tissue damage from hyperthermia. Int. J. Hyperth. 2003, 19, 267–294. [Google Scholar] [CrossRef] [PubMed]
  5. Sapareto, S.A.; Dewey, W.C. Thermal dose determination in cancer-therapy. Int. J. Radiat. Oncol. Biol. Phys. 1984, 10, 787–800. [Google Scholar] [CrossRef] [PubMed]
  6. Oleson, J.R.; Samulski, T.V.; Leopold, K.A.; Clegg, S.T.; Dewhirst, M.W.; Dodge, R.K.; George, S.L. Sensitivity of hyperthermia trial outcomes to temperature and time—Implications for thermal goals of treatment. Int. J. Radiat. Oncol. Biol. Phys. 1993, 25, 289–297. [Google Scholar] [CrossRef] [PubMed]
  7. Wust, P.; Hildebrandt, B.; Sreenivasa, G.; Rau, B.; Gellermann, J.; Riess, H.; Felix, R.; Schlag, P.M. Hyperthermia in combined treatment of cancer. Lancet Oncol. 2002, 3, 487–497. [Google Scholar] [CrossRef] [PubMed]
  8. van der Zee, J.; van Rhoon, G.C. Hyperthermia is effective in improving clinical radiotherapy results. Int. J. Radiat. Oncol. Biol. Phys. 2006, 66, 633–634. [Google Scholar] [CrossRef]
  9. Attaluri, A.; Kandala, S.K.; Wabler, M.; Zhou, H.; Cornejo, C.; Armour, M.; Hedayati, M.; Zhang, Y.; DeWeese, T.L.; Herman, C.; et al. Magnetic nanoparticle hyperthermia enhances radiation therapy: A study in mouse models of human prostate cancer. Int. J. Hyperth. 2015, 31, 359–374. [Google Scholar] [CrossRef]
  10. Tohnai, I.; Goto, Y.; Hayashi, Y.; Ueda, M.; Kobayashi, T.; Matsui, M. Preoperative thermochemotherapy of oral cancer using magnetic induction hyperthermia (implant heating system: IHS). Int. J. Hyperth. 1996, 12, 37–47. [Google Scholar] [CrossRef]
  11. Otsuka, T.; Yonezawa, M.; Kamiyama, F.; Matsushita, Y.; Matsui, N. Results of surgery and radio-hyperthermo-chemotherapy for patients with soft-tissue sarcoma. Int. J. Clin. Oncol. 2001, 6, 253–258. [Google Scholar] [CrossRef] [PubMed]
  12. Iwata, K.; Shakil, A.; Hur, W.J.; Makepeace, C.M.; Griffin, R.J.; Song, C.W. Tumour pO(2) can be increased markedly by mild hyperthermia. Br. J. Cancer 1996, 74, S217–S221. [Google Scholar]
  13. Brizel, D.M.; Scully, S.P.; Harrelson, J.M.; Layfield, L.J.; Dodge, R.K.; Charles, H.C.; Samulski, T.V.; Prosnitz, L.R.; Dewhirst, M.W. Radiation therapy and hyperthermia improve the oxygenation of human soft tissue sarcomas. Cancer Res. 1996, 56, 5347–5350. [Google Scholar] [PubMed]
  14. Elming, P.B.; Sorensen, B.S.; Oei, A.L.; Franken, N.A.P.; Crezee, J.; Overgaard, J.; Horsman, M.R. Hyperthermia: The Optimal Treatment to Overcome Radiation Resistant Hypoxia. Cancers 2019, 11, 60. [Google Scholar] [CrossRef] [PubMed]
  15. Zanker, K.S.; Lange, J. Whole-body hyperthermia and natural-killer cell-activity. Lancet 1982, 1, 1079–1080. [Google Scholar] [CrossRef] [PubMed]
  16. Yan, X.Y.; Xiu, F.M.; An, H.Z.; Wang, X.J.; Wang, J.L.; Cao, X.T. Fever range temperature promotes TLR4 expression and signaling in dendritic cells. Life Sci. 2007, 80, 307–313. [Google Scholar] [CrossRef] [PubMed]
  17. Dayanc, B.E.; Beachy, S.H.; Ostberg, J.R.; Repasky, E.A. Dissecting the role of hyperthermia in natural killer cell mediated anti-tumor responses. Int. J. Hyperth. 2008, 24, 41–56. [Google Scholar] [CrossRef] [PubMed]
  18. Roizintowle, L.; Pirro, J.P. The response of human and rodent cells to hyperthermia. Int. J. Radiat. Oncol. Biol. Phys. 1991, 20, 751–756. [Google Scholar] [CrossRef] [PubMed]
  19. Johannsen, M.; Gneveckow, U.; Eckelt, L.; Feussner, A.; Waldöfner, N.; Scholz, R.; Deger, S.; Wust, P.; Loening, S.A.; Jordan, A. Clinical hyperthermia of prostate cancer using magnetic nanoparticles:: Presentation of a new interstitial technique. Int. J. Hyperth. 2005, 21, 637–647. [Google Scholar] [CrossRef]
  20. LeBrun, A.; Manuchehrabadi, N.; Attaluri, A.; Wang, F.; Ma, R.H.; Zhu, L. MicroCT image-generated tumour geometry and SAR distribution for tumour temperature elevation simulations in magnetic nanoparticle hyperthermia. Int. J. Hyperth. 2013, 29, 730–738. [Google Scholar] [CrossRef]
  21. Gneveckow, U.; Jordan, A.; Scholz, R.; Brüss, V.; Waldöfner, N.; Ricke, J.; Feussner, A.; Hildebrandt, B.; Rau, B.; Wust, P. Description and characterization of the novel hyperthermia- and thermoablation-system MFH®300F for clinical magnetic fluid hyperthermia. Med. Phys. 2004, 31, 1444–1451. [Google Scholar] [CrossRef]
  22. Gleich, B.; Weizenecker, R. Tomographic imaging using the nonlinear response of magnetic particles. Nature 2005, 435, 1214–1217. [Google Scholar] [CrossRef]
  23. Murase, K.; Aoki, M.; Banura, N.; Nishimoto, K.; Mimura, A.; Kuboyabu, T.; Yabata, I. Usefulness of Magnetic Particle Imaging for Predicting the Therapeutic Effect of Magnetic Hyperthermia. Open J. Med. Imaging 2015, 5, 85–99. [Google Scholar] [CrossRef]
  24. Kuboyabu, T.; Yabata, I.; Aoki, M.; Banura, N.; Nishimoto, K.; Mimura, A.; Murase, K. Magnetic Particle Imaging for Magnetic Hyperthermia Treatment: Visualization and Quantification of the Intratumoral Distribution and Temporal Change of Magnetic Nanoparticles in Vivo. Open J. Med. Imaging 2016, 6, 1–15. [Google Scholar] [CrossRef]
  25. Tay, Z.W.; Chandrasekharan, P.; Chiu-Lam, A.; Hensley, D.W.; Dhavalikar, R.; Zhou, X.Y.; Conolly, S.M. Magnetic Particle Imaging-Guided Heating in Vivo Using Gradient Fields for Arbitrary Localization of Magnetic Hyperthermia Therapy. ACS Nano 2018, 12, 3699–3713. [Google Scholar] [CrossRef] [PubMed]
  26. Lei, S.; He, J.; Huang, X.; Hui, H.; An, Y.; Tian, J. A Novel Local Magnetic Fluid Hyperthermia Based on High Gradient Field Guided by Magnetic Particle Imaging. IEEE Trans. Biomed. Eng. 2024. ahead of print. [Google Scholar] [CrossRef]
  27. Le, T.A.; Hadadian, Y.; Yoon, J. A prediction model for magnetic particle imaging-based magnetic hyperthermia applied to a brain tumor model. Comput. Methods Programs Biomed. 2023, 235, 107546. [Google Scholar] [CrossRef]
  28. Carlton, H.; Weber, M.; Peters, M.; Arepally, N.; Lad, Y.S.; Jaswal, A.; Ivkov, R.; Attaluri, A.; Goodwill, P. HYPER: Pre-clinical device for spatially-confined magnetic particle hyperthermia. Int. J. Hyperth. 2023, 40, 2272067. [Google Scholar] [CrossRef]
  29. Buchholz, O.; Sajjamark, K.; Franke, J.; Wei, H.; Behrends, A.; Münkel, C.; Hofmann, U.G. In situ theranostic platform combining highly localized magnetic fluid hyperthermia, magnetic particle imaging, and thermometry in 3D. Theranostics 2024, 14, 324–340. [Google Scholar] [CrossRef]
  30. Behrends, A.; Wei, H.; Neumann, A.; Friedrich, T.; Bakenecker, A.C.; Franke, J.; Buzug, T.M. Integrable Magnetic Fluid Hyperthermia Systems for 3D Magnetic Particle Imaging. Nanotheranostics 2024, 8, 163–178. [Google Scholar] [CrossRef]
  31. Paulides, M.M.; Stauffer, P.R.; Neufeld, E.; Maccarini, P.F.; Kyriakou, A.; Canters, R.A.M.; Diederich, C.J.; Bakker, J.F.; Van Rhoon, G.C. Simulation techniques in hyperthermia treatment planning. Int. J. Hyperth. 2013, 29, 346–357. [Google Scholar] [CrossRef]
  32. Kok, H.P.; Crezee, J. Hyperthermia Treatment Planning: Clinical Application and Ongoing Developments. IEEE J. Electromagn. RF Microw. Med. Biol. 2021, 5, 214–222. [Google Scholar] [CrossRef]
  33. Stalling, D.; Seebass, M.; Zöckler, M.; Hege, H.-C. Hyperthermia Treatment Planning with HyperPlan—User’s Manual; Publication Server of Zuse Institute Berlin (ZIB): Berlin, Germany, 2000. [Google Scholar]
  34. Gellermann, J.; Wust, P.; Stalling, D.; Seebass, M.; Nadobny, J.; Beck, R.; Hege, H.C.; Deuflhard, P.; Felix, R. Clinical evaluation and verification of the hyperthermia treatment planning system HyperPlan. Int. J. Radiat. Oncol. Biol. Phys. 2000, 47, 1145–1156. [Google Scholar] [CrossRef]
  35. Sreenivasa, G.; Gellermann, J.; Rau, B.; Nadobny, J.; Schlag, P.; Deuflhard, P.; Felix, R.; Wust, P. Clinical use of the hyperthermia treatment planning system hyperplan to predict effectiveness and toxicity. Int. J. Radiat. Oncol. Biol. Phys. 2003, 55, 407–419. [Google Scholar] [CrossRef]
  36. Schmid, G.; Überbacher, R.; Samaras, T.; Tschabitscher, M.; Mazal, P.R. The dielectric properties of human pineal gland tissue and RF absorption due to wireless communication devices in the frequency range 400–1850 MHz. Phys. Med. Biol. 2007, 52, 5457–5468. [Google Scholar] [CrossRef] [PubMed]
  37. De Bruijne, M.; Wielheesen, D.H.M.; Van der Zee, J.; Chavannes, N.; Van Rhoon, G.C. Benefits of superficial hyperthermia treatment planning: Five case studies. Int. J. Hyperth. 2007, 23, 417–429. [Google Scholar] [CrossRef] [PubMed]
  38. Kok, H.P.; Kotte, A.; Crezee, J. Planning, optimisation and evaluation of hyperthermia treatments. Int. J. Hyperth. 2017, 33, 593–607. [Google Scholar] [CrossRef]
  39. Kok, H.P.; Crezee, J. Validation and practical use of Plan2Heat hyperthermia treatment planning for capacitive heating. Int. J. Hyperth. 2022, 39, 952–966. [Google Scholar] [CrossRef] [PubMed]
  40. Chen, X.; Diederich, C.J.; Wootton, J.H.; Pouliot, J.; Hsu, I.C. Optimisation-based thermal treatment planning for catheter-based ultrasound hyperthermia. Int. J. Hyperth. 2010, 26, 39–55. [Google Scholar] [CrossRef]
  41. Kandala, S.K.; Sharma, A.; Mirpour, S.; Liapi, E.; Ivkov, R.; Attaluri, A. Validation of a coupled electromagnetic and thermal model for estimating temperatures during magnetic nanoparticle hyperthermia. Int. J. Hyperth. 2021, 38, 611–622. [Google Scholar] [CrossRef]
  42. Kandala, S.K.; Liapi, E.; Whitcomb, L.L.; Attaluri, A.; Ivkov, R. Temperature-controlled power modulation compensates for heterogeneous nanoparticle distributions: A computational optimization analysis for magnetic hyperthermia. Int. J. Hyperth. 2019, 36, 115–129. [Google Scholar] [CrossRef]
  43. Banura, N.; Mimura, A.; Nishimoto, K.; Murase, K. Heat transfer simulation for optimization and treatment planning of magnetic hyperthermia using magnetic particle imaging. arXiv 2016, arXiv:1605.08139. [Google Scholar]
  44. Tang, Y.D.; Chen, M.; Flesch, R.C.C.; Jin, T. Extraction method of nanoparticles concentration distribution from magnetic particle image and its application in thermal damage of magnetic hyperthermia. Chin. Phys. B 2023, 32, 094401. [Google Scholar] [CrossRef]
  45. Bender, P.; Fock, J.; Frandsen, C.; Hansen, M.F.; Balceris, C.; Ludwig, F.; Johansson, C. Relating Magnetic Properties and High Hyperthermia Performance of Iron Oxide Nanoflowers. J. Phys. Chem. C 2018, 122, 3068–3077. [Google Scholar] [CrossRef]
  46. Hedayati, M.; Abubaker-Sharif, B.; Khattab, M.; Razavi, A.; Mohammed, I.; Nejad, A.; Ivkov, R. An optimised spectrophotometric assay for convenient and accurate quantitation of intracellular iron from iron oxide nanoparticles. Int. J. Hyperth. 2018, 34, 373–381. [Google Scholar] [CrossRef] [PubMed]
  47. Carlton, H.; Ivkov, R. A new method to measure magnetic nanoparticle heating efficiency in non-adiabatic systems using transient pulse analysis. J. Appl. Phys. 2023, 133, 044302. [Google Scholar] [CrossRef] [PubMed]
  48. Hergt, R.; Andra, W.; d’Ambly, C.G.; Hilger, I.; Kaiser, W.A.; Richter, U.; Schmidt, H.G. Physical limits of hyperthermia using magnetite fine particles. IEEE Trans. Magn. 1998, 34, 3745–3754. [Google Scholar] [CrossRef]
  49. Andrä, W.; d’Ambly, C.G.; Hergt, R.; Hilger, I.; Kaiser, W.A. Temperature distribution as function of time around a small spherical heat source of local magnetic hyperthermia. J. Magn. Magn. Mater. 1999, 194, 197–203. [Google Scholar] [CrossRef]
  50. Dutz, S.; Hergt, R. Magnetic particle hyperthermia-a promising tumour therapy? Nanotechnology 2014, 25, 28. [Google Scholar] [CrossRef]
  51. Hergt, R.; Dutz, S. Magnetic particle hyperthermia-biophysical limitations of a visionary tumour therapy. J. Magn. Magn. Mater. 2007, 311, 187–192. [Google Scholar] [CrossRef]
  52. Nagy, J.A.; DiDonato, C.J.; Rutkove, S.B.; Sanchez, B. Permittivity of ex vivo healthy and diseased murine skeletal muscle from 10 kHz to 1 MHz. Sci. Data 2019, 6, 37. [Google Scholar] [CrossRef] [PubMed]
  53. National Institutes of Standards and Technology. Guidelines for Evaluating and Expressing the Uncertainty of NIST Measurement Results. Available online: http://physics.nist.gov/TN1297 (accessed on 15 October 2023).
  54. Karpavicius, A.; Coene, A.; Bender, P.; Leliaert, J. Advanced analysis of magnetic nanoflower measurements to leverage their use in biomedicine. Nanoscale Adv. 2021, 3, 1633–1645. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic representation of workflow for magnetic particle imaging (MPI)-guided thermal simulations. Magnetic nanoparticles (MNPs) at a pre-determined dose were injected intratumorally. After MNP injection, the tumor was imaged with both MPI and microCT for anatomical reference. Using fiducial markers, the MPI data were co-registered with anatomical images to establish a common datum. The tumor geometry was segmented, converted to a 3D part, and volume-meshed and imported into FEA software (COMSOL Multiphysics 6.2); additionally, voxel intensities obtained from the MPI scan were converted into volumetric thermal output values from a calibration curve developed using measured MNP heating data. With the MPI values spatially registered and calibrated for heat output, the FEA simulation was performed.
Figure 1. Schematic representation of workflow for magnetic particle imaging (MPI)-guided thermal simulations. Magnetic nanoparticles (MNPs) at a pre-determined dose were injected intratumorally. After MNP injection, the tumor was imaged with both MPI and microCT for anatomical reference. Using fiducial markers, the MPI data were co-registered with anatomical images to establish a common datum. The tumor geometry was segmented, converted to a 3D part, and volume-meshed and imported into FEA software (COMSOL Multiphysics 6.2); additionally, voxel intensities obtained from the MPI scan were converted into volumetric thermal output values from a calibration curve developed using measured MNP heating data. With the MPI values spatially registered and calibrated for heat output, the FEA simulation was performed.
Nanomaterials 14 01059 g001
Figure 2. Co-registration is necessary to correlate MPI voxel intensity values with 3D coordinates within the segmented tumor (Mouse 1 is shown as an example). The micro(μ)CT and MPI of tumored mice after injection are shown, where the tumor is circled in white. The micro(μ)CT and MPI scans were co-registered within Mimics, and a transformation matrix was created. The transformation matrix obtained for each tumor was applied to the MPI scan and overlaid onto the micro(μ)CT scan.
Figure 2. Co-registration is necessary to correlate MPI voxel intensity values with 3D coordinates within the segmented tumor (Mouse 1 is shown as an example). The micro(μ)CT and MPI of tumored mice after injection are shown, where the tumor is circled in white. The micro(μ)CT and MPI scans were co-registered within Mimics, and a transformation matrix was created. The transformation matrix obtained for each tumor was applied to the MPI scan and overlaid onto the micro(μ)CT scan.
Nanomaterials 14 01059 g002
Figure 3. MNP calibration is essential to generate accurate MPI-guided simulations. The max grayscale and thermal generation values (calculated from SLP) were fitted to a linear equation using least squares regression to yield y-intercept of 1031 ± 4 GV and slope = (1.10 ± 0.04) × 10−5 W·m−3·GV−1 with R2 = 0.9746. Uncertainty is represented by 95% confidence bounds, which are shown as dashed lines in the figure.
Figure 3. MNP calibration is essential to generate accurate MPI-guided simulations. The max grayscale and thermal generation values (calculated from SLP) were fitted to a linear equation using least squares regression to yield y-intercept of 1031 ± 4 GV and slope = (1.10 ± 0.04) × 10−5 W·m−3·GV−1 with R2 = 0.9746. Uncertainty is represented by 95% confidence bounds, which are shown as dashed lines in the figure.
Nanomaterials 14 01059 g003
Figure 4. After co-registration and segmentation, the represented MNP distribution and tumor geometry (Mouse 3 shown as an example) were prepared for COMSOL. (a) After converting the tumor mask to a .stl file, it was imported into 3-matic, where it was (b) smoothed, re-meshed, and (c) converted to a COMSOL mesh file. (d) The MPI scan was segmented, and the resulting mask was used to represent the MNP distribution in the computational tumor phantom. (e) The voxel intensities were mapped onto the values of the MPI calibration curve using the Mimics® material assignment tab. (f) The calibrated voxel values were then imported into COMSOL as a 3-argument interpolation function.
Figure 4. After co-registration and segmentation, the represented MNP distribution and tumor geometry (Mouse 3 shown as an example) were prepared for COMSOL. (a) After converting the tumor mask to a .stl file, it was imported into 3-matic, where it was (b) smoothed, re-meshed, and (c) converted to a COMSOL mesh file. (d) The MPI scan was segmented, and the resulting mask was used to represent the MNP distribution in the computational tumor phantom. (e) The voxel intensities were mapped onto the values of the MPI calibration curve using the Mimics® material assignment tab. (f) The calibrated voxel values were then imported into COMSOL as a 3-argument interpolation function.
Nanomaterials 14 01059 g004
Figure 5. Schematic diagram of computational tumor and muscle phantom used to simulate MPH with COMSOL. The imported tumor geometry was embedded in a phantom body represented by a cuboid having properties of muscle, suspended in a phantom AMF coil with a diameter similar to that of the coil in the HYPER device. Constant temperature boundary conditions were assumed, having values obtained from experimental rectal temperatures, on the bottom and sides of the phantom body. We imposed a convective boundary condition on the top of the phantom tumor and muscle, with an assumed heat transfer coefficient and ambient temperature.
Figure 5. Schematic diagram of computational tumor and muscle phantom used to simulate MPH with COMSOL. The imported tumor geometry was embedded in a phantom body represented by a cuboid having properties of muscle, suspended in a phantom AMF coil with a diameter similar to that of the coil in the HYPER device. Constant temperature boundary conditions were assumed, having values obtained from experimental rectal temperatures, on the bottom and sides of the phantom body. We imposed a convective boundary condition on the top of the phantom tumor and muscle, with an assumed heat transfer coefficient and ambient temperature.
Nanomaterials 14 01059 g005
Figure 6. MPI-guided simulations of magnetic hyperthermia accurately predict intratumor temperatures measured experimentally in the center of murine 4T1 tumors. The left column shows the distribution of MNPs within the tumor from the MPI data that were used as inputs for simulated MPH in computational phantoms. The second column shows the results of the continuous heating trials, where the gray region encompasses the maximum and minimum predicted temperature values obtained from the computational phantoms, the dashed lines show the average simulation temperature, and the solid blue lines are the experimentally measured temperatures from the 4T1 tumor heating trials in mouse cadavers. The third column shows similar results, which were obtained from the pulsed heating trials.
Figure 6. MPI-guided simulations of magnetic hyperthermia accurately predict intratumor temperatures measured experimentally in the center of murine 4T1 tumors. The left column shows the distribution of MNPs within the tumor from the MPI data that were used as inputs for simulated MPH in computational phantoms. The second column shows the results of the continuous heating trials, where the gray region encompasses the maximum and minimum predicted temperature values obtained from the computational phantoms, the dashed lines show the average simulation temperature, and the solid blue lines are the experimentally measured temperatures from the 4T1 tumor heating trials in mouse cadavers. The third column shows similar results, which were obtained from the pulsed heating trials.
Nanomaterials 14 01059 g006
Table 1. Tumor sizes were measured with calipers and associated MNP injection volume.
Table 1. Tumor sizes were measured with calipers and associated MNP injection volume.
SampleCaliper Measured Tumor Volume (mm3)MNP Injection Volume (µL) @ 50 mg Fe/mL
Tumor 11445.6
Tumor 21937.7
Tumor 328011.2
Tumor 435514.2
Tumor 545518.2
Table 2. Material properties used for COMSOL electromagnetic heating simulations.
Table 2. Material properties used for COMSOL electromagnetic heating simulations.
PropertyMuscleTumorRef.
Specific Heat at Constant Pressure (cp)3421 J/(kg-K)3760 J/(kg-K)[43]
Density (ρ)1090 kg/m31045 kg/m3[43]
Thermal Conductivity (k)0.49 W/(m-K)0.51 W/(m-K)[43]
Relative Permeability (µr)11
Relative Permittivity (εr)20002000[52]
Electrical Conductivity (σ)0.23 S/m0.23 S/m[52]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Carlton, H.; Arepally, N.; Healy, S.; Sharma, A.; Ptashnik, S.; Schickel, M.; Newgren, M.; Goodwill, P.; Attaluri, A.; Ivkov, R. Magnetic Particle Imaging-Guided Thermal Simulations for Magnetic Particle Hyperthermia. Nanomaterials 2024, 14, 1059. https://doi.org/10.3390/nano14121059

AMA Style

Carlton H, Arepally N, Healy S, Sharma A, Ptashnik S, Schickel M, Newgren M, Goodwill P, Attaluri A, Ivkov R. Magnetic Particle Imaging-Guided Thermal Simulations for Magnetic Particle Hyperthermia. Nanomaterials. 2024; 14(12):1059. https://doi.org/10.3390/nano14121059

Chicago/Turabian Style

Carlton, Hayden, Nageshwar Arepally, Sean Healy, Anirudh Sharma, Sarah Ptashnik, Maureen Schickel, Matt Newgren, Patrick Goodwill, Anilchandra Attaluri, and Robert Ivkov. 2024. "Magnetic Particle Imaging-Guided Thermal Simulations for Magnetic Particle Hyperthermia" Nanomaterials 14, no. 12: 1059. https://doi.org/10.3390/nano14121059

APA Style

Carlton, H., Arepally, N., Healy, S., Sharma, A., Ptashnik, S., Schickel, M., Newgren, M., Goodwill, P., Attaluri, A., & Ivkov, R. (2024). Magnetic Particle Imaging-Guided Thermal Simulations for Magnetic Particle Hyperthermia. Nanomaterials, 14(12), 1059. https://doi.org/10.3390/nano14121059

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