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

Modeling for neurosurgical laser interstitial thermal therapy with and without intracranial recording electrodes

Laser thermal ablation has become a prominent neurosurgical treatment approach, but in epilepsy patients it cannot currently be safely implemented with intracranial recording electrodes that are used to study interictal or epileptiform activity. There is a pressing need for computational models of laser interstitial thermal therapy (LITT) with and without intracranial electrodes to enhance the efficacy and safety of optical neurotherapies. In this paper, we aimed to build a biophysical bioheat and ray optics model to study the effects of laser heating in the brain, with and without intracranial electrodes in the vicinity of the ablation zone during the LITT procedure. COMSOL Multiphysics finite element method (FEM) solver software was used to create a bioheat thermal model of brain tissue, with and without blood flow incorporation via Penne's model, to model neural tissue response to laser heating. We report that the close placement of intracranial electrodes can increase the maximum temperature of the brain tissue volume as well as impact the necrosis region volume if the electrodes are placed too closely to the laser coupled diffuse fiber tip. The model shows that an electrode displacement of 4 mm could be considered a safe distance of intracranial electrode placement away from the LITT probe treatment area. This work, for the first time, models the impact of intracranially implanted recording electrodes during LITT, which could improve the understanding of the LITT treatment procedure on the brain's neural networks a sufficient safe distance to the implanted intracranial recording electrodes. We recommend modeling safe distances for placing the electrodes with respect to the infrared laser coupled diffuse fiber tip.

Experimental Validation of Realistic Measurement Setup for Quantitative UWB-Guided Hyperthermia Temperature Monitoring

Hyperthermia induces slight temperature increase of 4-8 °C inside the tumor, making it more responsive to radiation and drugs, thereby improving the outcome of the oncological treatment. To verify the level of heat in the tumor and to avoid damage of the healthy tissue, methods for non-invasive temperature monitoring are needed. Temperature estimation by means of microwave imaging is of great interest among the scientific community. In this paper, we present the results of experiments based on ultra-wideband (UWB) M-sequence technology. Our temperature estimation approach uses temperature dependency of tissue dielectric properties and relation of UWB images to the reflection coefficient on the boundary between tissue types. The realistic measurement setup for neck cancer hyperthermia considers three antenna arrangements. Data are processed with Delay and Sum beamforming and Truncated Singular Value Decomposition. Two types of experiments are presented in this paper. In the first experiment, relative permittivity of subsequently replaced tumor mimicking material is estimated, and in the second experiment, real temperature change in the tumor imitate is monitored. The results showed that the presented approach allows for qualitative as well as quantitative permittivity and temperature estimation. The frequency range for temperature estimation, preferable antenna configurations, and limitations of the method are indicated.

In vivo electrical properties of the healthy liver and the hepatic tumor in a mouse model between 1 Hz and 1 MHz during a thermal treatment

Objective: Understansding the changing patterns of in vivo electrical properties for the target tissue is crucial for the accurate temperature monitoring and the treatment efficacy in thermal therapy. Our research aims to investigate the changing patterns and the reversibility of in vivo electrical properties for both healthy livers and liver tumors in a mouse model over a frequency range of 1 Hz to 1 MHz at temperatures between 30 °C to 90 °C.

Methods and materials: The mice were anesthetized and the target organ was exposed. An 808-nm near-infrared laser was employed as the heating source to heat the organ in vivo. The four-needle electrode, connected to an impedance analyzer, was utilized to obtain the impedance at varying temperatures, which were monitored by a thermocouple.

Results: The findings indicated a gradual decline in impedance with an increase in temperature. Furthermore, the impedance was normalized to that at 30 °C, and the real part of the normalized impedance was defined as the k-values, which range from 0 to 1. The results demonstrated a linear correlation between k-values and temperatures (R2 > 0.9 for livers and R2 > 0.8 for tumors). Significant differences were observed between livers and tumors at 1, 10 and 50 kHz (p < 0.05). Additionally, it was demonstrated that the electrical properties could be reversed when the temperature was below or equal to 45 °C.

Conclusion: We believe that these results will contribute to the advancement of radiofrequency ablation systems and the development of techniques for temperature monitoring during liver thermal treatment.

Importance of the enhanced cooling system for more spherical ablation zones: Numerical simulation, ex vivo and in vivo validation

INTRODUCTION

This study aimed to investigate the efficacy of a small-gauge microwave ablation antenna (MWA) with an enhanced cooling system (ECS) for generating more spherical ablation zones.

METHODS

A comparison was made between two types of microwave ablation antennas, one with ECS and the other with a conventional cooling system (CCS). The finite element method was used to simulate in vivo ablation. Two types of antennas were used to create MWA zones for 5, 8, 10 min at 50, 60, and 80 W in ex vivo bovine livers (n = 6) and 5 min at 60 W in vivo porcine livers (n = 16). The overtreatment ratio, ablation aspect ratio, carbonization area, and other characteristcs of antennas were measured and compared using numerical simulation and gross pathologic examination.

RESULTS

In numerical simulation, the ECS antenna demonstrated a lower overtreatment ratio than the CCS antenna (1.38 vs 1.43 at 50 W 5 min, 1.19 vs 1.35 at 50 W 8 min, 1.13 vs 1.32 at 50 W 10 min, 1.28 vs 1.38 at 60 W 5 min, 1.14 vs 1.32 at 60 W 8 min, 1.10 vs 1.30 at 60 W 10 min). The experiments revealed that the ECS antenna generated ablation zones with a more significant aspect ratio (0.92 ± 0.03 vs 0.72 ± 0.01 at 50 W 5 min, 0.95 ± 0.02 vs 0.70 ± 0.01 at 50 W 8 min, 0.96 ± 0.01 vs 0.71 ± 0.04 at 50 W 10 min, 0.96 ± 0.01 vs 0.73 ± 0.02 at 60 W 5 min, 0.94 ± 0.03 vs 0.71 ± 0.03 at 60 W 8 min, 0.96 ± 0.02 vs 0.69 ± 0.04 at 60 W 10 min) and a smaller carbonization area (0.00 ± 0.00 cm2 vs 0.54 ± 0.06 cm2 at 50 W 5 min, 0.13 ± 0.03 cm2 vs 0.61 ± 0.09 cm2 at 50 W 8 min, 0.23 ± 0.05 cm2 vs 0.73 ± 0.05 m2 at 50 W 10 min, 0.00 ± 0.00 cm2 vs 1.59 ± 0.41 cm2 at 60 W 5 min, 0.23 ± 0.22 cm2 vs 2.11 ± 0.63 cm2 at 60 W 8 min, 0.57 ± 0.09 cm2 vs 2.55 ± 0.51 cm2 at 60 W 10 min). Intraoperative ultrasound images revealed a hypoechoic area instead of a hyperechoic area near the antenna. Hematoxylin-eosin staining of the dissected tissue revealed a correlation between the edge of the ablation zone and that of the hypoechoic area.

CONCLUSIONS

The ECS antenna can produce more spherical ablation zones with less charring and a clearer intraoperative ultrasound image of the ablation area than the CCS antenna.

Accurate Temperature Reconstruction in Radiofrequency Ablation for Atherosclerotic Plaques Based on Inverse Heat Transfer Analysis

Radio frequency ablation has emerged as a widely accepted treatment for atherosclerotic plaques. However, monitoring the temperature field distribution in the blood vessel wall during this procedure presents challenges. This limitation increases the risk of endothelial cell damage and inflammatory responses, potentially leading to lumen restenosis. The aim of this study is to accurately reconstruct the transient temperature distribution by solving a stochastic heat transfer model with uncertain parameters using an inverse heat transfer algorithm and temperature measurement data. The nonlinear least squares optimization method, Levenberg-Marquardt (LM), was employed to solve the inverse heat transfer problem for parameter estimation. Then, to improve the convergence of the algorithm and reduce the computational resources, a method of parameter sensitivity analysis was proposed to select parameters mainly affecting the temperature field. Furthermore, the robustness and accuracy of the algorithm were verified by introducing random noise to the temperature measurements. Despite the high level of temperature measurement noise (ξ = 5%) and larger initial guess deviation, the parameter estimation results remained closely aligned with the actual values, with an overall ERMS consistently below 0.05. The absolute errors between the reconstruction temperature at the measurement points TC1, TC2, and TC3, and the actual temperature, remained within 0.33 °C, 2.4 °C, and 1.17 °C, respectively. The Levenberg-Marquardt algorithm employed in this study proficiently tackled the ill-posed issue of inversion process and obtained a strong consistency between the reconstructed temperature the actual temperature.

