First validation of a model-based hepatic percutaneous microwave ablation planning on a clinical dataset

A model-based planning tool, integrated in an imaging system, is envisioned for CT-guided percutaneous microwave ablation. This study aims to evaluate the biophysical model performance, by comparing its prediction retrospectively with the actual ablation ground truth from a clinical dataset in liver. The biophysical model uses a simplified formulation of heat deposition on the applicator and a heat sink related to vasculature to solve the bioheat equation. A performance metric is defined to assess how the planned ablation overlaps the actual ground truth. Results demonstrate superiority of this model prediction compared to manufacturer tabulated data and a significant influence of the vasculature cooling effect. Nevertheless, vasculature shortage due to branches occlusion and applicator misalignment due to registration error between scans affects the thermal prediction. With a more accurate vasculature segmentation, occlusion risk can be estimated, whereas branches can be used as liver landmarks to improve the registration accuracy. Overall, this study emphasizes the benefit of a model-based thermal ablation solution in better planning the ablation procedures. Contrast and registration protocols must be adapted to facilitate its integration into the clinical workflow.

Morphometric characterization and temporal temperature measurements during hepatic microwave ablation in swine

PURPOSE

Heat-induced destruction of cancer cells via microwave ablation (MWA) is emerging as a viable treatment of primary and metastatic liver cancer. Prediction of the impacted zone where cell death occurs, especially in the presence of vasculature, is challenging but may be achieved via biophysical modeling. To advance and characterize thermal MWA for focal cancer treatment, an in vivo method and experimental dataset were created for assessment of biophysical models designed to dynamically predict ablation zone parameters, given the delivery device, power, location, and proximity to vessels.

MATERIALS AND METHODS

MWA zone size, shape, and temperature were characterized and monitored in the absence of perfusion in ex vivo liver and a tissue-mimicking thermochromic phantom (TMTCP) at two power settings. Temperature was monitored over time using implanted thermocouples with their locations defined by CT. TMTCPs were used to identify the location of the ablation zone relative to the probe. In 6 swine, contrast-enhanced CTs were additionally acquired to visualize vasculature and absence of perfusion along with corresponding post-mortem gross pathology.

RESULTS

Bench studies demonstrated average ablation zone sizes of 4.13±1.56cm2 and 8.51±3.92cm2, solidity of 0.96±0.06 and 0.99±0.01, ablations centered 3.75cm and 3.5cm proximal to the probe tip, and temperatures of 50 ºC at 14.5±13.4s and 2.5±2.1s for 40W and 90W ablations, respectively. In vivo imaging showed average volumes of 9.8±4.8cm3 and 33.2±28.4cm3 and 3D solidity of 0.87±0.02 and 0.75±0.15, and gross pathology showed a hemorrhagic halo area of 3.1±1.2cm2 and 9.1±3.0cm2 for 40W and 90W ablations, respectfully. Temperatures reached 50ºC at 19.5±9.2s and 13.0±8.3s for 40W and 90W ablations, respectively.

CONCLUSION

MWA results are challenging to predict and are more variable than manufacturer-provided and bench predictions due to vascular stasis, heat-induced tissue changes, and probe operating conditions. Accurate prediction of MWA zones and temperature in vivo requires comprehensive thermal validation sets.

Finding an effective MRI sequence to visualise the electroporated area in plant-based models by quantitative mapping

Plant-based models can reduce the number of animal studies for electroporation research in medical cancer treatment modalities like irreversible electroporation. Magnetic resonance imaging (MRI) provides volumetric visualisation of electroporated animal or plant tissues; however, contrast behaviour is complex, depending on tissue and sequence parameters. This study numerically analysed contrast between electroporated and non-electroporated tissue at 1.5 T in various MRI sequences (DWI, T1W, T2W, T2*W, PDW, FLAIR) performed 4 h after electroporation in apples (N = 4) and potatoes (N = 8). Sequence parameters (inversion time [TI], echo time [TE], b-value) for optimal contrast and electroporation-mediated changes in T1 and T2 relaxation times and apparent diffusion coefficient (ADC) were determined for potato (N = 4) using quantitative parameter mapping. FLAIR showed the electroporated zone in potatoes with best contrast, whereas no sequence yielded clear visibility in apples. After electroporation, T1 and T2 in potato decreased by 29% ([1245 ± 54 to 886 ± 119] ms) and 12% ([249 ± 17 to 217 ± 12] ms), respectively. ADC increased by 11% ([1303 ± 25 to 1449 ± 28] × 10-6 mm²/s). Optimal contrast was found for TI = 1000 ms, low TE and high b-value. T1 was most sensitive to EP-mediated tissue changes. Future research could use this methodology and findings to obtain high-contrast MR images of electroporated and non-electroporated biological tissues.

