Comparison of thermal damage calculated using magnetic resonance thermometry, with magnetic resonance imaging post-treatment and histology, after interstitial microwave thermal therapy of rabbit brain

Gd2O3-mesoporous silica/gold nanoshells: A potential dual T1/T2 contrast agent for MRI-guided localized near-IR photothermal therapy

A promising clinical trial utilizing gold-silica core-shell nanostructures coated with polyethylene glycol (PEG) has been reported for near-infrared (NIR) photothermal therapy (PTT) of prostate cancer. The next critical step for PTT is the visualization of therapeutically relevant nanoshell (NS) concentrations at the tumor site. Here we report the synthesis of PEGylated Gd₂O₃-mesoporous silica/gold core/shell NSs (Gd₂O₃-MS NSs) with NIR photothermal properties that also supply sufficient MRI contrast to be visualized at therapeutic doses (≥10⁸ NSs per milliliter). The nanoparticles have r₁ relaxivities more than three times larger than those of conventional T₁ contrast agents, requiring less concentration of Gd³⁺ to observe an equivalent signal enhancement in T₁-weighted MR images. Furthermore, Gd₂O₃-MS NS nanoparticles have r₂ relaxivities comparable to those of existing T₂ contrast agents, observed in agarose phantoms. This highly unusual combination of simultaneous T₁ and T₂ contrast allows for MRI enhancement through different approaches. As a rudimentary example, we demonstrate T₁/T₂ ratio MR images with sixfold contrast signal enhancement relative to its T₁ MRI and induced temperature increases of 20 to 55 °C under clinical illumination conditions. These nanoparticles facilitate MRI-guided PTT while providing real-time temperature feedback through thermal MRI mapping.

The utility of diffusion tractography for speech preservation in laser ablation of the dominant insula: illustrative case

BACKGROUND

Open surgical treatment of insular epilepsy holds particular risk of injury to middle cerebral artery branches, the operculum (through retraction), and adjacent language-related white matter tracts in the language-dominant hemisphere. Magnetic resonance imaging (MRI)-guided laser interstitial thermal therapy (LITT) is a surgical alternative that allows precise lesioning with potentially less operative risk. The authors presented the case of a 13-year-old girl with intractable, MRI-negative, left (dominant hemisphere) insular epilepsy that was treated with LITT. Diffusion tensor imaging (DTI) tractography was used to aid full posterior insular lesioning in the region of stereo electroencephalography–determined seizure onset while avoiding thermal injury to the language-related superior longitudinal fasciculus (SLF)/arcuate fasciculus (AF) and inferior fronto-occipital fasciculus (IFOF).

OBSERVATIONS

DTI tractography was used successfully in planning insular LITT and facilitated a robust insular ablation with sharp margins at the interfaces with the SLF/AF and IFOF. These tracts were spared, and no neurological deficits were induced through LITT.

LESSONS

Although it is technically demanding and has important limitations that must be understood, clinically available DTI tractography adds precision and confidence to insular laser ablation when used to protect important language-related white matter tracts.

Novel Paradigm of Subdural Cortical Stimulation Does Not Cause Thermal Damage in Brain Tissue: A Simulation-Based Study

The thermal effect of a novel effective electrical stimulation mapping (ESM) technique using an Ojemann's stimulation electrode in open craniotomy areas causes a nondestructive local increase in temperature. Another type of stimulating electrode is a subdural strip, routinely used in intraoperative electrocorticography (ECoG), which applies ESM in a covered subdural area over the motor cortex. ECoG electrode geometry produces a different electrical field, causing a different Joule heat distribution in tissue, one that is impossible to measure in subdural space. Therefore, the previous safety control study of the novel ESM technique needed to be extended to include an assessment of the thermal effect of ECoG strip electrodes. We adapted a previously well-validated numerical model and performed coupled complex electro-thermal transient simulations for short-time (28.4 ms) high-frequency (500 Hz) and hyperintense (peak 100 mA) ESM paradigm. The risk of heat-induced cellular damage was assessed by applying the Arrhenius equation integral on the computed time-dependent spatial distribution of temperature in the brain tissue during ESM stimulation and during the cooldown period. The results showed increases in temperature in the proximity around ECoG electrode discs in a safe range without destructive effects. As opposed to open craniotomy, subdural space is not cooled by the air; hence a higher - but still safe - induced temperature was observed. The presented simulation agrees with the previously published histopathological examination of the stimulated brain tissue, and confirms the safety of the novel ESM technique when applied using ECoG strip electrodes.

Picosecond Laser-Induced Photothermal Skin Damage Evaluation by Computational Clinical Trial

Background and Objectives

Computational clinical trial (CCT) in the field of laser medicine promotes clinical application of novel laser devices, because this trial carried out based on numerical modeling of laser-tissue interactions and simulation of a series of treatment process. To confirm the feasibility of the computational clinical trial of skin treatment with a novel picosecond laser, this paper presents an evaluation method of the safety.

Study Design/Materials and Methods

In this method, the light propagation and thermal diffusion process after ultrashort light pulse irradiation to a numerical skin model is calculated and the safety based on the photothermal damage is evaluated by computational modeling and simulation. As an example, the safety of a novel picosecond laser device was examined by comparing with several laser devices approved for clinical use.

Results

The ratio of the maximum thermal damage induced by picosecond laser irradiation was 1.2 × 10^-2 % at the epidermis, while that caused by approved laser irradiation was 99 % at the capillary vessels. The numerical simulation demonstrated that less thermal damage was observed compared with the approved devices. The results show the safety simulated by photothermal damage calculation was consistent with the reported clinical trials.

Conclusions

This computational clinical trial shows the feasibility of applying computational clinical trials for the safety evaluation of novel medical laser devices. In contrast to preclinical and clinical tests, the proposed computational method offers regulatory science for appropriately and quickly predicting and evaluating the safety of a novel laser device.

Multicenter validation of automated trajectories for selective laser amygdalohippocampectomy

Objective

Laser interstitial thermal therapy (LITT) is a novel minimally invasive alternative to open mesial temporal resection in drug‐resistant mesial temporal lobe epilepsy (MTLE). The safety and efficacy of the procedure are dependent on the preplanned trajectory and the extent of the planned ablation achieved. Ablation of the mesial hippocampal head has been suggested to be an independent predictor of seizure freedom, whereas sparing of collateral structures is thought to result in improved neuropsychological outcomes. We aim to validate an automated trajectory planning platform against manually planned trajectories to objectively standardize the process.

Methods

Using the EpiNav platform, we compare automated trajectory planning parameters derived from expert opinion and machine learning to undertake a multicenter validation against manually planned and implemented trajectories in 95 patients with MTLE. We estimate ablation volumes of regions of interest and quantify the size of the avascular corridor through the use of a risk score as a marker of safety. We also undertake blinded external expert feasibility and preference ratings.

Results

Automated trajectory planning employs complex algorithms to maximize ablation of the mesial hippocampal head and amygdala, while sparing the parahippocampal gyrus. Automated trajectories resulted in significantly lower calculated risk scores and greater amygdala ablation percentage, whereas overall hippocampal ablation percentage did not differ significantly. In addition, estimated damage to collateral structures was reduced. Blinded external expert raters were significantly more likely to prefer automated to manually planned trajectories.

Significance

Retrospective studies of automated trajectory planning show much promise in improving safety parameters and ablation volumes during LITT for MTLE. Multicenter validation provides evidence that the algorithm is robust, and blinded external expert ratings indicate that the trajectories are clinically feasible. Prospective validation studies are now required to determine if automated trajectories translate into improved seizure freedom rates and reduced neuropsychological deficits.

Modeling of Brain Tissue Heating Caused by Direct Cortical Stimulation for Assessing the Risk of Thermal Damage

This paper aims to employ the numerical simulations to assess the risk of cellular damage during the application of a novel paradigm of electrical stimulation mapping (ESM) used in neurosurgery. The core principle of the paradigm is the use of short, high-intensity and high-frequency stimulation pulses. We developed a complex numerical model and performed coupled electro-thermal transient simulations. The model was optimized by incorporating ESM electrodes' resistance obtained during multiple intraoperative measurements and validated by comparing them with the results of temperature distribution measurement acquired by thermal imaging. The risk of heat-induced cellular damage was assessed by applying the Arrhenius equation integral on the computed time-dependent spatial distribution of temperature in the brain tissue. Our results suggest that the impact of the temperature increase during our novel ESM paradigm is thermally non-destructive. The presented simulation results match the previously published thermographic measurement and histopathological examination of the stimulated brain tissue and confirm the safety of the novel ESM.