Modeling pulsed dye laser treatment of psoriatic plaques by combining numerical methods and image-derived lesion morphologies

OBJECTIVES

Knowledge of the physical effects of pulsed dye laser (PDL) treatment of psoriatic lesions is essential in unraveling the remedial mechanisms of this treatment and hence also in maximizing in its disease-modifying potential. Therefore, the main objective of this study was to provide estimates of these physical effects (for laser wavelengths of 585 and 595 nm), with the aim of identifying pathogenic processes that may be affected by these conditions.

METHODS

We modeled the laser light propagation and subsequent photothermal heating by numerically solving the transient diffusion and heat equations simultaneously. To this end, we used the finite element method in conjunction with an image-derived psoriatic lesion morphology (which was defined by segmenting blood vessels from a confocal microscopy image of a fluorescently labeled section of a 3 mm punch biopsy of a psoriatic lesion). The resulting predictions of the generated temperature field within the lesion were then used to assess the possibility of stalling or arresting some suspected pathogenic processes.

RESULTS

According to our results, it is conceivable that perivascular nerves are thermally denatured, as almost all locations that reach 60°C were found to be within 18 µm (at 585 nm) and 11 µm (at 595 nm) of a blood vessel wall. Furthermore, activation of TRPV1 and TRPV2 channels in perivascular neuronal and immune cells is highly likely, since a critical temperature of 43°C is generated at locations within up to 350 µm of a vessel wall (at both wavelengths) and sustained for up to 700 ms (at 585 nm) and 40 ms (at 595 nm), while a critical temperature of 52°C is reached by locations within 80 µm (at 585 nm) and 30 µm (at 595 nm) of a vessel wall and sustained for up to 100  ms (at 585 nm) and 30 ms (at 595 nm). Finally, we found that the blood vessel coagulation-inducing temperature of 70°C is sustained in the vascular epithelium for up to 19 and 5 ms at 585 and 595 nm, respectively, rendering partial or total loss of vascular functionality a distinct possibility.

CONCLUSIONS

The presented approach constitutes a useful tool to provide realistic estimates of the photothermal effects of PDL treatment of psoriatic plaques (as well as other selective photothermolysis-based treatments), yielding information that is essential in guiding future experimental studies toward unraveling the remedial mechanisms of these treatments.

Patient-specific temperature distribution prediction in laser interstitial thermal therapy: single-irradiation data-driven method

Laser interstitial thermal therapy (LITT) is popular for treating brain tumours and epilepsy. The strict control of tissue thermal damage extent is crucial for LITT. Temperature prediction is useful for predicting thermal damage extent. Accurately predictingin vivobrain tissue temperature is challenging due to the temperature dependence and the individual variations in tissue properties. Considering these factors is essential for improving the temperature prediction accuracy.Objective. To present a method for predicting patient-specific tissue temperature distribution within a target lesion area in the brain during LITT.Approach. A magnetic resonance temperature imaging (MRTI) data-driven estimation model was constructed and combined with a modified Pennes bioheat transfer equation (PBHE) to predict patient-specific temperature distribution. In the PBHE for temperature prediction, the individual specificity and temperature dependence of thermal tissue properties and blood perfusion, as well as the individual specificity of optical tissue properties were considered. Only MRTI data during one laser irradiation were required in the method. This enables the prediction of patient-specific temperature distribution and the resulting thermal damage region for subsequent ablations.Main results. Patient-specific temperature prediction was evaluated based on clinical data acquired during LITT in the brain, using intraoperative MRTI data as the reference standard. Our method significantly improved the prediction performance of temperature distribution and thermal damage region. The average root mean square error was decreased by 69.54%, the average intraclass correlation coefficient was increased by 37.5%, the average Dice similarity coefficient was increased by 43.14% for thermal damage region prediction.Significance. The proposed method can predict temperature distribution and thermal damage region at an individual patient level during LITT, providing a promising approach to assist in patient-specific treatment planning for LITT in the brain.

The Use of Virtual Tissue Constructs That Include Morphological Variability to Assess the Potential of Electrical Impedance Spectroscopy to Differentiate between Thyroid and Parathyroid Tissues during Surgery

Electrical impedance spectroscopy (EIS) has been proposed as a promising noninvasive method to differentiate healthy thyroid from parathyroid tissues during thyroidectomy. However, previously reported similarities in the in vivo measured spectra of these tissues during a pilot study suggest that this separation may not be straightforward. We utilise computational modelling as a method to elucidate the distinguishing characteristics in the EIS signal and explore the features of the tissue that contribute to the observed electrical behaviour. Firstly, multiscale finite element models (or 'virtual tissue constructs') of thyroid and parathyroid tissues were developed and verified against in vivo tissue measurements. A global sensitivity analysis was performed to investigate the impact of physiological micro-, meso- and macroscale tissue morphological features of both tissue types on the computed macroscale EIS spectra and explore the separability of the two tissue types. Our results suggest that the presence of a surface fascia layer could obstruct tissue differentiation, but an analysis of the separability of simulated spectra without the surface fascia layer suggests that differentiation of the two tissue types should be possible if this layer is completely removed by the surgeon. Comprehensive in vivo measurements are required to fully determine the potential for EIS as a method in distinguishing between thyroid and parathyroid tissues.

A Validated Methodological Approach to Prove the Safety of Clinical Electromagnetic Induction Systems in Magnetic Hyperthermia

The present study focuses on the development of a methodology for evaluating the safety of MNH systems, through the numerical prediction of the induced temperature rise in superficial skin layers due to eddy currents heating under an alternating magnetic field (AMF). The methodology is supported and validated through experimental measurements of the AMF's distribution, as well as temperature data from the torsos of six patients who participated in a clinical trial study. The simulations involved a computational model of the actual coil, a computational model of the cooling system used for the cooling of the patients during treatment, and a detailed human anatomical model from the Virtual Population family. The numerical predictions exhibit strong agreement with the experimental measurements, and the deviations are below the estimated combined uncertainties, confirming the accuracy of computational modeling. This study highlights the crucial role of simulations for translational medicine and paves the way for personalized treatment planning.

The State-of-the-Art and Perspectives of Laser Ablation for Tumor Treatment

Tumors significantly impact individuals' physical well-being and quality of life. With the ongoing advancements in optical technology, information technology, robotic technology, etc., laser technology is being increasingly utilized in the field of tumor treatment, and laser ablation (LA) of tumors remains a prominent area of research interest. This paper presents an overview of the recent progress in tumor LA therapy, with a focus on the mechanisms and biological effects of LA, commonly used ablation lasers, image-guided LA, and robotic-assisted LA. Further insights and future prospects are discussed in relation to these aspects, and the paper proposed potential future directions for the development of tumor LA techniques.

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).

Dynamic immuno-nanomedicines in oncology

Anti-cancer therapeutics have achieved significant advances due to the emergence of immunotherapies that rely on the identification of tumors by the patients' immune system and subsequent tumor eradication. However, tumor cells often escape immunity, leading to poor responsiveness and easy tolerance to immunotherapy. Thus, the potentiated anti-tumor immunity in patients resistant to immunotherapies remains a challenge. Reactive oxygen species-based dynamic nanotherapeutics are not new in the anti-tumor field, but their potential as immunomodulators has only been demonstrated in recent years. Dynamic nanotherapeutics can distinctly enhance anti-tumor immune response, which derives the concept of the dynamic immuno-nanomedicines (DINMs). This review describes the pivotal role of DINMs in cancer immunotherapy and provides an overview of the clinical realities of DINMs. The preclinical development of emerging DINMs is also outlined. Moreover, strategies to synergize the antitumor immunity by DINMs in combination with other immunologic agents are summarized. Last but not least, the challenges and opportunities related to DINMs-mediated immune responses are also discussed.

