Thermophysical and mechanical properties of biological tissues as a function of temperature: a systematic literature review

Automatic volumetric temperature regulation during in vivo MRI-guided laser-induced thermotherapy (MRg-LITT) with multiple laser probes

BACKGROUND

Clinical Laser-Induced Thermotherapy (LITT) currently lacks precise control of tissue temperature increase during the procedure. This study presents a new method to automatically regulate the maximum temperature increase in vivo at different positions by adjusting LITT power delivered by multiple laser probes using real-time volumetric MR-thermometry.

METHODS

The regulation algorithm was evaluated in vivo on a pig leg muscle. Temperature regulation was performed in volumes surrounding each laser probe tip. The power delivered to each laser probe was automatically adjusted every second using a feedback control algorithm by processing on-the-fly MR-thermometry images (10 slices/second) on a 1.5 T clinical scanner (1.56 mm × 1.56 mm x 3 mm resolution), using the proton-resonance frequency (PRF) shift technique. Several experimental conditions were tested with predefined temperature-time profiles corresponding to conditions of thermal ablation (+30 °C above body temperature) or moderate hyperthermia (+10 and + 15 °C). Control images were acquired after injection of Gadolinium at the end of experiment and were compared with the thermal dose images calculated from the thermometry images.

RESULTS

The mean difference and root mean squared error between target temperatures and measured ones remained below 0.5 °C and 2 °C respectively, for 5 min duration. Lesion sizes observed on thermal dose and on images acquired after gadolinium injection were in good agreement.

CONCLUSION

Automatic regulation of in vivo temperature increase during LITT procedures with multiple laser emitters control is feasible. The method provides an adaptative solution to improve the safety and efficacity of such clinical procedures.

Comparisons of computer simulations and experimental data for capacitive hyperthermia using different split-phantoms

INTRODUCTION

Several positive clinical trials have demonstrated that capacitive hyperthermia (CHT) improves the effectiveness of radiation therapy for the treatment of various cancer entities. However, the ability of CHT to induce significant heating throughout the body is under debate.

OBJECTIVES

To perform a pilot study involving comparisons of computer simulations and experimental data using different split-phantoms to validate hyperthermia treatment modeling for pre-planning for a clinical CHT system and to investigate the feasibility of split-phantom measurements in capacitive hyperthermia.

MATERIALS AND METHODS

The CHT system EHY-2030 (Oncotherm, Budapest, Hungary) was used. The system provides two electrode sizes, but only the smaller electrode, indicated as D200 electrode, was investigated in this pilot study. Horizontally and vertically splittable, different multi-slice phantoms with dielectric material properties simulating muscle and electrically low conductive fat were produced and heated. During the heating procedure, temperature-time curves were measured, and thermal images were captured. Specific absorption rate values were derived from the temperature rise (TR) values. Concomitantly, computer field simulations utilizing a detailed CAD-based model of the CHT system were performed using the simulation platform Sim4Life and compared with measurements.

RESULTS

For the investigated electrode D200 the system power of 75 W was applied, which is half of the maximum power of 150 W and lies in the range of usual values for this electrode applied in patient treatments in our clinic. For 75 W, a heating of 3.6 °C in 6 min in a depth of 1 cm in an agar-based, muscle tissue-equivalent phantom was achieved. The addition of a 1 cm thick, synthetic, low dielectric fat layer reduced the TR up until a depth of 8.5 cm by on average around 38% (from 8.5 cm onwards the absolute local TR is similar, deviations are ≤0.1 °C). In terms of point-to-point absolute SAR comparison (without any normalization), up to a depth of 11 cm in the phantoms central vertical plot, the simulation differs from the measured TR points by on average 25% (ranging from 7% to 36%) for the homogeneous phantom and by on average 43% (ranging from 26% to 60%) for the inhomogeneous phantom.

CONCLUSION

Computer simulations and experimental data were compared for the CHT system EHY-2030 using the D200 electrode, applying a thermal imaging technique for different vertically splittable phantoms. This pilot study data can be used as a guidance regarding the expected heating for this commonly used electrode size but also to further elucidate the significance of non-thermal anticancer effects. Further studies are needed for different sizes and geometries of electrodes and phantoms.

Performance of an Anti-Phase Technology-Powered Microwave Ablation System on Ex Vivo Liver, Lung and Kidney: Analysis of Temperature Trend, Ablation Size and Sphericity

PURPOSE

Investigating the performance of the new Dophi™ M150E Microwave Ablation System, in terms of temperature distribution, ablation size and shape, reproducibility.

MATERIALS AND METHODS

The Dophi™ M150E Microwave Ablation System was tested on ex vivo liver, lung and kidney, at 6 different settings of time, power and number of MW antennas (single antenna: 50 and 100 W at 5 and 10 min; double antenna: 75 W at 5 and 10 min). The temperature distribution was recorded by Fiber Bragg Grating sensors, placed at different distances from the antennas. The ablation axes were measured and the sphericity index was calculated.

