Concentration and volume effects in thermochemical ablation in vivo: Results in a porcine model

Enlarging the thermal coagulation volume during thermochemical ablation with alternating acid-base injection by shortening the injection interval: A computational study

BACKGROUND AND OBJECTIVES

Thermochemical ablation (TCA) is a cancer treatment that utilises the heat released from the neutralisation of acid and base to raise tissue temperature to levels sufficient to induce thermal coagulation. Computational studies have demonstrated that the coagulation volume produced by sequential injection is smaller than that with simultaneous injection. By injecting the reagents in an ensuing manner, the region of contact between acid and base is limited to a thin contact layer sandwiched between the distribution of acid and base. It is hypothesised that increasing the frequency of acid-base injections into the tissue by shortening the injection interval for each reagent can increase the effective area of contact between acid and base, thereby intensifying neutralisation and the exothermic heat released into the tissue.

METHODS

To verify this hypothesis, a computational model was developed to simulate the thermochemical processes involved during TCA with sequential injection. Four major processes that take place during TCA were considered, i.e., the flow of acid and base, their neutralisation, the release of exothermic heat and the formation of thermal damage inside the tissue. Equimolar acid and base at 7.5 M was injected into the tissue intermittently. Six injection intervals, namely 3, 6, 15, 20, 30 and 60 s were investigated.

RESULTS

Shortening of the injection interval led to the enlargement of coagulation volume. If one considers only the coagulation volume as the determining factor, then a 15 s injection interval was found to be optimum. Conversely, if one places priority on safety, then a 3 s injection interval would result in the lowest amount of reagent residue inside the tissue after treatment. With a 3 s injection interval, the coagulation volume was found to be larger than that of simultaneous injection with the same treatment parameters. Not only that, the volume also surpassed that of radiofrequency ablation (RFA); a conventional thermal ablation technique commonly used for liver cancer treatment.

CONCLUSION

The numerical results verified the hypothesis that shortening the injection interval will lead to the formation of larger thermal coagulation zone during TCA with sequential injection. More importantly, a 3 s injection interval was found to be optimum for both efficacy (large coagulation volume) and safety (least amount of reagent residue).

Quantitative dual-energy computed tomography with cesium as a novel contrast agent for localization of thermochemical ablation in phantoms and ex vivo models

BACKGROUND

Thermochemical ablation (TCA) is a minimally invasive therapy under development for hepatocellular carcinoma. TCA simultaneously delivers an acid (acetic acid, AcOH) and base (sodium hydroxide, NaOH) directly into the tumor, where the acid/base chemical reaction produces an exotherm that induces local ablation. However, AcOH and NaOH are not radiopaque, making monitoring TCA delivery difficult.

PURPOSE

We address the issue of image guidance for TCA by utilizing cesium hydroxide (CsOH) as a novel theranostic component of TCA that is detectable and quantifiable with dual-energy CT (DECT).

MATERIALS AND METHODS

To quantify the minimum concentration of CsOH that can be positively identified by DECT, the limit of detection (LOD) was established in an elliptical phantom (Multi-Energy CT Quality Assurance Phantom, Kyoto Kagaku, Kyoto, Japan) with two DECT technologies: a dual-source system (SOMATOM Force, Siemens Healthineers, Forchheim, Germany) and a split-filter, single-source system (SOMATOM Edge, Siemens Healthineers). The dual-energy ratio (DER) and LOD of CsOH were determined for each system. Cesium concentration quantification accuracy was evaluated in a gelatin phantom before quantitative mapping was performed in ex vivo models.

RESULTS

On the dual-source system, the DER and LOD were 2.94 and 1.36-mM CsOH, respectively. For the split-filter system, the DER and LOD were 1.41- and 6.11-mM CsOH, respectively. The signal on cesium maps in phantoms tracked linearly with concentration (R2  = 0.99) on both systems with an RMSE of 2.56 and 6.72 on the dual-source and split-filter system, respectively. In ex vivo models, CsOH was detected following delivery of TCA at all concentrations.