Model-based hepatic percutaneous microwaveablation planning. First validation on a clinical dataset

A model-based planning tool, integrated in an imaging system, is envisioned for CT-guided percutaneous microwave ablation. This study aims to evaluate the biophysical model performance, by comparing its prediction retrospectively with the actualablation ground truth from a clinical data set in liver. The biophysical model uses a simplified formulation of heat depositionon the applicator and a heat sink related to vasculature to solve the bioheat equation. A performance metric is defined toassess how the planned ablation overlaps the actual ground truth. Results demonstrate superiority of this model predictioncompared to manufacturer tabulated data and a significant influence of the vasculature cooling effect. Nevertheless, vasculatureshortage due to branches occlusion and applicator misalignment due to registration error between scans affects the thermalprediction. With a more accurate vasculature segmentation, occlusion risk can be estimated, whereas branches can be usedas liver landmarks to improve the registration accuracy. Overall, this study emphasizes the benefit of a model-based thermalablation solution in better planning the ablation procedures. Contrast and registration protocols must be adapted to facilitate itsintegration into the clinical workflow.

Global sensitivity study for irreversible electroporation: Towards treatment planning under uncertainty

BACKGROUND

Electroporation-based cancer treatments are minimally invasive, nonthermal interventional techniques that leverage cell permeabilization to ablate the target tumor. However, the amount of permeabilization is susceptible to the numerous uncertainties during treatment, such as patient-specific variations in the tissue, type of the tumor, and the resolution of imaging equipment. These uncertainties can reduce the extent of ablation in the tissue, thereby affecting the effectiveness of the treatment.

PURPOSE

The aim of this work is to understand the effect of these treatment uncertainties on the treatment outcome for irreversible electroporation (IRE) in the case of colorectal liver metastasis (CRLM). Understanding the nature and extent of these effects can help us identify the influential treatment parameters and build better models for predicting the treatment outcome.

METHODS

This is an in silico study using a static computational model with a custom applicator design, spherical tissue, and tumor geometry. A nonlinear electrical conductivity, dependent on the local electric field, is considered. Morris analysis is used to identify the influential treatment parameters on the treatment outcome. Seven treatment parameters pertaining to the relative tumor location with respect to the applicator, the tumor growth pattern, and the electrical conductivity of tissue are analyzed. The treatment outcome is measured in terms of the relative tumor ablation with respect to the target ablation volume and total ablation volume.

RESULTS

The Morris analysis was performed with 800 model evaluations, sampled from the seven dimensional input parameter space. Electrical properties of the tissue, especially the electrical conductivity of the tumor before ablation, were found to be the most influential parameter for relative tumor ablation and total ablation volume. This parameter was found to be about 4-15 times more influential than the least influential parameter, depending on the tumor size. The tumor border configuration was identified as the least important parameter for treatment effectiveness. The most desired treatment outcome is obtained by a combination of high healthy liver conductivity and low tumor conductivity. This information can be used to tackle worst-case scenarios in treatment planning. Finally, when the safety margins used in the clinical applications are accounted for, the effects of uncertainties in the treatment parameters reduce drastically.

CONCLUSIONS

The results of this work can be used to create an efficient surrogate estimator for uncertainty quantification in the treatment outcome, that can be utilized in optimal real-time treatment planning solutions.

Pulse Parameters and Thresholds for (ir)Reversible Electroporation on Hepatocellular Carcinoma Cells in Vitro

Hepatocellular carcinoma is a leading cause of cancer-related death in many parts of the world. Traditional treatment options are not always effective. During the promising minimally invasive electroporation-based therapies, biological cell membranes are exposed to an external, sufficiently high, pulsed electric field which creates so-called nanopores into the lipid bilayer of the cell membrane. These pores can either be permanent (irreversible electroporation (IRE)), leading to apoptosis, or repairable (reversible electroporation (RE)), with continued cell function. In tumor therapy, RE is used to increase the diffusion of a chemotherapeutic drug during electrochemotherapy. For both IRE and RE, the success of the treatment is dependent on application of the appropriate electric field. Therefore, this study aims to define the pulse parameters and thresholds for IRE and RE on hepatocellular carcinoma (HepG2) cells in-vitro.In a custom-made in-vitro setup, HepG2 cell viability (0, 5, 10, and 15 min), and the peak temperature were measured after electroporation with the different IRE and RE pulsing protocols, to determine the most successful settings for IRE and RE. A CAM/PI flow cytometric assay was performed to confirm cell permeabilization for the RE pulsing protocols with the highest cell viability.The results indicated that an IRE pulsing protocol (70 pulses, 100 µs pulse length, and 100 ms interval) with an electric field strength of 4000 V/cm was needed as threshold for almost complete cell death of HepG2 cells. A RE pulsing protocol (8 pulses, 100 µs pulse length, and 1000 ms interval) with an electric field strength of 1000 V/cm was needed as threshold for viable and permeabilized HepG2 cells. The low peak temperatures (max 30.1°C for IRE, max 23.1°C for RE) within this study indicated that the reduction in HepG2 cell viability was caused by the applied electric field and was not a result of Joule heating.