RESULTS

The standard deviation of ablation axes was < 5 mm, except at the highest energy and time setting for the lung. A maximum temperature rise of ~ 80 °C was measured. The measured ablation axes are overall comparable with the manufacture's values, especially at lower power and with one MW antenna (average maximum difference is 7 mm). The mean sphericity index of 0.95, 0.79 and 0.9 was obtained for the liver, lung and kidney, respectively, with a single antenna. With double antenna setup, the sphericity index was closer to 1 when 75 W for 10 min were used.

CONCLUSIONS

Dophi™ M150E allows good reproducibility of ablation axes for all cases except in the lung at the highest energy level. With one antenna, an almost spherical ablation area for the liver and kidney was obtained. Using double antenna results in more homogeneous temperature distribution within the tissue compared to single antenna.

Numerical simulation in magnetic resonance imaging radiofrequency dosimetry

Magnetic Resonance Imaging (MRI) employs a radiofrequency electromagnetic field to create pictures on a computer. The prospective biological consequences of exposure to radiofrequency electromagnetic fields (RF EMFs) have not yet been demonstrated, and there is not enough evidence on biological hazards to offer a definite response concerning possible RF health dangers. Therefore, it is crucial to research the health concerns in reaction to RF EMFs, considering the entire exposure in terms of patients receiving MRI. Monitoring increases in temperaturein-vivothroughout MRI scan is extremely invasive and has resulted in a rise in the utilization of computational methods to estimate distributions of temperatures. The purpose of this study is to estimate the absorbed power of the brain exposed to RF in patients undergoing brain MRI scan. A three-dimensional Penne's bio-heat equation was modified to computationally analyze the temperature distributions and potential thermal effects within the brain during MRI scans in the 0.3 T to 1.5 T range (12.77 MHz to 63.87 MHz). The instantaneous temperature distributions of thein-vivotissue in the brain temperatures measured at a time, t = 20.62 s is 0.2 °C and t = 30.92 s is 0.4 °C, while the highest temperatures recorded at 1.03 min and 2.06 min were 0.4 °C and 0.6 °C accordingly. From the temperature distributions of thein-vivotissue in the brain temperatures measured, there is heat build-up in patients who are exposed to electromagnetic frequency ranges, and, consequently, temperature increases within patients are difficult to prevent. The study has, however, indicated that lengthier imaging duration appears to be related to increasing body temperature.

Passive elasticity properties of Octopus rubescens arms

In this report, passive elasticity properties of Octopus rubescens arm tissue are investigated using a multidisciplinary approach encompassing biomechanical experiments, computational modeling, and analyses. Tensile tests are conducted to obtain stress-strain relationships of the arm under axial stretch. Rheological tests are also performed to probe the dynamic shear response of the arm tissue. Based on these tests, comparisons against three different viscoelasticity models are reported.

Calcium-induced upregulation of energy metabolism heats neurons during neural activity

Cellular temperature affects every biochemical reaction, underscoring its critical role in cellular functions. In neurons, temperature not only modulates neurotransmission but is also a key determinant of neurodegenerative diseases. Considering that the brain consumes a disproportionately high amount of energy relative to its weight, neural circuits likely generate a lot of heat, which can increase cytosolic temperature. However, the changes in temperature within neurons and the mechanisms of heat generation during neural excitation remain unclear. In this study, we achieved simultaneous imaging of Ca2+ and temperature using the genetically encoded indicators, B-GECO and B-gTEMP. We then compared the spatiotemporal distributions of Ca2+ responses and temperature. Following neural excitation induced by veratridine, an activator of the voltage-gated Na+ channel, we observed an approximately 2 °C increase in cytosolic temperature occurring 30 s after the Ca2+ response. The temperature elevation was observed in the non-nuclear region, while Ca2+ increased throughout the cell body. Moreover, this temperature increase was suppressed under Ca2+-free conditions and by inhibitors of ATP synthesis. These results indicate that Ca2+-induced upregulation of energy metabolism serves as the heat source during neural excitation.

Bone density estimation using tissue heat capacity

Osteoporosis onset is relatively asymptomatic, the condition often being identified only once a significant fracture occurs, leading to a potentially serious prognosis. Currently, early identification of osteoporosis is complicated by the difficulty in measuring bone density without using x-ray absorptiometry or quantitative ultrasound, so a simpler method for estimating bone density is needed. Given that bone is reported to have a lower specific heat than other tissues, we investigated the possibility of estimating bone density using this difference in tissue thermal properties. The tibia medial surface (shin) and medial malleolus (ankle) of 68 healthy volunteers were cooled using an ice bag, and skin surface temperatures and heat flow were recorded. These measurements were then used to calculate the heat energy transferred per unit temperature. Bone density was estimated by quantitative ultrasound using the T score OSISD, which is the participant's osteo sono-assessment index (OSI) compared to the average OSI of young adults. The heat energy transfer per unit temperature at the shin, but not the ankle, showed a significant negative correlation with T score OSISD (r = -0.413, p = 0.001). Multiple regression analysis showed that heat energy transfer per unit temperature at the shin was a significant predictor of T score OSISD, along with age and height. These results show that tissue thermal property measurements are useful for estimating bone density.

