- 1College of Medicine and Biological Information Engineering, Northeastern University, Shenyang, China
- 2Department of Cardiology, General Hospital of Northern Theater Command, Shenyang, China
- 3National Key Laboratory of Frigid Zone Cardiovascular Diseases, General Hospital of Northern Theater Command, Shenyang, China
- 4Graduate School of China Medical University, China Medical University, Shenyang, China
Background: Pulsed-field ablation (PFA) technique is a nonthermal ablation technique. No study has yet evaluated the effect of the positional relationship between the ablation electrode (AE) and the coronary metal stent (CMS) on the electric field distribution and temperature distribution in epicardial ablation. Our study aimed to evaluate the effect of the CMS on the electric field as well as the temperature distribution in different models.
Methods: Multi-angle modeling of the CMS and AE was performed. The PFA ablation region was evaluated with a field strength contour of 1,000 V/cm, which was used to assess the validity of the two-dimensional (2D) model simulation data as well as the distribution of the multi-angle electric field and temperature in the three-dimensional (3D) model.
Results: The presence of the CMS had little effect on the width of the ablation area (0.2 mm). In the 3D model, the temperature of the ablation area was highest when the angle between the AE and the CMS was in the 90° position (43.4°C, 41.3°C); a change in the distance between the AE and the CMS affected the temperature of the ablation area (maximum 2.1°C) and the width of the ablation (maximum 0.32 mm).
Conclusion: The presence of the CMS distorts the distribution of the electric field, but does not produce a change in the extent of the ablation damage, nor does it bring thermal damage to the ablation region. Different simulation models give similar results in PFA calculations, and this study effectively reduces the complexity of modeling simulation.
1 Introduction
Atrial fibrillation (AF), one of the most common clinical arrhythmias, and pulsed electric field ablation (PFA) has recently been proposed by researchers for the treatment of AF (1–6). PFA was used in the treatment of cardiac arrhythmias as early as 1980, but due to the limitations of the technology at that time, it could cause other complications. With the development of this ablation device, this ablation technique has rekindled the interest of researchers (7–10). The researchers found that compared with radiofrequency ablation and cryoablation, PFA has a better safety profile, such as reducing the chance of complications such as phrenic nerve injury and pulmonary vein stenosis during endocardial ablation (11–13). Early animal studies demonstrated that arrhythmias can be induced by stimulation of the cardiac nervous system. These nervous systems are interconnected with ganglionic plexus of the epicardium, which are usually embedded in the fatty layer of the epicardium (14–16).One study demonstrated that PFA for ablation of epicardial ganglionic plexus could effectively treat AF, and also found that the presence of ganglionic plexus distorted the normal electric field in the computational model, which may be due to the higher electrical conductivity of the ganglionic plexus compared to its surrounding tissues (17, 18). Follow-up researchers hypothesized that metal stents inside coronary arteries might have a more pronounced distorting effect on the electric field and cause a local temperature increase, and Hogens found experimentally and computationally that during electrode ablation of tumors, a distortion of the electric field distribution was observed when the metal stent was close to the ablation electrode. In another study, it was found that although the heat of the metal stent itself does not increase dramatically, it can be observed that the temperature near the ablation electrode increases due to the presence of the metal stent (19–21). In a later study, ANN et al. found that the presence of the metal stent caused a distorted distribution of the electric field of the ablation electrodes and increased the temperature of the tissues between the ablation electrodes, but this temperature change did not cause thermal damage to the nearby tissues by calculating the two-dimensional computational model. Later, it was found that in the 3D model calculation, the difference between the simplified 3D model and the real complex 3D model ablation area and the depth of the ablation area is negligible, which simplifies the complexity of 3D model modeling (22, 23). Given that it has been shown that the presence of metal brackets induces a distorted distribution of the normal electric field and that the size of the ablation region between the simplified 3D model and the real complex 3D model is negligible, it is necessary to analyze and calculate both the 3D model and the 2D model in the computational domain. To this end, we conduct a computer modeling study to evaluate the effect of the presence of the metal bracket in the 3D model and the 2D model and to analyze and compare the gap between the electric field distributions as well as the temperature distributions of the two models, to reduce the complexity of the modeling simulation calculations.
2 Methods
2.1 Modeling
To facilitate the comparison of the ablation effects of the two models, we only consider the ablation device and the surrounding area of interest when building the model, and this modeling method has been proven to be effective in previous studies (23–25). According to previous studies and the need for model simplification, we simplified the different tissues in the ablation area, modeled the ablation region of interest in layers (saline, fat, myocardium, and blood), and placed the ablation device above the fat layer according to the real structure of the human organ, in which the ablation device dimensions were set to realistically simulate the size of the ablation area, and the dimensions were identical to the real size of the ablation device. The size of the ablation device is the same as the real size of the ablation device, and the model is established as shown in Figure 1 (20, 23, 26). The PFA ablation position is shown in Figure 1, with the left coronary artery usually located directly below the ablation electrode. The ablation electrode has a diameter of 3.98 mm, a length of 2.54 mm, a 0.76 mm flushing hole in the center, and 0.5 mm saline on the upper side of the epicardial fat, which is used to act as a virtual electrode, and the ablation electrode and the fat are completely close to each other, which ensures that the ablation energy can be effectively transmitted to the ablation area (25–27). The computer simulation model is shown in Figure 2. The fat thickness in the model was 4.3 mm, the myocardial thickness was 2.7 mm, and the blood thickness and width models were set to X = 80 and Y = 40, respectively, concerning previous studies (23, 28).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1445424/fcvm-11-1445424-HTML-r1/image_m/fcvm-11-1445424-g001.jpg)
Figure 1. Schematic diagram of left atrial decortication (26, 29). Figure (A) Anatomical picture; Figure (B) Relation between the AE and the left coronary artery (LCx).
2.2 Boundary conditions and pulse parameters
Boundary conditions were set for each of the two simulation models, and 1,000 V was set for the PFA catheter ablation metal electrode with a pulse duration of 100 us and a pulse interval of 100 us (23, 30). The voltage pulse parameters are shown in Figure 3. For the metal electrode spacing insulator set 0 V, and set the voltage to 0 V for all the boundaries of the model, to show that the simulation process is not subject to external interference, and the electric field has the metal electrode spreading outward (25, 31).
To evaluate the effect of ablating different metal stents on temperature, we set thermal boundary conditions, for which thermal convection coefficients were applied at the air-electrode interface (20 W/m2 K, ambient temperature 21°C) and at the myocardium-blood interface to simulate blood circulation (1,417 W/m2 K, 37°C blood temperature) under high blood flow (24.4 cm/s) conditions (32–34). This is because blood flow occurs inside the coronary arteries due to the contraction of the heart chambers, so in the simulation, we took into account the characteristics of blood flow, and the cooling effect of tissues such as the heart chambers on the ablation area, we set the myocardial temperature to 37°C in the model. In the absence of a stent, we simplified the situation by assuming a flow rate of 0.5 m/s, half of the peak reported in, corresponding to a convection coefficient of 63.19 /m2 K calculated in (34, 35).
2.3 Material properties
The conductivity of different biological tissues varies from person to person, and appropriate values of biological properties can improve the calculation accuracy of computer simulation models. After long-term exploration by researchers, the conductivity of biological tissues shows an S-shaped distribution with the growth of voltage. This may be because the electric field applied by PFA changes the permeability of biological tissues, resulting in the gap in the cell membrane of biological tissues becoming larger when receiving the stimulation of the electric field, and the permeability of the electric current is enhanced resulting. In the paper, we also considered the change in conductivity of saline, blood, and other tissues with temperature so that the conductivity increased by 2% with temperature (30, 36–38). In our experiments, we used the sigmoid function to simulate the conductivity change of biological tissues, Eq. The various parameters of electrical and thermal properties needed for the simulation are shown in Table 1 (36, 39, 40).
where 0 and 1 are the conductivity before and after electroporation, respectively, and the values of 0 and 1 are related to the presence and size of pores in the cell membrane. Before the appearance of PFA, i.e., when the pores have not yet formed, the current flows only through the extracellular material, which is the same as when the tissue is subjected to low-frequency electrical excitation, and the cell membrane acts as an electrical insulator. Under practical conditions, it was maintained without significant change between 1 and 10 kHz. At higher pulse frequencies, the conductivity on the myocardium decreases due to pore formation resulting in current flowing not only extracellularly, but also through the cytoplasm through the pores, but the regional effect of PFA ablation is negligible in other pulse frequency conditions (23, 37).
2.4 Control equations
To simulate the PFA ablation phenomenon in a real scenario, similar to the previous research methods, we chose to build a finite element simulation model in the simulation software COMSOL, and the geometry and finite element mesh used in the study are shown above, in which the catheter electrodes are located in the saline layer, which is used to simulate the real environment of the human body ablation. The various physical parameters required for the control equations are shown in Table 2. In the computational model, we added the biological heat transfer module, so the model belongs to the electro-thermal coupling problem, and in the current module, we introduced Laplace's system of equations to simulate the simulated electric field (30, 41, 42).
where is the conductivity, the exact value of which is given above, V is the voltage, E is the electric field strength, and J is the current density.
For considering the heat generated during the PFA ablation process, we used the bioheat equation to investigate any possible side effects due to the heat source during the ablation process.
Where is the density in (kg/m3); c is the specific heat in (J/kg·K); T is the temperature in (°C); t is the time in (s); k is the thermal conductivity of the tissue in (W/m·K); Q is the heat source generated by the electric field, which is related to the field strength and conductivity in (W/m3); Qp is the heat loss due to blood flow in (W/m3); and Qmet is the heat loss due to the biological metabolism in (W/m3).
2.5 Analysis of results
PFA ablation was simulated using finite element simulation to assess whether the presence of a metal stent during epicardial ablation affects the electric field distribution during ablation and also to simulate the change in temperature of the ablation region during ablation. In the paper, we simulated different scenarios: the position and orientation between the arterial stent and the ablation electrode, the distance between the ablation electrode and the metal stent, and the position of the different saline layers about the ablation electrode. We utilized the 1,000 V/cm contour to assess the maximum width and depth of the PFA ablation area, due to a study that recently reported an electric field threshold of 1,000 V/cm for PFA-induced irreversible myocardial damage (43). Temperature distributions were calculated to assess possible adverse thermal side effects, especially in the vicinity of ablation electrodes and metal stents.
3 Results
3.1 Electric field distribution
The distribution of the electric field in the ablation region is shown in Figure 4 for different conditions at four positions. Our prerequisite was to perform ablation simulation experiments under the condition of ten pulses utilizing the ablation catheter (pulse duration of 0.1 ms and pulse interval of 0.1 ms). In the ablation experiments, we set the coronary arteries to contain a metal stent inside, and at the same time, the distance between the metal stent and the ablation catheter was 1 mm. It can be seen from Figure 4 that the presence of the metal stent inside the coronary arteries affects the normal distribution state of the electric field. This finding is similar to the results of previous studies, while the angular relationship between the coronary artery and the ablation catheter affects the distribution of the electric field around the artery, and it can be concluded from Figure 4 that when the direct angle between the ablation catheter and the coronary artery is 0°, the distribution of the electric field will be higher on both sides of the upper and lower sides of the coronary artery. At the same time, this phenomenon also occurs between different angles of the coronary artery and the ablation catheter, and this finding can also be derived from the electric field distributions of the simulated regions at different positional angles in Figure 4.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1445424/fcvm-11-1445424-HTML-r1/image_m/fcvm-11-1445424-g004.jpg)
Figure 4. Schematic diagram of the electric field distribution in the ablation area of the catheter and stent at different positions of 1 mm, with an ablation voltage of 1,000 V, and the number of ablation pulses of 10, with an interval of 0.1 ms between each pulse. Figures (A,B) show 2D simulations of epicardial pulsed fields; Figures (C–F) show 3D simulations of epicardial pulsed fields.
Due to the presence of the coronary arteries, an abnormal distortion of the electric field distribution occurs during PFA ablation, but it can be concluded from Figure 4 that this distorted electric field exists only inside the fat layer and has no effect on the cardiomyocytes. In Figure 4, we labeled the ablation regions at different angles, and inside the green wireframe is the region where the electric field strength is greater than 1,000/cm, and it is derived from the literature that this region may cause damage to cardiomyocytes. We measured the width and height of different angles to the ablation catheter in the axial direction of the coronary arteries, and we found that the width of the ablation changes as the angle between the coronary arteries and the metal stent changes, but these ablation areas do not penetrate the fat layer. The width of the ablation area was 11.47 mm in Figure 4A, 11.67 mm in Figure 4B, 11.29 mm in Figure 4C, 12.74 mm in Figure 4D, 18.57 mm in Figure 4E, and 25.37 mm in Figure 4F.
As shown in Figure 4, in the ablation cross-section, we can see that the region with electric field strength greater than 1,000/cm is located in the fat layer, which does not affect the myocardial layer cells, and we find that in the axial cross-section, the region with stronger electric field distribution is the position of the left and right sides of the coronary artery, and the electric field distribution in this region is distorted, which suggests that the presence of the coronary artery and the metal stent affects the distribution of the normal electric field, but this distortion of the electric field is confined to the fat layer. As shown in Figures 4C–F, for the axial position of the coronary artery, the electric field generated by the ablation electrodes at different angles showed an angular difference with the coronary artery, and the distortion of the electric field by the coronary artery still existed, and the location and direction of the electric field distortion and the angular direction of the ablation electrodes had an effect.
From Figure 4, it can be obtained that the AE generates the ablation electric field only in the fat layer, and there is no damage to the myocardium under the fat layer. Figure 5 shows the electric field distribution in the axial direction of the CMS at different angles between the ablation catheter and the coronary artery at a distance of 15 mm, and it can be obtained from Figure 5 that the distorting influence of the coronary artery on the distribution of the electric field still exists, and the width of the ablation region is 11.47 mm in Figure 5A, 11.45 mm in Figure 5B, 10.97 mm in Figure 5C, 12.67 mm in Figure 5D, 18.38 mm in Figure 4E, and 25.27 mm in Figure 4F.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1445424/fcvm-11-1445424-HTML-r1/image_m/fcvm-11-1445424-g005.jpg)
Figure 5. Schematic diagram of the electric field distribution in the ablation area of the catheter and stent at different positions of 1.5 mm, with an ablation voltage of 1,000 V, and the number of ablation pulses of 10, with an interval of 0.1 ms between each pulse. Figures (A,B) show 2D simulations of epicardial pulsed fields; Figures (C–F) show 3D simulations of epicardial pulsed fields.
Similar to Figure 4, the focus of electric field distribution is still on the left and right sides of the coronary arteries, and the distribution of electric field distortion is correlated with the position of the ablation electrode, which is similar to that shown in Figure 4. Analysis of Figures 4, 5 showed that the distance between the metal stent and the ablation electrode had a weak effect on the width of the ablation, which decreased when the distance between the two increased, and we also found that the change in the position of the metal stent did not affect the depth of ablation, which was limited to the fat layer in our study. The intensity of the electric field inside the coronary arteries was much less than the threshold for myocardial injury (1,000 v/cm), which suggests that the electric field did not produce a threshold injury inside the coronary arteries.
3.2 Temperature distribution
The temperature distribution of the ablation region was another focus of our study, and we performed temperature measurements for different positions of the ablation region to obtain the maximum temperature region at different positions. Figure 6 shows the maximum temperature distribution in the ablation region when the ablation catheter was 10 mm from the stromal artery. From Figure 6, it can be obtained that the presence of the metal stent and the coronary artery increases the maximum temperature in the ablation region, while the angle between the ablation electrode and the metal stent effects the temperature in the ablation region because of the presence of the ablation electrode and the metal stent. In Figure 6A the maximum temperature of the ablation region is 38.5 °C, in Figure 6B the maximum temperature of the ablation region is 42.2 °C, in Figure 6C the maximum temperature of the ablation region is 42.4 °C, in Figure 6D the maximum temperature of the ablation region is 41.3 °C, in Figure 6E the maximum temperature of the ablation region is 42.9 °C and in Figure 6F the maximum temperature of the ablation region is 43.4 °C. As shown in Figure 6, we also obtain that for different computer simulation models, the maximum temperature in the ablation region is almost unchanged (42.2°C for the 2D model and 42.4°C for the 3D model), and the temperature difference between the two is only 0.2°C. In our study, we found that the maximum temperature in the ablation area was 43.4°C when the angle between the metal stent and the ablation electrode was 90°. From this, we can conclude that the different positions of the ablation catheter do not have much effect on the temperature inside the ablation area, and the small temperature change will not cause discomfort to the clinical patients due to the change of temperature, it also provides meaningful simulation data to prove that pulsed-field ablation of atrial fibrillation can be performed for post-stenting patients.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1445424/fcvm-11-1445424-HTML-r1/image_m/fcvm-11-1445424-g006.jpg)
Figure 6. Schematic diagram of the electric field distribution in the ablation area of the catheter and stent at different positions of 1 mm, with an ablation voltage of 1,000 V, and the number of ablation pulses of 10, with an interval of 0.1 ms between each pulse. Figures (A,B) show 2D simulations of epicardial pulsed fields; Figures (C–F) show 3D simulations of epicardial pulsed fields.
To study the temperature variation of the ablation region in depth, we made simulations for different distances between the coronary artery and the ablation catheter, which is because this situation may exist in different patients, and the distributions of the temperature when the distance between the coronary artery and the ablation catheter was 1.5 mm are shown in Figure 7. The maximum temperature in the ablation region in Figure 7A is 38.5°C; the maximum temperature in the ablation region in Figure 7B is 41.2°C; the maximum temperature in the ablation region in Figure 7C is 41.2°C; the maximum temperature in the ablation region in Figure 7D is 40.8°C; the maximum temperature in the ablation region in Figure 7E is 41.0°C; and the maximum temperature in the ablation region in Figure 7F is 41.3°C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1445424/fcvm-11-1445424-HTML-r1/image_m/fcvm-11-1445424-g007.jpg)
Figure 7. Schematic diagram of the electric field distribution in the ablation area of the catheter and stent at different positions of 1 mm, with an ablation voltage of 1,000 V, and the number of ablation pulses of 10, with an interval of 0.1 ms between each pulse. Figures (A,B) show 2D simulations of epicardial pulsed fields; Figures (C–F) show 3D simulations of epicardial pulsed fields.
Figure 7 shows that the maximum temperature in the ablation region is 41.3°C and the minimum temperature in the ablation region is 38.5°C. We also found that the temperature in the ablation region changes when the angular position of the ablation electrode and the metal bracket changes, and the maximum temperature in the ablation region is 40.8°C when the angle between the two is 30° in this study. In this study, the maximum temperature in the ablation region was 40.8°C at an angle of 30°. The maximum temperature in the ablation region was 41.3°C at an angle of 90°. The change of angle has less effect on the temperature of the ablation region (maximum temperature difference is 0.5°C), and this result is consistent with that shown in Figure 6, in which the maximum temperature difference in the ablation region is 1.2°C.
4 Discussion
4.1 Findings
At present, the safety and efficacy of PFA in the treatment of AF are still worth exploring, and a large number of studies have proved that endocardial ablation can effectively reduce the chance of damaging the tissues of the esophagus, but there are fewer studies on epicardial ablation regarding tissue damage (44, 45). Recently, it has been demonstrated that due to the presence of metal stents inside the coronary arteries, the ablation process in organisms can cause an increase in temperature in the local area, which may be due to the conductive nature of the metal (21, 46). It has also been demonstrated that the presence of coronary arteries as well as metal stents affects the normal electric field and temperature distribution acted by the ablation catheter (17, 23, 27). However, to date, we are not aware of any study that has evaluated the effect of metal stents inside the coronary arteries on the electric field and temperature during epicardial ablation from different perspectives, and to the best of our knowledge, our study is the first to utilize computational modeling to evaluate the effect of the coronary arteries and metal stents on the distribution of the electric field and the distribution of the temperature during the epicardial ablation process from different perspectives. In our study, three-dimensional modeling was performed using computers to simulate the spatial position of the ablation catheter and coronary arteries. The main results of this study are as follows:
1. The presence of coronary arteries with metallic stents distorts the normal electric field distribution in the ablation region, and this again appears to be consistent in the three-dimensional and two-dimensional models.
2. There is an effect of the angle of the ablation electrode to the metal bracket on the width of the ablation area, but there is no change in the position of the cold spots of the electric field present on both sides of the metal bracket (left and right sides of the metal bracket).
3. Comparing different calculation models, the ablation area of the three-dimensional model is wider than the ablation area of the two-dimensional model, with a difference of 0.38 mm in the ablation area width when the distance between the ablation electrode and the metal bracket is 10 mm, and a difference of 0.46 mm in the ablation area width when the distance is 10 mm.
4. There is a temperature increase between the ablation electrode and the metal support, and the highest point of temperature distribution is near the ablation electrode, with the highest temperature in the ablation area differing by 1.2℃ when the distance between the ablation electrode and the metal support is 10 mm, and the highest temperature in the ablation area differing by 0.1℃ when the distance is 15 mm.
5. Different angles between the ablation catheter and the coronary artery affect the distribution of temperature in the ablation area, with the highest temperature in the ablation area decreasing and then increasing as the angle changes, independent of the distance between the two.
6. The presence of the metal support has a distorting effect on the electric and temperature fields, and this distortion is present in the fat layer region.
We determined the validity of the three-dimensional modeling by comparing the data from the two-dimensional computer model with the three-dimensional computer model. In our computational model, the ablation region of the 2D computational model was wider than that of the 3D computational model (approximately 0.4 mm), which may be attributed to the small difference in organ thickness and model meshing that we considered when modeling the 3D model. Our results are similar to those of previous studies in that the ablation area was mainly distributed near the ablation electrodes, a phenomenon that does not change due to the distance between the ablation electrodes and the coronary arteries (see Figures 4, 5).
We found that the width of the ablation area changes with the distance between the ablation electrode and the coronary artery, and the width of the ablation area increases by 0.32 mm when the distance between the ablation electrode and the coronary artery is increased by 0.5 mm, and this phenomenon does not change due to a change in the angle between the two (see Figures 4, 5). A comparison of Figures 4, 5 reveals that the distance between the coronary artery and the myocardial layer does not result in damage to the myocardial layer from the electric field (0.5 mm decrease in the distance between the coronary artery and the myocardial layer).
For the first time, we assessed the temperature within the ablation region using a three-dimensional model, and in our model, the distance between the ablation catheter and the coronary artery effected the temperature in the ablation region, with the temperature in the ablation region decreasing by 2.1℃ (in terms of the maximum temperature) when the distance between the two was increased, but these temperature changes did not cause damage to other tissues. The temperature in the ablation region is shown in Figures 6, 7 to be the smallest when the angle between the coronary artery and the ablation catheter is 30°, which may be becausethe cross-section we took did not contain an ablation electrode, and the temperature may increase in other ablation regions, but this increase in temperature does not exceed the temperatures shown in Figures 6C, 7C.
Finally, there is now a large body of research demonstrating the safety and efficacy of proving PFA in endocardial treatment of AF, including the fact that it can be effective in reducing the chance of complications such as phrenic nerve injury. In this study of ours, ablation treatment of AF from the epicardium using PFA means that the ablation catheter is closer to the metal stent, which would make the ablation energy to the target area of ablation more affected by the metal stent, and catheter ablation near the metal stent was safe in our simulations. This may mean that the ablation effect is less likely to be affected by the metal stent during endocardial ablation and validates the safety and reliability of endocardial PFA ablation.
4.2 Limitations
In our 3D computational model, the cross-section for measuring electric field and temperature was the radial cross-section of the coronary artery and the metal stent, which was to facilitate the comparison of the electric field distribution and temperature distribution on the same cross-section in different models, and we did not compare the data on other cross-section, which might result in the incomplete acquisition of the ablation data in the model, resulting in the under-analysis of the pulsed ablation characteristics. Nevertheless, we analyzed the electric field and temperature distributions at different angles between the ablation electrode and the coronary artery at the same cross-section in different models, and we also analyzed the two extreme conditions during epicardial ablation, namely, the extreme conditions of 0° and 90° angle between the ablation electrode and the coronary artery, and therefore we believe that our calculations have relevant clinical significance from a physical point of view.
5 Conclusion
This paper compares the simulation data of the 3D finite domain model and 2D finite domain model. 1. The consistency of the simulation data between the 2D finite domain model and 3D finite domain model is obtained, which provides data reference for the simplified model establishment afterwards. 2. The width of the ablation area is investigated by using 3D computational model, and the positional relationship between the coronary artery and the ablation catheter is found to effect the width of the ablation area. 3. The temperature of the ablation region was investigated using a three-dimensional computational model, and the highest temperature in the ablation region was obtained at the extreme condition of a 90° angle between the ablation electrode and the coronary artery. This may be due to the presence of the metal stent in the ablation area causing a normal distribution of the electric field, and this distortion of the electric field raises the temperature in the ablation area (less than 44°C, with the maximum temperature between the ablation electrode and the coronary artery). However, this increase in temperature does not produce thermal damage to the corresponding tissues.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
ZW: Writing – original draft, Conceptualization, Software, Validation. ML: Conceptualization, Investigation, Writing – review & editing. JS: Methodology, Software, Writing – review & editing. JZ: Formal Analysis, Investigation, Writing – review & editing. YL: Conceptualization, Investigation, Writing – review & editing. LX: Methodology, Software, Writing – review & editing. YH: Methodology, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by The Open Project of the Central Logistics Support Department grant number BLB20J012; Science and Technology Program of Liaoning Provincial grant number 2022JH2/101500017 and Science and Technology Program of Liaoning Provincial grant number 2020JH1/10300002.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Turagam MK, Neuzil P, Schmidt B, Reichlin T, Neven K, Metzner A, et al. Clinical outcomes by sex after pulsed field ablation of atrial fibrillation. JAMA Cardiol. (2023) 8(12):1142–51. doi: 10.1001/jamacardio.2023.3752
2. Zhang W, Chen Y, Huang Q, Wang JG. Rationale and design of a randomized controlled trial on intensive management of blood PRESSure and cholesterol in elderly Chinese with hypertension and atrial FibrillatION (IMPRESSION). Cardiol Discov. (2021) 1(3):173–8. doi: 10.1097/CD9.0000000000000026
3. Martin RCG II, McFarland K, Ellis S, Velanovich V. Irreversible electroporation therapy in the management of locally advanced pancreatic adenocarcinoma. J Am Coll Surg. (2012) 215(3):361–9. doi: 10.1016/j.jamcollsurg.2012.05.021
4. Garcia PA, Rossmeisl JH Jr, Ellis TL, Davalos RV. Nonthermal Irreversible Electroporation as a Focal Ablation Treatment for Brain Cancer[M]//Tumors of the Central Nervous System, Volume 12: Molecular Mechanisms, Children’s Cancer, Treatments, and Radiosurgery. Dordrecht: Springer Netherlands (2013). 171–82.
5. Reddy VY, Gerstenfeld EP, Natale A, Whang W, Cuoco FA, Patel C, et al. Pulsed field or conventional thermal ablation for paroxysmal atrial fibrillation. N Engl J Med. (2023) 389(18):1660–71. doi: 10.1056/NEJMoa2307291
6. Verma A, Haines DE, Boersma LV, Sood N, Natale A, Marchlinski FE, et al. Pulsed field ablation for the treatment of atrial fibrillation: PULSED AF pivotal trial. Circulation. (2023) 147(19):1422–32. doi: 10.1161/CIRCULATIONAHA.123.063988
7. Fisher JD, Brodman R, Kim SG, Matos JA, Brodman LE, Wallerson D, et al. Attempted nonsurgical electrical ablation of accessory pathways via the coronary sinus in the Wolff-Parkinson-white syndrome. J Am Coll Cardiol. (1984) 4(4):685–94. doi: 10.1016/S0735-1097(84)80394-0
8. Morady F, Scheinman MM, Kou WH, Griffin JC, Dick M, Herre J, et al. Long-term results of catheter ablation of a posteroseptal accessory atrioventricular connection in 48 patients. Circulation. (1989) 79(6):1160–70. doi: 10.1161/01.CIR.79.6.1160
9. Bardy GH, Coltorti F, Stewart RB, Greene HL, Ivey TD. Catheter-mediated electrical ablation: the relation between current and pulse width on voltage breakdown and shock-wave generation. Circ Res. (1988) 63(2):409–14. doi: 10.1161/01.RES.63.2.409
10. Bertacchini C, Margotti PM, Bergamini E, Lodi A, Ronchetti M, Cadossi R. Design of an irreversible electroporation system for clinical use. Technol Cancer Res Treat. (2007) 6(4):313–20. doi: 10.1177/153303460700600408
11. Parikh V, Kowalski M. Comparison of phrenic nerve injury during atrial fibrillation ablation between different modalities, pathophysiology and management. J Atr Fibrillation. (2015) 8(4):1314. doi: 10.4022/jafib.1314
12. Saad EB, Rossillo A, Saad CP, Martin DO, Bhargava M, Erciyes D, et al. Pulmonary vein stenosis after radiofrequency ablation of atrial fibrillation: functional characterization, evolution, and influence of the ablation strategy. Circulation. (2003) 108(25):3102–7. doi: 10.1161/01.CIR.0000104569.96907.7F
13. Cochet H, Nakatani Y, Sridi-Cheniti S, Cheniti G, Ramirez FD, Nakashima T, et al. Pulsed field ablation selectively spares the oesophagus during pulmonary vein isolation for atrial fibrillation. Europace. (2021) 23(9):1391–9. doi: 10.1093/europace/euab090
14. Armour JA, Hageman GR, Randall WC. Arrhythmias induced by local cardiac nerve stimulation. Am J Physiol Legacy Cont. (1972) 223(5):1068–75. doi: 10.1152/ajplegacy.1972.223.5.1068
15. Armour JA, Murphy DA, Yuan BX, Donald SM, Hopkins DA. Gross and microscopic anatomy of the human intrinsic cardiac nervous system. Anat Rec. (1997) 247(2):289–98. doi: 10.1002/(SICI)1097-0185(199702)247:2%3C289::AID-AR15%3E3.0.CO;2-L
16. Scherlag BJ, Yamanashi W, Patel U, Lazzara R, Jackman WM. Autonomically induced conversion of pulmonary vein focal firing into atrial fibrillation. J Am Coll Cardiol. (2005) 45(11):1878–86. doi: 10.1016/j.jacc.2005.01.057
17. Madhavan M, Venkatachalam KL, Swale MJ, Swale MJ, Desimone CV, Gard JJ, et al. Novel percutaneous epicardial autonomic modulation in the canine for atrial fibrillation: results of an efficacy and safety study. Pacing Clin Electrophysiol. (2016) 39(5):407–17. doi: 10.1111/pace.12824
18. Verma A, Asivatham SJ, Deneke T, Castellvi Q, Neal RE. Primer on pulsed electrical field ablation: understanding the benefits and limitations. Circ Arrhythm Electrophysiol. (2021) 14(9):e010086. doi: 10.1161/CIRCEP.121.010086
19. Månsson C, Nilsson A, Karlson BM. Severe complications with irreversible electroporation of the pancreas in the presence of a metallic stent: a warning of a procedure that never should be performed. Acta Radiol Short Rep. (2014) 3(11):2047981614556409. doi: 10.1177/2047981614556409
20. Hogenes AM, Slump CH, Te Riet OG, Scholten GA, Meijerink MR, Fütterer JJ, et al. Effect of irreversible electroporation parameters and the presence of a metal stent on the electric field line pattern. Sci Rep. (2020) 10(1):13517. doi: 10.1038/s41598-020-70308-3
21. Scheffer HJ, Vogel JA, Van Den Bos W, Neal RE, van Lienden KP, Besselink MGH, et al. The influence of a metal stent on the distribution of thermal energy during irreversible electroporation. PLoS One. (2016) 11(2):e0148457. doi: 10.1371/journal.pone.0148457
22. González-Suárez A, O’Brien B, O’Halloran M, Elahi A. Pulsed electric field ablation of epicardial autonomic ganglia: computer analysis of monopolar electric field across the tissues involved. Bioengineering. (2022) 9(12):731. doi: 10.3390/bioengineering9120731
23. González-Suárez A, Irastorza RM, Deane S, O’Brien B, O’Halloran M, Elahi A. Full torso and limited-domain computer models for epicardial pulsed electric field ablation. Comput Methods Programs Biomed. (2022) 221:106886. doi: 10.1016/j.cmpb.2022.106886
24. Russo AM. Pulsed field ablation: is it better than conventional thermal ablation for treatment of atrial fibrillation? Circulation. (2023) 147(19):1433–5. doi: 10.1161/CIRCULATIONAHA.123.064329
25. Caluori G, Odehnalova E, Jadczyk T, Pesl M, Pavlova I, Valikova L, et al. AC Pulsed field ablation is feasible and safe in atrial and ventricular settings: a proof-of-concept chronic animal study. Front Bioeng Biotechnol. (2020) 8:552357. doi: 10.3389/fbioe.2020.552357
26. Padmanabhan D, Naksuk N, Killu AK, Kapa S, Witt C, Sugrue A, et al. Electroporation of epicardial autonomic ganglia: safety and efficacy in medium-term canine models. J Cardiovasc Electrophysiol. (2019) 30(4):607–15. doi: 10.1111/jce.13860
27. Sánchez-Quintana D, López-Mínguez JR, Macías Y, Cabrera JA, Saremi F. Left atrial anatomy relevant to catheter ablation. Cardiol Res Pract. (2014) 2014. doi: 10.1155/2014/289720
28. González-Suárez A, Pérez JJ, O’Brien B, Elahi A. In silico modelling to assess the electrical and thermal disturbance provoked by a metal intracoronary stent during epicardial pulsed electric field ablation. J Cardiovasc Dev Dis. (2022) 9(12):458. doi: 10.3390/jcdd9120458
29. Yokokawa M, Sundaram B, Garg A, Stojanovska J, Oral H, Morady F, et al. Impact of mitral isthmus anatomy on the likelihood of achieving linear block in patients undergoing catheter ablation of persistent atrial fibrillation. Heart Rhythm. (2011) 8(9):1404–10. doi: 10.1016/j.hrthm.2011.04.030
30. Ding L, Moser M, Luo Y, Zhang WJ, Zhang B. Treatment planning optimization in irreversible electroporation for complete ablation of variously sized cervical tumors: a numerical study. J Biomech Eng. (2021) 143(1):014503. doi: 10.1115/1.4047551
31. Trainito CI, Français O, Le Pioufle B. Analysis of pulsed electric field effects on cellular tissue with cole–cole model: monitoring permeabilization under inhomogeneous electrical field with bioimpedance parameter variations. Innov Food Sci Emerg Technol. (2015) 29:193–200. doi: 10.1016/j.ifset.2015.02.004
32. Suárez AG, Hornero F, Berjano EJ. Mathematical modeling of epicardial RF ablation of atrial tissue with overlying epicardial fat. Open Biomed Eng J. (2010) 4:47–55. doi: 10.2174/1874120701004020047
33. Tungjitkusolmun S, Vorperian VR, Bhavaraju N, Cao H, Tsai JZ, Webster JG. Guidelines for predicting lesion size at common endocardial locations during radio-frequency ablation. IEEE Trans Biomed Eng. (2001) 48(2):194–201. doi: 10.1109/10.909640
34. Anjaneyulu A, Raghu K, Chandramukhi S, Satyajit GM, Arramraja S, Raghavaraju P, et al. Evaluation of left main coronary artery stenosis by transthoracic echocardiography. J Am Soc Echocardiogr. (2008) 21(7):855–60. doi: 10.1016/j.echo.2007.12.016
35. González-Suárez A, Berjano E. Comparative analysis of different methods of modeling the thermal effect of circulating blood flow during RF cardiac ablation. IEEE Trans Biomed Engs. (2015) 63(2):250–9. doi: 10.1109/TBME.2015.2451178
36. Neal RE II, Garcia PA, Robertson JL, Davalos RV. Experimental characterization and numerical modeling of tissue electrical conductivity during pulsed electric fields for irreversible electroporation treatment planning. IEEE Trans Biomed Eng. (2012) 59(4):1076–85. doi: 10.1109/TBME.2012.2182994
37. Zhao Y, Bhonsle S, Dong S, Lv YP, Liu HM, Safaai-Jazi A, et al. Characterization of conductivity changes during high-frequency irreversible electroporation for treatment planning. IEEE Trans Biomed Eng. (2017) 65(8):1810–9. doi: 10.1109/TBME.2017.2778101
38. Fang Z, Mao H, Moser MAJ, Zhang WJ, Qian ZQ, Zhang B. Irreversible electroporation enhanced by radiofrequency ablation: an in vitro and computational study in a 3D liver tumor model. Ann Biomed Eng. (2021) 49(9):2126–38. doi: 10.1007/s10439-021-02734-x
39. Pérez JJ, Ewertowska E, Berjano E. Computer modeling for radiofrequency bipolar ablation inside ducts and vessels: relation between pullback speed and impedance progress. Lasers Surg Med. (2020) 52(9):897–906. doi: 10.1002/lsm.23230
40. Gabriel C, Peyman A, Grant EH. Electrical conductivity of tissue at frequencies below 1 MHz. Phys Med Biol. (2009) 54(16):4863. doi: 10.1088/00319155/54/16/002
41. Sheehan MC, Srimathveeravalli G. Pulsed electric fields. Principles and technologies for electromagnetic based therapies. In: Prakash P, Srimathveeravalli G, editors. Principles and Technologies for Electromagnetic Energy Based Therapies. Cambridge, MA: Elsevier Academic Press (2021). p. 71–106. ISBN: 9780128205945.
42. Berjano EJ. Theoretical modeling for radiofrequency ablation: state-of-the-art and challenges for the future. Biomed Eng Online. (2006) 5(1):24. doi: 10.1186/1475-925X-5-24
43. Avazzadeh S, O’Brien B, Coffey K, O’Halloran M, Keane D, Quinlan LR, et al. Establishing irreversible electroporation electric field potential threshold in a suspension in vitro model for cardiac and neuronal cells. J Clin Med. (2021) 10(22):5443. doi: 10.3390/jcm10225443
44. Xie F, Chen Y, Chen X, Zhao ZH. Irreversible electroporation ablation for atrial fibrillation: status and challenges. Cardiol Discov. (2022) 2(01):41–50. doi: 10.1097/CD9.0000000000000045
45. Howard B, Haines DE, Verma A, Packer D, Kirchhof N, Barka N, et al. Reduction in pulmonary vein stenosis and collateral damage with pulsed field ablation compared with radiofrequency ablation in a canine model. Circ Arrhythm Electrophysiol. (2020) 13(9):e008337. doi: 10.1161/CIRCEP.120.008337
Keywords: pulsed electric field ablation technique, atrial fibrillation therapy, computer modeling, metal stent, ablation region
Citation: Wang Z, Liang M, Sun J, Zhang J, Li Y, Xu L and Han Y (2024) Epicardial pulsed-field ablation-impact of electric field and heat distribution induced by coronary metallic stents. Front. Cardiovasc. Med. 11:1445424. doi: 10.3389/fcvm.2024.1445424
Received: 7 June 2024; Accepted: 13 August 2024;
Published: 29 August 2024.
Edited by:
Marco Schiavone, Monzino Cardiology Center (IRCCS), ItalyReviewed by:
Kohei Sawasaki, Hamamatsu Medical Center, JapanMark Gallagher, St George’s University Hospitals NHS Foundation Trust, United Kingdom
Copyright: © 2024 Wang, Liang, Sun, Zhang, Li, Xu and Han. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yaling Han, hanyaling@163.net; Ming Liang, lming000919@sina.com