Quantitative MR thermometry based on phase-drift correction PRF shift method at 0.35 T

Background

Noninvasive magnetic resonance thermometry (MRT) at low-field using proton resonance frequency shift (PRFS) is a promising technique for monitoring ablation temperature, since low-field MR scanners with open-configuration are more suitable for interventional procedures than closed systems. In this study, phase-drift correction PRFS with first-order polynomial fitting method was proposed to investigate the feasibility and accuracy of quantitative MR thermography during hyperthermia procedures in a 0.35 T open MR scanner.

Methods

Unheated phantom and ex vivo porcine liver experiments were performed to evaluate the optimal polynomial order for phase-drift correction PRFS. The temperature estimation approach was tested in brain temperature experiments of three healthy volunteers at room temperature, and in ex vivo porcine liver microwave ablation experiments. The output power of the microwave generator was set at 40 W for 330 s. In the unheated experiments, the temperature root mean square error (RMSE) in the inner region of interest was calculated to assess the best-fitting order for polynomial fit. For ablation experiments, relative temperature difference profile measured by the phase-drift correction PRFS was compared with the temperature changes recorded by fiber optic temperature probe around the microwave ablation antenna within the target thermal region.

Results

The phase-drift correction PRFS using first-order polynomial fitting could achieve the smallest temperature RMSE in unheated phantom, ex vivo porcine liver and in vivo human brain experiments. In the ex vivo porcine liver microwave ablation procedure, the temperature error between MRT and fiber optic probe of all but six temperature points were less than 2 °C. Overall, the RMSE of all temperature points was 1.49 °C.

Conclusions

Both in vivo and ex vivo experiments showed that MR thermometry based on the phase-drift correction PRFS with first-order polynomial fitting could be applied to monitor temperature changes during microwave ablation in a low-field open-configuration whole-body MR scanner.

A methodology for thermal dose model parameter development using perioperative MRI

Post-treatment imaging is the principal method for evaluating thermal lesions following image-guided thermal ablation procedures. While real-time temperature feedback using magnetic resonance temperature imaging (MRTI) is a complementary tool that can be used to optimise lesion size throughout the procedure, a thermal dose model is needed to convert temperature-time histories to estimates of thermal damage. However, existing models rely on empirical parameters derived from laboratory experiments that are not direct indicators of post-treatment radiologic appearance. In this work, we investigate a technique that uses perioperative MR data to find novel thermal dose model parameters that are tailored to the appearance of the thermal lesion on post-treatment contrast-enhanced imaging. Perioperative MR data were analysed for five patients receiving magnetic resonance-guided laser-induced thermal therapy (MRgLITT) for brain metastases. The characteristic enhancing ring was manually segmented on post-treatment T1-weighted imaging and registered into the MRTI geometry. Post-treatment appearance was modelled using a coupled Arrhenius-logistic model and non-linear optimisation techniques were used to find the maximum-likelihood kinetic parameters and dose thresholds that characterise the inner and outer boundary of the enhancing ring. The parameter values and thresholds were consistent with previous investigations, while the average difference between the predicted and segmented boundaries was on the order of one pixel (1 mm). The areas predicted using the optimised model parameters were also within 1 mm of those predicted by clinically utilised dose models. This technique makes clinically acquired data available for investigating new thermal dose model parameters driven by clinically relevant endpoints.

Evaluation of MR thermometry with proton resonance frequency method at 7T

Quantitative and non-invasive temperature mapping using magnetic resonance imaging (MRI) provides a unique way to measure temperature evolution inside biological tissues. The method is widely used in thermal ablation procedures with magnetic fields at or below 3T. In this paper, the sensitivity of the MRI thermometry at 7T was studied using a proton resonance frequency (PRF)-based technique. We first used an agarose gel phantom with MR-compatible thermometry to calibrate the temperature coefficient, and then this temperature coefficient was employed to measure the internal temperature in both ex vivo (beef muscle) and in vivo (rat) experiments using focused ultrasound heating. The temperature coefficient calibrated by the phantom was 0.0095 ppm/°C, and both the ex vivo and in vivo experiments exhibited clear temperature evolution. This quantitative study confirmed the sensitivity (<1 °C) of MR temperature mapping at 7T.

ANALYTICAL SOLUTIONS IN THE MODELING OF THE LOCAL RF ABLATION

Coupled mathematical models for the radiofrequency (RF) ablation performed in biomedical sciences have been developed based on the bioheat transfer theory. The heat exchange problem is important to be analytically studied in order to control the size of the necrosis zone caused by RF ablation. This lesion size in the tissue may be predicted by the knowledge of the internal tissue temperature. We propose an analytical solution for the Pennes heat transfer equation in bi- and tri-region domains, applicable to the RF ablation of cancerigeneous tissue — a clinical relevant problem. The model consists of two partial differential equations describing the spatio-temporal interactions between the electric and thermic effects. The aim is to find simple algebraic expressions of analytical solutions that may allow to generate quantitative results which in turn may be interpreted (including uncertainties). The dependence of the temperature as function of the electrothermal parameters in both diseased and surrounding healthy tissues is pointed out. Two cases, namely the tumor–tissue and tumor–tissue–skin systems, are graphically computed, and important findings include the fact that the presence of tissue with smaller value parameters protects somehow healthy cells. Moreover, the graphical representations are conducted to highlight the link of the profile of current density distribution in the physiological problem with the (neither oval nor circular) shape of the temperature isoclinic lines.

A statistical model describing combined irreversible electroporation and electroporation-induced blood-brain barrier disruption

Background

Electroporation-based therapies such as electrochemotherapy (ECT) and irreversible electroporation (IRE) are emerging as promising tools for treatment of tumors. When applied to the brain, electroporation can also induce transient blood-brain-barrier (BBB) disruption in volumes extending beyond IRE, thus enabling efficient drug penetration. The main objective of this study was to develop a statistical model predicting cell death and BBB disruption induced by electroporation. This model can be used for individual treatment planning.

Material and methods

Cell death and BBB disruption models were developed based on the Peleg-Fermi model in combination with numerical models of the electric field. The model calculates the electric field thresholds for cell kill and BBB disruption and describes the dependence on the number of treatment pulses. The model was validated using in vivo experimental data consisting of rats brains MRIs post electroporation treatments.

Results

Linear regression analysis confirmed that the model described the IRE and BBB disruption volumes as a function of treatment pulses number (r² = 0.79; p < 0.008, r² = 0.91; p < 0.001). The results presented a strong plateau effect as the pulse number increased. The ratio between complete cell death and no cell death thresholds was relatively narrow (between 0.88-0.91) even for small numbers of pulses and depended weakly on the number of pulses. For BBB disruption, the ratio increased with the number of pulses. BBB disruption radii were on average 67% ± 11% larger than IRE volumes.

Conclusions

The statistical model can be used to describe the dependence of treatment-effects on the number of pulses independent of the experimental setup.

Tissue Ablation Dynamics During Magnetic Resonance-Guided, Laser-Induced Thermal Therapy

BACKGROUND

Magnetic resonance-guided, laser-induced thermal therapy is a real-time magnetic resonance thermometry-guided, minimally invasive procedure used in the treatment of intracranial tumors, epilepsy, and pain. Little is known about its dynamics and the effects of various pathologies on overall ablation.

OBJECTIVE

To determine the relationship between thermal energy delivery and the time to maximal estimated thermal damage and whether differences exist between various intracranial pathologies.

METHODS

We used real-time ablation data from 28 patients across 5 unique intracranial pathologies. All ablations were performed using the Visualase Thermal Therapy System (Medtronic, Inc, Minneapolis, Minnesota), which uses a 980-nm diffusing tip diode laser. The thermal damage area was plotted against time for each ablation. We then estimated the duration of time required to reach 50% (t50) and 97% (t97) of maximal damage. Comparisons were then made between different intracranial pathologies.

RESULTS

The duration required to reach maximal thermal damage estimate (TDE) among all ablations was 159 ± 62 seconds, and the t50 and t97 were 43 ± 21 and 136 ± 57 seconds, respectively, where t97 was reached at an average of 23 seconds before the maximal TDE. The t97 was shorter in the recurrent metastasis/radiation necrosis and epilepsy groups compared with the previously untreated glioblastoma multiforme group.