Spatially modulated ablation driven by chaotic attractors in human lung epithelial cancer cells

A significant modification in photoinduced energy transfer in cancer cells is reported by the assistance of a dynamic modulation of the beam size of laser irradiation. Human lung epithelial cancer cells in monolayer form were studied. In contrast to the quantum and thermal ablation effect promoted by a standard focused Gaussian beam, a spatially modulated beam can caused around 15% of decrease in the ablation threshold and formation of a ring-shaped distribution of the photothermal transfer effect. Optical irradiation was conducted in A549 cells by a 532 nm single-beam emerging from a Nd:YVO4 system. Ablation effects derived from spatially modulated convergent waves were controlled by an electrically focus-tunable lens. The proposed chaotic behavior of the spatial modulation followed an Arneodo chaotic oscillator. Fractional dynamic thermal transport was analyzed in order to describe photoenergy in propagation through the samples. Immediate applications of chaos theory for developing phototechnology devices driving biological functions or phototherapy treatments can be considered.

Non-Fourier Bioheat Transfer Analysis in Brain Tissue During Interstitial Laser Ablation: Analysis of Multiple Influential Factors

This work presents the dual-phase lag-based non-Fourier bioheat transfer model of brain tissue subjected to interstitial laser ablation. The finite element method has been utilized to predict the brain tissue's temperature distributions and ablation volumes. A sensitivity analysis has been conducted to quantify the effect of variations in the input laser power, treatment time, laser fiber diameter, laser wavelength, and non-Fourier phase lags. Notably, in this work, the temperature-dependent thermal properties of brain tissue have been considered. The developed model has been validated by comparing the temperature obtained from the numerical and ex vivo brain tissue during interstitial laser ablation. The ex vivo brain model has been further extended to in vivo settings by incorporating the blood perfusion effects. The results of the systematic analysis highlight the importance of considering temperature-dependent thermal properties of the brain tissue, non-Fourier behavior, and microvascular perfusion effects in the computational models for accurate predictions of the treatment outcomes during interstitial laser ablation, thereby minimizing the damage to surrounding healthy tissue. The developed model and parametric analysis reported in this study would assist in a more accurate and precise prediction of the temperature distribution, thus allowing to optimize the thermal dosage during laser therapy in the brain.

Magnetic resonance imaging of blood perfusion rate based on Helmholtz decomposition of heat flux

Objective.Thermal property (TP) maps of human tissues are useful for tumor treatment and diagnosis. In particular, the blood perfusion rate is significantly different for tumors and healthy tissues. Noninvasive techniques that reconstruct TPs from the temperature measured via magnetic resonance imaging (MRI) by solving an inverse bioheat transfer problem have been developed. A few conventional methods can reconstruct spatially varying TP distributions, but they have several limitations. First, most methods require the numerical Laplacian computation of the temperature, and hence they are sensitive to noise. In addition, some methods require the division of a region of interest (ROI) into sub-regions with homogeneous TPs using prior anatomical information, and they assume an unmeasurable initial temperature distribution. We propose a novel robust reconstruction method without the division of an ROI or the assumption of an initial temperature distribution.Approach.The proposed method estimates blood perfusion rate maps from relative temperature changes. This method avoids the computation of the Laplacian by using integral representations of the Helmholtz decomposition of the heat flux.Main Result.We compare the reconstruction results of the conventional and proposed methods using numerical simulations. The results indicate the robustness of the proposed method.Significance.This study suggests the feasibility of thermal property mapping with MRI using the robust proposed method.

Key Biophysical and Physiological Properties Impacting the Oxygenation Status of Breast Cancers During Thermo-radiotherapy