CONCLUSIONS

DECT can be used to detect and quantify the concentration of cesium in phantom and ex vivo tissue models. When incorporated in TCA, CsOH performs as a theranostic agent for quantitative DECT image-guidance.

Pangolin-inspired untethered magnetic robot for on-demand biomedical heating applications

Untethered magnetic miniature soft robots capable of accessing hard-to-reach regions can enable safe, disruptive, and minimally invasive medical procedures. However, the soft body limits the integration of non-magnetic external stimuli sources on the robot, thereby restricting the functionalities of such robots. One such functionality is localised heat generation, which requires solid metallic materials for increased efficiency. Yet, using these materials compromises the compliance and safety of using soft robots. To overcome these competing requirements, we propose a pangolin-inspired bi-layered soft robot design. We show that the reported design achieves heating > 70 °C at large distances > 5 cm within a short period of time <30 s, allowing users to realise on-demand localised heating in tandem with shape-morphing capabilities. We demonstrate advanced robotic functionalities, such as selective cargo release, in situ demagnetisation, hyperthermia and mitigation of bleeding, on tissue phantoms and ex vivo tissues.

Quantitative Dual-Energy CT Image Guidance for Thermochemical Ablation: In Vivo Results in the Rabbit VX2 Model

PURPOSE

To evaluate the feasibility of using dual-energy computed tomography (CT) and theranostic cesium hydroxide (CsOH) for image guidance of thermochemical ablation (TCA) in a rabbit VX2 tumor model.

MATERIALS AND METHODS

In vivo experiments were performed on New Zealand white rabbits, where VX2 tumor fragments (0.3 mL) were inoculated into the right and left flanks (n = 16 rabbits, 32 tumors). Catheters were placed in the approximate center of 1- to 2-cm diameter tumors under ultrasound guidance. TCA was delivered in 1 of 3 treatment groups: untreated control, 5-M TCA, or 10-M TCA. The TCA base reagent was doped with 250-mM CsOH. Dual-energy CT was performed before and after TCA. Cesium (CS)-specific images were postprocessed on the basis of previous phantom calibrations to determine Cs concentration. Line profiles were drawn through the ablation center. Twenty-four hours after TCA, subjects were euthanized, and the resulting damage was evaluated with histopathology.

RESULTS

Cs was detected in 100% of treated tumors (n = 21). Line profiles indicated highest concentrations at the injection site and decreased concentrations at the tumor margins, with no Cs detected beyond the ablation zone. The maximum detected Cs concentration ranged from 14.39 to 137.33 mM. A dose-dependent trend in tissue necrosis was demonstrated between the 10-M TCA and 5-M TCA treatment groups (P = .0005) and untreated controls (P = .0089).

CONCLUSIONS

Dual-energy CT provided image guidance for delivery, localization, and quantification of TCA in the rabbit VX2 model.

An in silico derived dosage and administration guide for effective thermochemical ablation of biological tissues with simultaneous injection of acid and base

BACKGROUND AND OBJECTIVES

Thermochemical ablation (TCA) is a thermal ablation technique involving the injection of acid and base, either sequentially or simultaneously, into the target tissue. TCA remains at the conceptual stage with existing studies unable to provide recommendations on the optimum injection rate, and reagent concentration and volume. Limitations in current experimental methodology have prevented proper elucidation of the thermochemical processes inside the tissue during TCA. Nevertheless, the computational TCA framework developed recently by Mak et al. [Mak et al., Computers in Biology and Medicine, 2022, 145:105494] has opened new avenues in the development of TCA. Specifically, a recommended safe dosage is imperative in driving TCA research beyond the conceptual stage.

METHODS

