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

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.

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.

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.