Mild hyperthermia at 39-43 °C for 30-60 min is applied locoregionally to improve the oxygenation status of recurrent breast cancers, thus enhancing the efficacy of radio-, chemo-, and immunotherapy. In this context, estimated (or even conflicting) data are often used in computational modelling of tumour oxygenation and simulation of O2 transport. In this chapter, we present information that may help to improve adjuvant thermotherapy delivered immediately prior to radiotherapy of recurrent breast cancers. Data are preferentially derived from clinical investigations; in some cases, measurements in experimental breast cancers are included.The biophysical properties presented for healthy, mostly postmenopausal, human breast (composite glandular-adipose-fibrous tissue) measured under normothermic (NT) conditions and in therapeutically heated breast cancers include tissue water content and tissue density. In general, averaged values of parameters reported for NT conditions are higher in breast cancers than in normal breast tissue, i.e., all ratios breast cancer/normal breast are >1. Mean values observed in breast cancers during mild hyperthermia (mHT) are consistently higher than those in NT tumours. Parameters determining convective transports in healthy breast tissue and breast cancer include blood flow rates, blood volume, exchanging water space, arterio-venous shunt flow, interstitial fluid flow rate, interstitial fluid pressure, microvascular permeability, interstitial hydraulic conductivity, and interstitial flow velocity. In general, averaged values of parameters measured under NT conditions are higher in breast cancers than in healthy breast. Except for interstitial fluid pressure, these values increase upon mHT treatment of cancers. Prime factors determining and describing the oxygenation status of the healthy breast, and in NT- versus mHT-treated breast cancers, include: oxygen (O2) delivery rates, O2- extractions, O2- consumption rates, subepidermal microvascular HbO2, tissue oxygen solubility, oxygen diffusion coefficients, mean O2 partial pressures pO2, hypoxic fractions HF <5 mmHg, oxygen enhancement ratio, and mitochondrial ROS production. With the exception of the mean pO2, O2 extraction rate and tissue O2 saturation all parameters listed are distinctly higher in breast cancers under NT conditions compared to normal breast. Mild hyperthermia results in therapeutically relevant improvements of the oxygenation status of cancers and enhances mitochondrial ROS production, thus improving radiosensitivity. Note: The oxygenation status of the healthy (postmenopausal) breast is very similar to that of the normal human subcutis.

Research on evolution process of full-layer incision of skin tissue under different laser incidences

Considering difficulties of achieving vertical incidence of beam in different positions of skin, it is significant to study potential effects of incidence angles of laser on incisions. Surgical platform with a 1064 nm continuous fiber laser was established. Incident angle was adopted and real-time temperature fluctuations in laser operating area could be monitored. The rats were treated with laser at day 0 and day 3 after incision modeling, and H&E, Masson, Sirius Red, and Immuno-histochemical staining and enzyme-linked immunosorbent assay were adopted at day 3, 7, 14 to analyze the performance of healing. Laser with energy density of 67.54 J/mm2 can effectively accelerate wound healing in vivo, in which a laser with incident angle around 60° can effectively avoid scar hyperplasia. Therefore, the use of low energy laser with a small deflection angle has a good clinical application prospect in promoting wound healing.

Thermal modulation of skin friction at the finger pad

Preliminary human studies show that reduced skin temperature minimises the risk of mechanically induced skin damage. However, the mechanisms by which cooling enhances skin tolerance to pressure and shear remain poorly understood. We hypothesized that skin cooling below thermo-neutral conditions will decrease kinetic friction at the skin-material interface. To test our hypothesis, we measured the friction coefficient of a thermally pre-conditioned index finger pad sliding at a normal load (5N) across a plate maintained at three different temperatures (38, 24, and 16 °C) in 8 healthy young adults (29±5y). To quantify the temperature distribution of the skin tissue, we used 3D surface scanning and Optical Coherence Tomography to develop an anatomically representative thermal model of the finger. Our group-level data indicated that the sliding finger with thermally affected tissues (up to 8 mm depth) experienced significantly lower frictional forces (p<0.01) at plate temperatures of 16 °C (i.e. 32% decrease) and 24 °C (i.e. 13% decrease) than at 38 °C, respectively. This phenomenon occurred consistently across participants (i.e. N = 6/8, 75%) and without large changes in skin hydration during sliding. Our complementary experimental and theoretical results provide new insights into thermal modulation of skin friction that can be employed for developing thermal technologies to maintain skin integrity under mechanical loading and shearing.

Protection of laryngeal mucosa and function in laryngeal burns by heat absorption of perilaryngeal tissue

OBJECTIVE

The laryngeal tissue carries most of the heat during inhalation injury. This study aims to explore the heat transfer process and the severity of injury inside laryngeal tissue by horizontally studying the temperature rise process at various anatomical layers of the larynx and observing the thermal damage in various parts of the upper respiratory tract.

METHODS

The 12 healthy adult beagles were randomly divided into four groups, and inhaled room temperature air (control group), dry hot air of 80 °C (group I), 160 °C (group II), and 320 °C (group III) for 20 min, respectively. The temperature changes of the glottic mucosal surface, the inner surface of the thyroid cartilage, the external surface of the thyroid cartilage, and subcutaneous tissue were measured every minute. All animals were immediately sacrificed after injury, and pathological changes in various parts of laryngeal tissue were observed and evaluated under a microscope.

RESULTS