The aforesaid computational TCA framework for sequential injection was applied and adapted to simulate TCA with simultaneous injection of acid and base at equimolar and equivolume. The developed framework, which describes the flow of acid and base, their neutralisation, the rise in tissue temperature and the formation of thermal damage, was solved numerically using the finite element method. The framework will be used to investigate the effects of injection rate, reagent concentration, volume and type (weak/strong acid-base combination) on temperature rise and thermal coagulation formation.

RESULTS

A higher injection rate resulted in higher temperature rise and larger thermal coagulation. Reagent concentration of 7500 mol/m³ was found to be optimum in producing considerable thermal coagulation without the risk of tissue overheating. Thermal coagulation volume was found to be consistently larger than the total volume of acid and base injected into the tissue, which is beneficial as it reduces the risk of chemical burn injury. Three multivariate second-order polynomials that express the targeted coagulation volume as functions of injection rate and reagent volume, for the weak-weak, weak-strong and strong-strong acid-base combinations were also derived based on the simulated data.

CONCLUSIONS

A guideline for a safe and effective implementation of TCA with simultaneous injection of acid and base was recommended based on the numerical results of the computational model developed. The guideline correlates the coagulation volume with the reagent volume and injection rate, and may be used by clinicians in determining the safe dosage of reagents and optimum injection rate to achieve a desired thermal coagulation volume during TCA.

A computational framework to simulate the thermochemical process during thermochemical ablation of biological tissues

Thermochemical ablation (TCA) is a thermal ablation therapy that utilises heat released from acid-base neutralisation reaction to destroy tumours. This procedure is a promising low-cost solution to existing thermal ablation treatments such as radiofrequency ablation (RFA) and microwave ablation (MWA). Studies have demonstrated that TCA can produce thermal damage that is on par with RFA and MWA when employed properly. Nevertheless, TCA remains a concept that is tested only in a few animal trials due to the risks involved as the result of uncontrolled infusion and incomplete acid-base reaction. In this study, a computational framework that simulates the thermochemical process of TCA is developed. The proposed framework consists of three physics, namely chemical flow, neutralisation reaction and heat transfer. An important parameter in the TCA framework is the neutralisation reaction rate constant, which has values in the order of 10⁸ m³/(mol⋅s). The present study will demonstrate that since the rate constant impacts only the rate and direction of the reaction but has little influence on the extent of reaction, it is possible to replicate the thermochemical process of TCA by employing significantly smaller values of rate constant that are numerically tractable. Comparisons of the numerical results against experimental studies from the literature supports this. The aim of this framework is for researchers to advance and develop TCA to gain an in-depth understanding of the fundamental mechanisms of TCA and to develop a safe treatment protocol of TCA in the hope of advancing TCA into clinical trials.

Nanoparticle transport phenomena in confined flows

Nanoparticles submerged in confined flow fields occur in several technological applications involving heat and mass transfer in nanoscale systems. Describing the transport with nanoparticles in confined flows poses additional challenges due to the coupling between the thermal effects and fluid forces. Here, we focus on the relevant literature related to Brownian motion, hydrodynamic interactions and transport associated with nanoparticles in confined flows. We review the literature on the several techniques that are based on the principles of non-equilibrium statistical mechanics and computational fluid dynamics in order to simultaneously preserve the fluctuation-dissipation relationship and the prevailing hydrodynamic correlations. Through a review of select examples, we discuss the treatments of the temporal dynamics from the colloidal scales to the molecular scales pertaining to nanoscale fluid dynamics and heat transfer. As evident from this review, there, indeed has been little progress made in regard to the accurate modeling of heat transport in nanofluids flowing in confined geometries such as tubes. Therefore the associated mechanisms with such processes remain unexplained. This review has revealed that the information available in open literature on the transport properties of nanofluids is often contradictory and confusing. It has been very difficult to draw definitive conclusions. The quality of work reported on this topic is non-uniform. A significant portion of this review pertains to the treatment of the fluid dynamic aspects of the nanoparticle transport problem. By simultaneously treating the energy transport in ways discussed in this review as related to momentum transport, the ultimate goal of understanding nanoscale heat transport in confined flows may be achieved.