CONCLUSION

The optimal duration can be estimated by the t97, which can be achieved in less than 3 minutes and differs across ablation targets. TDE expansion decelerates with prolonged ablation. Future studies are needed to examine the radiographic and clinical outcomes as well as the effects of ablation power, irrigation speed, and the effect of previous therapies on ablation dynamics.

Enhancing laser therapy using PEGylated gold nanoparticles combined with ultrasound and microbubbles

BACKGROUND

Gold nanorod (AuNR) laser therapy (LT) is a non-invasive method of increasing the temperature of a target tissue using near infrared light. In this study, the effects of ultrasound and microbubbles (USMB) with AuNR and LT were investigated on cell viability.

METHODS

MDA-MB-231 cells in suspension were treated with three different treatment combinations of AuNR, LT and USMB (Pneg=0.6 or 1.0 MPa): (1) AuNR with USMB followed by LT, (2) AuNR and LT followed by USMB, and (3) USMB followed by AuNR and LT. Cells were also exposed to USMB and LT without AuNR. The USMB conditions were: 500 kHz frequency, 16 cycles, 1kHz pulse repetition frequency for 1 min in the presence of Definity microbubbles (1.7% v/v). AuNR and LT conditions were: mPEG coated AuNR at 3×10(11) np/mL and 1.9 W/cm(2) for 3 min. Following the treatment, cell viability was assessed using propidium iodide (PI) fluorescent marker and flow cytometry (VPI), and colony assay (VCA). Cell viabilities were compared using a non-parametric Mann-Whitney U-test and synergism was assessed using the Bliss Independence Model.

RESULTS AND DISCUSSION

USMB improved cell death when combined with AuNR and LT. VPI of 17±2% (at 0.6 MPa) and 11±4% (at 1.0 MPa) were observed with combined treatment of AuNR and USMB followed by LT compared to VPI of 60±2% (at 0.6 MPa) and 42±3% (at 1.0 MPa) with USMB alone and VPI of 22±3% for AuNR and LT. The combined effect of AuNR and LT with USMB was additive regardless of treatment order. VCA results agreed with the additive effect caused by combining AuNR and LT with USMB for all treatment orders. In the absence of AuNR, samples exposed to LT prior to USMB at 0.6 MPa increased VPI by 13% (p<0.01) showing a protective effect.

CONCLUSION

Combining AuNR and LT with USMB resulted in an additive effect on cell viability compared to AuNR and LT, or USMB. In addition, cells exposed to low intensity NIR light appear to be protected against USMB exposure.

A numerical investigation of the electric and thermal cell kill distributions in electroporation-based therapies in tissue

Electroporation-based therapies are powerful biotechnological tools for enhancing the delivery of exogeneous agents or killing tissue with pulsed electric fields (PEFs). Electrochemotherapy (ECT) and gene therapy based on gene electrotransfer (EGT) both use reversible electroporation to deliver chemotherapeutics or plasmid DNA into cells, respectively. In both ECT and EGT, the goal is to permeabilize the cell membrane while maintaining high cell viability in order to facilitate drug or gene transport into the cell cytoplasm and induce a therapeutic response. Irreversible electroporation (IRE) results in cell kill due to exposure to PEFs without drugs and is under clinical evaluation for treating otherwise unresectable tumors. These PEF therapies rely mainly on the electric field distributions and do not require changes in tissue temperature for their effectiveness. However, in immediate vicinity of the electrodes the treatment may results in cell kill due to thermal damage because of the inhomogeneous electric field distribution and high current density during the electroporation-based therapies. Therefore, the main objective of this numerical study is to evaluate the influence of pulse number and electrical conductivity in the predicted cell kill zone due to irreversible electroporation and thermal damage. Specifically, we simulated a typical IRE protocol that employs ninety 100-µs PEFs. Our results confirm that it is possible to achieve predominant cell kill due to electroporation if the PEF parameters are chosen carefully. However, if either the pulse number and/or the tissue conductivity are too high, there is also potential to achieve cell kill due to thermal damage in the immediate vicinity of the electrodes. Therefore, it is critical for physicians to be mindful of placement of electrodes with respect to critical tissue structures and treatment parameters in order to maintain the non-thermal benefits of electroporation and prevent unnecessary damage to surrounding healthy tissue, critical vascular structures, and/or adjacent organs.

Accumulation of phase-shift nanoemulsions to enhance MR-guided ultrasound-mediated tumor ablation in vivo

Magnetic resonance-guided high intensity focused ultrasound (MRgHIFU) is being explored as a non-invasive technology to treat solid tumors. However, the clinical use of HIFU for tumor ablation applications is currently limited by the long treatment times required. Phase-shift nanoemulsions (PSNE), consisting of liquid perfluorocarbon droplets that can be vaporized into microbubbles, are being developed to accelerate HIFU-mediated heating. The purpose of this study was to examine accumulation of PSNE in intramuscular rabbit tumors in vivo. MR images were acquired before and after intravenous injection of gadolinium-containing PSNE. MR signal enhancement was observed in rabbit tumors up to six hours after injection, indicating that PSNE accumulated in the tumors. In addition, PSNE vaporization was detected in the tumor with B-mode ultrasound imaging, and MR thermometry measurements indicated that PSNE accelerated the rate of HIFU-mediated heating. These results suggest that PSNE could dramatically improve the efficiency and clinical feasibility of MRgHIFU.

Heat localization for targeted tumor treatment with nanoscale near-infrared radiation absorbers

Focusing heat delivery while minimizing collateral damage to normal tissues is essential for successful nanoparticle-mediated laser-induced thermal cancer therapy. We present thermal maps obtained via magnetic resonance imaging characterizing laser heating of a phantom tissue containing a multiwalled carbon nanotube inclusion. The data demonstrate that heating continuously over tens of seconds leads to poor localization (∼ 0.5 cm) of the elevated temperature region. By contrast, for the same energy input, heat localization can be reduced to the millimeter rather than centimeter range by increasing the laser power and shortening the pulse duration. The experimental data can be well understood within a simple diffusive heat conduction model. Analysis of the model indicates that to achieve 1 mm or better resolution, heating pulses of ∼2 s or less need to be used with appropriately higher heating power. Modeling these data using a diffusive heat conduction analysis predicts parameters for optimal targeted delivery of heat for ablative therapy.

High-frequency irreversible electroporation (H-FIRE) for non-thermal ablation without muscle contraction

BACKGROUND

Therapeutic irreversible electroporation (IRE) is an emerging technology for the non-thermal ablation of tumors. The technique involves delivering a series of unipolar electric pulses to permanently destabilize the plasma membrane of cancer cells through an increase in transmembrane potential, which leads to the development of a tissue lesion. Clinically, IRE requires the administration of paralytic agents to prevent muscle contractions during treatment that are associated with the delivery of electric pulses. This study shows that by applying high-frequency, bipolar bursts, muscle contractions can be eliminated during IRE without compromising the non-thermal mechanism of cell death.

METHODS

A combination of analytical, numerical, and experimental techniques were performed to investigate high-frequency irreversible electroporation (H-FIRE). A theoretical model for determining transmembrane potential in response to arbitrary electric fields was used to identify optimal burst frequencies and amplitudes for in vivo treatments. A finite element model for predicting thermal damage based on the electric field distribution was used to design non-thermal protocols for in vivo experiments. H-FIRE was applied to the brain of rats, and muscle contractions were quantified via accelerometers placed at the cervicothoracic junction. MRI and histological evaluation was performed post-operatively to assess ablation.

RESULTS

No visual or tactile evidence of muscle contraction was seen during H-FIRE at 250 kHz or 500 kHz, while all IRE protocols resulted in detectable muscle contractions at the cervicothoracic junction. H-FIRE produced ablative lesions in brain tissue that were characteristic in cellular morphology of non-thermal IRE treatments. Specifically, there was complete uniformity of tissue death within targeted areas, and a sharp transition zone was present between lesioned and normal brain.

CONCLUSIONS

H-FIRE is a feasible technique for non-thermal tissue ablation that eliminates muscle contractions seen in IRE treatments performed with unipolar electric pulses. Therefore, it has the potential to be performed clinically without the administration of paralytic agents.