Genetically encoded mediators for sonogenetics and their applications in neuromodulation

Sonogenetics is an emerging approach that harnesses ultrasound for the manipulation of genetically modified cells. The great penetrability of ultrasound waves enables the non-invasive application of external stimuli to deep tissues, particularly advantageous for brain stimulation. Genetically encoded ultrasound mediators, a set of proteins that respond to ultrasound-induced bio-effects, play a critical role in determining the effectiveness and applications of sonogenetics. In this context, we will provide an overview of these ultrasound-responsive mediators, delve into the molecular mechanisms governing their response to ultrasound stimulation, and summarize their applications in neuromodulation.

A review on radiofrequency, laser, and microwave ablations and their thermal monitoring through fiber Bragg gratings

Thermal ablation of tumors aims to apply extreme temperatures inside the target tissue to achieve substantial tumor destruction in a minimally invasive manner. Several techniques are comprised, classified according to the type of energy source. However, the lack of treatment selectivity still needs to be addressed, potentially causing two risks: i) incomplete tumor destruction and recurrence, or conversely, ii) damage of the surrounding healthy tissue. Therefore, the research herein reviewed seeks to develop sensing systems based on fiber Bragg gratings (FBGs) for thermal monitoring inside the lesion during radiofrequency, laser, and microwave ablation. This review shows that, mainly thanks to multiplexing and minimal invasiveness, FBGs provide an optimal sensing solution. Their temperature measurements are the feedback to control the ablation process and allow to investigate different treatments, compare their outcomes, and quantify the impact of factors such as proximity to thermal probe and blood vessels, perfusion, and tissue type.

Model-Optimizing Radiofrequency Parameters of 3D Finite Element Analysis for Ablation of Benign Thyroid Nodules

Radiofrequency (RF) ablation represents an efficient strategy to reduce the volume of thyroid nodules. In this study, a finite element model was developed with the aim of optimizing RF parameters, e.g., input power and treatment duration, in order to achieve the target volume reduction rate (VRR) for a thyroid nodule. RF ablation is modelled as a coupled electro-thermal problem wherein the electric field is applied to induce tissue heating. The electric problem is solved with the Laplace equation, the temperature distribution is estimated with the Pennes bioheat equation, and the thermal damage is evaluated using the Arrhenius equation. The optimization model is applied to RF electrode with different active tip lengths in the interval from 5 mm to 40 mm at the 5 mm step. For each case, we also explored the influence of tumour blood perfusion rate on RF ablation outcomes. The model highlights that longer active tips are more efficient as they require lesser power and shorter treatment time to reach the target VRR. Moreover, this condition is characterized by a reduced transversal ablation zone. In addition, a higher blood perfusion increases the heat dispersion, requiring a different combination of RF power and time treatment to achieve the target VRR. The model may contribute to an improvement in patient-specific RF ablation treatment.

Full-thickness skin rejuvenation by a novel dual-length microneedle radiofrequency device: A proof-of-concept study using human skin

BACKGROUND

A novel dual-length microneedle radiofrequency (DLMR) device has been developed to achieve full-thickness skin rejuvenation by stimulating the papillary and reticular dermis simultaneously. This device's dual-level targeting concept need to be validated on human skin, although its clinical efficacy has been demonstrated in a previous study.

OBJECTIVES

This study evaluated the dual-depth targeting capability and the ability to induce rejuvenation in each layer of vertical skin anatomy, that is, the epidermis, papillary dermis, and reticular dermis, using full-thickness human facial skin samples.

METHODS

Human facial skin samples were obtained from 13 Asian patients who had facelift surgery. To validate the dual-depth targeting concept, DMLR-treated skin samples were analyzed using a digital microscope, thermal imaging, and hematoloxylin and eosin (H&E) staining immediately after DLMR application. On samples stained with H&E, Masson's tricrome, and Verhoeff-Van Gieson, histological observation and morphometric analysis were performed. Total collagen assay (TCA) and quantitative real-time polymerase chain reaction (qPCR) were used to assess changes in total collagen content and mRNA expression levels of collagen types I/III and vimentin, respectively.

RESULTS

The DLMR device successfully induced thermal stimulation in the papillary and reticular dermis. The thickness, stacks, and dermal-epidermal junction convolution of the epidermis treated with DLMR were significantly increased. Collagen bundles in the dermis treated with DLMR exhibited a notable increase in thickness, density, and horizontal alignment. Dermal collagen levels were significantly higher in the morphometric and TCA data, as well as in the qPCR data for dermal matrix proteins.

CONCLUSIONS

Our DLMR device independently and precisely targeted the papillary and reticular dermis, and it appears to be an effective modality for implementing full-thickness rejuvenation.

Measurement of Thermal Conductivity and Thermal Diffusivity of Porcine and Bovine Kidney Tissues at Supraphysiological Temperatures up to 93 °C

This experimental study aimed to characterize the thermal properties of ex vivo porcine and bovine kidney tissues in steady-state heat transfer conditions in a wider thermal interval (23.2-92.8 °C) compared to previous investigations limited to 45 °C. Thermal properties, namely thermal conductivity (k) and thermal diffusivity (α), were measured in a temperature-controlled environment using a dual-needle probe connected to a commercial thermal property analyzer, using the transient hot-wire technique. The estimation of measurement uncertainty was performed along with the assessment of regression models describing the trend of measured quantities as a function of temperature to be used in simulations involving heat transfer in kidney tissue. A direct comparison of the thermal properties of the same tissue from two different species, i.e., porcine and bovine kidney tissues, with the same experimental transient hot-wire technique, was conducted to provide indications on the possible inter-species variabilities of k and α at different selected temperatures. Exponential fitting curves were selected to interpolate the measured values for both porcine and bovine kidney tissues, for both k and α. The results show that the k and α values of the tissues remained rather constant from room temperature up to the onset of water evaporation, and a more marked increase was observed afterward. Indeed, at the highest investigated temperatures, i.e., 90.0-92.8 °C, the average k values were subject to 1.2- and 1.3-fold increases, compared to their nominal values at room temperature, in porcine and bovine kidney tissue, respectively. Moreover, at 90.0-92.8 °C, 1.4- and 1.2-fold increases in the average values of α, compared to baseline values, were observed for porcine and bovine kidney tissue, respectively. No statistically significant differences were found between the thermal properties of porcine and bovine kidney tissues at the same selected tissue temperatures despite their anatomical and structural differences. The provided quantitative values and best-fit regression models can be used to enhance the accuracy of the prediction capability of numerical models of thermal therapies. Furthermore, this study may provide insights into the refinement of protocols for the realization of tissue-mimicking phantoms and the choice of tissue models for bioheat transfer studies in experimental laboratories.

Mathematical Modelling with the Exact Solution of Three Different Bioheat Conduction Models of a Skin Tissue Shocked by Thermoelectrical Effect

This research deals with the temperature increment and responsiveness of skin tissue to a continuous flow of surface heat induced by a constant-voltage electrical current. The exact analytical solution for the dual-phase-lag (DPL) of bioheat transfer has been obtained. It is used to confine the variables to a limited domain to solve the governing equations. The transition temperature reactions have been measured and investigated. The figures provide a comparison of the Pennes, Tzou models, and Vernotte-Cattaneo models. The numerical results demonstrate the values of the voltage, resistance, electric shock time, and dual-phase-lag time parameters which have significant influences on the distributions of the dynamic and conductive temperature rise through the skin tissue.

Characterization of the Optical and Thermal Properties of Cardiac Tissue as a Function of Temperature