A molecular dynamics approach towards evaluating osmotic and thermal stress in the extracellular environment

OBJECTIVE

A molecular dynamics approach to understanding fundamental mechanisms of combined thermal and osmotic stress induced by thermochemical ablation (TCA) is presented.

METHODS

Structural models of fibronectin and fibronectin bound to its integrin receptor provide idealized models for the effects of thermal and osmotic stress in the extracellular matrix. Fibronectin binding to integrin is known to facilitate cell survival. The extracellular environment produced by TCA at the lesion boundary was modelled at 37 °C and 43 °C with added sodium chloride (NaCl) concentrations (0, 40, 80, 160, and 320 mM). Atomistic simulations of solvated proteins were performed using the GROMOS96 force field and TIP3P water model. Computational results were compared with the results of viability studies of human hepatocellular carcinoma (HCC) cell lines HepG2 and Hep3B under matching thermal and osmotic experimental conditions.

RESULTS

Cell viability was inversely correlated with hyperthermal and hyperosmotic stresses. Added NaCl concentrations were correlated with a root mean square fluctuation increase of the fibronectin arginylglycylaspartic acid (RGD) binding domain. Computed interaction coefficients demonstrate preferential hydration of the protein model and are correlated with salt-induced strengthening of hydrophobic interactions. Under the combined hyperthermal and hyperosmotic stress conditions (43 °C and 320 mM added NaCl), the free energy change required for fibronectin binding to integrin was less favorable than that for binding under control conditions (37 °C and 0 mM added NaCl).

CONCLUSION

Results quantify multiple measures of structural changes as a function of temperature increase and addition of NaCl to the solution. Correlations between cell viability and stability measures suggest that protein aggregates, non-functional proteins, and less favorable cell attachment conditions have a role in TCA-induced cell stress.

Image-guided chemistry altering biology: An in vivo study of thermoembolization

RATIONALE

Advances in image-guided drug delivery for liver cancer have shown a significant survival benefit. However, incomplete treatment is common and residual disease is often found in explanted liver specimens. In addition, the need to treat a malignancy from multiple mechanisms at the same time for optimal outcomes is becoming more widely appreciated. To address this, we hypothesized that an exothermic chemical reaction could be performed in situ. Such a strategy could in principle combine several angles of attack, including ischemia, hyperthermia, acidic protein denaturation, and metabolic modulation of the local environment.

METHODS

The University of Texas MD Anderson Cancer Center Institutional Animal Care and Use Committee approved this study. Outbred swine (25-35 kg, 5 control and 5 experimental) were treated under general anesthesia. Embolization was performed with coaxial microcatheter technique in a segmental hepatic arterial branch using either ethiodized oil as control or with thermoembolic solutionBlood samples were obtained before, immediately after, and the day following the procedure just before CT scans and euthanasia. Livers were explanted and samples were obtained for histologic analysis.

RESULTS

All animals survived the procedure and laboratory values of the control and experimental groups remained within normal limits. The control group had a diffuse or cloudy pattern of attenuation on follow-up CT scan the day after, consistent with gradual antegrade sinusoidal transit of the embolic material. The experimental group had clearly defined vascular casts with some degree of peripheral involvement. At histology, the control group samples had the appearance of normal liver, whereas the experimental group had coagulative necrosis in small pale, punctate areas extending several hundred microns away from the treated vessels and a brisk inflammatory response just outside the margins.

CONCLUSION

In situ chemistry via thermoembolization shows early promise as a fundamentally new tactic for image-guided therapy of solid tumors.

First In Vivo Test of Thermoembolization: Turning Tissue Against Itself Using Transcatheter Chemistry in a Porcine Model

PURPOSE