A parametric study delineating irreversible electroporation from thermal damage based on a minimally invasive intracranial procedure

BACKGROUND

Irreversible electroporation (IRE) is a new minimally invasive technique to kill undesirable tissue in a non-thermal manner. In order to maximize the benefits from an IRE procedure, the pulse parameters and electrode configuration must be optimized to achieve complete coverage of the targeted tissue while preventing thermal damage due to excessive Joule heating.

METHODS

We developed numerical simulations of typical protocols based on a previously published computed tomographic (CT) guided in vivo procedure. These models were adapted to assess the effects of temperature, electroporation, pulse duration, and repetition rate on the volumes of tissue undergoing IRE alone or in superposition with thermal damage.

RESULTS

Nine different combinations of voltage and pulse frequency were investigated, five of which resulted in IRE alone while four produced IRE in superposition with thermal damage.

CONCLUSIONS

The parametric study evaluated the influence of pulse frequency and applied voltage on treatment volumes, and refined a proposed method to delineate IRE from thermal damage. We confirm that determining an IRE treatment protocol requires incorporating all the physical effects of electroporation, and that these effects may have significant implications in treatment planning and outcome assessment. The goal of the manuscript is to provide the reader with the numerical methods to assess multiple-pulse electroporation treatment protocols in order to isolate IRE from thermal damage and capitalize on the benefits of a non-thermal mode of tissue ablation.

Quantitative comparison of thermal dose models in normal canine brain

PURPOSE

Minimally invasive thermal ablative therapies as alternatives to conventional surgical management of solid tumors and other pathologies is increasing owing to the potential benefits of performing these procedures in an outpatient setting with reduced complications and comorbidity. Magnetic resonance temperature imaging (MRTI) measurement allows existing thermal dose models to use the spatiotemporal temperature history to estimate the thermal damage to tissue. However, the various thermal dose models presented in the literature employ different parameters and thresholds, affecting the reliability of thermal dosimetry. In this study, the authors quantitatively compared three thermal dose models (Arrhenius rate process, CEM43, and threshold temperature) using the dice similarity coefficient (DSC).

METHODS

The DSC was used to compare the spatial overlap between the region of thermal damage as predicted by the models for in vivo normal canine brain during thermal therapy to the region of thermal damage as revealed by contrast-enhanced T1-weighted images acquired immediately after therapy (< 20 min). The outer edge of the hyperintense rim of the ablation region was used as the surrogate marker for the limits of thermal coagulation. The DSC was also used to investigate the impact of varying the thresholds on each models' ability to predict the zone of thermal necrosis.

RESULTS

At previously reported thresholds, the authors found that all three models showed good agreement (defined as DSC > 0.7) with post-treatment imaging. All three models examined across the range of commonly applied thresholds consistently showed highly accurate spatial overlap, low variability, and little dependence on temperature uncertainty. DSC values corresponding to cited thresholds were not significantly different from peak DSC values.

CONCLUSIONS

Thus, the authors conclude that the all three thermal dose models can be used as a reliable surrogate for postcontrast tissue damage verification imaging in rapid ablation procedures and can also be used to enhance the capability of MRTI to control thermal therapy in real time.

An Electrode Array for Limiting Blood Loss During Liver Resection: Optimization via Mathematical Modeling

Liver resection is the current standard treatment for patients with both primary and metastatic liver cancer. The principal causes of morbidity and mortality after liver resection are related to blood loss (typically between 0.5 and 1 L), especially in cases where transfusion is required. Blood transfusions have been correlated with decreased long-term survival, increased risk of perioperative mortality and complications. The goal of this study was to evaluate different designs of a radiofrequency (RF) electrode array for use during liver resection. The purpose of this electrode array is to coagulate a slice of tissue including large vessels before resecting along that plane, thereby significantly reducing blood loss. Finite Element Method models were created to evaluate monopolar and bipolar power application, needle and blade shaped electrodes, as well as different electrode distances. Electric current density, temperature distribution, and coagulation zone sizes were measured. The best performance was achieved with a design of blade shaped electrodes (5 x 0.1 mm cross section) spaced 1.5 cm apart. The electrodes have power applied in bipolar mode to two adjacent electrodes, then switched sequentially in short intervals between electrode pairs to rapidly heat the tissue slice. This device produces a ~1.5 cm wide coagulation zone, with temperatures over 97 masculineC throughout the tissue slice within 3 min, and may facilitate coagulation of large vessels.

Young's modulus reconstruction for radio-frequency ablation electrode-induced displacement fields: a feasibility study

Radio-frequency (RF) ablation is a minimally invasive treatment for tumors in various abdominal organs. It is effective if good tumor localization and intraprocedural monitoring can be done. In this paper, we investigate the feasibility of using an ultrasound-based Young's modulus reconstruction algorithm to image an ablated region whose stiffness is elevated due to tissue coagulation. To obtain controllable tissue deformations for abdominal organs during and/or intermediately after the RF ablation, the proposed modulus imaging method is specifically designed for using tissue deformation fields induced by the RF electrode. We have developed a new scheme under which the reconstruction problem is simplified to a 2-D problem. Based on this scheme, an iterative Young's modulus reconstruction technique with edge-preserving regularization was developed to estimate the Young's modulus distribution. The method was tested in experiments using a tissue-mimicking phantom and on ex vivo bovine liver tissues. Our preliminary results suggest that high contrast modulus images can be successfully reconstructed. In both experiments, the geometries of the reconstructed modulus images of thermal ablation zones match well with the phantom design and the gross pathology image, respectively.

Temperature mapping of thermal ablation using MRI

MRI is a unique tool for minimally invasive thermal ablation in that it can provide both targeting, monitoring and control during the procedure. Monitoring is achieved by using MRI temperature mapping. In this review the relevant physics is explained as a background to the state-of-the-art methods for computing temperature maps as well as the more cutting edge methods. The review covers both methods to monitor heating and cooling of tissue and explains temperature mapping using Proton Resonance Frequency shift, T1 mapping, diffusion mapping, R2* mapping and thermal models.

Quantitative MRI-based temperature mapping based on the proton resonant frequency shift: Review of validation studies

MRI-based temperature imaging that exploits the temperature-sensitive water proton resonant frequency shift is currently the only available method for reliable quantification of temperature changes in vivo. Extensive pre-clinical work has been performed to validate this method for guiding thermal therapies. That work has shown the method to be useful for all stages of the thermal therapy, from resolving heating below the threshold for damage to ensuring that the thermal exposure is sufficient within the target volume and protecting surrounding critical structures and to accurately predicting the extent of the ablated volume. In this paper, these validation studies will be reviewed. In addition, clinical studies that have shown this method feasible in human treatments will be overviewed.

Optimal power deposition patterns for ideal high temperature therapy/hyperthermia treatments

If it were possible to achieve, an ideal high temperature therapy or hyperthermia treatment would involve a single heating session and yield a desired thermal dose distribution in the tumour that would be attained in the shortest possible treatment time without heating critical normal tissues excessively. Simultaneously achieving all of these goals is impossible in practice, thus requiring trade-offs that allow clinicians to approach more closely some of these ideal goals at the expense of others. To study the basic nature of a subset of these trade-offs, the present simulation study looked at a simple, ideal case in which the tumour is heated by a single, optimized (with respect to space) power pulse, with no power deposition in the normal tissue. Results were obtained for two different clinical strategies (i.e. trade-off approaches), including: (1) an 'aggressive' approach, wherein the desired, uniform thermal dose is completely delivered to the tumour during the power-on period. This approach gives the clinician the satisfaction of knowing that the tumour was treated completely while power was being delivered, and yields the shortest attainable tumour dose delivery time. However, that benefit is attained at the cost of both 'overdosing' the tumour during the subsequent cool down period and, paradoxically, requiring a longer, overall treatment time. Here, the treatment time is considered as that time interval from the initiation of the heating pulse to the time at which the entire tumour has decayed to a specified 'safe' temperature--below 43 degrees C for our calculations. And, (2) a 'conservative' approach is considered, wherein the desired uniform dose is attained at the post-heating time at which the complete tumour cools back down to 'basal' conditions, taken as 4 h in this study. This conservative approach requires less applied power and energy and avoids the 'overdosing' problem, but at the cost of having a tumour dose delivery time that can be significantly longer than the heating pulse duration. This approach can require that clinicians wait a significant time after the power has been turned off before being able to confirm that the desired tumour thermal dose was reached. The present findings show that: (1) for both clinical strategies, an optimal power deposition shape (with respect to position in the tumour) can always be found that provides the desired uniform thermal dose in the tumour, regardless of the heating pulse duration chosen or the tumour perfusion pattern; and (2) shorter heating pulses are preferable to longer ones in that they require less total energy, take less total time to treat the patients, and have optimal power deposition patterns less influenced by perfusion. On the other hand, shorter pulses always require higher temperatures, and for the 'aggressive' clinical approach, they give significantly larger excess thermal doses in the tumour. The aggressive approach always requires longer treatment times than comparable conservative treatments. The optimal power patterns for both strategies involve a high-power density at the tumour boundary, which frequently creates a 'thermal wave' that contributes significantly to the final thermal dose distribution attained.