In this work, we devised the first characterization of the optical and thermal properties of ex vivo cardiac tissue as a function of different selected temperatures, ranging from room temperature to hyperthermic and ablative temperatures. The broadband (i.e., from 650 nm to 1100 nm) estimation of the optical properties, i.e., absorption coefficient (μa) and reduced scattering coefficient $({\mu ^{\prime}}_s)$, was performed by means of time-domain diffuse optics. Besides, the measurement of the thermal properties was based on the transient hot-wire technique, employing a dual-needle probe to estimate the tissue thermal conductivity (k), thermal diffusivity (α), and volumetric heat capacity (Cv). Increasing the tissue temperature led to variations in the spectral characteristics of μa (e.g., the redshift of the 780 nm peak, the rise of a new peak at 840 nm, and the formation of a valley at 900 nm). Moreover, an increase in the values of ${\mu ^{\prime}}_s$ was assessed as tissue temperature raised (e.g., for 800 nm, at 25 °C ${\mu ^{\prime}}_s = 9.8{\text{ c}}{{\text{m}}^{{\text{ - 1}}}}$, while at 77 °C ${\mu ^{\prime}}_s = 29.1{\text{ c}}{{\text{m}}^{{\text{ - 1}}}}$). Concerning the thermal properties characterization, k was almost constant in the selected temperature interval. Conversely, α and Cv were subjected to an increase and a decrease with temperature, respectively; thus, they registered values of 0.190 mm2/s and 3.03 MJ/(m3•K) at the maximum investigated temperature (79 °C), accordingly.Clinical Relevance- The experimentally obtained optical and thermal properties of cardiac tissue are useful to improve the accuracy of simulation-based tools for thermal therapy planning. Furthermore, the measured properties can serve as a reference for the realization of tissue-mimicking phantoms for medical training and testing of medical instruments.

Sonomechanobiology: Vibrational stimulation of cells and its therapeutic implications

All cells possess an innate ability to respond to a range of mechanical stimuli through their complex internal machinery. This comprises various mechanosensory elements that detect these mechanical cues and diverse cytoskeletal structures that transmit the force to different parts of the cell, where they are transcribed into complex transcriptomic and signaling events that determine their response and fate. In contrast to static (or steady) mechanostimuli primarily involving constant-force loading such as compression, tension, and shear (or forces applied at very low oscillatory frequencies (≤ 1 Hz) that essentially render their effects quasi-static), dynamic mechanostimuli comprising more complex vibrational forms (e.g., time-dependent, i.e., periodic, forcing) at higher frequencies are less well understood in comparison. We review the mechanotransductive processes associated with such acoustic forcing, typically at ultrasonic frequencies (> 20 kHz), and discuss the various applications that arise from the cellular responses that are generated, particularly for regenerative therapeutics, such as exosome biogenesis, stem cell differentiation, and endothelial barrier modulation. Finally, we offer perspectives on the possible existence of a universal mechanism that is common across all forms of acoustically driven mechanostimuli that underscores the central role of the cell membrane as the key effector, and calcium as the dominant second messenger, in the mechanotransduction process.

Sonoprinting nanoparticles on cellular spheroids via surface acoustic waves for enhanced nanotherapeutics delivery

Nanotherapeutics, on their path to the target tissues, face numerous physicochemical hindrances that affect their therapeutic efficacy. Physical barriers become more pronounced in pathological tissues, such as solid tumors, where they limit the penetration of nanocarriers into deeper regions, thereby preventing the efficient delivery of drug cargo. To address this challenge, we introduce a novel approach that employs surface acoustic wave (SAW) technology to sonoprint and enhance the delivery of nanoparticles onto and into cell spheroids. Our SAW platform is designed to generate focused and unidirectional acoustic waves for creating vigorous acoustic streaming while promoting Bjerknes forces. The effect of SAW excitation on cell viability, as well as the accumulation and penetration of nanoparticles on human breast cancer (MCF 7) and mouse melanoma (YUMM 1.7) cell spheroids were investigated. The high frequency, low input voltage, and contact-free nature of the proposed SAW system ensured over 92% cell viability for both cell lines after SAW exposure. SAW sonoprinting enhanced the accumulation of 100 nm polystyrene particles on the periphery of the spheroids to near four-fold, while the penetration of nanoparticles into the core regions of the spheroids was improved up to three times. To demonstrate the effectiveness of our SAW platform on the efficacy of nanotherapeutics, the platform was used to deliver nanoliposomes encapsulated with the anti-cancer metal compound copper diethyldithiocarbamate (CuET) to MCF 7 and YUMM 1.7 cell spheroids. A three-fold increase in the cytotoxic activity of the drug was observed in spheroids under the effect of SAW, compared to controls. The capacity of SAW-based devices to be manufactured as minuscule wearable patches can offer highly controllable, localized, and continuous acoustic waves to enhance drug delivery efficiency to target tissues.

Use of Cold Atmospheric Plasma in the Treatment of Squamous Cell Carcinoma: in vitro Effects and Clinical Application in Feline Tumors: A Pilot Study

Cold atmospheric plasma (CAP) has shown promising results against squamous cell carcinoma (SCC) in both in vivo and in vitro assays, mainly in humans and mice. Its applicability for treatment of feline tumors, however, remains unknown. This study aimed to evaluate the anticancer effects of CAP on a head and neck squamous cell carcinoma (HNSCC) cell lineage and against a clinical case of cutaneous SCC in a cat. Control and treatment groups employing the HNSCC cell line (SCC-25) were used, the latter exposed to CAP for 60 seconds, 90 seconds, or 120 seconds. The cells were subjected to the MTT assay nitric oxidation assay and thermographic in vitro analyses. The clinical application was performed in one cat with cutaneous SCC (3 sites). The lesions were treated and evaluated by thermographic, histopathological, and immunohistochemical examinations (caspase-3 and TNF-alpha). Treatment of the SCC-25 cells for 90 seconds and 120 seconds resulted in a significant nitrite concentration increase. Decreased cell viability was observed after 24 hours and 48 hours, regardless of exposure time. However, the cell viability reduction observed at 72 hours was significant only in the 120 seconds treatment. In vitro, the temperature decreased for all treatment times, while the plasma induced a slight increase in mean temperature (0.7°C) in the in vivo assay. Two of the 3 clinical tumors responded to the treatment: one with a complete response and the other, partial, while the third (lower lip SCC) remained stable. Both remaining tumors displayed apoptotic areas and increased expression of caspase-3 and TNF-alpha. Adverse effects were mild and limited to erythema and crusting. The CAP exhibited an in vitro anticancer effect on the HNSCC cell line, demonstrated by a dose-dependent cell viability reduction. In vivo, the therapy appears safe and effective against feline cutaneous SCC. The treatment did not result in a clinical response for 1 of 3 lesions (proliferative lower lip tumor), however, a biological effect was still demonstrated by the higher expression of apoptosis indicators.

Gold nanorods assisted photothermal therapy of bladder cancer in mice: A computational study on the effects of gold nanorods distribution at the centre, periphery, and surface of bladder cancer

BACKGROUND AND OBJECTIVES

Gold nanorod-assisted photothermal therapy (GNR-PTT) is a cancer treatment whereby GNRs incorporated into the tumour act as photo-absorbers to elevate the thermal destruction effect. In the case of bladder, there are few possible routes to target the tumour with GNRs, namely peri/intra-tumoural injection and intravesical instillation of GNRs. These two approaches lead to different GNR distribution inside the tumour and can affect the treatment outcome.

METHODOLOGY

The present study investigates the effects of heterogeneous GNR distribution in a typical setup of GNR-PTT. Three cases were considered. Case 1 considered the GNRs at the tumour centre, while Case 2 represents a hypothetical scenario where GNRs are distributed at the tumour periphery; these two cases represent intratumoural accumulation with different degree of GNR spread inside the tumour. Case 3 is achieved when GNRs target the exposed tumoural surface that is invading the bladder wall, when they are delivered by intravesical instillation.

RESULTS