Plant-based model for the visual evaluation of electroporated area after irreversible electroporation and its comparison to in-vivo animal data

Electroporation (EP) is widely used in medicine, such as cancer treatment, in form of electrochemotherapy or irreversible electroporation (IRE). For EP device testing, living cells or tissue inside a living organism (including animals) are needed. Plant-based models seem to be a promising alternative to substitute animal models in research. The aim of this study is to find a suitable plant-based model for visual evaluation of IRE, and to compare the geometry of electroporated areas with in-vivo animal data.For this purpose, a variety of fruit and vegetables were selected and visually evaluated after 0/1/2/4/6/8/12/16/24 h post-EP. Apple and potato were found to be suitable models as they enabled a visual evaluation of the electroporated area. For these models, the size of the electroporated area was determined after 0/1/2/4/6/8/12/16/24 h. For apples, a well-defined electroporated area was visual within two hours, while in potatoes it reached a plateau after eight hours only. The electroporated area of apple, which showed the fastest visual results was then compared to a retrospectively evaluated swine liver IRE dataset which had been obtained for similar conditions. The electroporated area of the apple and swine liver both showed a spherical geometry of comparable size. For all experiments, the standard protocol for human liver IRE was followed. To conclude, potato and apple were found to be suitable plant-based models for the visual evaluation of electroporated area after irreversible EP, with apple being the best choice for fast visual results. Given the comparable range, the size of the electroporated area of the apple may be promising as a quantitative predictor in animal tissue. Even if plant-based models cannot completely replace animal experiments, they can be used in the early stages of EP device development and testing, decreasing animal experiments to the necessary minimum.

Tuning the Pennes Perfusion Rate to Model Large Vessel Cooling Effects in Hepatic Radiofrequency Ablation

Radio-frequency ablation (RFA) has become a popular method for the minimally invasive treatment of liver cancer. However, the success rate of these treatments depends heavily on the amount of experience the clinician possesses. Mathematical modelling can help mitigate this problem by providing an indication of the treatment outcome. Thermal lesions in RFA are affected by the cooling effect of both fine-scale and large-scale blood vessels. The exact model for large-scale blood vessels is advection-diffusion, i.e. a model capable of producing directional effects, which are known to occur in certain cases. In previous research, in situations where directional effects do not occur, the advection term in the blood vessel model has been typically replaced with the Pennes perfusion term, albeit with a higher-than usual perfusion rate. Whether these values of the perfusion rate appearing in literature are optimal for the particular vessel radii in question, has not been investigated so far. The present work aims to address this issue. An attempt has been made to determine, for values of vessel radius between 0.55 mm and 5 mm, best estimates for the perfusion rate which minimize the error in thermal lesion volumes between the perfusion-based model and the advection-based model. The results for the best estimate of the perfusion rate presented may be used in existing methods for fast estimation of RFA outcomes. Furthermore, the possible improvements to the presented methodology have been highlighted.

Effective models of microwave antennae for ablation treatment planning

The level of detail of typical numerical models of microwave tumor ablations poses a challenge to the development of generic, model based treatment planning tools aiming at real time performance. The present contribution describes a flexible and accurate approximation of the microwave heat absorption that aims at mitigating these issues.

Study of flow effects on temperature-controlled radiofrequency ablation using phantom experiments and forward simulations

PURPOSE

Blood flow is known to add variability to hepatic radiofrequency ablation (RFA) treatment outcomes. However, few studies exist on its impact on temperature-controlled RFA. Hence, we investigate large-scale blood flow effects on temperature-controlled RFA in flow channel experiments and numerical simulations.

METHODS