Conformal thermal therapy using planar ultrasound transducers and adaptive closed-loop MR temperature control: demonstration in gel phantoms andex vivotissues

MRI-guided transurethral ultrasound therapy offers a minimally invasive approach for the treatment of localized prostate cancer. Integrating a multi-element planar transducer with active MR temperature feedback can enable three-dimensional conformal thermal therapy of a target region within the prostate gland while sparing surrounding normal tissues. Continuous measurement of the temperature distribution in tissue enables dynamic compensation for unknown changes in blood flow and tissue properties during treatment. The main goal of this study was to evaluate the feasibility of using active temperature feedback on a clinical 1.5 T MR imager for conformal thermal therapy. MR thermometry was performed during heating in both gel phantoms and excised tissue with a transurethral heating applicator, and the rotation rate and power were varied based on the thermal measurements. The capability to produce a region of thermal damage that matched a target boundary was evaluated. The influence of a cooling gradient (to simulate cooling of the rectum or urethra) on the desired pattern of thermal damage was also investigated in gel phantoms. Results showed high correlation between the desired target boundary and the 55 degrees C isotherm generated during heating with an average distance error of 0.9 mm +/- 0.4 mm (n = 6) in turkey breasts, 1.4 mm +/- 0.6 mm (n = 4) in gel phantoms without rectal cooling and 1.4 mm +/- 0.6 mm (n = 3) in gel phantoms with rectal cooling. The results were obtained using a temporal update rate of 5 s, a spatial resolution of 3 x 3 x 10 mm for the control point, and a temperature uncertainty of approximately 1 degrees C. The performance of the control algorithm under these conditions was comparable to that of simulations conducted previously by our group. Overall, the feasibility of generating targeted regions of thermal damage with a transurethral heating applicator and active MR temperature feedback has been demonstrated experimentally. This method of treatment appears capable of accounting for unpredictable and varying tissue properties during the treatment.

MRI-guided Thermal Ablation Therapy: Model and Parameter Estimates to Predict Cell Death from MR Thermometry Images

Solid tumors and other pathologies can be treated using laser thermal ablation under interventional magnetic resonance imaging (iMRI) guidance. A model was developed to predict cell death from magnetic resonance (MR) thermometry measurements based on the temperature-time history, and validated using in vivo rabbit brain data. To align post-ablation T2-weighted spin-echo MR lesion images to gradient-echo MR images, from which temperature is derived, a registration method was used that aligned fiducials placed near the thermal lesion. The outer boundary of the hyperintense rim in the post-ablation MR lesion image was used as the boundary for cell death, as verified from histology. Model parameters were simultaneously estimated using an iterative optimization algorithm applied to every interesting voxel in 328 images from multiple experiments having various temperature histories. For a necrotic region of 766 voxels across all lesions, the model provided a voxel specificity and sensitivity of 98.1 and 78.5%, respectively. Mislabeled voxels were typically within one voxel from the segmented necrotic boundary with median distances of 0.77 and 0.22 mm for false positives (FP) and false negatives (FN), respectively. As compared to the critical temperature cell death model and the generalized Arrhenius model, our model typically predicted fewer FP and FN. This is good evidence that iMRI temperature maps can be used with our model to predict therapeutic regions in real-time during treatment.

Evaluation of endorectal and urethral cooling devices during MR-guided ultrasound thermal ablation in canine prostate

High-temperature thermal therapy for the treatment of prostate cancer is currently being applied as a minimally-invasive alternative over traditional forms of treatment. Catheter-based interstitial and transurethral ultrasound applicators are being developed for controlled and selective thermal ablation of prostaric tissues with concurrent MR thermal imaging. As part of this treatment strategy we have devised a transurethral cooling catheter and a cooling jacket to be placed over the endorectal MR imaging coil to protect the urethral mucosa and rectal wall from thermal damage during treatment. The cooling efficiencies and protective abilities of these devices were evaluated in vivo within three canine prostate glands. Invasive and MR derived temperature measurements within the prostate and rectal wall indicate that the protective influence of the endorectal cooling extends 5-10 mm from the rectal wall into the dorsal prostate. The urethral cooling extends approximately 5 mm from the cooling balloon. The protective capabilities were further verified with subsequent histological analysis with TTC stained tissue sections and contrast enhanced T1-weighted MR images post treatment. Both of these cooling devices are compatible with the MR thermometry and can be used to protect the urethral mucosa and rectal wall during prostate thermal ablation with interstitial and transurethral ultrasound devices.

MRI-guided laser thermal ablation: model to predict cell death from MR thermometry images for real-time therapy monitoring

Solid tumors and other pathologies can be treated using laser thermal ablation under interventional magnetic resonance imaging (iMRI) guidance. We developed a model to predict cell death from MR thermometry measurements and applied it to in vivo rabbit brain data. To align post-ablation T2-weighted spin-echo MR lesion images to gradient echo MR images, from which temperature is derived, we used a registration method that aligned fiducials placed near the thermal lesion. We used the outer boundary of the hyperintense rim in the post-ablation MR lesion image as the boundary for cell death, as verified from histology. Model parameters were simultaneously estimated using an iterative optimization algorithm applied to every interesting pixel in 328 images from multiple experiments having various temperature histories. For a necrotic region of 766 voxels across all lesions, the model gave a voxel specificity and sensitivity of 98.1% and 78.4%, respectively. Median distance between the segmented necrotic boundary and the mislabeled voxels was within one MR voxel. Furthermore, our model predicted fewer errors as compared to the critical temperature cell death model. This is good evidence that iMRI temperature maps can be used with our model to predict therapeutic regions in real-time.

An investigation of the use of transmission ultrasound to measure acoustic attenuation changes in thermal therapy

The potential of using a commercial ultrasound transmission imaging system to quantitatively monitor tissue attenuation changes after thermal therapy was investigated. The ultrasound transmission imaging system used, the AcoustoCam (Imperium Inc., MD) allows ultrasonic images to be captured using principles similar to that of a CCD-type camera that collects light. Ultrasound energy is focused onto a piezoelectric array by an acoustic lens system, creating a gray scale acoustic image. In this work, the pixel values from the acoustic images were assigned acoustic attenuation values by imaging polyacrylamide phantoms of varying known attenuation. After the calibration procedure, data from heated polyacrylamide/bovine serum albumin (BSA) based tissue-mimicking (TM) phantoms and porcine livers were acquired. Samples were heated in water at temperatures of 35, 45, 55, 65, and 75 degrees C for 1 h. Regions of interest were chosen in the images and acoustic attenuation values before and after heating were compared. An increase in ultrasound attenuation was found in phantoms containing BSA and in porcine liver. In the presence of BSA, attenuation in the TM phantom increased by a factor of 1.5, while without BSA no significant changes were observed. The attenuation of the porcine liver increased by up to a factor of 2.4, consistent with previously reported studies. The study demonstrates the feasibility of using a quantitative ultrasound transmission imaging system for monitoring thermal therapy.

Closed-form solution for the thermal dose delivered during single pulse thermal therapies