Embolotherapies are commonly used for management of primary liver cancer. Explant studies of treated livers, however, reveal an untreated tumor in a high fraction of cases. To improve on this, we propose a new concept referred to as thermoembolization. In this technique, the embolic material reacts in local tissues. Highly localized heat energy is released simultaneously with the generation of acid in the target vascular bed. Combined with ischemia, this should provide a multiplexed attack. We report herein our initial results testing the feasibility of this method in vivo.

MATERIALS AND METHODS

Institutional approval was obtained, and three outbred swine were treated in a segmental hepatic artery branch (right or left medial lobe) with thermoembolic material (100, 400, or 500 µL). Solutions (2 or 4 mol/L) of an acid chloride were made using ethiodized oil as the vehicle. Animals were housed overnight, scanned by CT, and euthanized. Necropsy samples of treated tissue were obtained for histologic analysis.

RESULTS

All animals survived the procedure. Vascular stasis occurred rapidly in all cases despite the small volumes used. The lower concentration (2 mol/L) penetrated more distally than the 4 mol/L solution. At CT the following day, vascular casts of ethiodized oil were observed, indicating recanalization had not occurred. Histology specimens demonstrated coagulative necrosis centered on the vessel lumen extending for several hundred microns with a peripheral inflammatory infiltrate.

CONCLUSIONS

Thermoembolization is a new technique for embolization with initial promise. However, results indicate much work must be done to optimize the technique.

Microwave ablation at 915 MHz vs 2.45 GHz: A theoretical and experimental investigation

PURPOSE

The relationship between microwave ablation system operating frequency and ablation performance is not currently well understood. The objective of this study was to comparatively assess the differences in microwave ablation at 915 MHz and 2.45 GHz.

METHODS

Analytical expressions for electromagnetic radiation from point sources were used to compare power deposition at the two frequencies of interest. A 3D electromagnetic-thermal bioheat transfer solver was implemented with the finite element method to characterize power deposition and thermal ablation with asymmetrical insulated dipole antennas (single-antenna and dual-antenna synchronous arrays). Simulation results were validated against experiments in ex vivo tissue.

RESULTS

Theoretical, computational, and experimental results indicated greater power deposition and larger diameter ablation zones when using a single insulated microwave antenna at 2.45 GHz; experimentally, 32±4.1 mm and 36.3±1.0 mm for 5 and 10 min, respectively, at 2.45 GHz, compared to 24±1.7 mm and 29.5±0.6 mm at 915 MHz, with 30 W forward power at the antenna input port. In experiments, faster heating was observed at locations 5 mm (0.91 vs 0.49 °C/s) and 10 mm (0.28 vs 0.15 °C/s) from the antenna operating at 2.45 GHz. Larger ablation zones were observed with dual-antenna arrays at 2.45 GHz; however, the differences were less pronounced than for single antennas.

CONCLUSIONS

Single- and dual-antenna arrays systems operating at 2.45 GHz yield larger ablation zone due to greater power deposition in proximity to the antenna, as well as greater role of thermal conduction.

Development of a Simple Miniature Thermochemical Ablation Device Suitable for Tumor Ablation Research in Rodent Models

Thermochemical ablation is a recently developed minimally invasive method with potential for solid tumor treatment such as in liver cancer. A recently described prototype device, however, is too large for use in the more common rodent models of cancer. In this report we describe a simple, low-cost variant of the device that is easy to assemble, small enough to be readily applicable to small animal models, and then demonstrate its use in an ex vivo model for ablation. It should therefore enable study of the method without requiring specialized equipment or access to a machine shop for device manufacture.

Chemothermal therapy for localized heating and ablation of tumor

Chemothermal therapy is a new hyperthermia treatment on tumor using heat released from exothermic chemical reaction between the injected reactants and the diseased tissues. With the highly minimally invasive feature and localized heating performance, this method is expected to overcome the ubiquitous shortcomings encountered by many existing hyperthermia approaches in ablating irregular tumor. This review provides a relatively comprehensive review on the latest advancements and state of the art in chemothermal therapy. The basic principles and features of two typical chemothermal ablation strategies (acid-base neutralization-reaction-enabled thermal ablation and alkali-metal-enabled thermal/chemical ablation) are illustrated. The prospects and possible challenges facing chemothermal ablation are analyzed. The chemothermal therapy is expected to open many clinical possibilities for precise tumor treatment in a minimally invasive way.