After inhaling hot air of 80 °C, 160 °C and 320 °C, the increase of laryngeal temperature in each group was ΔT = 3.57 ± 0.25 °C, 7.83 ± 0.15 °C, 11.93 ± 0.21 °C. The tissue temperature was approximately uniformly distributed, and the difference was not statistically significant. The average laryngeal temperature-time curve showed that the laryngeal tissue temperature in group I and group II showed a trend of "first decrease and then increase", except that the temperature of group III directly increased with time. The prominent pathological changes after thermal burns mainly concluded necrosis of epithelial cells, loss of the mucosal layer, atrophy of submucosal glands, vasodilatation, erythrocytes exudation, and degeneration of chondrocytes. Mild degeneration of cartilage and muscle layers was also observed in mild thermal injury. Pathological scores indicated that the pathological severity of laryngeal burns increased significantly with the increase of temperature, and all layers of laryngeal tissue were seriously damaged by 320 °C hot air.

CONCLUSIONS

The high efficiency of tissue heat conduction enabled the larynx to quickly transfer heat to the laryngeal periphery, and the heat-bearing capacity of perilaryngeal tissue has a certain degree of protective effect on laryngeal mucosa and function in mild to moderate inhalation injury. The laryngeal temperature distribution was in accordance with the pathological severity, and the pathological changes of laryngeal burns provided a theoretical basis for the early clinical manifestations and treatment of inhalation injury.

Hydrogel Micropillar Array for Temperature Sensing in Fluid

Localized temperature sensing and control on a micron-scale have diverse applications in biological systems. We present a micron-sized hydrogel pillar array as potential temperature probes and actuators by exploiting sensitive temperature dependence of their volume change. Soft lithography-based molding processes were presented to fabricate poly N-isopropyl acrylamide (p-NIPAAm)-based hydrogel pillar array on a glass substrate. Au nanorods as light-induced heating elements were embedded within the hydrogel pillars, and near-infrared (NIR) light was used to modulate temperature in a local area. First, static responses of the micro-pillar array were characterized as a function of its temperature. It was shown that the hydrogel had a sensitive volume transition near its low critical solution temperature (LCST). Furthermore, we showed that LCST could be readily adjusted by utilizing copolymerizing with acrylamide (AAM). To demonstrate the feasibility of spatiotemporal temperature mapping and modulation using the presented pillar array, pulsed NIR light was illuminated on a local area of the hydrogel pillar array, and its responses were recorded. Dynamic temperature change in water was mapped based on the abrupt volume change characteristics of the hydrogel pillar, and its potential actuation using NIR light was successfully demonstrated. Considering that the structure can be arrayed in a two-dimensional pixel format with high spatial resolution and high sensitive temperature characteristics, the presented method and the device structure can have diverse applications to change and sense local temperatures in liquid. This is particularly useful in biological systems, where their physiological temperature can be modulated and mapped with high spatial resolution.

Energy Transfer Characteristics of Lipid Bilayer Membranes of Liposomes Examined with Picosecond Time-Resolved Raman Spectroscopy

A number of biochemical reactions proceed inside biomembranes. Since the rate of a chemical reaction is influenced by chemical properties of the surrounding environment, it is important to examine the chemical environment inside the biomembranes. Although the energy transfer characteristics are a basic and important property of a reaction medium, experimental investigation of the thermal conducting capabilities of the biomembranes is a challenging task. We have examined the energy transfer characteristics of lipid bilayer membranes of liposomes, a good model system for the biomembrane, with picosecond time-resolved Raman spectroscopy. The cooling kinetics of the first excited singlet (S1) state of trans-stilbene solubilized within the lipid bilayer membranes is observed as a peak shift of the 1570 cm-1 Raman band of S1 trans-stilbene. The cooling rate constant of S1 trans-stilbene is obtained in six lipid bilayer membranes formed by phospholipids with different hydrocarbon chains, DSPC, DPPC, DMPC, DLPC, DOPC, and egg-PC. We estimate the thermal diffusivity of the lipid bilayer membranes with a known correlation between the cooling rate constant and the thermal diffusivity of the solvent. The thermal diffusivity estimated for the liquid-crystal-phase lipid bilayer membranes is 8.9 × 10-8 to 9.4 × 10-8 m2 s-1, while that for the gel-phase lipid bilayer membranes is 8.4 × 10-8 to 8.5 × 10-8 m2 s-1. The difference in thermal diffusivity between the two phases is explained by a one-dimensional diffusion equation of heat.

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.

Optimization of laser dosimetry based on patient-specific anatomical models for the ablation of pancreatic ductal adenocarcinoma tumor