This study provides a closed form, analytical expression for the thermal dose delivered by a single heating pulse. The solution is derived using the effective cooling method and the non-linear Sapareto-Dewey equation to determine the thermal dose delivered by the time-temperature history of a treatment. The analytical solutions are used to determine the optimal treatment conditions, i.e. those that exactly deliver the desired thermal dose at a specified time. For purposes of illustration, this study focuses on a 'conservative' clinical approach in which the desired thermal dose is delivered at the end of the 'cool down' period. The analytical results show that, after a clinical strategy has been chosen (e.g. conservative, aggressive or intermediate), the user can only specify two free variables for such an optimal treatment. Results are presented which suggest that a practical approach would be to specify both (1) the desired thermal dose to be delivered to the target (the clinically relevant outcome) and (2) the peak temperature to be reached (a measurable, clinically useful, patient dependent response variable that can be employed in feedback control systems); and then determine the associated, optimal heating magnitude and duration that need to be used to reach that dose and temperature. The results also reveal that, with a given patient condition and power deposition distribution (together specifying an effective cooling time constant for the treatment) and a specified thermal dose, there is a maximum allowable peak temperature that, if exceeded, will result in 'over-dosing' the heated tissue. The results also show that avoiding such non-optimal 'over-dosing' will be difficult in most high temperature therapies since, when high temperatures are produced in tissues, the temperature decay must be very fast in order to avoid over-dosing during the cooling period. Such rapid cooling can only occur if short effective cooling time constants are present-either as a result of large tissue blood flows in the patient or due to large conduction effects induced by the use of highly localized power deposition sources.

MRI guided and monitored focused ultrasound thermal ablation methods: a review of progress

This paper reviews the current status in using magnetic resonance imaging (MRI) to guide and monitor thermal coagulation of tumours using focused ultrasound. The patient treatment procedure with a second generation phased array system will be described. Several clinical trials have found that patient treatments are feasible and that MRI thermometry allows noninvasive monitoring of clinical treatments. Overall, this emerging modality holds significant potential for non-invasive tumour treatment of both benign and malignant tumours.

Interstitial microwave thermal therapy and its application to the treatment of recurrent prostate cancer

Interstitial microwave thermal therapy may be an effective alternative to surgery for the treatment of some solid tumours. Arrays of helical antennae can produce complex heating patterns which when combined with active cooling of normal tissue structures can provide conformal heating for thermal coagulation of tumours. The development of a clinical protocol involving phantom and animal model studies, treatment planning, tissue property measurement and methods for on-line treatment monitoring is reviewed. The technology developed has been applied to the problem of recurrent prostate cancer following failed radiation treatment where available curative options are associated with high normal tissue morbidity. The purpose was to develop a treatment option for this group of patients with a very low side-effect profile that would not preclude further treatment if the disease progressed. Results of a Phase I/II trial demonstrate safety, promising efficacy and a low complication rate. As the technology for delivering this treatment matures, larger multi-institutional trials should be considered.

MRI-guided interstitial ultrasound thermal therapy of the prostate: A feasibility study in the canine model

The feasibility of MRI-guided interstitial ultrasound thermal therapy of the prostate was evaluated in an in vivo canine prostate model. MRI compatible, multielement interstitial ultrasound applicators were developed using 1.5 mm diameter cylindrical piezoceramic transducers (7 to 8 MHz) sectored to provide 180 degrees of angular directional heating. Two in vivo experiments were performed in canine prostate. The first using two interstitial ultrasound applicators, the second using three ultrasound applicators in conjunction with rectal and urethral cooling. In both experiments, the applicators were inserted transperineally into the prostate with the energy directed ventrally, away from the rectum. Electrical power levels of 5-17 W per element (approximately 1.6-5.4 W acoustic output power) were applied for heating periods of 18 and 48 min. Phase-sensitive gradient-echo MR imaging was used to monitor the thermal treatment in real-time on a 0.5 T interventional MRI system. Contrast-enhanced T1-weighted images and vital-stained serial tissue sections were obtained to assess thermal damage and correlate to real-time thermal contour plots and calculated thermal doses. Results from these studies indicated a large volume of ablated (nonstained) tissue within the prostate, extending 1.2 to 2.0 cm from the applicators to the periphery of the gland, with the dorsal margin of coagulation well-defined by the applicator placement and directionality. The shape of the lesions correlated well to the hypointense regions visible in the contrast-enhanced T1-weighted images, and were also in good agreement with the contours of the 52 degrees C threshold temperature and t43 > 240 min. This study demonstrates the feasibility of using directional interstitial ultrasound in conjunction with MRI thermal imaging to monitor and possibly control thermal coagulation within a targeted tissue volume while potentially protecting surrounding tissue, such as rectum, from thermal damage.

Changes in the dielectric properties of rat prostate ex vivo at 915 MHz during heating

Changes in the tissue dielectric properties at 915 MHz of rat prostate tissue due to heating at temperatures in the range 45-75 degrees C were measured. The changes were found to be reversible with temperature and independent of the time of heating. The dielectric properties at 23 +/- 1 degrees C were measured as epsilon' = 62.8 +/- 2.7 and sigma = 1.17 +/- 0.07 S/m, while the linear temperature coefficients for reversible changes were 1.10 +/- 0.11%/ degree C for conductivity and -0.31 +/- 0.05% /degree C for relative permittivity. These properties and their temperature coefficients can be utilized in microwave treatment planning programmes to provide insight into the effects of dielectric changes that arise during microwave thermal therapy of prostate cancer.

Noninvasive thermal ablation of hepatocellular carcinoma by using magnetic resonance imaging-guided focused ultrasound

A number of minimally invasive methods have been tested for the thermal ablation of liver tumors as an alternative to surgical resection. The use of focused ultrasound transducers to ablate deep tumors offers the first completely noninvasive alternative to these techniques. By increasing the flexibility of this technology with modern phased-array transducer design and by combining it with magnetic resonance imaging for targeting and online guidance, a powerful tool results with the potential to offer treatment to a larger population of patients, to reduce trauma to the patient, and to reduce the cost of treatment. In this article, we review previous work with focused ultrasound in the liver and recent experimental results with magnetic resonance imaging guidance.

Technical developments for cerebral thermal treatment: water-cooled diffusing laser fibre tips and temperature-sensitive MRI using intersecting image planes

The aim was to determine if water-cooled diffusing tips could produce larger and safer (better controlled) thermal lesions than non-cooled diffusing tips at 980 nm. Thermal lesions were induced in beef myocardium in vitro with and without water cooling using a 980 nm diode laser at various power levels. Seven intracerebral treatments were performed in six canines using water-cooled diffusing tips with four animals having intracerebral transmissible venereal tumours grown from inoculate. Magnetic resonance thermal imaging (MRTI)-based feedback software using a fast, radio frequency-spoiled gradient echo acquisition with two intersecting image planes was used for on-line monitoring and control of treatment and for the evaluation of in vivo laser lesion production. In cases where two-plane MRTI was employed, the maximum calculated temperature was compared in each plane. Using water-cooled tips and 400 micro m core diameter laser diffusing fibres in in vitro beef myocardium, power of up to 9.5 W was applied for 8 min without tip failure. Without cooling, tip failure occurred in under 4 min at 6 W, in under 2 min at 7 W and instantaneously at 8 W. Additionally, char accompanied lesions made with uncooled tips while cooled application resulted in only minimal char at only the highest thermal dose. Achieved lesion cross-sectional diameters in in vitro samples were up to 26.5 x 23.3 mm when water cooling was used. In canine brain and transmissible venereal tumours, up to 18.1 x 21.4 mm lesions were achieved. It is concluded that water cooling allows safe application of higher power to small core diameter diffusing tip fibres, which results in larger thermal lesions than can be achieved without cooling. Two-plane MRTI enhances on-line monitoring and feedback of thermal treatment.

Interstitial microwave thermal therapy for prostate cancer

Thermal therapy is used to kill tumors by heating them to temperatures >50 degrees C for an extended period of time. Cell death results from thermal coagulation. The energy sources available for this approach include radiofrequency electrodes, microwave antennas, laser fiberoptics, and ultrasound transducers. Each of these modalities has the potential to be delivered in a minimally invasive manner, and many theoretical and experimental investigations of these devices have been performed. This review describes current knowledge of interstitial microwave thermal therapy for prostate cancer. Examples are given from an ongoing trial in patients who have recurrent or persistent disease following radiation therapy. Future directions for pretreatment planning and real-time monitoring and control are discussed. These techniques have the potential to optimize treatments on a patient-specific basis and will be instrumental in planned future trials of this therapy as first line for prostate cancer.

Prediction of subtle thermal histopathological change using a novel analysis of Gd-DTPA kinetics

PURPOSE

To investigate Gd-DTPA kinetics as predictors of histopathological changes following focused ultrasound (FUS) thermal ablation for improved planning and assessment.

