Interstitial ultrasound ablation of vertebral and paraspinal tumours: parametric and patient-specific simulations

Safety Review of Therapeutic Ultrasound for Spinal Cord Neuromodulation and Blood-Spinal Cord Barrier Opening

New focused ultrasound spinal cord applications have emerged, particularly those improving therapeutic agent delivery to the spinal cord via blood-spinal cord barrier opening and the neuromodulation of spinal cord tracts. One hurdle in the development of these applications is safety. It may be possible to use safety trends from seminal and subsequent works in focused ultrasound to guide the development of safety guidelines for spinal cord applications. We collated data from decades of pre-clinical studies and illustrate a clear relationship between damage, time-averaged spatial peak intensity and exposure duration. This relationship suggests a thermal mechanism underlies ultrasound-induced spinal cord damage. We developed minimum and mean thresholds for damage from these pre-clinical studies. When these thresholds were plotted against the parameters used in recent pre-clinical ultrasonic spinal cord neuromodulation studies, the majority of the neuromodulation studies were near or above the minimum threshold. This suggests that a thermal neuromodulatory effect may exist for ultrasonic spinal cord neuromodulation, and that the thermal dose must be carefully controlled to avoid damage to the spinal cord. By contrast, the intensity-exposure duration threshold had no predictive value when applied to blood-spinal cord barrier opening studies that employed injected contrast agents. Most blood-spinal cord barrier opening studies observed slight to severe damage, except for small animal studies that employed an active feedback control method to limit pressures based on measured bubble oscillation behavior. The development of new focused ultrasound spinal cord applications perhaps reflects the recent success in the development of focused ultrasound brain applications, and recent work has begun on the translation of these technologies from brain to spinal cord. However, a great deal of work remains to be done, particularly with respect to developing and accepting safety standards for these applications.

Image-based computer modeling assessment of microwave ablation for treatment of adrenal tumors

PURPOSE

To assess the feasibility of delivering microwave ablation for targeted treatment of aldosterone producing adenomas using image-based computational models.

METHODS

We curated an anonymized dataset of diagnostic ¹¹C-metomidate PET/CT images of 14 patients with aldosterone producing adenomas (APA). A semi-automated approach was developed to segment the APA, adrenal gland, and adjacent organs within 2 cm of the APA boundary. The segmented volumes were used to implement patient-specific 3D electromagnetic-bioheat transfer models of microwave ablation with a 2.45 GHz directional microwave ablation applicator. Ablation profiles were quantitatively assessed based on the extent of the APA target encompassed by an ablative thermal dose, while limiting thermal damage to the adjacent normal adrenal tissue and sensitive critical structures.

RESULTS

Across the 14 patients, adrenal tumor volumes ranged between 393 mm³ and 2,395 mm³. On average, 70% of the adrenal tumor volumes received an ablative thermal dose of 240CEM43, while limiting thermal damage to non-target structures, and thermally sparing 83.5-96.4% of normal adrenal gland. Average ablation duration was 293 s (range: 60-600 s). Simulations indicated coverage of the APA with an ablative dose was limited when the axis of the ablation applicator was not well aligned with the major axis of the targeted APA.

CONCLUSIONS

Image-based computational models demonstrate the potential for delivering microwave ablation to APA targets within the adrenal gland, while limiting thermal damage to surrounding non-target structures.

Modeling of Interstitial Ultrasound Ablation for Continuous Applicator Rotation With MR Validation

The primary objective of cancer intervention is the selective removal of malignant cells while conserving surrounding healthy tissues. However, the accessibility, size and shape of the cancer can make achieving appropriate margins a challenge. One minimally invasive treatment option for these clinical cases is interstitial needle based therapeutic ultrasound (NBTU). In this work, we develop a finite element model (FEM) capable of simulating continuous rotation of a directional NBTU applicator. The developed model was used to simulate the thermal deposition for different rotation trajectories. The actual thermal deposition patterns for the simulated trajectories were then evaluated using magnetic resonance thermal imaging (MRTI) in a porcine skin gelatin phantom. An MRI-compatible robot was used to control the rotation motion profile of the physical NBTU applicator to match the simulated trajectory. The model showed agreement when compared to experimental measurements with Pearson correlation coefficients greater than 0.839 when comparing temperature fields within an area of 12.6 mm radius from the ultrasound applicator. The average temperature error along a 6.3 mm radius profile from the applicator was 1.27 °C. The model was able to compute 1 s of thermal deposition by the applicator in 0.2 s on average with a 0.1 mm spatial resolution and 0.5 s time steps. The developed simulation demonstrates performance suitable for real-time control which may enable robotically-actuated closed-loop conformal tumor ablation.