Results indicate that for a laser power of 0.6 W and GNR volume fraction of 0.01%, Case 2 and 3 were successful in achieving complete tumour eradication after 330 and 470 s of laser irradiation, respectively. Case 1 failed to form complete tumour damage when the GNRs are concentrated at the tumour centre but managed to produce complete tumour damage if the spread of GNRs is wider. Results from Case 2 also demonstrated a different heating profile from Case 1, suggesting that thermal ablation during GNR-PTT is dependant on the GNRs distribution inside the tumour. Case 3 shows similar results to Case 2 whereby gradual but uniform heating is observed. Cases 2 and 3 show that uniformly heating the tumour can reduce damage to the surrounding tissues.

CONCLUSIONS

Different GNR distribution associated with the different methods of introducing GNRs to the bladder during GNR-PTT affect the treatment outcome of bladder cancer in mice. Insufficient spreading during intratumoural injection of GNRs can render the treatment ineffective, while administered via intravesical instillation. GNR distribution achieved through intravesical instillation present some advantages over intratumoural injection and is worthy of further exploration.

An in silico assessment on the potential of using saline infusion to overcome non-confluent coagulation zone during two-probe, no-touch bipolar radiofrequency ablation of liver cancer

No-touch bipolar radiofrequency ablation (bRFA) is known to produce incomplete tumour ablation with a 'butterfly-shaped' coagulation zone when the interelectrode distance exceeds a certain threshold. Although non-confluent coagulation zone can be avoided by not implementing the no-touch mode, doing so exposes the patient to the risk of tumour track seeding. The present study investigates if prior infusion of saline into the tissue can overcome the issues of non-confluent or butterfly-shaped coagulation. A computational modelling approach based on the finite element method was carried out. A two-compartment model comprising the tumour that is surrounded by healthy liver tissue was developed. Three cases were considered; i) saline infusion into the tumour centre; ii) one-sided saline infusion outside the tumour; and iii) two-sided saline infusion outside the tumour. For each case, three different saline volumes were considered, i.e. 6, 14 and 22 ml. Saline concentration was set to 15% w/v. Numerical results showed that saline infusion into the tumour centre can overcome the butterfly-shaped coagulation only if the infusion volume is sufficient. On the other hand, one-sided infusion outside the tumour did not overcome this. Two-sided infusion outside the tumour produced confluent coagulation zone with the largest volume. Results obtained from the present study suggest that saline infusion, when carried out correctly, can be used to effectively eradicate liver cancer. This presents a practical solution to address non-confluent coagulation zone typical of that during two-probe bRFA treatment.

Temperature Dependence of Thermal Properties of Ex Vivo Porcine Heart and Lung in Hyperthermia and Ablative Temperature Ranges

This work proposes the characterization of the temperature dependence of the thermal properties of heart and lung tissues from room temperature up to > 90 °C. The thermal diffusivity (α), thermal conductivity (k), and volumetric heat capacity (C_v) of ex vivo porcine hearts and deflated lungs were measured with a dual-needle sensor technique. α and k associated with heart tissue remained almost constant until ~ 70 and ~ 80 °C, accordingly. Above ~ 80 °C, a more substantial variation in these thermal properties was registered: at 94 °C, α and k respectively experienced a 2.3- and 1.5- fold increase compared to their nominal values, showing average values of 0.346 mm²/s and 0.828 W/(m·K), accordingly. Conversely, C_v was almost constant until 55 °C and decreased afterward (e.g., C_v = 2.42 MJ/(m³·K) at 94 °C). Concerning the lung tissue, both its α and k were characterized by an exponential increase with temperature, showing a marked increment at supraphysiological and ablative temperatures (at 91 °C, α and k were equal to 2.120 mm²/s and 2.721 W/(m·K), respectively, i.e., 13.7- and 13.1-fold higher compared to their baseline values). Regression analysis was performed to attain the best-fit curves interpolating the measured data, thus providing models of the temperature dependence of the investigated properties. These models can be useful for increasing the accuracy of simulation-based preplanning frameworks of interventional thermal procedures, and the realization of tissue-mimicking materials.

A coupled thermal-electrical-structural model for balloon-based thermoplasty treatment of atherosclerosis

OBJECTIVES

The outcome of balloon-based atherosclerosis thermoplasty is closely related to the temperature/stress distribution during the treatment. For precise prediction of a required thermal lesion in the heterogeneous and thin atherosclerotic vessel, a numerical model incorporating heat-induced tissue expansion or shrinkage and the strain caused by balloon dilation is necessary.

METHODS

A fully coupled thermal-electrical-structural new model was established. The model features a heterogeneous structure including eccentric plaque, healthy artery and surrounding tissue. Tissue expansion/shrinkage and hyperelasticity material model were taken into consideration. Different heating strategies and plaque mechanical properties were investigated. The temperature distribution was compared with the traditional thermal-electrical coupled model. The possibility of thermoplasty treatment using balloons with different sizes was also explored.

RESULTS

The temperature, the electrical intensity and the stress during the thermoplasty were obtained. Lower stress was found in the heating region where tissue shrinkage occurred. The ablation depth was predicted to be ∼0.42 mm larger without coupling the biomechanical influence. The mechanical properties and input condition significantly affect the temperature and stress distribution considering the small dimensions of the tissue. Besides, with a 12.5% reduction of balloon diameter, the largest Von Mises stress decreases by 25.4%.

CONCLUSIONS

It is confirmed that a coupled thermal-electrical-structural model is needed for precise temperature prediction in the balloon-based thermoplasty of the heterogeneous and thin tissue. The model presented may help with future development of optimized treatment planning considering both ablation depth and minimum stress.

Review of the Delivery Kinetics of Thermosensitive Liposomes

Thermosensitive liposomes (TSL) are triggered nanoparticles that release the encapsulated drug in response to hyperthermia. Combined with localized hyperthermia, TSL enabled loco-regional drug delivery to tumors with reduced systemic toxicities. More recent TSL formulations are based on intravascular triggered release, where drug release occurs within the microvasculature. Thus, this delivery strategy does not require enhanced permeability and retention (EPR). Compared to traditional nanoparticle drug delivery systems based on EPR with passive or active tumor targeting (typically <5%ID/g tumor), TSL can achieve superior tumor drug uptake (>10%ID/g tumor). Numerous TSL formulations have been combined with various drugs and hyperthermia devices in preclinical and clinical studies over the last four decades. Here, we review how the properties of TSL dictate delivery and discuss the advantages of rapid drug release from TSL. We show the benefits of selecting a drug with rapid extraction by tissue, and with quick cellular uptake. Furthermore, the optimal characteristics of hyperthermia devices are reviewed, and impact of tumor biology and cancer cell characteristics are discussed. Thus, this review provides guidelines on how to improve drug delivery with TSL by optimizing the combination of TSL, drug, and hyperthermia method. Many of the concepts discussed are applicable to a variety of other triggered drug delivery systems.

Transient changes during microwave ablation simulation : a comparative shape analysis

Microwave ablation therapy is a hyperthermic treatment for killing cancerous tumours whereby microwave energy is dispersed into a target tissue region. Modelling can provide a prediction for the outcome of ablation, this paper explores changes in size and shape of temperature and Specific absorption rate fields throughout the course of simulated treatment with different probe concepts. Here, an axisymmetric geometry of a probe embedded within a tissue material is created, solving coupled electromagnetic and bioheat equations using the finite element method, utilizing hp discretisation with the NGSolve library. Results show dynamic changes across all metrics, with different responses from different probe concepts. The sleeve probe yielded the most circular specific absorption rate pattern with circularity of 0.81 initially but suffered the largest reduction throughout ablation. Similarly, reflection coefficients differ drastically from their initial values, with the sleeve probe again experiencing the largest change, suggesting that it is the most sensitive the changes in the tissue dielectric properties in these select probe designs. These collective characteristic observations highlight the need to consider dielectric property changes and probe specific responses during the design cycle.