Laser-induced thermotherapy has shown promising potential for the treatment of unresectable primary pancreatic ductal adenocarcinoma tumors. Nevertheless, heterogeneous tumor environment and complex thermal interaction phenomena that are established under hyperthermic conditions can lead to under/over estimation of laser thermotherapy efficacy. Using numerical modeling, this paper presents an optimized laser setting for Nd:YAG laser delivered by a bare optical fiber (300 µm in diameter) at 1064 nm working in continuous mode within a power range of 2-10 W. For the thermal analysis, patient-specific 3D models were used, consisting of tumors in different portions of the pancreas. The optimized laser power and time for ablating the tumor completely and producing thermal toxic effects on the possible residual tumor cells beyond the tumor margins were found to be 5 W for 550 s, 7 W for 550 s, and 8 W for 550 s for the pancreatic tail, body, and head tumors, respectively. Based on the results, during the laser irradiation at the optimized doses, thermal injury was not evident either in the 15 mm lateral distances from the optical fiber or in the nearby healthy organs. The present computational-based predictions are also in line with the previous ex vivo and in vivo studies, hence, they can assist in the estimation of the therapeutic outcome of laser ablation for pancreatic neoplasms prior to clinical trials.

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.

Theoretical Estimation of Tissue Thermal Response and Associated Thermal Damage During Gold Nanorod-enhanced Photothermal Therapy of Tumors

In the present work, we implemented a computational framework of in vivo gold nanorod (GNR)-enhanced photothermal therapy (PTT) for tumor treatment. The temperature-dependent thermophysical properties of biological tissue and the optical properties of both GNRs and the biological media were included. The latter were modulated during the treatment simulation to account for their variation, from the native to the coagulated state. The contribution of tissue injury-dependent blood perfusion was also considered. The developed model allowed for the estimation of temperature distribution during the photothermal procedure at different procedural settings and amounts of GNRs embedded in the tumor region (i.e., 12.5 μg, 25 μg, and 50 μg). Furthermore, the influence of GNRs on thermal injury, estimated with different damage models, was assessed. The inclusion of GNRs in the tumor entailed an increment of maximum tissue temperature, and faster heating kinetics, as witnessed by the lower time needed to reach complete thermal damage at the tumor center. The percentage of tumor thermal damage evaluated at the end of the simulated treatment was 48%, 69%, and 90%, for PTT in the presence of 12.5 μg, 25 μg, and 50 μg of GNRs, respectively.Clinical Relevance-This establishes that simulation-based tools, modeling the tissue properties variation during the photothermal treatment, can serve as promising preplanning platforms for nanoparticle-assisted light therapies.

Effect of temperature on dynamic compressive behavior of periodontal ligament

Periodontal ligament (PDL) attaches tooth root to the surrounding bone. Its existence between tooth and jaw bone is of utmost importance due to its significant role in absorbing and distributing physiological and para-physiological loading. According to the previous studies, various mechanical tests have been performed to characterize the mechanical properties of the PDL; however, all of them have been done at room temperature. To the best of our knowledge, this is the first study in which the testing was performed at body temperature. The present research was planned to measure the dependency of PDL's viscoelastic behavior on temperature and frequency. Three different temperatures, including body and room temperature, were opted to perform the dynamic compressive tests of the bovine PDL. In addition, a Generalized Maxwell model (GMM) was presented based on empirical outcomes. At 37 °C, amounts of loss factor were found to be greater than those in 25 °C, which demonstrates that the viscous phase of the PDL in higher temperatures plays a critical role. Likewise, by raising the temperature from 25 °C to 37 °C, the model parameters show an enlargement in the viscous part and lessening in the elastic part. It was concluded that the PDL's viscosity in body temperature is much higher than that in room temperature. This model would be functional for a more accurate computational analysis of the PDL at the body temperature (37 °C) in various loading conditions such as orthodontic simulations, mastication, and impact.

Real-time automatic temperature regulation during in vivo MRI-guided laser-induced thermotherapy (MR-LITT)

Precise control of tissue temperature during Laser-Induced Thermotherapy (LITT) procedures has the potential to improve the clinical efficiency and safety of such minimally invasive therapies. We present a method to automatically regulate in vivo the temperature increase during LITT using real-time rapid volumetric Magnetic Resonance thermometry (8 slices acquired every second, with an in-plane resolution of 1.4 mmx1.4 mm and a slice thickness of 3 mm) using the proton-resonance frequency (PRF) shift technique. The laser output power is adjusted every second using a feedback control algorithm (proportional-integral-derivative controller) to force maximal tissue temperature in the targeted region to follow a predefined temperature-time profile. The root-mean-square of the difference between the target temperature and the measured temperature ranged between 0.5 °C and 1.4 °C, for temperature increases between + 5 °C to + 30 °C above body temperature and a long heating duration (up to 15 min), showing excellent accuracy and stability of the method. These results were obtained on a 1.5 T clinical MRI scanner, showing a potential immediate clinical application of such a temperature controller during MR-guided LITT.

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.

Hyperhydration of Cancers: A Characteristic Biophysical Trait Strongly Increasing O2, CO2, Glucose and Lactate Diffusivities, and Improving Thermophysical Properties of Solid Malignancies