In silico feasibility assessment of extracorporeal delivery of low-intensity pulsed ultrasound to intervertebral discs within the lumbar spine

Low intensity pulsed ultrasound (LIPUS) may have utility for non-invasive treatment of discogenic lower back pain through stimulating, remodeling and accelerating healing of injured or degenerated intervertebral disc (IVD) tissues. This study investigates the feasibility of delivering LIPUS to lumbar IVDs between L2 and S1 spine vertebra using a planar extracorporeal phased array (8 × 8 cm, 1024 elements, 500 kHz). Three 3D anatomical models with heterogenous tissues were generated from patient CT image sets and used in the simulation-based analysis. Time-reversal acoustic modeling techniques were applied to optimize posterior-lateral placement of the array with respect to the body to facilitate energy deposition in discrete target regions spanning the annulus fibrosus and central nucleus of each IVD. Forward acoustic and biothermal simulations were performed with time-reversal optimized array placements and driving amplitude/phase settings to predict LIPUS intensity distributions at target sites and to investigate off-target energy deposition and heating potential. Simulation results demonstrate focal intensity gain of 5-168 across all IVD targets and anatomical models, with greater average intensity gain (>50) and energy localization in posterior, posterolateral, and lateral target sites of IVDs. Localized LIPUS delivery was enhanced in thinner patient anatomies and in the high lumbar levels (L2-L3 and L3-L4). Multiple amplitude/phasing illumination patterns could be sequenced at a fixed array position for larger regional energy coverage in the IVD. Biothermal simulations demonstrated that LIPUS-appropriate exposures of 100 mW cm^-2 I_SPTA to the target disc region would result in <1 °C global peak temperature elevation for all cases. Hence, simulations suggest that spatially-precise extracorporeal delivery of therapeutically relevant LIPUS doses to discrete regions of lumbar IVDs is feasible and may be useful in clinical management of discogenic back pain.

An incoherent HIFU transducer for treatment of the medial branch nerve: Numerical study and in vivo validation

BACKGROUND

Chronic back pain due to facet related degenerative changes affects 4-6 million patients a year in the United States. Patients refractory to conservative therapy may warrant targeted injections of steroids into the joint or percutaneous medial branch nerve denervation with radiofrequency ablation. We numerically tested a novel noninvasive high intensity focused ultrasound transducer to optimize nerve ablation near a bone-soft tissue interface.

METHODS

A transducer with 4 elements operating in an incoherent mode was modeled numerically and tested pre-clinically under fluoroscopic guidance. After 6 lumbar medial branch nerve ablations were performed in 2 pigs, they were followed clinically for 1 week and then sacrificed for pathological evaluation.

RESULTS

Simulations show that the acoustic spot size in water at 6 dB was 14mm axial x 1.6mm lateral and 52mm axial x 1.6mm lateral for coherent and incoherent modes, respectively. We measured the size of N = 6 lesions induced in vivo in a pig model and compared them to the size of the simulated thermal dose. The best match between the simulated and measured lesion size was found with a maximum absorption coefficient in the cortical bone adjusted to 30 dB/cm/MHz. This absorption was used to simulate clinical scenarios in humans to generate lesions with no potential side effects at 1000 and 1500 J.

CONCLUSION

The elongated spot obtained with the incoherent mode facilitates the targeting during fluoroscopic-guided medial branch nerve ablation.

Thermal ablation of biological tissues in disease treatment: A review of computational models and future directions