Effects of human tissue acoustic properties, abdominal wall shape, and respiratory motion on ultrasound-mediated hyperthermia for targeted drug delivery to pancreatic tumors

Background

PanDox is a Phase-1 trial of chemotherapeutic drug delivery to pancreatic tumors using ultrasound-mediated hyperthermia to release doxorubicin from thermally sensitive liposomes. This report describes trial-related hyperthermia simulations featuring: (i) new ultrasonic properties of human pancreatic tissues, (ii) abdomen deflections imposed by a water balloon, and (iii) respiration-driven organ motion.

Methods

Pancreas heating simulations were carried out using three patient body models. Pancreas acoustic properties were varied between values found in the literature and those determined from our human tissue study. Acoustic beam distortion was assessed with and without balloon-induced abdomen deformation. Target heating was assessed for static, normal respiratory, and jet-ventilation-controlled pancreas motion.

Results

Human pancreatic tumor attenuation is 63% of the literature values, so that pancreas treatments require commensurately higher input intensity to achieve adequate hyperthermia. Abdominal wall deformation decreased the peak field pressure by as much as 3.5 dB and refracted the focal spot by as much as 4.5 mm. These effects were thermally counteracted by sidelobe power deposition, so the net impact on achieving mild hyperthermia was small. Respiratory motion during moving beam hyperthermia produced localized regions overheated by more than 8.0 °C above the 4.0 °C volumetric goal. The use of jet ventilation reduced this excess to 0.7 °C and yielded temperature field uniformity that was nearly identical to having no respiratory motion.

Conclusion

Realistic modeling of the ultrasonic propagation environment is critical to achieving adequate mild hyperthermia without the use of real time thermometry for targeted drug delivery in pancreatic cancer patients.

Laser-induced heating of polydimethylsiloxane-magnetite nanocomposites for hyperthermic inhibition of triple-negative breast cancer cell proliferation

This paper presents the results of an experimental and computational study of the effects of laser-induced heating provided by magnetite nanocomposite structures that are being developed for the localized hyperthermic treatment of triple-negative breast cancer. Magnetite nanoparticle-reinforced polydimethylsiloxane (PDMS) nanocomposites were fabricated with weight percentages of 1%, 5%, and 10% magnetite nanoparticles. The nanocomposites were exposed to incident Near Infrared (NIR) laser beams with well-controlled powers. The laser-induced heating is explored in: (i) heating liquid media (deionized water and cell growth media [Leibovitz L15+]) to characterize the photothermal properties of the nanocomposites, (ii) in vitro experiments that explore the effects of localized heating on triple-negative breast cancer cells, and (iii) experiments in which the laser beams penetrate through chicken tissue to heat up nanocomposite samples embedded at different depths beneath the chicken skin. The resulting plasmonic laser-induced heating is explained using composite theories and heat transport models. The results show that the laser/nanocomposite interactions decrease the viability of triple-negative breast cancer cells (MDA-MB-231) at temperatures in the hyperthermia domain between 41 and 44°C. Laser irradiation did not cause any observed physical damage to the chicken tissue. The potential in vivo performance of the PDMS nanocomposites was also investigated using computational finite element models of the effects of laser/magnetite nanocomposite interactions on the temperatures and thermal doses experienced by tissues that surround the nanocomposite devices. The implications of the results are then discussed for the development of implantable nanocomposite devices for localized treatment of triple-negative breast cancer tissue via hyperthermia.

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.

Gold Nanoparticles-Mediated Photothermal Therapy of Pancreas Using GATE: A New Simulation Platform

This work presents the first investigation of gold nanorods (GNRs)-based photothermal therapy of the pancreas tumor using the Monte Carlo-based code implemented with Geant4 Application for Emission Tomography (GATE). The model of a human pancreas was obtained by segmenting an abdominal computed tomography (CT) scan, and its physical and chemical properties, were obtained from experimental and theoretical data. In GATE, GNRs-mediated hyperthermal therapy, simple heat diffusion as well as interstitial laser ablation were then modeled in the pancreas tumor by defining the optical parameters of this tissue when it is loaded with GNRs. Two different experimental setups on ex vivo pancreas tissue and GNRs-embedded water were devised to benchmark the developed Monte Carlo-based model for the hyperthermia in the pancreas alone and with GNRs, respectively. The influence of GNRs on heat distribution and temperature increase within the pancreas tumor was compared for two different power values (1.2 W and 2.1 W) when the tumor was exposed to 808 nm laser irradiation and with two different laser applicator diameters. Benchmark tests demonstrated the possibility of the accurate simulating of NPs-assisted thermal therapy and reproducing the experimental data with GATE software. Then, the output of the simulated GNR-mediated hyperthermia emphasized the importance of the precise evaluation of all of the parameters for optimizing the preplanning of cancer thermal therapy. Simulation results on temperature distribution in the pancreas tumor showed that the temperature enhancement caused by raising the power was increased with time in both the tumor with and without GNRs, but it was higher for the GNR-load tumor compared to tumor alone.

Radiofrequency impedance monitoring of the biological tissues under local heating by optical radiation

Development of methods for simultaneous control of state of biological tissues during optical treatment is the important tasks in laser surgery. We introduce a novel approach for the monitoring of the state of biological tissues in the process of its local heating by optical radiation. It is based on measurements of the electrical radiofrequency impedance kinetics of the sample during irradiation. The obtained data are processed using interconnected mathematical modeling of corresponding thermodynamic, optical and electrical problems. Experimental applications of this approach, represented in the paper, reveal its high sensitivity, repeatability and consistency with the model. The introduced method can be used for the selection and optimization of radiation parameters of medical laser sources as an alternative or an addition to histological techniques. Radiofrequency impedance measurement can be used directly in the course of surgical operations for monitoring the treated tissues state, including its temperature and degree of damage.

Effects of neoadjuvant hyperthermic intraperitoneal chemotherapy on chemotherapy response score and recurrence in high-grade serous ovarian cancer patients with advanced disease: A multicentre retrospective cohort study

OBJECTIVE

To investigate whether the combination of neoadjuvant hyperthermic intraperitoneal chemotherapy (NHIPEC) plus intravenous neoadjuvant chemotherapy (IV NACT) has superior efficacy to IV NACT alone.

DESIGN

Retrospective cohort study.

SETTING

Two tertiary referral university hospitals.

POPULATION

Patients with ovarian cancer who received NACT-interval debulking surgery (IDS) between 2012 and 2020.

METHODS

The tumour response to NACT was evaluated with the chemotherapy response score (CRS) system. Survival outcomes were compared.

MAIN OUTCOME MEASURES

CRS 3, progression-free survival (PFS), and overall survival (OS).

RESULTS

In total, 127 patients were included, and 46 received NHIPEC plus IV NACT. The addition of NHIPEC was independently associated with an increased likelihood of CRS 3 (p = 0.033). Patients who received NHIPEC + IV NACT had significantly improved PFS compared with those who received IV NACT alone (median PFS: 22 versus 16 months, p < 0.001). The use of NHIPEC was identified as an independent predictor of PFS (p < 0.0001). OS did not differ significantly between treatment groups (p = 0.062), although a trend favouring NHIPEC was noted. Incidence of grade 3-4 adverse events and the surgical complexity score of IDS were similar between the two groups.

CONCLUSIONS

Compared with IV NACT alone, the combination of NHIPEC and IV NACT resulted in improved tumour response and longer PFS. The addition of NHIPEC did not increase the risk of adverse effects or affect the complexity of IDS.

Temperature dependence of dielectric properties of blood at 10 Hz-100 MHz