MATERIALS AND METHODS

Twenty-nine FUS lesions were created in the thigh muscle of eight rabbits under MR-guidance at 1.5 Tesla. Three rabbits were killed at four hours; and 11 lesions were analyzed with histopathology. Temperature-sensitive MRI using proton-resonant frequency-shift was used for time-dependent temperature measurements. Analysis of the uptake kinetics of Gd-DTPA was performed after Gd-DTPA injection, within 20 minutes after heating and again at two hours after heating. The resulting kinetic maps, permeability (K(trans)) and leakage space (v(e)), were correlated to peak temperatures, T(2)-weighted MR, and histopathology.

RESULTS

Images of K(trans) and v(e) reveal regions of histopathological change not visible on conventional post-therapy MR. At early times after heating, v(e) predicts the area of injury more accurately than T(2) (7 +/- 2% vs. 25 +/- 6% underestimation). A circular region of extensive structural/vascular disruption is indicated only on K(trans) maps. The sharp decrease in K(trans) at the boundary of this region occurs at 47.5 +/- 0.5 degrees C, and may be a better estimate of cell death than the conventional method of temperature threshold (55 degrees C for coagulation) used in therapy planning.

CONCLUSION

Our results suggest Gd-DTPA kinetics can predict different histopathological changes following FUS ablation and may be valuable for early prediction.

MR thermometry-based feedback control of laser interstitial thermal therapy at 980 nm

BACKGROUND AND OBJECTIVES

The goal of this study was to explore the feasibility of magnetic resonance thermal imaging (MRTI)-based feedback control of intracerebral laser interstitial thermal therapy (LITT), using a computer workstation and 980-nm diode laser interfaced to an MR scanner.

STUDY DESIGN/MATERIALS AND METHODS

A computer-controlled laser thermal therapy system was used to produce 12 ex vivo lesions in 3 canine and porcine brains and 16 in vivo lesions in 6 canines with diffusing tip fiberoptic applicators and energies from 54 to 900 J. MRTI predictions of thermal damage were correlated with histopathologic analysis.

RESULTS

Under feedback control, no carbonization, vaporization, or applicator damage was observed. MRTI-based prediction of thermal dose was not significantly different from histological evaluation of achieved thermal necrosis.

CONCLUSIONS

The computer-controlled thermal therapy system was effective at regulating heating, eliminating carbonization and vaporization, and protecting fiberoptic applicators. MRTI estimation of thermal dose accurately predicted achieved thermal necrosis.

Highly directional transurethral ultrasound applicators with rotational control for MRI-guided prostatic thermal therapy

Transurethral ultrasound applicators with highly directional energy deposition and rotational control were investigated for precise treatment of benign prostatic hyperplasia (BPH) and adenocarcinoma of the prostate (CaP). Two types of catheter-based applicators were fabricated, using either 90 degrees sectored tubular (3.5 mm OD x 10 mm) or planar transducers (3.5 mm x 10 mm). They were constructed to be MRI compatible, minimally invasive and allow for manual rotation of the transducer array within a 10 mm cooling balloon. In vivo evaluations of the applicators were performed in canine prostates (n = 3) using MRI guidance (0.5 T interventional magnet). MR temperature imaging (MRTI) utilizing the proton resonance frequency shift method was used to acquire multiple-slice temperature overlays in real time for monitoring and guiding the thermal treatments. Post-treatment T1-weighted contrast-enhanced imaging and triphenyl tetrazolium chloride stained tissue sections were used to define regions of tissue coagulation. Single sonications with the 90 degrees tubular applicator (9-15 W, 12 min, 8 MHz) produced coagulated zones covering an 80 degrees wedge of the prostate extending from 1-2 mm outside the urethra to the outer boundary of the gland (16 mm radial coagulation). Single sonications with the planar applicator (15-20 W, 10 min, approximately 8 MHz) generated thermal lesions of approximately 30 degrees extending to the prostate boundary. Multiple sequential sonications (sweeping) of a planar applicator (12 W with eight rotations of 30 degrees each) demonstrated controllable coagulation of a 270 degrees contiguous section of the prostate extending to the capsule boundary. The feasibility of using highly directional transurethral ultrasound applicators with rotational capabilities to selectively coagulate regions of the prostate while monitoring and controlling the treatments with MRTI was demonstrated in this study.

Usefulness of contrast kinetics for predicting and monitoring tissue changes in muscle following thermal therapy in long survival studies

PURPOSE

To investigate Gd-DTPA kinetics as indicators of subacute and subchronic histopathological changes following focused ultrasound (FUS) thermal therapy for improved evaluation.

MATERIALS AND METHODS

A total of 18 FUS lesions were created in the thigh muscle of five rabbits under magnetic resonance (MR) guidance at 1.5 Tesla. The rabbits were killed at different times: 40 hours, three days, and seven days. All lesions were analyzed histologically. An analysis of the uptake kinetics of Gd-DTPA, injected within two hours postheating and before sacrifice, was performed. The resulting kinetic maps, permeability (K(trans)) and leakage space (v(e)), were correlated to T(2)-weighted MR and histology.

RESULTS

Images of K(trans) and v(e) better differentiate subacute and subchronic changes not visible on conventional MR in the days following therapy and are consistent with the histopathology observed. In particular, the border between nonviable and viable tissue is well demarcated. The extent of damage is best indicated on v(e), whereas the borders of inflammation are shown on K(trans). The total lesion extent is relatively stable over the 7 days posttherapy and can be predicted by v(e) or T(2)-weighted MR at early times after heating.

CONCLUSION

Our results suggest that Gd-DTPA kinetics can complement conventional MR for improved evaluation of FUS thermal therapy by providing finer differentiation of necrotic states, inflammation, and repair processes.

MRI investigation of the threshold for thermally induced blood-brain barrier disruption and brain tissue damage in the rabbit brain

The ability of MRI-derived thermometry to predict thermally induced tissue changes in the brain was tested, and the thermal thresholds for blood-brain barrier (BBB) disruption and brain tissue damage were estimated. In addition, the ability of standard MRI to detect threshold-level effects was confirmed. These safety thresholds are being investigated to provide guidelines for clinical thermal ablation studies in the brain. MRI-monitored focused ultrasound heating was delivered to 63 locations in 26 rabbits. Tissue changes were detected in T(2)-weighted imaging and T(1)-weighted imaging (with and without contrast) and with light microscopy. The probability for tissue damage as a function of the accumulated thermal dose, the peak temperature achieved, the applied acoustic energy, and the peak acoustic power was estimated with probit regression. The discriminative abilities of these parameters were compared using the areas under the receiver operator characteristic (ROC) curves. In MRI, BBB disruption was observed in contrast-enhanced T(1)-weighted imaging shortly after the ultrasound exposures, sometimes accompanied by changes in T(2)-weighted imaging. Two days later, changes in T(2)-weighted imaging were observed, sometimes accompanied by changes in T(1)-weighted imaging. In histology, tissue damage was seen at every location where MRI changes were observed, ranging from small (diameter <1.0 mm) areas of tissue necrosis to severe vascular damage and associated hemorrhagic infarct. In one location, small (diameter: 0.8 mm) damage was not detected in MRI. The thermal dose and peak temperature thresholds were between 12.3-40.1 equivalent min at 43 degrees C and 48.0-50.8 degrees C, respectively, and values of 17.5 equivalent min at 43 degrees C and 48.4 degrees C were estimated to result in tissue damage with 50% probability. Thermal dose and peak temperature were significantly better predictors than the applied acoustic energy and peak acoustic power (P < 0.01). BBB disruption was always accompanied by tissue damage. The temperature information was better than the applied acoustic power or energy for predicting the damage than the ultrasound parameters. MRI was sensitive in detecting threshold-level damage.

Thermal therapy of canine cerebral tumors using a 980 nm diode laser with MR temperature-sensitive imaging feedback

BACKGROUND AND PURPOSE

The laser-induced thermal therapy (LITT) of cerebral tumors has conventionally been performed using Nd:YAG lasers and is associated with a risk of high focal temperatures potentially followed by cavitation that could result in boiling and/or explosive char. We have developed small diffusing laser fiber tips to better distribute the energy deposition and a computer controlled feedback system to monitor therapy and prevent excess temperature buildup. In this study, we evaluated the feasibility of using magnetic resonance temperature imaging (MRTI)-based feedback system for the thermal treatment of experimental intracerebral tumors using 980 nm laser irradiation delivered through these diffusing tips.