Cancers are complex, heterogeneous, dynamic and aggressive diseases exhibiting a series of characteristic biophysical traits which complement the original biological hallmarks of cancers favouring progressive growth, metastasis, and contributing to immune evasion and treatment resistance. One of the prevalent differences between most solid tumors and their corresponding, healthy tissues is a significantly higher water content (hyperhydration) in cancers. As a consequence, cancers have distinctly higher (Fick's) diffusion coefficients D [cm2 s-1] for the respiratory gases O2 and CO2, the key substrate glucose, and for the oncometabolite lactate. In addition, cancers have (a) clearly increased specific heat capacities cp [J g-1 K-1], thus representing high-capacity-tissues upon therapeutic heating induced by electromagnetic irradiation, and (b) higher thermal conductivities k [W m-1 K-1], i.e., increased abilities to conduct heat. Therefore, in diffusion analyses (e.g., when describing critical O2 and glucose supplies or CO2 removal, and the development of hypoxic subvolumes) and for modeling temperature distributions in hyperthermia treatment planning, these specific cancer-related data must be considered in order to reliably reflect oncologic thermo-radiotherapy settings.

Thermo-Visco-Elastometry of RF-Wave-Heated and Ablated Flesh Tissues Containing Au Nanoparticles

We report non-contact laser-based Brillouin light-scattering (BLS) spectroscopy measurements of the viscoelastic properties of hyperthermally radiofrequency (RF)-heated and ablated bovine liver and chicken flesh tissues with embedded gold nanoparticles (AuNPs). The spatial lateral profile of the local surface temperature in the flesh samples during their hyperthermia was measured through optical backscattering reflectometry (OBR) using Mg−silica-NP-doped sensing fibers distributed with an RF applicator and correlated with viscoelastic variations in heat-affected and ablated tissues. Substantial changes in the tissue stiffness after heating and ablation were directly related to their heat-induced structural modifications. The main proteins responsible for muscle elasticity were denatured and irreversibly aggregated during the RF ablation. At T > 100 °C, the proteins constituting the flesh further shrank and became disorganized, leading to substantial plastic deformation of biotissues. Their uniform destruction with larger thermal lesions and a more viscoelastic network was attained via AuNP-mediated RF hyperthermal ablation. The results demonstrated here pave the way for simultaneous real-time hybrid optical sensing of viscoelasticity and local temperature in biotissues during their denaturation and gelation during hyperthermia for future applications that involve mechanical- and thermal-property-controlled theranostics.

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.

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.

Survival benefit of thermal ablation therapy for patients with stage II-III non-small cell lung cancer: A propensity-matched analysis

Background

Thermal ablation (TA) is considered a safe alternative to surgical resection for the treatment of non-small cell lung cancer (NSCLC). While previous studies have shown that TA is beneficial for stage I NSCLC patients, however, few have reported on TA efficacy in patients with stage II-III NSCLC. The current study investigated the impact of TA on the overall survival (OS) and cancer-specific survival (CSS) of patients with stage II-III NSCLC.

Methods

Data on patients with stage II-III NSCLC who did not undergo surgical resection between 2004 and 2015 were extracted from the Surveillance, Epidemiology, and End Results (SEER) database. Propensity score matching (PSM), Kaplan-Meier survival curves, and Cox regression were used for statistical analyses.

Results

A total of 57,959 stage II-III NSCLC patients who did not undergo surgical resection were included in this study, 261 of whom received TA. Overall, TA was associated with a longer OS (p = 0.035) and CSS (p = 0.005) than non-ablation. After 1:3 PSM, 252 patients receiving TA and 732 patients not receiving ablation were enrolled in the matched cohort. The OS (p = 0.047) and CSS (p = 0.029) remained higher in the TA group than in the non-ablation group after PSM. Cox regression analysis showed that age, sex, primary tumor site, pathological type, tumor size, radiotherapy, chemotherapy, and thermal ablation were independently associated with OS and CSS (p <0.05). Subgroup analysis found that the advantages of TA were more pronounced among individuals ≥70 years of age, with tumor size ≤3.0 cm, or who did not receive radiotherapy.

Conclusion

TA could be an effective alternative treatment for stage II-III NSCLC patients unsuitable for surgical resection, particularly those ≥70 years of age, with tumor size ≤3.0 cm, or who have not received radiotherapy.

Estimation of porcine pancreas optical properties in the 600-1100 nm wavelength range for light-based therapies

This work reports the optical properties of porcine pancreatic tissue in the broad wavelength range of 600–1100 nm. Absorption and reduced scattering coefficients (µ_a and µ_s′) of the ex vivo pancreas were obtained by means of Time-domain Diffuse Optical Spectroscopy. We have investigated different experimental conditions—including compression, repositioning, spatial sampling, temporal stability—the effect of the freezing procedure (fresh vs frozen-thawed pancreas), and finally inter-sample variability. Good repeatability under different experimental conditions was obtained (median coefficient of variation less than 8% and ~ 16% for µ_a and µ_s′, respectively). Freezing–thawing the samples caused an irreversible threefold reduction of µ_s′ and no effect on µ_a. The absorption and reduced scattering spectra averaged over different samples were in the range of 0.12–0.74 cm⁻¹ and 12–21 cm⁻¹ with an inter-sample variation of ~ 10% and ~ 40% for µ_a and µ_s′, respectively. The calculated effective transport coefficient (µ_eff) for fresh pancreatic tissue shows that regions between 800–900 nm and 1050–1100 nm are similar and offer the lowest tissue attenuation in the considered range (i.e., µ_eff ranging from 2.4 to 2.7 cm⁻¹). These data, describing specific light-pancreas interactions in the therapeutic optical window for the first time, provide pivotal information for planning of light-based thermotherapies (e.g., laser ablation) and instruction of light transport models for biophotonic applications involving this organ.

Intracellular Heat Transfer and Thermal Property Revealed by Kilohertz Temperature Imaging with a Genetically Encoded Nanothermometer

Despite improved sensitivity of nanothermometers, direct observation of heat transport inside single cells has remained challenging for the lack of high-speed temperature imaging techniques. Here, we identified insufficient temperature resolution under short signal integration time and slow sensor kinetics as two major bottlenecks. To overcome the limitations, we developed B-gTEMP, a nanothermometer based on the tandem fusion of mNeonGreen and tdTomato fluorescent proteins. We visualized the propagation of heat inside intracellular space by tracking the temporal variation of local temperature at a time resolution of 155 μs and a temperature resolution 0.042 °C. By comparing the fast in situ temperature dynamics with computer-simulated heat diffusion, we estimated the thermal diffusivity of live HeLa cells. The present thermal diffusivity in cells was about 1/5.3 of that of water and much smaller than the values reported for bulk tissues, which may account for observations of heterogeneous intracellular temperature distributions.

Unscented Kalman Filtering for Real Time Thermometry During Laser Ablation Interventions

We present a data-assimilation Bayesian framework in the context of laser ablation for the treatment of cancer. For solving the nonlinear estimation of the tissue temperature evolving during the therapy, the Unscented Kalman Filter (UKF) predicts the next thermal status and controls the ablation process, based on sparse temperature information. The purpose of this paper is to study the outcome of the prediction model based on UKF and to assess the influence of different model settings on the framework performances. In particular, we analyze the effects of the time resolution of the filter and the number and the location of the observations. Clinical Relevance - The application of a data-assimilation approach based on limited temperature information allows to monitor and predict in real-time the thermal effects induced by thermal therapy for tumors.

Three-Phase-Lag Bio-Heat Transfer Model of Cardiac Ablation

Significant research efforts have been devoted in the past decades to accurately modelling the complex heat transfer phenomena within biological tissues. These modeling efforts and analysis have assisted in a better understanding of the intricacies of associated biological phenomena and factors that affect the treatment outcomes of hyperthermic therapeutic procedures. In this contribution, we report a three-dimensional non-Fourier bio-heat transfer model of cardiac ablation that accounts for the three-phase-lags (TPL) in the heat propagation, viz., lags due to heat flux, temperature gradient, and thermal displacement gradient. Finite element-based COMSOL Multiphysics software has been utilized to predict the temperature distributions and ablation volumes. A comparative analysis has been conducted to report the variation in the treatment outcomes of cardiac ablation considering different bio-heat transfer models. The effect of variations in the magnitude of different phase lags has been systematically investigated. The fidelity and integrity of the developed model have been evaluated by comparing the results of the developed model with the analytical results of the recent studies available in the literature. This study demonstrates the importance of considering non-Fourier lags within biological tissue for predicting more accurately the characteristics important for the efficient application of thermal therapies.

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.

Model-Based Thermometry for Laser Ablation Procedure Using Kalman Filters and Sparse Temperature Measurements

In this work, we implement a data-assimilation Bayesian framework for the reconstruction of the spatiotemporal profile of the tissue temperature during laser irradiation. The predictions of a physical model simulating the heat transfer in the tissue are associated with sparse temperature measurements, using an Unscented Kalman Filter. We compare a standard state-estimation filtering procedure with a joint-estimation (state and parameters) approach: whereas in the state-estimation only the temperature is evaluated, in the joint-estimation the filter corrects also uncertain model parameters (i.e., the medium thermal diffusivity, and laser beam properties). We have tested the method on synthetic temperature data, and on the temperature measured on agar-gel phantom and porcine liver with fiber optic sensors. The joint-estimation allows retrieving an accurate estimate of the temperature distribution with a maximal error < 1.5 C in both synthetic and liver 1D data, and < 2 C in phantom 2D data. Our approach allows also suggesting a strategy for optimizing the temperature estimation based on the positions of the sensors. Under the constraint of using only two sensors, optimal temperature estimations are obtained when one sensor is placed in proximity of the source, and the other one is in a non-symmetrical position.