A spatiotemporal computational model of focused ultrasound heat-induced nano-sized drug delivery system in solid tumors

Focused Ultrasound (FUS)-triggered nano-sized drug delivery, as a smart stimuli-responsive system for treating solid tumors, is computationally investigated to enhance localized delivery of drug and treatment efficacy. Integration of thermosensitive liposome (TSL), as a doxorubicin (DOX)-loaded nanocarrier, and FUS, provides a promising drug delivery system. A fully coupled partial differential system of equations, including the Helmholtz equation for FUS propagation, bio-heat transfer, interstitial fluid flow, drug transport in tissue and cellular spaces, and a pharmacodynamic model is first presented for this treatment approach. Equations are then solved by finite element methods to calculate intracellular drug concentration and treatment efficacy. The main objective of this study is to present a multi-physics and multi-scale model to simulate drug release, transport, and delivery to solid tumors, followed by an analysis of how FUS exposure time and drug release rate affect these processes. Our findings not only show the capability of model to replicate this therapeutic approach, but also confirm the benefits of this treatment with an improvement of drug aggregation in tumor and reduction of drug delivery in healthy tissue. For instance, the survival fraction of tumor cells after this treatment dropped to 62.4%, because of a large amount of delivered drugs to cancer cells. Next, a combination of three release rates (ultrafast, fast, and slow) and FUS exposure times (10, 30, and 60 min) was examined. Area under curve (AUC) results show that the combination of 30 min FUS exposure and rapid drug release leads to a practical and effective therapeutic response.

Numerical Modeling and Simulation of Non-Invasive Acupuncture Therapy Utilizing Near-Infrared Light-Emitting Diode

Acupuncture is one of the most extensively used complementary and alternative medicine therapies worldwide. In this study, we explore the use of near-infrared light-emitting diodes (LEDs) to provide acupuncture-like physical stimulus to the skin tissue, but in a completely non-invasive way. A computational modeling framework has been developed to investigate the light-tissue interaction within a three-dimensional multi-layer model of skin tissue. Finite element-based analysis has been conducted, to obtain the spatiotemporal temperature distribution within the skin tissue, by solving Pennes' bioheat transfer equation, coupled with the Beer-Lambert law. The irradiation profile of the LED has been experimentally characterized and imposed in the numerical model. The experimental validation of the developed model has been conducted through comparing the numerical model predictions with those obtained experimentally on the agar phantom. The effects of the LED power, treatment duration, LED distance from the skin surface, and usage of multiple LEDs on the temperature distribution attained within the skin tissue have been systematically investigated, highlighting the safe operating power of the selected LEDs. The presented information about the spatiotemporal temperature distribution, and critical factors affecting it, would assist in better optimizing the desired thermal dosage, thereby enabling a safe and effective LED-based photothermal therapy.

hBN Nanoparticle-Assisted Rapid Thermal Cycling for the Detection of Acanthamoeba

Acanthamoeba are widely distributed in the environment and are known to cause blinding keratitis and brain infections with greater than 90% mortality rate. Currently, polymerase chain reaction (PCR) is a highly sensitive and promising technique in Acanthamoeba detection. Remarkably, the rate of heating-cooling and convective heat transfer of the PCR tube is limited by low thermal conductivity of the reagents mixture. The addition of nanoparticles to the reaction has been an interesting approach that could augment the thermal conductivity of the mixture and subsequently enhance heat transfer through the PCR tube. Here, we have developed hexagonal boron nitride (hBN) nanoparticle-based PCR assay for the rapid detection of Acanthamoeba to amplify DNA from low amoeba cell density. As low as 1 × 10^-4 wt % was determined as the optimum concentration of hBN nanoparticles, which increased Acanthamoeba DNA yield up to ~16%. Further, it was able to reduce PCR temperature that led to a ~2.0-fold increase in Acanthamoeba DNA yield at an improved PCR specificity at 46.2 °C low annealing temperature. hBN nanoparticles further reduced standard PCR step time by 10 min and cycles by eight; thus, enhancing Acanthamoeba detection rapidly. Enhancement of Acanthamoeba PCR DNA yield is possibly due to the high adsorption affinity of hBN nanoparticles to purine (Guanine-G) due to the higher thermal conductivity achieved in the PCR mixture due to the addition of hBN. Although further research is needed to demonstrate these findings in clinical application, we propose that the interfacial layers, Brownian motion, and percolation network contribute to the enhanced thermal conductivity effect.

Analysis of thyroid thermographic images for detection of thyroid tumor: An experimental-numerical study

Thermography is a developing and noninvasive medical imaging technique that can be used for diagnosis of body disorders based on temperature deviation from normal body temperature. This research investigates the feasibility of thermography method in conjunction with artificial neural networks (ANNs) for detection of thyroid tumors. For this purpose, first, a 3-D model of the healthy human neck is constructed based on patient-specific computed tomography (CT) images. This model is used for analyzing bio-heat transfer in the human neck. The healthy thyroid gland is considered as a heat source and generates heat according to its temporal temperature. Finite element results verify the thermography potential for detection of thyroid gland location and estimation of its butterfly shape on the neck thermogram. The numerical analysis is carried out on 35 models with varying thermo-physical parameters of the healthy thyroid gland, including heat generation and blood perfusion. The acquired thermograms are used to develop an ANN for correlating the thermo-physical parameters of the gland and temperature profile on the neck surface. In the next stage, dynamic thermal images are captured from 10 healthy and three cancerous human cases. The experimental thermal images are analyzed by the developed ANN and the corresponding thermo-physical parameters are obtained. Results show that the estimated heat generation values for the healthy cases are about 3000 Wm3 while it increases to more than 12 000 Wm3 for the cases with tumors. This significant variation confirms the potential of dynamic thermography in diagnosis of thyroid tumors.