STUDY DESIGN/MATERIALS AND METHODS

Transmissible venereal tumors (TVTs) were grown via inoculation in the right cerebral hemisphere of seven canines. The laser fiber tips were inserted into a total of 10 independent TVT-suspected regions in the seven animals. Margins for the target area in each animal were prescribed on the basis of pretreatment MR images. MRTI-based feedback software was used to measure and regulate both temperature and the delivered thermal dose to achieve the desired thermal ablation and prevent excess heating. The effects of treatment were verified by results of histologic analyses.

RESULTS

Treatments resulted in contiguous areas of thermal necrosis in tumors and adjacent brain margin. The feedback software successfully cut off the laser power once the desired treatment volume was achieved, and prevented focal temperatures from exceeding predefined thresholds. Follow-up MRI studies showed 1.4- to 2.9-fold LITT-induced lesion expansion within 1-6 days after treatment.

CONCLUSIONS

Targeted thermal coagulation of small intracerebral tumors is feasible using MRTI-based feedback and diffused 980 nm diode laser light.

Models and measurements of light intensity changes during laser interstitial thermal therapy: implications for optical monitoring of the coagulation boundary location

We have developed a multi-region spherical Monte Carlo (MC) model to simulate the dynamic changes in light intensity measured during laser interstitial thermal therapy (LITT). Model predictions were validated experimentally in tissue-simulating albumen phantoms with well-characterized optical properties that vary dynamically with LITT in a way similar to tissue. For long treatments (2.5 W, approximately 1800 s), the transient light intensity changes demonstrated better qualitative agreement with a three-region MC model (with an inner layer of fully coagulated optical properties, a middle layer of partially coagulated properties and an outer region of native properties); for short treatments (4 W, approximately 240 s), better qualitative agreement was seen with a two-region MC model (with an inner layer of fully coagulated properties and outer region of native properties). These differences were attributed to differences in coagulation formation during low- and high-powered heating regimes, respectively. At the end of heating, a three-region coagulation zone was observed for both heating schemes. Quantitatively, final light intensity changes at the end of heating were compared with changes predicted by both two- and three-region MC for the same experimentally measured coagulation size and found to agree within approximately 30% for both models. The developed MC model helps lend insight into the nature of thermal coagulation events occurring for low and high power LITT irradiation schemes.

Variation correction algorithm: analysis of phase suppression and thermal profile fidelity for proton resonance frequency magnetic resonance thermometry at 0.2 T

PURPOSE

To develop and analyze the performance of the variation correction algorithm (VCA), a phase correction technique that mitigates the contribution of background phase variations by combining accurate alignment of echoes, K-space-based phase correction (as opposed to spatial polynomials), and extraction of alias-free phase difference images.

MATERIALS AND METHODS

A series of echo-shifted gradient-recalled echo (GRE) images was processed with K-space alignment and phase corrected with increasing sizes of M x M masks of central K-space coefficients. The extent of background phase variation suppression due to magnet field drift was assessed. Further, a simulated thermal profile was superimposed on the same data in a related experiment. Residual errors in reconstructed simulated thermal profiles were quantitatively characterized to estimate algorithm performance.

RESULTS

Using a 3 x 3 K-space mask, the VCA was able to 1) maintain the typical mean background error in a 35 x 35 pixel region of interest (ROI) at -0.1 degrees C; and 2) reconstruct, relative to the applied thermal profile, a phase-corrected profile that typically contains a 1.7 degrees C underestimation of peak temperature difference and a mean error along the 60 degrees C line of -0.8 degrees C.

CONCLUSION

The results suggest that thermal profiles can be accurately reconstructed at 0.2 T using the VCA, even in the presence of over 1 ppm spatially and temporally dependent field drift over a 1-hour time frame.

Multiplanar MR temperature-sensitive imaging of cerebral thermal treatment using interstitial ultrasound applicators in a canine model

PURPOSE

To study the feasibility of an interleaved gradient-echo, echo-planar imaging (iGE-EPI) sequence for multiplanar magnetic resonance temperature imaging (MRTI) to monitor intracerebral thermal treatment three-dimensionally using multielement ultrasound applicators.

MATERIALS AND METHODS

Transmissible venereal tumor (TVT) fragments were injected into the right cerebral hemisphere of five dogs. Guided by MRI, an interstitial ultrasound applicator was inserted into the tumor or normal brain tissue. The iGE-EPI sequence was used to estimate temperature changes by computing the complex phase-difference induced by temperature-dependent shifts in the proton resonance frequency of water. The thermal dose maps were updated every 6-8 seconds for five to seven image planes during treatment. The results of MRTI were compared with those of post-treatment MRI and histologic analysis.

RESULTS

The multiplanar MRTI monitored temperature and thermal dose distributions in tumor and normal brain tissue over the entire user-defined treatment volume. The ultrasound applicators produced contiguous areas of coagulative necrosis, resulting in 1.5-4.0 cm(3) volumes of tissue necrosis. MRTI-based assessments of thermal-dose distributions were consistent with the results of post-treatment MRI and histologic analysis.

CONCLUSION

Multiplanar MRTI is feasible for measuring necrosing thermal doses during intracerebral thermal delivery by interstitial ultrasound applicators.

MRI-guided thermal therapy of transplanted tumors in the canine prostate using a directional transurethral ultrasound applicator

PURPOSE

To evaluate MRI-based techniques for visual guidance, thermal monitoring, and assessment during transurethral ultrasound thermal therapy of implanted tumors in an in vivo canine prostate model.

MATERIALS AND METHODS

Transmissible venereal tumors (TVT) were grown in the right lobe of the prostate in four dogs. High-temperature thermal therapy was selectively applied to the tumor-bearing lobe using a transurethral ultrasound applicator with a 180 degrees directional heating pattern. Temperature-sensitive MRI (MRTI) using a fast interleaved gradient-echo echo-planar (iGE-EPI) imaging sequence was used for cumulative thermal dose calculations in multiple image planes during the treatment. The results from MRTI-based dose maps and post-treatment MRI were compared to those from histologic analysis.

RESULTS

MRTI monitoring in multiple planes across the prostate guided the use and control of a directive ultrasound applicator for the selective ablation of the sections of the prostate that contained implanted tumors. Findings in gadolinium enhanced MRI obtained immediately after thermal therapy slightly underestimated the size of tissue necrosis after treatment, as verified by histopathologic analysis.

CONCLUSION

The use of multiplanar MRTI with a transurethral ultrasound applicator shows significant potential for selective thermal ablation of prostate tumor and tissue.

The use of quantitative temperature images to predict the optimal power for focused ultrasound surgery: in vivo verification in rabbit muscle and brain

In this study, we investigated the use of MRI-derived thermal imaging for determining the exposure parameters for focused ultrasound (FUS) surgery. Since the temperature rise induced by a FUS beam scales linearly with power, the temperature maps acquired during subthreshold sonications can be used to determine the power necessary to produce thermal tissue damage with a desired size. Thermal images acquired during multiple sonications delivered at different locations in rabbit thigh muscle and brain tissue in vivo were analyzed to test this hypothesis. First, the linearity of the induced temperature rise with the acoustic power was tested. Next, the temperature maps acquired during preliminary low power sonications were scaled up until the estimated size of the tissue damage was equal to the tissue damage size of subsequent high power sonications. A threshold thermal dose was used to estimate the onset of thermal damage. The predicted power (based on amount of scaling required to reach the target size) was then compared to the true high power value. Overall, the temperature rise varied linearly with power (slope of deltaThigh/deltaTlow vs Power(high)/Power(low) = 0.97, 0.93 for pairs of sonications at each location in brain, muscle). The predicted power matched the true high power in the brain sonications (slope = 1.04). The predicted power underestimated the true high power in the muscle sonications (slope = 0.87). This under-prediction was due to a deviation from linearity in those cases where tissue damage was detected in subsequent MR images (slope of deltaThigh/deltaTlow vs Power(high)/Power(low) = 1.02, 0.84 for no tissue damage, tissue damage). The source of this deviation was not clear from these experiments. Even with this underestimation of the power, this method will be useful because it will allow an estimate of the proper power to use during FUS surgery without exact knowledge of the tissue parameters.