The temperature dependence of the dielectric properties of blood is important for studying the biological effects of electromagnetic fields, electromagnetic protection, disease diagnosis, and treatment. However, owing to the limitations of measurement methods, there are still some uncertainties regarding the temperature characteristics of the dielectric properties of blood at low and medium frequencies. In this study, we designed a composite impedance measurement box with high heat transfer efficiency that allowed for a four/two-electrode measurement method. Four-electrode measurements were carried out at 10 Hz-1 MHz to overcome the influence of electrode polarization, and two-electrode measurements were carried out at 100 Hz-100 MHz to avoid the influence of distribution parameters, and the data was integrated to achieve dielectric measurements at 10 Hz-100 MHz. At the same time, the temperature of fresh blood from rabbits was controlled at 17-39°C in combination with a temperature-controlled water sink. The results showed that the temperature coefficient for the real part of the resistivity of blood remained constant from 10 Hz to 100 kHz (-2.42%/°C) and then gradually decreased to -0.26%/°C. The temperature coefficient of the imaginary part was positive and bimodal from 6.31 kHz to 100 MHz, with peaks of 5.22%/°C and 4.14%/°C at 126 kHz and 39.8 MHz, respectively. Finally, a third-order function model was developed to describe the dielectric spectra at these temperatures, in which the resistivity parameter in each dispersion zone decreased linearly with temperature and each characteristic frequency increased linearly with temperature. The model could estimate the dielectric properties at any frequency and temperature in this range, and the maximum error was less than 1.39%, thus laying the foundation for subsequent studies.

Complex Permittivity of Ex-Vivo Human, Bovine and Porcine Brain Tissues in the Microwave Frequency Range

Accurate knowledge about the dielectric properties of biological tissues in the microwave frequency range may lead to advancement of biomedical applications based on microwave technology. However, the published data are very scarce, especially for human brain tissues. The aim of this work was to measure and report the complex permittivity of brain white matter, grey matter and cerebellum. Complex permittivity was measured on human, bovine and porcine brain tissues in the microwave frequency range from 0.5 to 18 GHz using an open-ended coaxial probe. The results present a valuable addition to the available data on the brain tissue complex permittivity. Some noticeable variations between the results lead to several conclusions. Complex permittivity variation within the same tissue type of the individual species was comparable to interspecies variation. The difference was prominent between human brains obtained from autopsies, while bovine brains obtained from healthy animals showed very similar complex permittivity. We hypothesize that the difference might have been caused by the basic pathologies of the patients, where the associated therapies could have affected the brain water content. We also examined the effect of excised tissue degradation on its complex permittivity over the course of three days, and the results suggest the gradual dehydration of the samples.

[Material design and temperature field simulation analysis of tumor radiofrequency ablation needle]

To solve the problems of small one-time ablation range and easy charring of the tissue around the electrode associated with the tumor radiofrequency ablation needle, based on the multiphysical field coupling analysis software COMSOL, the effects of needle material, the number of sub needles and the bending angle of sub needles on the ablation effect of radiofrequency ablation electrode needle were studied. The results show that compared with titanium alloy and stainless steel, nickel titanium alloy has better radiofrequency energy transmission efficiency and it is the best material for electrode needle. The number of sub needles has a great influence on the average necrosis depth and the maximum necrosis diameter. Under the same conditions, the more the number of sub needles, the larger the volume of coagulation necrosis area. The bending angle of the needle has a great effect on the maximum diameter of the coagulated necrotic area, but has little effect on the average necrotic depth. Under the same other conditions, the coagulation necrosis area formed by ablation increased with the increase of the bending angle of the sub needle. For the three needles with bending angles of 60 °, 90 ° and 120 ° analyzed in this paper, the one with bending angle of 120 ° can obtain the largest coagulation necrosis area. In general, the design of nickel titanium alloy with 120 ° bending 8-pin is the optimal. The average depth of radiofrequency ablation necrosis area is 32.40 mm, and the maximum necrosis diameter is 52.65 mm. The above optimized design parameters can provide guidance for the structure and material design of tumor radiofrequency ablation needle.

Validating a simulation model for laser-induced thermotherapy using MR thermometry

OBJECTIVES

We want to investigate whether temperature measurements obtained from MR thermometry are accurate and reliable enough to aid the development and validation of simulation models for Laser-induced interstitial thermotherapy (LITT).

METHODS

Laser-induced interstitial thermotherapy (LITT) is applied to ex-vivo porcine livers. An artificial blood vessel is used to study the cooling effect of large blood vessels in proximity to the ablation zone. The experimental setting is simulated using a model based on partial differential equations (PDEs) for temperature, radiation, and tissue damage. The simulated temperature distributions are compared to temperature data obtained from MR thermometry.

RESULTS

The overall agreement between measurement and simulation is good for two of our four test cases, while for the remaining cases drift problems with the thermometry data have been an issue. At higher temperatures local deviations between simulation and measurement occur in close proximity to the laser applicator and the vessel. This suggests that certain aspects of the model may need some refinement.

CONCLUSION

Thermometry data is well-suited for aiding the development of simulations models since it shows where refinements are necessary and enables the validation of such models.

Probe Contact Force Monitoring during Conductivity Measurements of the Left Atrial Appendage to Support the Design of Novel Diagnostic and Therapeutic Procedures

The electrical properties of many biological tissues are freely available from the INRC and the IT'IS databases. However, particularly in lower frequency ranges, few studies have investigated the optimal measurement protocol or the key confounders that need to be controlled, monitored, and reported. However, preliminary work suggests that the contact force of the measurement probe on the tissue sample can affect the measurements. The aim of this paper is to investigate the conductivity change due to the probe contact force in detail. Twenty ex vivo bovine heart samples are used, and conductivity measurements are taken in the Left Atrial Appendage, a common target for medical device developments. The conductivity measurements reported in this work (between 0.14 S/m and 0.24 S/m) align with the literature. The average conductivity is observed to change by -21% as the contact force increases from 2 N to 10 N. In contrast, in conditions where the fluid concentration in the measurement area is expected to be lower, very small changes are observed (less than 2.5%). These results suggest that the LAA conductivity is affected by the contact force due to the fluid concentration in the tissue. This work suggests that contact force should be controlled for in all future experiments.

Radiofrequency impedance spectroscopy of biological tissues under heating by homogeneous laser radiation

We have developed an original method to measure the temperature dependence of the biological tissue electrical properties based on its heating by homogeneous optical radiation.The conventional approach in hyperthermia research involves heating due to a thermal conductivity through the surface of the sample. The novel technique is based on the emission of heat sources in the sample volume caused by the absorption of optical radiation.The method was verified using chicken liver and aloe parenchyma samples, which were uniformly irradiated in a special chamber with an optically scattering inner coating. The electrical impedance of the samples was measured using a 4-electrode technique in the frequency range 100 Hz-1 MHz. In order to approximate and analyze the electrical impedance module, an equivalent electrical circuit based on the Cole-Cole function was used and the dependences of the approximation parameters on time and temperature were obtained. Applying the Arrhenius formulation to the kinetics of low-frequency resistance, we obtained the parameters of the kinetics of degradation of the biological tissues (critical temperatureTcrand activation energyEa):Ea=(16±4)·105J·mol^-1,Tcr=63±1°C for the aloe parenchyma andEa=(4.5±2)·105J·mol^-1,Tcr=83±1°Cfor the chicken liver.

Temperature dependence studies of tissue-mimicking phantoms for ultra-wideband microwave breast tumor detection

Microwave imaging (MWI) systems are being investigated for breast cancer diagnostics as an alternative to conventional X-ray mammography and breast ultrasound. This work aims at a next generation of tissue-mimicking phantoms modelling the temperature-dependent dielectric properties of breast tissue over a large frequency bandwidth. Such phantoms can be used to develop a novel kind of MWI systems that exploit the temperature-dependent permittivity of tissue as a natural contrast agent. Due to the higher water content in tumor tissue, a temperature increase leads to a different change in the complex permittivity compared to surrounding tissue. This will generate a tumor dominated scattering response when the overall tissue temperature increases by a few degrees, e.g. through the use of microwave hyperthermia systems. In that case a differential diagnostic image can be calculated between microwave measurements at reference (around 37◦C) and elevated temperature conditions. This work proposes the design and characterization of agar-oil-glycerin phantoms for fatty, glandular, skin and tumor tissue. The characterization includes measurements with an open-ended coaxial probe and a network analyzer for the frequency range from 50 MHz to 20 GHz in a temperature-controlled environment covering the temperature range from 25◦C to 46◦C. The phantoms show an unique temperature response over the considered frequency bandwidth leading to significant changes in the real and imaginary part of the complex permittivity. Comparative studies with porcine skin and fat tissue show a qualitative agreement.