Percutaneous thermal ablation has proven to be an effective modality for treating both benign and malignant tumours in various tissues. Among these modalities, radiofrequency ablation (RFA) is the most promising and widely adopted approach that has been extensively studied in the past decades. Microwave ablation (MWA) is a newly emerging modality that is gaining rapid momentum due to its capability of inducing rapid heating and attaining larger ablation volumes, and its lesser susceptibility to the heat sink effects as compared to RFA. Although the goal of both these therapies is to attain cell death in the target tissue by virtue of heating above 50°C, their underlying mechanism of action and principles greatly differs. Computational modelling is a powerful tool for studying the effect of electromagnetic interactions within the biological tissues and predicting the treatment outcomes during thermal ablative therapies. Such a priori estimation can assist the clinical practitioners during treatment planning with the goal of attaining successful tumour destruction and preservation of the surrounding healthy tissue and critical structures. This review provides current state-of-the-art developments and associated challenges in the computational modelling of thermal ablative techniques, viz., RFA and MWA, as well as touch upon several promising avenues in the modelling of laser ablation, nanoparticles assisted magnetic hyperthermia and non-invasive RFA. The application of RFA in pain relief has been extensively reviewed from modelling point of view. Additionally, future directions have also been provided to improve these models for their successful translation and integration into the hospital work flow.

Delivering Focused Ultrasound to Intervertebral Discs Using Time-Reversal

Chronic low back pain causes more disability worldwide than any other condition and is thought to arise in part through loss of biomechanical function of degenerate intervertebral discs (IVDs). Current treatments can involve replacing part or all of the degenerate IVDs by invasive surgery. Our vision is to develop a minimally invasive approach in which high intensity focused ultrasound (HIFU) is used to mechanically fractionate degenerate tissue in an IVD; a fine needle is then used to first remove the fractionated tissue and then inject a biomaterial able to restore normal physiologic function. The goal of this manuscript is to demonstrate the feasibility of trans-spinal HIFU delivery using simulations of 3-D ultrasound propagation in models derived from patient computed tomography (CT) scans. The CT data were segmented into bone, fat and other soft tissue for three patients. Ultrasound arrays were placed around the waist of each patient model, and time-reversal was used to determine the source signals necessary to create a focus in the center of the disc. The simulations showed that for 0.5 MHz ultrasound, a focus could be created in most of the lumbar IVDs, with the pressure focal gain ranging from 3.2-13.7. In conclusion, it is shown that with patient-specific planning, focusing ultrasound into an IVD is possible in the majority of patients despite the complex acoustic path introduced by the bony structures of the spine.

Theoretical investigation of transgastric and intraductal approaches for ultrasound-based thermal therapy of the pancreas

Background

The goal of this study was to theoretically investigate the feasibility of intraductal and transgastric approaches to ultrasound-based thermal therapy of pancreatic tumors, and to evaluate possible treatment strategies.

Methods

This study considered ultrasound applicators with 1.2 mm outer diameter tubular transducers, which are inserted into the tissue to be treated by an endoscopic approach, either via insertion through the gastric wall (transgastric) or within the pancreatic duct lumen (intraductal). 8 patient-specific, 3D, transient, biothermal and acoustic finite element models were generated to model hyperthermia (n = 2) and ablation (n = 6), using sectored (210°–270°, n = 4) and 360° (n = 4) transducers for treatment of 3.3–17.0 cm³ tumors in the head (n = 5), body (n = 2), and tail (n = 1) of the pancreas. A parametric study was performed to determine appropriate treatment parameters as a function of tissue attenuation, blood perfusion rates, and distance to sensitive anatomy.

Results

Parametric studies indicated that pancreatic tumors up to 2.5 or 2.7 cm diameter can be ablated within 10 min with the transgastric and intraductal approaches, respectively. Patient-specific simulations demonstrated that 67.1–83.3% of the volumes of four sample 3.3–11.4 cm³ tumors could be ablated within 3–10 min using transgastric or intraductal approaches. 55.3–60.0% of the volume of a large 17.0 cm³ tumor could be ablated using multiple applicator positions within 20–30 min with either transgastric or intraductal approaches. 89.9–94.7% of the volume of two 4.4–11.4 cm³ tumors could be treated with intraductal hyperthermia. Sectored applicators are effective in directing acoustic output away from and preserving sensitive structures. When acoustic energy is directed towards sensitive structures, applicators should be placed at least 13.9–14.8 mm from major vessels like the aorta, 9.4–12.0 mm from other vessels, depending on the vessel size and flow rate, and 14 mm from the duodenum.