Non-Fourier thermal transport induced structural hierarchy and damage to collagen ultrastructure subjected to laser irradiation

Comprehending the mechanism of thermal transport through biological tissues is an important factor for optimal ablation of cancerous tissues and minimising collateral tissue damage. The present study reports detailed mapping of the rise in internal temperature within the tissue mimics due to NIR (1064 nm) laser irradiation, both for bare mimics and with gold nanostructures infused. Gold nanostructures such as mesoflowers and nanospheres have been synthesised and used as photothermal converters to enhance the temperature rise, resulting in achieving the desired degradation of malignant tissue in targeted region. Thermal history was observed experimentally and simulated considering non-Fourier dual phase lag (DPL) model incorporated Pennes bio-heat transfer equation using COMSOL Multiphysics software. The gross deviation in temperature i.e. rise from the classical Fourier model for bio-heat conduction suggests additional effects of temperature rise on the secondary structures and morphological and physico-chemical changes to the collagen ultrastructures building the tissue mass. The observed thermal denaturation in the collagen fibril morphologies have been explained based on the physico-chemical structure of collagen and its response to thermal radiation. The large shift in frequency of amides A and B is pronounced at a depth of maximum temperature rise compared with other positions in tissue phantom. Observations for change in band of amide I, amide II, and amide III are found to be responsible for damage to collagen ultra-structure. Variation in the concentration of gold nanostructures shows the potentiality of localised hyperthermia treatment subjected to NIR radiation through a proposed free radical mechanism.

Non-Fourier based thermal-mechanical tissue damage prediction for thermal ablation

Prediction of tissue damage under thermal loads plays important role for thermal ablation planning. A new methodology is presented in this paper by combing non-Fourier bio-heat transfer, constitutive elastic mechanics as well as non-rigid motion of dynamics to predict and analyze thermal distribution, thermal-induced mechanical deformation and thermal-mechanical damage of soft tissues under thermal loads. Simulations and comparison analysis demonstrate that the proposed methodology based on the non-Fourier bio-heat transfer can account for the thermal-induced mechanical behaviors of soft tissues and predict tissue thermal damage more accurately than classical Fourier bio-heat transfer based model.

Experimental measurement of microwave ablation heating pattern and comparison to computer simulations

INTRODUCTION

For computational models of microwave ablation (MWA), knowledge of the antenna design is necessary, but the proprietary design of clinical applicators is often unknown. We characterised the specific absorption rate (SAR) during MWA experimentally and compared to a multi-physics simulation.

METHODS

An infrared (IR) camera was used to measure SAR during MWA within a split ex vivo liver model. Perseon Medical's short-tip (ST) or long-tip (LT) MWA antenna were placed on top of a tissue sample (n = 6), and microwave power (15 W) was applied for 6 min, while intermittently interrupting power. Tissue surface temperature was recorded via IR camera (3.3 fps, 320 × 240 resolution). SAR was calculated intermittently based on temperature slope before and after power interruption. Temperature and SAR data were compared to simulation results.

RESULTS

Experimentally measured SAR changed considerably once tissue temperatures exceeded 100 °C, contrary to simulation results. The ablation zone diameters were 1.28 cm and 1.30 ± 0.03 cm (transverse), and 2.10 cm and 2.66 ± -0.22 cm (axial), for simulation and experiment, respectively. The average difference in temperature between the simulation and experiment were 5.6 °C (ST) and 6.2 °C (LT). Dice coefficients for 1000 W/kg SAR iso-contour were 0.74 ± 0.01 (ST) and 0.77 (± 0.03) (LT), suggesting good agreement of SAR contours.

CONCLUSION

We experimentally demonstrated changes in SAR during MWA ablation, which were not present in simulation, suggesting inaccuracies in dielectric properties. The measured SAR may be used in simplified computer simulations to predict tissue temperature when the antenna geometry is unknown.

Prediction of tissue thermal damage

This paper presents a method to characterize tissue thermal damage by taking into account the thermal-mechanical effect of soft tissues for thermal ablation. This method integrates the bio-heating conduction and non-rigid motion dynamics to describe thermal-mechanical behaviors of soft tissues and further extends the traditional tissue damage model to characterize thermal-mechanical damage of soft tissues. Simulations and comparison analysis demonstrate that the proposed method can effectively predict tissue thermal damage and it also provides reliable guidelines for control of the thermal ablation procedure.

Bio-heat transfer simulation of retinal laser irradiation

Retinopathy is a surgical process in which maladies of the human eye are treated by laser irradiation. A two-dimensional numerical model of the human eye geometry has been developed to investigate transient thermal effects due to laser radiation. In particular, the influence of choroidal pigmentation and that of choroidal blood convection-parameterized as a function of choroidal blood perfusion-are investigated in detail. The Pennes bio-heat transfer equation is invoked as the governing equation, and finite volume formulation is employed in the numerical method. For a 500-μm diameter spot size, laser power of 0.2 W, and 100% absorption of laser radiation in the retinal pigmented epithelium (RPE) region, the peak RPE temperature is observed to be 103 °C at 100 ms of the transient simulation of the laser surgical period. Because of the participation of pigmented layer of choroid in laser absorption, peak temperature is reduced to 94 °C after 100 ms of the laser surgery period. The effect of choroidal blood perfusion on retinal cooling is found to be negligible during transient simulation of retinopathy. A truncated three-dimensional model incorporating multiple laser irradiation of spots is also developed to observe the spatial effect of choroidal blood perfusion and choroidal pigmentation. For a circular array of seven uniformly distributed spots of identical diameter and laser power of 0.2 W, transient temperature evolution using simultaneous and sequential mode of laser surgical process is presented with analysis.