Numerical mesoscale tissue model of electrochemotherapy in liver based on histological findings

Electrochemotherapy (ECT) and irreversible electroporation (IRE) are being investigated for treatment of hepatic tumours. The liver is a highly heterogeneous organ, permeated with a network of macro- and microvasculature, biliary tracts and connective tissue. The success of ECT and IRE depends on sufficient electric field established in whole target tissue; therefore, tissue heterogeneity may affect the treatment outcome. In this study, we investigate electroporation in the liver using a numerical mesoscale tissue model. We numerically reconstructed four ECT experiments in healthy porcine liver and computed the electric field distribution using our treatment planning framework. We compared the computed results with histopathological changes identified on microscopic images after treatment. The mean electric field threshold that best fitted the zone of coagulation necrosis was 1225 V/cm, while the mean threshold that best fitted the zone of partially damaged liver parenchyma attributed to IRE was 805 V/cm. We evaluated how the liver macro- and microstructures affect the electric field distribution. Our results show that the liver microstructure does not significantly affect the electric field distribution on the level needed for treatment planning. However, major hepatic vessels and portal spaces significantly affect the electric field distribution, and should be considered when planning treatments.