Conclusions

This study demonstrated the feasibility of generating shaped or conformal ablative or hyperthermic temperature distributions within pancreatic tumors using transgastric or intraductal ultrasound.

A review of radiofrequency ablation: Large target tissue necrosis and mathematical modelling

Radiofrequency ablation (RFA) is an effective clinical method for tumour ablation with minimum intrusiveness. However, the use of RFA is mostly restricted to small tumours, especially those <3cm in diameter. This paper discusses the state-of-the-art of RFA, drawn from experimental and clinical results, for large tumours (i.e. ⩾3cm in diameter). In particular, the paper analyses clinical results related to target tissue necrosis (TTN) and mathematical modelling of the RFA procedure to understand the mechanism whereby the TTN is limited to under 3cm with RFA. This paper also discusses a strategy of controlling of the temperature of target tissue in the RFA procedure with the state-of-art device, which has the potential to increase the size of TTN. This paper ends with a discussion of some future ideas to solve the so-called 3-cm problem with RFA.

MRI-guided focused ultrasound ablation of lumbar medial branch nerve: Feasibility and safety study in a swine model

PURPOSE

About 10-40% of chronic low back pain cases involve facet joints, which are commonly treated with lumbar medial branch (MB) radiofrequency neurotomy. Magnetic resonance imaging-guided focused ultrasound (MRgFUS), a non-invasive, non-ionising ablation modality used to treat tumours, neuropathic pain and painful bone metastasis can also be used to disrupt nerve conduction. This work's purpose was to study the feasibility and safety of direct MRgFUS ablation of the lumbar MB nerve in acute and subacute swine models.

MATERIALS AND METHODS

In vivo MRgFUS ablation was performed in six swine (three acute and three subacute) using a clinical MRgFUS system and a 3-T MRI scanner combination. Behavioural assessment was performed, and imaging and histology were used to assess the treatment.

RESULTS AND CONCLUSIONS

Histological analysis of the in vivo studies confirmed thermal necrosis of the MB nerve could be achieved without damaging the spinal cord or adjacent nerve roots. MRgFUS did not cause changes in the animals' behaviour or ambulation.

Modelling the temperature evolution of bone under high intensity focused ultrasound

Magnetic resonance-guided high intensity focused ultrasound (MR-HIFU) has been clinically shown to be effective for palliative pain management in patients suffering from skeletal metastasis. The underlying mechanism is supposed to be periosteal denervation caused by ablative temperatures reached through ultrasound heating of the cortex. The challenge is exact temperature control during sonication as MR-based thermometry approaches for bone tissue are currently not available. Thus, in contrast to the MR-HIFU ablation of soft tissue, a thermometry feedback to the HIFU is lacking, and the treatment of bone metastasis is entirely based on temperature information acquired in the soft tissue adjacent to the bone surface. However, heating of the adjacent tissue depends on the exact sonication protocol and requires extensive modelling to estimate the actual temperature of the cortex. Here we develop a computational model to calculate the spatial temperature evolution in bone and the adjacent tissue during sonication. First, a ray-tracing technique is used to compute the heat production in each spatial point serving as a source term for the second part, where the actual temperature is calculated as a function of space and time by solving the Pennes bio-heat equation. Importantly, our model includes shear waves that arise at the bone interface as well as all geometrical considerations of transducer and bone geometry. The model was compared with a theoretical approach based on the far field approximation and an MR-HIFU experiment using a bone phantom. Furthermore, we investigated the contribution of shear waves to the heat production and resulting temperatures in bone. The temperature evolution predicted by our model was in accordance with the far field approximation and agreed well with the experimental data obtained in phantoms. Our model allows the simulation of the HIFU treatments of bone metastasis in patients and can be extended to a planning tool prior to MR-HIFU treatments.

Thermal therapy of pancreatic tumours using endoluminal ultrasound: Parametric and patient-specific modelling

PURPOSE

The aim of this study is to investigate endoluminal ultrasound applicator configurations for volumetric thermal ablation and hyperthermia of pancreatic tumours using 3D acoustic and biothermal finite element models.

MATERIALS AND METHODS