Temperature and Impedance Variations During Tumor Treating Fields (TTFields) Treatment

Tumor Treating Fields (TTFields) is an FDA-approved cancer treatment technique used for glioblastoma multiforme (GBM). It consists in the application of alternating (100–500 kHz) and low-intensity (1–3 V/cm) electric fields (EFs) to interfere with the mitotic process of tumoral cells. In patients, these fields are applied via transducer arrays strategically positioned on the scalp using the NovoTAL™ system. It is recommended that the patient stays under the application of these fields for as long as possible. Inevitably, the temperature of the scalp increases because of the Joule effect, and it will remain above basal values for most part of the day. Furthermore, it is also known that the impedance of the head changes throughout treatment and that it might also play a role in the temperature variations. The goals of this work were to investigate how to realistically account for these increases and to quantify their impact in the choice of optimal arrays positions using a realistic head model with arrays positions obtained through NovoTAL™. We also studied the impedance variations based on the log files of patients who participated in the EF-14 clinical trial. Our computational results indicated that the layouts in which the arrays were very close to each other led to the appearance of a temperature hotspot that limited how much current could be injected which could consequently reduce treatment efficacy. Based on these data, we suggest that the arrays should be placed at least 1 cm apart from each other. The analysis of the impedance showed that the variations seen during treatment could be explained by three main factors: slow and long-term variations, array placement, and circadian rhythm. Our work indicates that both the temperature and impedance variations should be accounted for to improve the accuracy of computational results when investigating TTFields.

Review of Thermal and Physiological Properties of Human Breast Tissue

Electromagnetic thermal therapies for cancer treatment, such as microwave hyperthermia, aim to heat up a targeted tumour site to temperatures within 40 and 44 °C. Computational simulations used to investigate such heating systems employ the Pennes' bioheat equation to model the heat exchange within the tissue, which accounts for several tissue properties: density, specific heat capacity, thermal conductivity, metabolic heat generation rate, and blood perfusion rate. We present a review of these thermal and physiological properties relevant for hyperthermia treatments of breast including fibroglandular breast, fatty breast, and breast tumours. The data included in this review were obtained from both experimental measurement studies and estimated properties of human breast tissues. The latter were used in computational studies of breast thermal treatments. The measurement methods, where available, are discussed together with the estimations and approximations considered for values where measurements were unavailable. The review concludes that measurement data for the thermal and physiological properties of breast and tumour tissue are limited. Fibroglandular and fatty breast tissue properties are often approximated from those of generic muscle or fat tissue. Tumour tissue properties are mostly obtained from approximating equations or assumed to be the same as those of glandular tissue. We also present a set of reliable data, which can be used for more accurate modelling and simulation studies to better treat breast cancer using thermal therapies.

On the measurement of skeletal muscle anisotropic permittivity property with a single cross-shaped needle insertion

Application of minimally invasive methods to enable the measurement of tissue permittivity in the neuromuscular clinic remain elusive. This paper provides a theoretical and modeling study on the measurement of the permittivity of two-dimensional anisotropic tissues such as skeletal muscle with a multi-electrode cross-shaped needle. For this, we design a novel cross-shaped needle with multiple-electrodes and analyse apparent impedance corresponding to the measured impedance. In addition, we propose three methods of estimate anisotropic muscle permittivity. Compared to existing electrical impedance-based needle methods that we have developed, the new needle design and numerical methods associated enable estimating in vivo muscle permittivity values with only a single needle insertion. Being able to measure muscle permittivity directly with a single needle insertion could open up an entirely new area of research with direct clinical application, including using these values to assist in neuromuscular diagnosis and to assess subtle effects of therapeutic intervention on muscle health.

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.

Experimental Evaluation of Radiation Response and Thermal Properties of NPs-Loaded Tissues-Mimicking Phantoms

Many efforts have recently concentrated on constructing and developing nanoparticles (NPs) as promising thermal agent for optical hyperthermia and photothermal therapy. However, thermal energy transfer in biological tissue is a complex process involving different mechanisms such as conduction, convection, radiation. Therefore, having information about thermal properties of tissue especially when NPs are embedded in is a necessity for predicting the heat transfer during hyperthermia. In this work, the thermal properties of solid phantom based on agar in the presence of three different nanoparticles (BPSi, tNAs, GNRs) and alone were measured and reported as a function of temperature (ranging from 22 to 62 °C). The thermal response of these NPs to an 808 nm laser beam with three different powers were studied in the water comparatively. Agar and tNAs have almost constant thermal properties in the considered range. Among the three NPs, gold has the highest conductivity and diffusivity. At 62 °C BPSi NPs have the similar amount of increase for the diffusivity. The thermal parameters reported in this paper can be useful for the mathematical modeling. Irradiation of the NPs-loaded water phantom displayed the highest radiosensitivity of gold among the three mentioned NPs. However, for the higher power of irradiation, BPSi and tNAs NPs showed the increased absorption of heat during shorter time and the increased temperature gradient slope for the initial 15 s after the irradiation started. The three NPs showed different thermal and irradiation response behavior; however, this comparison study notes the worth of having information about thermal parameters of NPs-loaded tissue for pre-clinical planning.

Computer modeling of radiofrequency cardiac ablation: 30 years of bioengineering research

This review begins with a rationale of the importance of theoretical, mathematical and computational models for radiofrequency (RF) catheter ablation (RFCA). We then describe the historical context in which each model was developed, its contribution to the knowledge of the physics of RFCA and its implications for clinical practice. Next, we review the computer modeling studies intended to improve our knowledge of the biophysics of RFCA and those intended to explore new technologies. We describe the most important technical details of the implementation of mathematical models, including governing equations, tissue properties, boundary conditions, etc. We discuss the utility of lumped element models, which despite their simplicity are widely used by clinical researchers to provide a physical explanation of how RF power is absorbed in different tissues. Computer model verification and validation are also discussed in the context of RFCA. The article ends with a section on the current limitations, i.e. aspects not yet included in state-of-the-art RFCA computer modeling and on future work aimed at covering the current gaps.

Systematic Characterization of High-Power Short-Duration Ablation: Insight From an Advanced Virtual Model

Background: High-power short-duration (HPSD) recently emerged as a new approach to radiofrequency (RF) catheter ablation. However, basic and clinical data supporting its effectiveness and safety is still scarce.

Objective: We aim to characterize HPSD with an advanced virtual model, able to assess lesion dimensions and complications in multiple conditions and compare it to standard protocols.

Methods: We evaluate, on both atrium and ventricle, three HPSD protocols (70 W/8 s, 80 W/6 s, and 90 W/4 s) through a realistic 3D computational model of power-controlled RF ablation, varying catheter tip design (spherical/cylindrical), contact force (CF), blood flow, and saline irrigation. Lesions are defined by the 50°C isotherm contour. Ablations are deemed safe or complicated by pop (tissue temperature >97°C) or charring (blood temperature >80°C). We compared HPSD with standards protocols (30–40 W/30 s). We analyzed the effect of a second HPSD application.

Results: We simulated 432 applications. Most (79%) associated a complication, especially in the atrium. The three HPSD protocols performed similarly in the atrium, while 90 W/4 s appeared the safest in the ventricle. Low irrigation rate led frequently to charring (72%). High-power short-duration lesions were 40–60% shallower and smaller in volume compared to standards, although featuring similar width. A second HPSD application increased lesions to a size comparable to standards.

Conclusion: High-power short-duration lesions are smaller in volume and more superficial than standards but comparable in width, which can be advantageous in the atrium. A second application can produce lesions similar to standards in a shorter time. Despite its narrow safety margin, HPSD seems a valuable new clinical approach.