Dual mode single agent thermochemical ablation by simultaneous release of heat energy and acid: hydrolysis of electrophiles

PURPOSE

This study aimed to investigate two readily available electrophilic reagents, acetyl chloride (AcCl), and acetic anhydride (Ac(2)O), for their potential in tissue ablation.

MATERIALS AND METHODS

Reagents were diluted in diglyme as solutions up to 8 mol/L and tested in a gel phantom with NaOH solutions and ex vivo in porcine liver. Temperature, pH, and volume measurements were obtained. Infrared and gross pathological images were obtained in bisected specimens immediately after injection.

RESULTS

AcCl was much more reactive than Ac(2)O and AcCl was therefore used in the tissue studies. Temperature increases of up to 37°C were noted in vitro and 30°C in ex vivo tissues using 4 mol/L AcCl solutions. Experiments at 8 mol/L were abandoned due to the extreme reactivity at this higher concentration. A change in pH of up to 4 log units was noted with 4 mol/L solutions of AcCl with slight recovery over time. Ablated volumes were consistently higher than injected volumes.

CONCLUSIONS

Reaction of electrophiles in tissues shows promise as a new thermochemical ablation technique by means of only a single reagent. Further studies in this area are warranted.

Multiple applicator hepatic ablation with interstitial ultrasound devices: theoretical and experimental investigation

PURPOSE

To evaluate multiple applicator implant configurations of interstitial ultrasound devices for large volume ablation of liver tumors.

METHODS

A 3D bioacoustic-thermal model using the finite element method was implemented to assess multiple applicator implant configurations for thermal ablation with interstitial ultrasound energy. Interstitial applicators consist of linear arrays of up to four 10 mm-long tubular ultrasound transducers, each under separate and dynamic power control, enclosed within a water-cooled delivery catheter (2.4 mm OD). The authors considered parallel implants with two and three applicators (clustered configuration), spaced 2-3 cm apart, to simulate open surgical placement. In addition, the authors considered two applicator implants with applicators converging and diverging at angles of ∼20°, 30°, and 45° to simulate percutaneous placement. Heating experiments (10-15 min) were performed and compared against simulations employing the same experimental parameters. To estimate the performance of parallel, multiple applicator configurations in an in vivo setting, simulations were performed taking into account a range of blood perfusion levels (0, 5, 12, and 15 kg m(-3) s(-1)) that may occur in tumors of varying vascularity. The impact of tailoring the power supplied to individual transducer elements along the length of applicators is explored for applicators inserted in non-parallel (converging and diverging) configurations. Thermal dose (t(43) > 240 min) and temperature thresholds (T > 52 °C) were used to define the ablation zones, with dynamic changes to tissue acoustic and thermal properties incorporated within the model.

RESULTS

Experiments in ex vivo bovine liver yielded ablation zones ranging between 4.0-5.6 cm × 3.2-4.9 cm, in cross section. Ablation zone dimensions predicted by simulations with similar parameters to the experiments were in close agreement (within 5 mm). Simulations of in vivo heating showed that 15 min heating and interapplicator spacing less than 3 cm are required to obtain contiguous, complete ablation zones. The ability to create complete ablation zone profiles for nonparallel implants was illustrated by tailoring applied power levels along the length of applicators.

CONCLUSIONS

Parallel implants consisting of three interstitial ultrasound applicators in a triangular configuration yield complete ablation zones measuring up to 6.2 cm × 5.7 cm after 15 min heating. At larger interapplicator spacing, the level of blood perfusion in the tumor may yield indentations along the periphery of the ablation zone. Tailoring applied power along the length of the applicator can accommodate for nonparallel implants, without compromising safety.