Subscriber access provided by University of South Florida
Communication
Assessment of Gold Nano-particle mediated Enhanced Hyperthermia using MR Guided High-Intensity Focused Ultrasound Ablation Procedure Surendra Balaji Devarakonda, Matthew R Myers, Matthew Lanier, Charles Dumoulin, and Rupak Kumar Banerjee Nano Lett., Just Accepted Manuscript • DOI: 10.1021/acs.nanolett.7b00272 • Publication Date (Web): 13 Mar 2017 Downloaded from http://pubs.acs.org on March 13, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Nano Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1 2 3
Nano Letters
Assessment of Gold Nano-particle mediated Enhanced Hyperthermia using MR Guided High-Intensity Focused Ultrasound Ablation Procedure
4 5 6 7
Surendra B. Devarakonda1, Matthew R. Myers2, Mathew Lanier3, Charles Dumoulin3, and Rupak K. Banerjee1*
8 9
Department of Mechanical, Materials Engineering1
10
College of Engineering and Applied Science,
11
University of Cincinnati,
12
Cincinnati, OH 45221
13
Tel: 513-556-2124; Fax: 513-556-3390; email:
[email protected] 14
Division of Solid and Fluid Mechanics2
15
Center for Devices and Radiological Health,
16
U.S. Food and Drug Administration,
17
Silver Spring, MD 20993
18
Department of Radiology3
19
Cincinnati Children's Hospital Medical Center,
20
Cincinnati, Ohio, USA
21 22 23 24
Keywords: High Intensity Focused Ultrasound; HIFU; Hyperthermia; Gold Nano-particles; Cancer therapy, MRI
25 1 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
26 27
Abstract High-intensity focused ultrasound (HIFU) has gained increasing popularity as a non-invasive
28
therapeutic procedure to treat solid tumors. However, collateral damage due to the use of high
29
acoustic powers during HIFU procedures remains a challenge. The objective of this
30
study is to assess the utility of using gold nano-particles (gNPs) during HIFU procedures to
31
locally enhance heating at low powers, thereby reducing the likelihood of collateral damage.
32
Phantoms containing tissue-mimicking material (TMM) and physiologically relevant
33
concentrations (0%, 0.0625%, and 0.125%) of gNPs were fabricated. Sonications at acoustic
34
powers of 10 W, 15 W, and 20 W were performed for a duration of 16 sec using an MR-HIFU
35
system. Temperature rises and lesion volumes were calculated and compared for phantoms with
36
and without gNPs. For an acoustic power of 10 W, the maximum temperature rise increased by
37
32% and 43% for gNPs concentrations of 0.0625% and 0.125%, respectively, when compared to
38
the 0% gNPs concentration. For the power of 15 W, a lesion volume of 0 mm3, 44.5 ± 7 mm3,
39
and 63.4 ± 32 mm3 was calculated for the gNPs concentration of 0%, 0.0625%, and 0.125%,
40
respectively. For a power of 20 W, it was found that the lesion volume doubled and tripled for
41
concentrations of 0.0625% and 0.125% gNPs, respectively, when compared to the concentration
42
of 0% gNPs. We conclude that gNPs have the potential to locally enhance the heating and reduce
43
damage to healthy tissue during tumor ablation using HIFU.
44
Keywords: High Intensity Focused Ultrasound; HIFU; Hyperthermia; Gold Nano-particles;
45
Cancer therapy; MRI
46
2 ACS Paragon Plus Environment
Page 2 of 24
Page 3 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
47
In high-intensity focused ultrasound (HIFU) ablation procedures, ultrasound energy from an
48
external transducer is focused onto target tissue within the body. HIFU therapy has increasingly
49
gained clinical interest as a novel tool for ablation procedures. Using focused ultrasound energy,
50
temperatures greater than 60°C can be induced in tissues essentially instantaneously, achieving
51
cellular necrosis and vascular cauterization1-3. A large tumor area can be treated by placing
52
multiple focal zones sequentially and in close proximity until the entire area is ablated3.
53
Characterization and testing of HIFU procedures has been performed using computational
54
methods4-8, in-vitro tissue phantoms, and animal models9-14.
55
Furusawa et al.15, performed MR-guided HIFU on patients with breast cancer and were able
56
to successfully ablate tumors using acoustic powers of greater than 60 W. However, they
57
reported third degree skin burns in 3% of the patients. Li et al.16, performed HIFU ablations on
58
17 human patients with recurrent and metastatic abdominal tumors to determine the
59
complications arising during and after the HIFU treatment. They used focal peak intensities of
60
2000 to 5000 Wcm-2and the treatment time was between 3120 sec to 8950 sec. Skin burns were
61
found in all 17 patients and mild enteroparalysis was found in 15 patients. These adverse events
62
are possibly due to the usage of higher acoustic power. Although HIFU is experiencing higher
63
clinical application, issues related to high acoustic intensity, such as skin burns and damage to
64
neighboring healthy tissue17,18 remain.
65
To reduce the required acoustic intensity and consequently the likelihood of adverse events
66
such as skin burns, absorption-enhancing agents can be used. One promising agent is
67
nanoparticles, which can be injected to the tumor site to achieve targeted and intense heating of
68
the tumor cells 17. Sun et al.18, 19, reported the effects of magnetic microcapsules (PLGA-coated
69
Fe3O4) during HIFU therapy, and found enhancement of hyperthermia due to the microcapsules. 3 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
70
However, the authors themselves pointed out that the microcapsules size was very high (500 –
71
800 nm). The power used was extremely high (180 W - 250 W) which can lead to skin burns
72
and collateral damage to neighboring tissue17. Also, the MR thermometry images reported in this
73
study have limited clarity regarding the temperature increase in the tumor.
74
Using 1% and 3% concentrations of magnetic nanoparticles (mNPs) embedded in tissue
75
mimicking material (TMM), Dibaji et al.20, measured the HIFU-induced temperature rise using
76
embedded thermocouples (TCs). They determined that the peak temperature rise increased by 1.6
77
and 2 times when mNPs concentration of 1% and 3% were used, respectively, for an acoustic
78
power of 14.2 W. Temperature changes in HIFU can also be measured using temperature
79
sensitive MR. However, the magnetism and magnetic susceptibility of mNPs have been
80
reported to interfere with Magnetic Resonance Imaging (MRI). This results in decreased signal
81
to noise ratio, geometric distortions and low quality images21. It may be noted that the mNPs
82
concentration of 1% and 3% for the pilot study was relatively high from the perspective of
83
physiological relevance.
84
To help assess the efficacy of nanoparticle introduction pre-clinically, it is useful to make
85
thermal measurements in tissue phantoms. HIFU temperature measurements can be performed
86
using invasive modalities such as TCs, or by non-invasive modalities such as magnetic resonance
87
(MR) thermometry. Direct measurement of temperature using TCs can lead to errors such as
88
viscous-heating22, 23 artifacts and positioning errors24 caused by misalignment between the HIFU
89
beam and the TCs. Dasgupta et al.23, have conducted HIFU experiments in phantoms using
90
powers of 5 W, 10.3 W, 17.3 W, and 24.8 W and measured the temperature rise using TCs.
91
They have reported that the thermal artifact can be up to 2 times the local tissue temperature rise
92
when direct sonications were performed on TCs. Alternatively, with the use of MR thermometry, 4 ACS Paragon Plus Environment
Page 4 of 24
Page 5 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
93
thermocouple artifacts can be avoided25-27. By optimizing the MR thermometry parameters such
94
as T1 and T2 relaxation times, proton resonance frequency (PRF), and the temperature-dependent
95
diffusion coefficient, temperature measurement can be achieved during HIFU sonication with
96
reasonable accuracy28-32.
97
In the past, researchers from our lab have performed MR thermometry on in-vivo porcine
98
livers using acoustic powers of 10 W, 30 W, and 40 W20. This study reported that MR
99
thermometry can be used during HIFU sonication for acquiring the maximum temperature rise of
100
57 °C with localized cavitation. On a similar note, MR-guided HIFU procedures33 were
101
performed in-vitro in turkey breast samples and in-vivo in transcranial rat brains with acoustic
102
powers of 40 W. In this study, tissue displacements due to acoustic force were detected.
103
In the present study, the effect of gold nanoparticles (gNPs) on the HIFU-induced
104
temperature rise and lesion volume have been assessed using MR thermometry. To the authors’
105
knowledge, quantification of temperature rise and lesion volume using MR thermometry with
106
nano-particles, for example, gNPs has not been previously reported. The TMM phantoms with
107
physiologically relevant concentrations of 0% (control), 0.0625%, and 0.125% gNPs by volume
108
were fabricated. Low acoustic powers - 10 W, 15 W, and 20 W - were used for the HIFU
109
sonications. The temperature rises and the lesion volumes were calculated from the MR
110
thermometry and were compared for different concentrations of gNPs.
111 112
Three cylindrical fixtures with a length of 3 cm and inner diameter of 2.5 cm (volume ~ 15
113
cm3) were developed (Fig. 1A). The protocol reported by King et al.34 was used to fabricate the
114
tissue phantom with 0% gNPs concentration. To prepare the tissue phantoms with 0.0625% and
115
0.125% gNPs, the method explained in Dibaji et al.20 was used. The size of the gNPs embedded 5 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
116
in the TMM was 15 nm (US Research Nanomaterials, Inc., Houston, TX). Thorough mixing
117
using a magnetic stir bar was performed for 2 hours to achieve uniform distribution of the gNPs
118
inside the TMM during the fabrication of the phantoms21. All the fabricated phantoms were kept
119
for 12 hours at room temperature to ensure complete solidification of TMM.
120
A clinical MR-HIFU system (SonalleveTM V2, Philips Medical Systems, Vantaa, Finland),
121
integrated into a 1.5 Tesla (T) whole-body scanner (Philips IngeniaTM, Healthcare, Best, The
122
Netherlands) was used for scanning the phantom (Fig. 1B). The MR-HIFU system contains a
123
256-element phased array HIFU transducer which can be used to focus the ultrasound energy
124
into small volumes within the phantom. Radius of curvature and the aperture of the transducer
125
array are 120 mm and 130 mm, respectively35. The operating frequency was 1.2 MHz. Focal
126
spot of the HIFU beam was approximately 2 mm in the radial direction.
127
For targeting and sonicating inside the TMM, phantoms were placed above the HIFU
128
transducer on the MR-HIFU tabletop in a vertical position. Scans were performed to ensure
129
proper acoustic coupling and for detecting air bubbles. Any air bubbles detected near the surface
130
of the TMM were removed. The phantoms were placed inside a 75 mm x 75 mm single loop
131
tunnel and matched imaging coil to achieve an optimized resolution. A resolution of 0.7 mm x
132
0.7 mm in the radial direction and a slice thickness of 6 mm were achieved.
133
The electrical powers used by the MR-HIFU system to generate the desired acoustic powers
134
were derived from radiation force balance measurements. In order to optimize the conversion of
135
electrical field to acoustic waves, an automated electrical matching routine was performed by the
136
system after placing the phantom inside the imaging coil.
137 138
After aligning the phantoms inside the imaging coil, MR images using T1-weighted and T2weighted sequences were obtained. These images were obtained in both sagittal and coronal 6 ACS Paragon Plus Environment
Page 6 of 24
Page 7 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
139
planes. Using the T1-weighted images, it was possible to avoid any obstacles in the path of the
140
HIFU beam such as the plastic fixture of the phantoms. The imaging technique was two
141
dimensional (2D) “fast field echo” (FFE) and the flip angle was 20°. The repetition time (TR)
142
and the echo time (TE) were 53 msec and 23 msec, respectively. Additionally, for the T2-
143
weighted sequence, the TR, TE, and flip angle were 1300 msec, 130 msec, and 90°, respectively.
144
The HIFU beam was oriented inside the tissue phantom in such a way that the outer plastic
145
fixture was not in the path of incoming or outgoing beam. Three acoustic powers of 10 W, 15 W,
146
and 20 W were used for the HIFU sonication. A period of 16 sec was chosen for the HIFU
147
sonication of the phantoms.
148
During the sonication and cooling periods of 16 sec and 15 sec, respectively, MR
149
thermometry was used to acquire the temperature data. The initial temperature of the phantoms
150
and the surrounding medium was fixed at 37 °C. Using a fast field echo (FFE) and echo planar
151
technique (EPI), the temperature maps at the focus in both radial and axial direction were
152
captured. The temperature changes inside the phantom were calculated using the proton
153
resonance frequency-shift (PRFS) method. The optimized temporal resolution of the
154
temperature measurements was 6 sec. The maximum temperature values discussed in the
155
manuscript are the volumetric average of all the temperature values inside a pixel. Temperature
156
measurements were recorded from three trials (n=3) for each acoustic power (m=3).
157
For cases where temperature values were desired outside the 6-second intervals, the
158
temperature values obtained from the MR thermometry were curve fit to a function representing
159
a solution to the heat equation. The functional form of the fitting function is exponential-
160
integral36. The exponential-integral solution to the heat equation was derived assuming a heat
161
source with a Gaussian profile in radial direction37. The heat source is assumed to be constant in 7 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 24
162
the axial direction over the distances of interest. The solution is applicable to cases where the
163
radial dimension of the heat source is much smaller than the axial dimension, as is the case with
164
the HIFU beam. The temperature distribution as a function of time at a distance r from the beam
165
axis is:
166
167
, =
−
!
(1a)
where Ei(x) is the exponential integral: ( $ %&
" = − #)
168
'
*+
(1b)
169
where r is the radial coordinate, r0 is the beam radius width of the Gaussian intensity distribution,
170
α is the absorption coefficient (dB/m), ,- is the intensity on the beam axis (W/mm3), κ is the
171
thermal diffusivity (m2/s), ρ0 is the density (kg/m3), and cp is the specific heat (J/kg.K) of TMM.
172
The properties of the TMM are presented in Table 1.
173
In order to calculate the lesion volume, the MR-HIFU algorithm calculates the thermal dose
174
using the method developed by Sapareto and Dewey 38. The thermal dose parameter is given by
175
3
. ", /, 0 = #3- 45678 1 .2 *
(2)
176
where . is the thermal dose at the reference temperature of 43 ºC, 9:; M 0.25 IℎKL+K
(3)
180
Thermal dose contours calculated by the treatment planning algorithm for the temperature maps
181
were used in determining the region that has been subjected to a thermal dose of at least 240
182
equivalent minutes. This volume is defined as the lesion volume. Line contours in the middle 8 ACS Paragon Plus Environment
Page 9 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
183
zone of the temperature maps were measured using the MR software across the focal diameter
184
(d) and along the beam axial path (l) at the focus. Symmetry was assumed in the radial direction
185
as the properties of the TMM did not vary in the radial direction. These dimensions were then
186
used in equation 2, which calculates the lesion volume by assuming an ellipsoid shape for the
187
lesions inside the phantom:
N = P* Q O
188
(4)
189
Here d (mm) is the focal diameter of the lesion in radial direction, l (mm) is the length of the
190
lesion in axial direction. Lesions were measured from three trials (n=3) for each acoustic power
191
(m=3). Uncertainties in lesion volume represent the standard deviation of the volumes computed
192
from the three measurements.
193
In addition to temperature rise and lesion volume, absorbed ultrasound energy due to the
194
presence of gNPs was also calculated. This was computed from the initial temperature values
195
calculated from the MR thermometry, as heat conduction can be neglected for initial duration.
196
Thus, temperature rise is related to the absorbed energy Q by:
197
U2
R- ST U = V = 2 W ,
(5)
198
The ratio of the absorbed ultrasound energy with and without gNPs was computed for the
199
different powers considered. Statistical analysis has been performed on the maximum
200
temperature rise and lesion volume data using a t-test to determine if the enhancement of heating
201
by gNPs is statistically significant.
202
At each power level (10 W, 15 W, 20 W), the measured temperatures were averaged over
203
three trials (n=3). Results are presented as mean ± SD. The temperature maps (Fig. 2) show the 9 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 24
204
variation of temperature for the concentrations of 0%, 0.0625%, and 0.125% gNPs and for the
205
powers of 10 W, 15 W, and 20 W at the plane of maximum pixel temperature. A trending
206
increase in temperature (increase in number of red pixels) can be seen in moving from left to
207
right (increasing gNP) concentration, particularly at the two lower powers. Figure 3 shows the
208
temperature traces derived using the exponential-integral function (1) fitted using the MR
209
thermometry values. The temperature rise at 7 seconds (the second MR data value) for the
210
0.0625% gNP concentration was about 1.1 times that for no gNP’s, and for 0.125% gNP
211
concentration the ratio was about 1.2. The acoustic power was 10 W.
212
increase in temperature rise at the end of 7 seconds of sonication was nearly 1.2 times and 1.6
213
times for the concentration of 0.0625% and 0.125% gNPs, respectively, when compared to 0%
214
gNPs. Finally, for a power of 20 W, the increase in temperature rise for the initial 7 seconds of
215
sonication was nearly 1.2 times and 1.5 times for the concentration of 0.0625% and 0.125%
216
gNPs, respectively, when compared to 0% gNPs.
For a power of 15 W, the
217
Figure 4A shows the comparison of maximum temperature rise, which occurs at the end of
218
the 13-second sonication time. For a power of 10 W, the maximum temperature rises of 16±1
219
°C, 21±3 °C, and 23±3 °C were observed for the concentrations of 0%, 0.0625%, and 0.125%
220
gNPs, respectively. The changes relative to 0% gNP’s were 32% and 43% for the two phantoms
221
with gNP’s. For a power of 15 W, the maximum temperature rises of 25±0.7 °C, 29±4 °C, and
222
36±5 °C were observed for the concentrations of 0%, 0.0625%, and 0.125% gNPs. The
223
increases relative to baseline (0% gNP) were 16% and 43% for the 0.0625% and 0.125% gNP
224
values. For a power of 20 W, the maximum temperature rises of 34±3 °C (0% gNP), 43±0.4 °C
225
(0.0625% gNP), and 46±0.7 °C (0.125% gNP). The increase relative to baseline for the 0.0625%
10 ACS Paragon Plus Environment
Page 11 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
226
concentration was 24% and for the 0.125% concentration was 34%. It is also evident that the
227
enhancement of heating by gNPs is statistically significant (p < 0.05).
228
The amount of absorbed ultrasound energy, determined using an initial time of 1 sec for
229
slope calculation, is shown in Fig. 4B. For the power of 10 W, the relative energy absorption
230
was similar for all the concentrations. It is possible that the temperature rise for the time period
231
of 1 sec is very small in relation to 6 sec, particularly for the lower power of 10 W. Therefore,
232
the small temperature rise may have remained undetected by the MRI sensor. Consequently, the
233
change in the slope of temperature rise between no gNPs (0%) and with gNPs (00625% and
234
0.125%) for the power of 10 W is not apparent. For the power of 15 W, approximately 2 times
235
more energy was absorbed with the 0.0625% gNPs concentration than with no gNPs, and 2.7
236
times more energy relative to control for the 0.125% concentration. For the 20 W power level,
237
the relative energy absorption was 4 for the 0.0625% gNPs concentration and 5.7 for the 0.125%
238
gNPs concentration.
239
Figure 4C presents the lesion volumes obtained for the 0%, 0.0625%, and 0.125% gNPs
240
concentrations for powers of 10 W, 15 W, and 20 W. The white line contours in the middle zone
241
of the temperature maps (Fig 2) indicate the predicted lesion area at that cross-section. The
242
lesion area has been exposed to at least a thermal dose of 240 equivalent minutes41. To calculate
243
the predicted lesion volumes, the white lines from the temperature maps in both radial and axial
244
directions have been measured.
245
For the concentration of 0% gNPs, at powers of 10 W (Fig 2A1) and 15 W (Fig 2B1),
246
insufficient heat was generated to produce a predicted lesion. For the power of 15 W, a predicted
247
lesion volume of 44.5 ± 7 mm3 and 63.4 ± 32 mm3 was calculated for the gNPs concentration of
248
0.0625% and 0.125%, respectively. For 20 W, the lesion volume increased roughly two and three 11 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
249
times when gNPs concentration of 0.0625% and 0.125% were used compared to 0% gNPs
250
concentration.
251 252
For the same power, the lesion volume increased nearly 3 fold with the addition of gNP’s
253
(Fig. 4C). Thus, it can be observed from Fig. 4A and 4C that the enhancement of heating by
254
gNPs (0.0625% and 0.125%) is statistically significant (p < 0.05) when compared to 0% gNPs at
255
13 sec time point. Through the use of lower acoustic powers (e.g. 15 W) and shorter sonication
256
times afforded by optimized concentrations of gNP’s, undesired damage can potentially be
257
reduced in a clinical setting. Currently, researchers have used acoustic powers in excess of 60 W
258
- 250 W to achieve cell necrosis during HIFU procedures15.
259
The enhancement of ablation efficiency does not appear to be linear with gNP concentration.
260
For example, for a power of 15 W, the relative absorbed energy (Fig. 4B) is 2.0 for a gNP
261
concentration of 0.0625%, but doubling this concentration increases the absorbed energy to only
262
2.7. A similar trend holds at 20 W. Likewise, a 0.0625% gNP concentration produces a lesion
263
volume of 44 mm3, but doubling the gNP concentration increases the lesion volume by a factor
264
of 1.4 (power = 15 W). At 20 W, doubling the gNP concentration also leads to a factor of 1.4
265
increase in lesion volume. Based upon the current data of our study, enhancement in ablation
266
efficiency scales roughly as the square root of the gNP concentration.
267
Temperature rise in the presence of gNP’s is a linear function of power, as it is for no gNP’s
268
(Fig. 3A). In that sense, the effect of incorporating gNP’s in tissue volume can be modeled as a
269
local increase in acoustic absorption.
12 ACS Paragon Plus Environment
Page 12 of 24
Page 13 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
270
The assumption of an ellipsoid shape for a lesion at higher powers may lead to some error in
271
calculation of lesion volume. This is because the ellipsoid shape is known to get distorted due to
272
cavitation at higher energy deposition, although no cavitation was detected in the present
273
experiments.
274
The temperatures obtained from the MR thermometry are an average of all the temperatures
275
with in a pixel. To estimate the effect of averaging, the temperature distribution was assumed to
276
possess the same functional form as the intensity distribution. Such an approximation would be
277
most accurate at the beginning of sonication, prior to significant heat conduction. For the
278
calculations, the intensity profile provided by Wu and Du39 (based on a Gaussian-beam
279
assumption) was used. Averaging over the 6 mm pixel thickness in the axial direction results in
280
a temperature error of approximately 4.2% at the focus. Away from the focus the error is less, as
281
a change in temperature in the +z direction (relative to the center of the pixel) is partially
282
compensated by a change of the opposite sign in the –z direction. In the radial direction, the 0.7
283
mm resolution results in an averaging error of approximately 1%.
284
The HIFU absorption in media embedded with gNPs depend on the thermal processes within
285
the viscous and phonon layers at the interface of gNPs as well as on the intrinsic absorption
286
properties of the media40-43. Propagating HIFU waves in a medium interact with the thermal
287
phonons and consequently, a part of the wave is absorbed. Due to the wave absorption, there is
288
an increase in the momentum of thermal phonons inside the medium leading to temperature rise.
289
It is expected that the attenuation due to phonon layer is the dominating mechanism for the gNPs
290
size of 15 nm leading to temperature rise, but further analysis is required. Such analysis may
291
need an independent theoretical-experimental characterization study that we have embarked on.
13 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Gold nanoparticles can be made to accumulate at a desired location in the tumor tissue by
292 293
injecting the nanoparticles directly into the tumor tissue intravenously. The gNPs are expected
294
to accumulate in the tumor tissue due to the enhanced permeability and retention effect (EPR)44-
295
46
296
performed using nanoparticles injected into the tumor site intravenously. Also, localized
297
injection of nanoparticles could be performed in superficial tumors. MR thermometry can then
298
be performed similar to the work presented in this study to evaluate the influence of gNPs.
. To demonstrate this technique, future in-vivo experiments, preferably in animals, can be
299 300
Corresponding Author
301
Rupak K. Banerjee
302
Department of Mechanical, Materials Engineering1
303
College of Engineering and Applied Science,
304
University of Cincinnati,
305
Cincinnati, OH 45221
306
Tel: 513-556-2124; Fax: 513-556-3390; email:
[email protected] 307
Author contributions
308
Surendra B. Devarakonda performed the experiments and wrote the manuscript. Rupak Banerjee
309
supervised the project. Mathew R. Myers, Mathew Lanier, Charles Dumoulin and Rupak K.
310
Banerjee reviewed and revised the manuscript.
311
Competing financial interests
312
The authors declare no competing financial interests.
313
Funding Sources
14 ACS Paragon Plus Environment
Page 14 of 24
Page 15 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
314
This work was supported by NSF grant 1403356.
315
Acknowledgements
316
This work was supported by NSF grant 1403356.
15 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
317
Page 16 of 24
Table 1: Properties of the TMM. Property
Value
Density, R-
1040 kg/m3
Absorption Coefficient, α Specific Heat, Cp
45 dB/m 4064 J/kg.K 1.4 x 10-7 m2/s
Thermal Diffusivity, κ 318 319
16 ACS Paragon Plus Environment
Page 17 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
A
320
B
321 322
Figure 1: A) MR image of the phantom filled with tissue-mimicking material (TMM) (diameter
323
of 25 mm) and B) MR-HIFU setup - vertical cross-sectional view of the TMM within MR coil
324
aligned with HIFU transducer. 17 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
325 326
Figure 2: Representative temperature maps of TMM phantoms at maximum temperature with
327
0%, 0.0625%, and 0.125% gNPs concentration for the acoustic powers of 10 W, 15 W, and 20
328
W. The line contours in the center shows calculated lesion volume.
329
330
18 ACS Paragon Plus Environment
Page 18 of 24
Page 19 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
331 332 333 334 335
Figure 23: Comparison of temporal variation of temperature rise (°C) using MR thermometry for the concentrations of 0%, 0.0625%, and 0.125% gNPs and powers of A) 10 W, B) 15 W, and C) 20 W. The sonication period was 16 sec. Measurements for each time point was conducted in triplicate (n=3). 19 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
336 337 338 339 340 341 342
Figure 4: A) Comparison of maximum temperature rise (°C) for 0%, 0.0625%, and 0.125% gNPs concentration for powers of 10 W, 15 W, and 20 W, B) The ratio of the initial slope of the temperature trace, normalized by the slope for the case of 0% gNPs at each selected (10 W, 15 W, and 20 W) power, and C) Comparison of lesion volume (mm3) for 0%, 0.0625%, and 0.125% gNPs concentration for the acoustic powers of 15 W and 20 W. 20 ACS Paragon Plus Environment
Page 20 of 24
Page 21 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Nano Letters
343 344
Table of contents Graphic
345
21 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
346
References:
347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391
1. ter Haar, G. Ultrasound. Med. Biol. 1995, 21, 1089-100. 2. Hill, C. R.; ter Haar, G. R. Br. J. Radiol. 1995, 68, 1296-1303. 3. Gelet, A.; Chapelon, J. Y.; Bouvier, R.; Pangaud, C.; Lasne, Y. J. Urol. 1999, 161, 156-62. 4. Wu, J. R.; Du, G. H. Ultrasound. Med. Biol. 1990, 16, 489-98. 5. Curra, F. P.; Mourad, P. D.; Khokhlova, V. A.; Cleveland, R. O.; Crum, L. A. IEEE Trans. Ultrason., Ferroelectr., Freq. Control 2000, 47, 1077-89. 6. Soneson, J. E. AIP Conf. Proc. 2009, 1113, 165-169. 7. Myers, M. R.; Soneson, J. E. J. Acoust. Soc. Am. 2009, 126, 425-33. 8. Dibaji, S. A. R.; Banerjee, R. K.; Soneson, J. E.; Myers, M. R. J. Acoust. Soc. Am. 2013, 134, 34353445. 9. Bailey, M. R.; Couret, L. N.; Sapozhnikov, O. A.; Khokhlova, V. A.; ter Haar, G.; Vaezy, S.; Shi, X.; Martin, R.; Crum, L. A. Ultrasound. Med. Biol. 2001, 27, 695-708. 10. McLaughlan, J.; Rivens, I.; Leighton, T.; Ter Haar, G. Ultrasound. Med. Biol. 2010, 36, 1327-44. 11. Kyriakou, Z.; Corral-Baques, M. I.; Amat, A.; Coussios, C. C. Ultrasound. Med. Biol. 2011, 37, 56879. 12. Dasgupta, S.; Das, P.; Wansapura, J.; Hariharan, P.; Pratt, R.; Witte, D.; Myers, M. R.; Banerjee, R. K. J. Nanotechnol. Eng. Med 2011, 2, 024501-024501. 13. Mesiwala, A. H.; Farrell, L.; Wenzel, H. J.; Silbergeld, D. L.; Crum, L. A.; Winn, H. R.; Mourad, P. D. Ultrasound. Med. Biol. 2002, 28, 389-400. 14. Solomon, S. B.; Nicol, T. L.; Chan, D. Y.; Fjield, T.; Fried, N.; Kavoussi, L. R. Invest. Radiol. 2003, 38, 293-301. 15. Furusawa, H.; Namba, K.; Thomsen, S.; Akiyama, F.; Bendet, A.; Tanaka, C.; Yasuda, Y.; Nakahara, H. J. Am. Coll. Surg. 2006, 203, 54-63. 16. Li, J.-J.; Xu, G.-L.; Gu, M.-F.; Luo, G.-Y.; Rong, Z.; Wu, P.-H.; Xia, J.-C. World J. Gastroenterol. 2007, 13, 2747-2751. 17. Day, E. S.; Morton, J. G.; West, J. L. J. Biomech. Eng. 2009, 131, 074001. 18. Sun, Y.; Zheng, Y.; Li, P.; Wang, D.; Niu, C.; Gong, Y.; Huang, R.; Wang, Z.; Ran, H. BMC cancer 2014, 14, 800. 19. Sun, Y.; Zheng, Y.; Ran, H.; Zhou, Y.; Shen, H.; Chen, Y.; Chen, H.; Krupka, T. M.; Li, A.; Li, P.; Wang, Z. Biomaterials 2012, 33, 5854-64. 20. Dibaji, S. A. R., Al-Rjoub, M. F., Myers, M. R., Banerjee, R. K. J. Nanotechnol. Eng. Med 2014, 4, 040902. 21. Etheridge, M. L.; Hurley, K. R.; Zhang, J.; Jeon, S.; Ring, H. L.; Hogan, C.; Haynes, C. L.; Garwood, M.; Bischof, J. C. Technology 2014, 2, 214-228. 22. Morris, H.; Rivens, I.; Shaw, A.; Haar, G. T. Phys. Med. Biol. 2008, 53, 4759-76. 23. Dasgupta, S.; Banerjee, R. K.; Hariharan, P.; Myers, M. R. Ultrasonics 2011, 51, 171-80. 24. Kohler, M. O.; Mougenot, C.; Quesson, B.; Enholm, J.; Le Bail, B.; Laurent, C.; Moonen, C. T.; Ehnholm, G. J. Med. Phys. 2009, 36, 3521-35. 25. Quesson, B.; Laurent, C.; Maclair, G.; de Senneville, B. D.; Mougenot, C.; Ries, M.; Carteret, T.; Rullier, A.; Moonen, C. T. NMR Biomed. 2011, 24, 145-53. 26. Hynynen, K.; Pomeroy, O.; Smith, D. N.; Huber, P. E.; McDannold, N. J.; Kettenbach, J.; Baum, J.; Singer, S.; Jolesz, F. A. Radiology 2001, 219, 176-85. 27. McDannold, N. J.; King, R. L.; Jolesz, F. A.; Hynynen, K. H. Radiology 2000, 216, 517-23. 28. Canney, M. S.; Bailey, M. R.; Crum, L. A.; Khokhlova, V. A.; Sapozhnikov, O. A. J. Acoust. Soc. Am. 2008, 124, 2406-20.
22 ACS Paragon Plus Environment
Page 22 of 24
Page 23 of 24
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417
Nano Letters
29. Salomir, R.; Vimeux, F. C.; de Zwart, J. A.; Grenier, N.; Moonen, C. T. W. Magn. Reson. Med. 2000, 43, 342-347. 30. Chen, L.; Bouley, D.; Yuh, E.; D'Arceuil, H.; Butts, K. J. Magn. Reson. Imaging 1999, 10, 146-53. 31. Bohris, C.; Jenne, J. W.; Rastert, R.; Simiantonakis, I.; Brix, G.; Spoo, J.; Hlavac, M.; Nemeth, R.; Huber, P. E.; Debus, J. Magn. Reson. Imaging 2001, 19, 167-75. 32. Bohris, C.; Schreiber, W. G.; Jenne, J.; Simiantonakis, I.; Rastert, R.; Zabel, H. J.; Huber, P.; Bader, R.; Brix, G. Magn. Reson. Imaging 1999, 17, 603-10. 33. Larrat, B.; Pernot, M.; Aubry, J. F.; Dervishi, E.; Sinkus, R.; Seilhean, D.; Marie, Y.; Boch, A. L.; Fink, M.; Tanter, M. Phys. Med. Biol. 2010, 55, 365-88. 34. King, R. L.; Liu, Y.; Maruvada, S.; Herman, B. A.; Wear, K. A.; Harris, G. R. IEEE Trans. Ultrason., Ferroelectr., Freq. Control 2011, 58, 1397-405. 35. Partanen, A.; Tillander, M.; Yarmolenko, P. S.; Wood, B. J.; Dreher, M. R.; Köhler, M. O. Med. Phys. 2013, 40, 013301-n/a. 36. Hariharan, P.; Dibaji, S. A.; Banerjee, R. K.; Nagaraja, S.; Myers, M. R. IEEE Trans. Ultrason., Ferroelectr., Freq. Control 2014, 61, 2019-31. 37. Dillon, C. R.; Vyas, U.; Payne, A.; Christensen, D. A.; Roemer, R. B. Phys. Med. Biol. 2012, 57, 4527-44. 38. Sapareto, S. A.; Dewey, W. C. Int. J. Radiat. Oncol. Biol. Phys. 1984, 10, 787-800. 39. Wu, J.; Du, G. Ultrasound. Med. Biol. 1993, 19, 167-76. 40. Allegra, J. R.; Hawley, S. A. J. Acoust. Soc. Am. 1972, 51, 1545-1564. 41. Simons, S. Proc. Phys. Soc. 1964, 83, 749. 42. Brawer, S. Phys. Rev. B 1973, 7, 1712-1717. 43. Pinkerton, J. M. M. Proc. Phys. Soc., B. 1949, 62, 129. 44. Maeda, H.; Matsumura, Y. Adv. Drug Deliv. Rev. 2011, 63, 129-30. 45. Iyer, A. K.; Khaled, G.; Fang, J.; Maeda, H. Drug Discov. Today 2006, 11, 812-8. 46. Byrne, J. D.; Betancourt, T.; Brannon-Peppas, L. Adv. Drug Deliv. Rev. 2008, 60, 1615-26.
418 419
23 ACS Paragon Plus Environment
Nano Letters
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
365x129mm (96 x 96 DPI)
ACS Paragon Plus Environment
Page 24 of 24