Parametric studies compared endoluminal heating performance for varying applicator transducer configurations (planar, curvilinear-focused, or radial-diverging), frequencies (1-5 MHz), and anatomical conditions. Patient-specific pancreatic head and body tumour models were used to evaluate feasibility of generating hyperthermia and thermal ablation using an applicator positioned in the duodenal or stomach lumen. Temperature and thermal dose were calculated to define ablation (> 240 EM(43 °C)) and moderate hyperthermia (40-45 °C) boundaries, and to assess sparing of sensitive tissues. Proportional-integral control was incorporated to regulate maximum temperature to 70-80 °C for ablation and 45 °C for hyperthermia in target regions.

RESULTS

Parametric studies indicated that 1-3 MHz planar transducers are the most suitable for volumetric ablation, producing 5-8 cm(3) lesion volumes for a stationary 5-min sonication. Curvilinear-focused geometries produce more localised ablation to 20-45 mm depth from the GI tract and enhance thermal sparing (T(max) < 42 °C) of the luminal wall. Patient anatomy simulations show feasibility in ablating 60.1-92.9% of head/body tumour volumes (4.3-37.2 cm(3)) with dose < 15 EM(43 °C) in the luminal wall for 18-48 min treatment durations, using 1-3 applicator placements in GI lumen. For hyperthermia, planar and radial-diverging transducers could maintain up to 8 cm(3) and 15 cm(3) of tissue, respectively, between 40-45 °C for a single applicator placement.

CONCLUSIONS

Modelling studies indicate the feasibility of endoluminal ultrasound for volumetric thermal ablation or hyperthermia treatment of pancreatic tumour tissue.

Catheter-based ultrasound technology for image-guided thermal therapy: current technology and applications

Catheter-based ultrasound (CBUS) is applied to deliver minimally invasive thermal therapy to solid cancer tumours, benign tissue growth, vascular disease, and tissue remodelling. Compared to other energy modalities used in catheter-based surgical interventions, unique features of ultrasound result in conformable and precise energy delivery with high selectivity, fast treatment times, and larger treatment volumes. We present a concise review of CBUS technology being currently utilized in animal and clinical studies or being developed for future applications. CBUS devices have been categorised into interstitial, endoluminal and endovascular/cardiac applications. Basic applicator designs, site-specific evaluations and possible treatment applications have been discussed in brief. Particular emphasis has been given to ablation studies that incorporate image guidance for applicator placement, therapy monitoring, feedback control, and post-procedure assessment. Examples of devices included here span the entire spectrum of the development cycle from preliminary simulation-based design studies to implementation in clinical investigations. The use of CBUS under image guidance has the potential for significantly improving precision and applicability of thermal therapy delivery.

Review of temperature dependence of thermal properties, dielectric properties, and perfusion of biological tissues at hyperthermic and ablation temperatures

The application of supraphysiological temperatures (>40°C) to biological tissues causes changes at the molecular, cellular, and structural level, with corresponding changes in tissue function and in thermal, mechanical and dielectric tissue properties. This is particularly relevant for image-guided thermal treatments (e.g. hyperthermia and thermal ablation) delivering heat via focused ultrasound (FUS), radiofrequency (RF), microwave (MW), or laser energy; temperature induced changes in tissue properties are of relevance in relation to predicting tissue temperature profile, monitoring during treatment, and evaluation of treatment results. This paper presents a literature survey of temperature dependence of electrical (electrical conductivity, resistivity, permittivity) and thermal tissue properties (thermal conductivity, specific heat, diffusivity). Data of soft tissues (liver, prostate, muscle, kidney, uterus, collagen, myocardium and spleen) for temperatures between 5 to 90°C, and dielectric properties in the frequency range between 460 kHz and 3 GHz are reported. Furthermore, perfusion changes in tumors including carcinomas, sarcomas, rhabdomyosarcoma, adenocarcinoma and ependymoblastoma in response to hyperthmic temperatures up to 46°C are presented. Where appropriate, mathematical models to describe temperature dependence of properties are presented. The presented data is valuable for mathematical models that predict tissue temperature during thermal therapies (e.g. hyperthermia or thermal ablation), as well as for applications related to prediction and monitoring of temperature induced tissue changes.