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.

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.

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 mm2/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.

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.

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.

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 modeling 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. This 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.