Ablation zones were induced in tissue-mimicking, thermochromic phantoms with a single flow channel, using an RF generator with temperature-controlled power delivery and a monopolar needle electrode. Channels were generated by molding the phantom around a removable rod. Channel radius and saline flow rate were varied to study the impact of flow on (i) the ablated cross-sectional area, (ii) the delivered generator power, and (iii) the occurrence of directional effects on the thermal lesion. Finite volume simulations reproducing the experimental geometry, flow conditions, and generator power input were conducted and compared to the experimental ablation outcomes.

RESULTS

Vessels of different channel radii r affected the ablation outcome in different ways. For r = 0.275  mm, the ablated area decreased with increasing flow rate while the energy input was hardly affected. For r = 0.9  mm and r = 2.3  mm, the energy input increased toward larger flow rates; for these radii, the ablated area decreased and increased toward larger flow rates, respectively, while still being reduced overall as compared to the reference experiment without flow. Directional effects, that is, local shrinking of the lesion upstream of the needle and an extension thereof downstream, were observed only for the smallest radius. The simulations qualitatively confirmed these observations. As compared to performing the simulations without flow, including flow effects in the simulations reduced the mean absolute error between experimental and simulated ablated areas from 0.23 to 0.12.

CONCLUSION

While the temperature control mechanism did not detect the heat sink effect in the case of the smallest channel radius, it counteracted the heat sink effect in the case of the larger channel radii with an increased energy input; this explains the increase in ablated area toward high flow rates (for r = 2.3  mm). The experiments in a simple phantom setup, thus, contribute to a good understanding of the phenomenon and are suitable for model validation.

Simulation study of the cooling effect of blood vessels and blood coagulation in hepatic radio-frequency ablation

PURPOSE

Computer simulations of hepatic radio-frequency ablation (RFA) were performed to: (i) determine the dependence of the vessel wall heat transfer coefficient on geometrical parameters; (ii) study the conditions required for the occurrence of the directional effect of blood; and (iii) classify blood vessels according to their effect on the thermal lesion while considering blood coagulation. The information thus obtained supports the development of a multi-scale bio-heat model tailored for more accurate prediction of hepatic RFA outcomes in the vicinity of blood vessels.

MATERIALS AND METHODS

The simulation geometry consisted of healthy tissue, tumor tissue, a mono-polar RF-needle, and a single cylindrical blood vessel. The geometrical parameters of interest were the RF-needle active length and those describing blood vessel configuration. A simple, novel method to incorporate the effects of blood coagulation into the simulation was developed and tested.

RESULTS

A closed form expression giving the dependence of the vessel wall heat transfer coefficient on geometrical parameters was obtained. Directional effects on the thermal lesion were found to occur for blood vessel radii between 0.4 mm and 0.5 mm. Below 0.4 mm blood coagulation blocked the flow.

CONCLUSIONS

The closed form expression for the heat transfer coefficient can be used in models of RFA to speed up computation. The conditions on vessel radii required for the occurrence of directional effects on the thermal lesion were determined. These conditions allow the classification of blood vessels. Different approximations to the thermal equation can thus be used for these vessel classes.

Modeling the interference between shear and longitudinal waves under high intensity focused ultrasound propagation in bone

Magnetic resonance-guided high intensity focused ultrasound (MR-HIFU) is a noninvasive thermal technique that enables rapid heating of a specific area in the human body. Its clinical relevance has been proven for the treatments of soft tissue tumors, like uterine fibroids, and for the treatments of solid tumors in bone. In MR-HIFU treatment, MR-thermometry is used to monitor the temperature evolution in soft tissue. However, this technique is currently unavailable for bone tissue. Computer models can play a key role in the accurate prediction and monitoring of temperature. Here, we present a computer ray tracing model that calculates the heat production density in the focal region. This model accounts for both the propagation of shear waves and the interference between longitudinal and shear waves. The model was first compared with a finite element approach which solves the Helmholtz equation in soft tissue and the frequency-domain wave equation in bone. To obtain the temperature evolution in the focal region, the heat equation was solved using the heat production density generated by the raytracer as a heat source. Then, we investigated the role of the interaction between shear and longitudinal waves in terms of dissipated power and temperature output. The results of our model were in agreement with the results obtained by solving the Helmholtz equation and the frequency-domain wave equation, both in soft tissue and bone. Our results suggest that it is imperative to include both shear waves and their interference with longitudinal waves in the model when simulating high intensity focused ultrasound propagation in solids. In fact, when modeling HIFU treatments, omitting the interference between shear and longitudinal waves leads to an over-estimation of the temperature increase in the tissues.