Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 17 March 2022
Sec. Computational Physiology and Medicine
This article is part of the Research Topic Exploring Mechanisms of Cardiac Rhythm Disturbances Using Novel Computational Methods: Prediction, Classification, and Therapy View all 19 articles

Anti-atrial Fibrillation Effects of Pulmonary Vein Isolation With or Without Ablation Gaps: A Computational Modeling Study

Ze JinZe JinInseok HwangInseok HwangByounghyun LimByounghyun LimOh-Seok KwonOh-Seok KwonJe-Wook ParkJe-Wook ParkHee-Tae YuHee-Tae YuTae-Hoon KimTae-Hoon KimBoyoung JoungBoyoung JoungMoon-Hyoung LeeMoon-Hyoung LeeHui-Nam Pak
Hui-Nam Pak*
  • Yonsei University College of Medicine, Yonsei University Health System, Seoul, South Korea

Background: Although pulmonary vein isolation (PVI) gaps contribute to recurrence after atrial fibrillation (AF) catheter ablation, the mechanism is unclear. We used realistic computational human AF modeling to explore the AF wave-dynamic changes of PVI with gaps (PVI-gaps).

Methods: We included 40 patients (80% male, 61.0 ± 9.8 years old, 92.5% persistent AF) who underwent AF catheter ablation to develop our realistic computational AF model. We compared the effects of a complete PVI (CPVI) and PVI-gap (2-mm × 4) on the AF wave-dynamics by evaluating the dominant frequency (DF), spatial change of DF, maximal slope of the action potential duration restitution curve (Smax), and AF defragmentation rate (termination or change to atrial tachycardia), and tested the effects of additional virtual interventions and flecainide on ongoing AF with PVI-gaps.

Results: Compared with the baseline AF, CPVIs significantly reduced extra-PV DFs (p < 0.001), but PVI-gaps did not. COV-DFs were greater after CPVIs than PVI-gaps (p < 0.001). Neither CPVIs nor PVI-gaps changed the mean Smax. CPVIs resulted in higher AF defragmentation rates (80%) than PVI-gaps (12.5%, p < 0.001). In ongoing AF after PVI-gaps, the AF defragmentation rates after a wave-breaking gap ablation, extra-PV DF ablation, or flecainide were 60.0, 34.3, and 25.7%, respectively (p = 0.010).

Conclusion: CPVIs effectively reduced the DF, increased its spatial heterogeneity in extra-PV areas, and offered better anti-AF effects than extra-PV DF ablation or additional flecainide in PVI-gap conditions.

Introduction

Catheter ablation is the most effective rhythm control method for atrial fibrillation (AF; Turagam et al., 2019). The cornerstone of AF catheter ablation (AFCA) is pulmonary vein isolation (PVI; Oral et al., 2006). However, even after AFCA with an adequate PVI, the AF recurrence rate is 40–50% within 5 years (Saguner et al., 2018). In patients with a post-AFCA recurrence during a repeat procedure, the pulmonary vein (PV) reconnection rate, which is the leading cause of recurrence, reaches 37–95% (Lin et al., 2015). PV reconnections are due to the technical limitations of the PVI durability (Arujuna et al., 2012). On the other hand, extra-PV triggers are also an important AF mechanism in AF patients with significant atrial remodeling (Velagapudi et al., 2013). Because AF is a progressive, degenerative disease, an empirical extra-PV ablation has been attempted, especially in patients with persistent or long-standing AF (Haissaguerre et al., 2005). Nevertheless, there is no evidence that an empirical extra-PV linear, electrogram-guided, or rotor ablation has a rhythm outcome superior to a complete PVI (CPVI) alone for ablating persistent AF (Verma et al., 2015). Indeed, the mechanism by which a CPVI and additional extra-PV ablation provide equivalent rhythm control is unknown. Little is known about the optimal procedure for AF patients with recurrence and an incomplete PVI (Park et al., 2019).

Therefore, this study used realistic human AF computational modeling to explore how a CPVI affects the extra-PV AF wave-dynamics and how a PVI-gap affects the recurrence mechanism. Current AF computational modeling takes advantage of recent improvements in the computational power to precisely simulate human AF by applying a personalized anatomy, histology, and electrophysiology (Lim et al., 2020a). Various virtual interventions and virtual anti-arrhythmic drug tests have become possible using controlled in silico conditions (Shim et al., 2017; Kim et al., 2019). In this study, we hypothesized that a CPVI would affect the extra-PV AF wave-dynamics. In addition, we evaluated the effects of a PVI with-gaps on the AF maintenance and compared the interventions and drug challenges to find the optimal anti-AF outcomes.

Materials and Methods

Ethical Approval

The study protocol followed the Declaration of Helsinki and was approved by the Institutional Review Board of Severance Cardiovascular Hospital, Yonsei University Health System. All participants were included in the Yonsei AF Ablation Cohort Database (ClinicalTrials.gov Identifier: NCT02138695) and provided written informed consent for us to use their cardiac CT images and clinical electrophysiological mapping data for computational modeling studies.

3D Computational Model of the Left Atrium

Figure 1 illustrates the protocol for the computational atrial modeling and AF simulations used in this study. We developed an ionic currents model according to the human atrial action potential model proposed by Courtemanche et al. (1998). For AF atrial ionic remodeling, the sodium current (INa), transient outward potassium current (Ito), L-type calcium current (ICaL), and ultrarapid outward current (IKur) decreased by 10, 70, 50, and 50%, respectively, and the inwardly rectifying potassium current (IK1) and Na+/Ca2+ exchanger (INCX) increased by 100 and 40%, respectively (Lee et al., 2016). The surface of the left atrial (LA) 3D model was composed of triangular meshes containing 400,000–500,000 geometric elements, and the mean distance between adjacent elements was 235.1 ± 32.1 μm. Interpolated voltage data were generated from bipolar electrograms recorded from >500 points on the atrial surface using a circular mapping catheter and CT images (Figure 1A). Artifact caused by the patient’s breathing was removed by trimming the ostial positions on the PVs and mitral valve. The coordinates of the electroanatomical map (NavX, Abbott, Inc., Chicago, IL, United States; CARTO System, Biosense Webster, Diamond Bar, CA, United States) were precisely aligned with the patients’ clinical heart CT images, and then the registration between the electroanatomical maps and clinical CT data was completed. We used the inverse distance weighting method (Ugarte et al., 2015) to represent the interpolation of the electroanatomical map values during the simulation procedures. Our graphical user interface software (Model: SH01, CUVIA ver. 2.5; Laonmed Inc., Seoul, Korea) integrated the fibrosis formation and fiber orientation into the LA surface and enabled a virtual AF induction and AF wave-dynamic changes by using the dominant frequency (DF) and maximal slope of the action potential duration restitution curve (Smax; Lim et al., 2017). The fiber orientations were defined in the meshes of each patient’s geometry and adjusted based on the clinical local activation time map (Pashakhanloo et al., 2016). Bipolar voltage data obtained from catheter ablation mapping were matched onto the computational nodes of the LA 3D model, and the locations of the fibrotic areas were determined using the map (Figure 1B). The fibrosis status of each node was numerically defined (Hwang et al., 2019). We simulated the clinical local activation data by using the model, which reflected the cardiac structural orientation and fiber orientation (Figure 1C). The conductivity of the model was based on the status and shape of the fibrosis (Zahid et al., 2016). For the ion currents of fibrotic cells, the inward rectifier potassium current (IK1), L-type calcium current (ICaL), and sodium current (INa) were decreased by 50, 50, and 40%, respectively, as compared to normal cells (Zahid et al., 2016). The reaction–diffusion equation for the cardiac wave propagation was solved numerically and adjusted based on the specific conduction velocity in each case to represent personalized AF simulations (Lee et al., 2016).

FIGURE 1
www.frontiersin.org

Figure 1. Study protocol of the computational atrial modeling and AF simulation. (A) CT merged 3D-clinical electroanatomical map. (B) Clinical map integrated computational modeling. (C) Baseline AF induction in activation map. (D) Application of CPVI and PVI-gap. The yellow arrows indicate gaps in the PVI. (E) DF-based wave-dynamics analyses. (F) Study protocol. CT indicates computed tomography; EP, electrophysiology; LAT, local activation time; PVI-gap, pulmonary vein isolation with gap; CPVI, complete pulmonary vein isolation; DF, dominant frequency; AF, atrial fibrillation; AT, atrial tachycardia.

AF Simulation

Figure 1F shows the process used in the study protocol. We induced AF in each case using AF pacing from 200 to 120 ms with eight beats per cycle lasting a total of 11,520 ms, based on appropriate ion current settings. Each virtual pacing location corresponded to the clinical activation time map for a realistic LA modeling, and the pacing sites were matched precisely to reflect the personalized LA models. AF maintenance was observed for 32 s. We defined a successful AF induction according to the electrograms in each LA model. AF defragmentation involved AF termination and AF conversion to atrial tachycardia.

Virtual Interventions

We applied a virtual ablation and virtual anti-arrhythmia drug to our realistic AF model. For the virtual ablation, the membrane potential of the ablated regions was set as zero to produce a permanent conduction block interrupting the cardiac wave propagation. First, we performed a virtual PV isolation with four gaps located on the anterior and posterior sides of each pulmonary vein isolation line (Figure 1D). Then we ablated all those gaps to perform a circumferential pulmonary vein isolation. The size of each gap was 2 mm. We initiated AF induction after performing each ablation. We applied high dose virtual flecainide (15 μm) to the failed AF defragmentation cases. With the Courtemanche-Ramirez-Nattel model (Courtemanche et al., 1998; Sossalla et al., 2010; Grandi et al., 2011) as the baseline, the effects of flecainide were implemented by applying the percent changes for the specific ion channels. Compared with baseline AF, the Ito, ICaL, IKr, and INa were decreased by 5, 5, 30 and 45%, respectively (Wang et al., 1993; Hilliard et al., 2010; Kramer et al., 2013; Crumb et al., 2016; Geng et al., 2018; Sutanto et al., 2019).

Analysis of the Spatial Changes in the AF Wave-Dynamics

We analyzed the wave-dynamics of the DF and Smax from 17 to 23 s. The DF was defined as the frequency with the highest power. It was analyzed using a Fourier transform for 6 s of action potentials at each node and the power spectra density function (Lim et al., 2017). We calculated the DF values for all nodes in the 3D LA model (Figure 1E). To examine the stability of the DF after each intervention, the coefficient of variation (COV) of the high-DF was calculated as the standard deviation divided by the mean. The Smax values were defined at every node in the LA regions of each patient.

Statistical Analysis

The continuous variables are represented as the median and range. The proportion of categorical variables was compared among the groups using a Fisher’s exact test. Comparisons of the DF, Smax, and COV-DF were conducted using t-testing or Mann–Whitney testing, depending on the distribution. A value of p <0.05 was considered statistically significant. Statistical analyses were conducted using SPSS (Statistical Package for Social Sciences, Chicago, IL, United States) software for Windows (version 26.0).

Results

Wave-Dynamic Changes After the PVI

We induced virtual AF in realistic computational models obtained from 40 patients (Supplementary Table S1). PV interventions (40 CPVI and 40 PVI-gap) significantly increased the mean AF cycle-lengths (140.9 ms [129.9, 153.3] to 147.1 ms [131.6, 200.3], p = 0.045), but they did not change the mean DF (7.69 Hz [7.31, 8.47] to 7.65 Hz [6.27, 8.30], p = 0.150) or Smax (0.97 [0.81, 1.26] to 0.88 [0.66, 1.13], p = 0.144). The AF termination and defragmentation rates after the PV interventions were 20 and 46.3%, respectively. Any episode terminated before 17 s was excluded from the AF cycle-length, DF, and Smax analyses.

Effects of the CPVI vs. PVI-Gaps on the Extra-PV Area

Table 1 compares the AF wave-dynamic changes between the CPVI and PVI-gap procedures. Compared with baseline AF, the CPVI significantly increased the mean AF cycle-lengths (p < 0.001), but the PVI-gap did not (p = 0.581). The CPVI also reduced the mean DF (p < 0.001), but the PVI-gap did not (p = 0.354, Figure 2). The amount of the DF reduction was significantly greater in the CPVI group than PVI-gap group (p < 0.001). The COV-DF, which reflected the spatial instability of the DF, was significantly greater in the CPVI group than PVI-gap group (p < 0.001). However, neither the CPVI (p = 0.445) nor PVI-gap (p = 0.078) changed the mean Smax (Table 1).

TABLE 1
www.frontiersin.org

Table 1. Wave-dynamic changes and defragmentation after virtual PV interventions.

FIGURE 2
www.frontiersin.org

Figure 2. Wave-dynamic changes after the virtual PVI-gap and CPVI. (A) The ECGs were obtained at the blue * sites in the DF maps. (B) The red arrows indicate gaps in the PVI. (C) The 3D DF maps and ECGs indicate that the CPVI reduced the mean DF, increased the AF cycle lengths, and defragmented the AF, but the PVI-gap did not. However, neither the CPVI nor PVI-gap changed the Smax. DF indicates dominant frequency; Smax, the maximal slope of the restitution curves; PVI-gap, pulmonary vein isolation with gaps; CPVI, complete pulmonary vein isolation; AF, atrial fibrillation; AT, atrial tachycardia; ECG, electrocardiogram.

Defragmentation Rates After the CPVI vs. PVI-Gap and Flecainide

The AF defragmentation and termination rates after the PVI are listed in Table 1. There was no AF termination or defragmentation during the 32 s waiting period during the baseline AF. The AF termination rate after the CPVI was 32.5%, and it was 7.5% after the PVI-gap (p = 0.010). The AF defragmentation rates were also higher in the CPVI group (80.0%) than PVI-gap group (5%, p < 0.001). This finding was consistent after changing the location and number of PVI-gaps (Supplementary Table S2).

We added virtual flecainide to the failed defragmented episodes after the CPVI (n = 8) and PVI-gap (n = 35). The post-flecainide termination rate in the group after the CPVI (25.0%) tended to be higher than that in the PVI-gap group without statistical significance (2.9%, p = 0.084, Table 2).

TABLE 2
www.frontiersin.org

Table 2. Flecainide-induced defragmentation rates of ongoing AF after the CPVI vs. PVI-gap.

Mechanism of the AF Maintenance After the PVI-Gap

We evaluated the PV and extra-PV wave-dynamic interactions in 35 ongoing AF episodes after the PVI-gap intervention (Figure 3). The extra-PV DF was greater than the PV-DF in 28 of those episodes (80%), and the inside PV-DF was greater than the extra-PV DF in 7 episodes (20%). The wavelet interactions were maintained through the PVI-gaps for as long as the AF was maintained. The wave-breaks generally appeared at the wavelet exit of the gaps by a source-sink mismatch (Figure 4). The rate of a wave-break did not differ between the episodes in which the wavelet moved from the higher DF site to the lower DF site and that in which it moved from the lower DF site to the higher DF site (48.8% vs. 50%, p = 1.000).

FIGURE 3
www.frontiersin.org

Figure 3. Wave-dynamic interactions between the PVs and extra-PVs in 35 ongoing AF episodes after the PVI-gap intervention. DF indicates dominant frequency; PV, pulmonary vein.

FIGURE 4
www.frontiersin.org

Figure 4. Wave-break at the wavelet exit of the gaps by a source-sink mismatch. (A) The 3D DF map indicates that the extra-PV DFs was greater than the PV-DFs. The activation map and ECG indicate that a regular wavelet turned into a wave-break when moving from the posterior wall of the LA to inside the PV through the gap. (B) The 3D DF map indicates that the inside PV-DF was greater than the extra-PV DF. The activation map and ECG indicate that a regular wavelet turned into a wave-break when moving from inside the PV to the lateral wall of the LA through the gap. DF indicates dominant frequency; PV, pulmonary vein; LA, left atrium; LSPV, left superior pulmonary vein; LIPV, left inferior pulmonary vein; ECG, Electrocardiogram.

Among the 35 episodes of ongoing AF after the PVI-gap, filling the wave-breaking PVI-gaps (total 91 gaps), ablation of the highest DF site without touching the PVI-gaps (112 extra-PV and 51 intra-PV sites), and additional virtual flecainide defragmented the AF in 60.0, 34.3, and 25.7% of cases, respectively (p = 0.010, Figure 5). The AF termination rate was significantly greater after filling the PVI-gaps than after ablating the highest DF or administering flecainide (p = 0.010, Table 3).

FIGURE 5
www.frontiersin.org

Figure 5. Wave-dynamic changes after filling the wave-breaking PVI-gaps, ablating the highest DF sites, and adding virtual flecainide. ECGs were obtained at the green * sites in the DF map, and the green arrows indicate the gaps in the PVI. (A) AF terminated after ablating the two gaps on the posterior wall of the LA. (B) AF terminated after ablating the high DF site (white areas indicated by black arrows). (C) AF converted to AT after adding virtual flecainide without additional ablation. PVI-gap indicates pulmonary vein isolation with gaps; PV, pulmonary vein; DF, dominant frequency; AF, atrial fibrillation; AT, atrial tachycardia; ECG, Electrocardiogram.

TABLE 3
www.frontiersin.org

Table 3. Defragmentation rates after interventions for 35 AF PVI-gap episodes.

Discussion

Main Findings

In this study, we used realistic human AF computational modeling to explore how a CPVI affected the extra-PV wave-dynamics and PVI-gaps as a mechanism for AF recurrence. The CPVI significantly reduced the mean extra-PV DF and its spatial stability (increased COV-DF), but it had no effect on the Smax, an index of wave-breakability. In the episodes of ongoing AF after the PVI-gap, wave-breaks commonly occurred at the wavelet exit of the gaps. Additional ablation of wave-breaking PVI-gaps had a greater defragmentation effect than extra-PV DF ablation or virtual flecainide. Therefore, the CPVI effectively reduced the mean DF in the extra-PV area, and under the PVI-gap condition, filling the PVI-gaps had anti-AF effects superior to those of an extra-PV DF ablation or additional flecainide.

Extra-PV Effects of the CPVI

The CPVI has traditionally been the cornerstone of AFCA (Oral et al., 2006), and it has several potential anti-AF mechanisms. First, the CPVI blocks the PV triggers (Oral et al., 2006). Second, the CPVI partially denervates cardiac autonomic nerves, including the ganglionated plexi around the PV antrum (Po et al., 2009). Third, a wide CPVI significantly reduces the atrial critical mass (Hwang et al., 2012). We have here identified another novel mechanism by which the CPVI significantly reduces the mean DF and spatial heterogeneity of the DF at extra-PV sites. Among the 40 patients who were included in this study, 93% had persistent AF, and we integrated their voltage-activation maps into our realistic computational models. We found that the CPVI still plays an important role in the extra-PV AF maintenance mechanism, even in AF patients with atrial substrate remodeling. To date, no empirical extra-PV LA ablation, including linear, electrogram-guided, or rotor ablation, has shown equivalent rhythm outcomes in patients with a CPVI, despite multiple randomized clinical trials (Verma et al., 2015). Our results reconfirm the importance of a CPVI in the LA, even though extra-PV areas, including the right atrium, contribute to the maintenance mechanisms of AF (Lim et al., 2020b).

Role of PV-Gaps as a Mechanism of Recurrence After a PVI

Irrespective of the uncertain unifying mechanisms of AF, a PV electrical isolation is an objective, standard, and widely accepted minimal requirement for AFCA (Calkins et al., 2012). However, despite advances in the catheter efficiency, a long-term durable PVI still has technical limitations. The rate of PV reconnections during a repeat ablation has been reported to range from 36% to more than 95% in repeat procedures (Lin et al., 2015). Many studies have reported that PV reconnections are the leading cause of arrhythmia recurrences. However, PV reconnections have also been observed in patients without AF recurrence (Nilsson et al., 2006). In this study, we evaluated how PV gaps affected the AF maintenance mechanisms. Of the four 2-mm PV gaps that remained in our model, 65% (91/140) were accompanied by wave-breaks and 19.3% (27/140) contributed to wave-breaks in the PV to LA direction. The PV gap size also had an essential influence on the AF maintenance. Following Herweg’s report, we applied 2-mm gaps, which were suitable for generating wave-breaks by a source-sink mismatch (Herweg et al., 2021).

Appropriate Management of Recurrent AF After AFCA

The most effective repeat procedures for patients with recurrence after AFCA is not yet known. Depending on the recurrent AF burden or the operators’ discretion, gap fillings for the reconnected PVs or an empirical extra-PV ablation, such as a posterior wall isolation, have been performed (Lee et al., 2019). In patients with a well-maintained PVI during redo-mapping, extra-PV triggers play a significant role in AF recurrence, and the outcomes of a repeat ablation are worse than in patients with PV reconnections (Kim et al., 2021). In this study, we examined the AF maintenance mechanism during the PVI-gap state by using computational models made with the electroanatomical maps of 40 patients, most of whom had persistent AF. Under those controlled conditions, filling the wave-breaking gap produced a more effective AF defragmentation than ablation of the extra-PV DF sites or using antiarrhythmic drugs. That was consistent with the results of a recent randomized clinical trial, suggesting that a durable PVI is more effective than an empirical extra-PV ablation as a repeat AF ablation procedure (Kim et al., 2022).

Limitations

This study had some limitations in its computational simulations. First, bi-atrial modeling manifesting interatrial conduction was not applied in this study. Second, we did not consider the LA wall thickness in our LA model. Third, because we applied fibrosis based on voltage-map, not a magnetic resonance image (Boyle et al., 2019; Baek et al., 2021), it is not clear whether the detailed structure of cardiac fibrosis with microchannels is adequately reflected in our morel. Fourth, our personalized LA model consisted of a monolayer, not bilayers representing the endocardial and epicardial layers. Fifth, the rate-dependent effect of flecainide on INa may not be reflected sufficiently in this study (Moreno et al., 2011).

Conclusion

The CPVI effectively reduced the mean DF and increased its spatial heterogeneity in the extra-PV areas. Filling the PVI-gaps had anti-AF effects superior to those of an extra-PV DF ablation or additional flecainide under the PVI-gap condition.

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 author.

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Board of Severance Cardiovascular Hospital, Yonsei University Health System. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

ZJ contributed to the data, statistical analyses, and writing of the manuscript. IH contributed to the statistical analyses and data acquisition. O-SK contributed to the software programming and data acquisition. BL confirmed the data acquisition and references. J-WP contributed the clinical data acquisition. H-TY, T-HK, BJ, and M-HL contributed to the clinical data acquisition and interpretation of clinical data. H-NP contributed to the study design, clinical data acquisition, data interpretation, and revision of manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants (HI19C0114) and (H21C0011) from the Ministry of Health and Welfare and grants (NRF-2020R1A2B01001695) and (NRF-2019R1C1C1009075) from the Basic Science Research Program run by the National Research Foundation of Korea (NRF), which is funded by the Ministry of Science, ICT, and Future Planning (MSIP).

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.

Acknowledgments

We would like to thank Mr. John Martin for his linguistic assistance.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2022.846620/full#supplementary-material

References

Arujuna, A., Karim, R., Caulfield, D., Knowles, B., Rhode, K., Schaeffter, T., et al. (2012). Acute pulmonary vein isolation is achieved by a combination of reversible and irreversible atrial injury after catheter ablation: evidence from magnetic resonance imaging. Circ. Arrhythm. Electrophysiol. 5, 691–700. doi: 10.1161/CIRCEP.111.966523

PubMed Abstract | CrossRef Full Text | Google Scholar

Baek, Y. S., Kwon, O. S., Lim, B., Yang, S. Y., Park, J. W., Yu, H. T., et al. (2021). Clinical outcomes of computational virtual mapping-guided catheter ablation in patients With persistent atrial fibrillation: A multicenter prospective randomized clinical trial. Front. Cardiovasc. Med. 8:772665. doi: 10.3389/fcvm.2021.772665

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyle, P. M., Zghaib, T., Zahid, S., Ali, R. L., Deng, D., Franceschi, W. H., et al. (2019). Computationally guided personalized targeted ablation of persistent atrial fibrillation. Nat. Biomed. Eng. 3, 870–879. doi: 10.1038/s41551-019-0437-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Calkins, H., Kuck, K. H., Cappato, R., Brugada, J., Camm, A. J., Chen, S. A., et al. (2012). 2012 HRS/EHRA/ECAS expert consensus statement on catheter and surgical ablation of atrial fibrillation: recommendations for patient selection, procedural techniques, patient management and follow-up, definitions, endpoints, and research trial design. J. Interv. Card. Electrophysiol. 33, 171–257. doi: 10.1007/s10840-012-9672-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Courtemanche, M., Ramirez, R. J., and Nattel, S. (1998). Ionic mechanisms underlying human atrial action potential properties: insights from a mathematical model. Am. J. Phys. 275, H301–H321. doi: 10.1152/ajpheart.1998.275.1.H301

PubMed Abstract | CrossRef Full Text | Google Scholar

Crumb, W. J. Jr., Vicente, J., Johannesen, L., and Strauss, D. G. (2016). An evaluation of 30 clinical drugs against the comprehensive in vitro proarrhythmia assay (CiPA) proposed ion channel panel. J. Pharmacol. Toxicol. Methods 81, 251–262. doi: 10.1016/j.vascn.2016.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Geng, L., Kong, C. W., Wong, A. O. T., Shum, A. M., Chow, M. Z. Y., Che, H., et al. (2018). Probing flecainide block of INa using human pluripotent stem cell-derived ventricular cardiomyocytes adapted to automated patch-clamping and 2D monolayers. Toxicol. Lett. 294, 61–72. doi: 10.1016/j.toxlet.2018.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Grandi, E., Pandit, S. V., Voigt, N., Workman, A. J., Dobrev, D., Jalife, J., et al. (2011). Human atrial action potential and Ca2+ model: sinus rhythm and chronic atrial fibrillation. Circ. Res. 109, 1055–1066. doi: 10.1161/CIRCRESAHA.111.253955

PubMed Abstract | CrossRef Full Text | Google Scholar

Haissaguerre, M., Hocini, M., Sanders, P., Sacher, F., Rotter, M., Takahashi, Y., et al. (2005). Catheter ablation of long-lasting persistent atrial fibrillation: clinical outcome and mechanisms of subsequent arrhythmias. J. Cardiovasc. Electrophysiol. 16, 1138–1147. doi: 10.1111/j.1540-8167.2005.00308.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Herweg, B., Nellaiyappan, M., Welter-Frost, A. M., Tran, T., Mabry, G., Weston, K., et al. (2021). Immuno-electrophysiological mechanisms of functional electrical connections Between recipient and donor heart in patients With Orthotopic heart transplantation presenting With atrial arrhythmias. Circ. Arrhythm. Electrophysiol. 14:e008751. doi: 10.1161/CIRCEP.120.008751

CrossRef Full Text | Google Scholar

Hilliard, F. A., Steele, D. S., Laver, D., Yang, Z., Le Marchand, S. J., Chopra, N., et al. (2010). Flecainide inhibits arrhythmogenic Ca2+ waves by open state block of ryanodine receptor Ca2+ release channels and reduction of Ca2+ spark mass. J. Mol. Cell. Cardiol. 48, 293–301. doi: 10.1016/j.yjmcc.2009.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, M., Kim, J., Lim, B., Song, J. S., Joung, B., Shim, E. B., et al. (2019). Multiple factors influence the morphology of the bipolar electrogram: an in silico modeling study. PLoS Comput. Biol. 15:e1006765. doi: 10.1371/journal.pcbi.1006765

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, E. S., Nam, G. B., Joung, B., Park, J., Lee, J. S., Shim, J., et al. (2012). Significant reduction of atrial defibrillation threshold and inducibility by catheter ablation of atrial fibrillation. Pacing Clin. Electrophysiol. 35, 1428–1435. doi: 10.1111/j.1540-8159.2012.03517.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Shim, J., Kim, Y. G., Yu, H. T., Kim, T. H., Uhm, J. S., et al. (2021). Malnutrition and risk of procedural complications in patients With atrial fibrillation undergoing catheter ablation. Front. Cardiovasc. Med. 8:736042. doi: 10.3389/fcvm.2021.736042

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D., Yu, H. T., Kim, T.-H., Uhm, J.-S., Joung, B., Lee, M.-H., et al. (2022). Electrical posterior box isolation in repeat ablation for atrial fibrillation: a prospective randomized clinical study. J. Am. Coll. Cardiol. Clin. Electrophysiol. doi: 10.1016/j.jacep.2022.01.003 (in press)

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, I. S., Lim, B., Shim, J., Hwang, M., Yu, H. T., Kim, T. H., et al. (2019). Clinical usefulness of computational modeling-guided persistent atrial fibrillation ablation: updated outcome of multicenter randomized study. Front. Physiol. 10:1512. doi: 10.3389/fphys.2019.01512

PubMed Abstract | CrossRef Full Text | Google Scholar

Kramer, J., Obejero-Paz, C. A., Myatt, G., Kuryshev, Y. A., Bruening-Wright, A., Verducci, J. S., et al. (2013). MICE models: superior to the HERG model in predicting torsade de pointes. Sci. Rep. 3:2100. doi: 10.1038/srep02100

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, Y. S., Hwang, M., Song, J. S., Li, C., Joung, B., Sobie, E. A., et al. (2016). The contribution of ionic currents to rate-dependent action potential duration and pattern of reentry in a mathematical model of human atrial fibrillation. PLoS One 11:e0150779. doi: 10.1371/journal.pone.0150779

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. M., Shim, J., Park, J., Yu, H. T., Kim, T. H., Park, J. K., et al. (2019). The electrical isolation of the left atrial Posterior Wall in catheter ablation of persistent atrial fibrillation. JACC Clin. Elect. 5, 1253–1261. doi: 10.1016/j.jacep.2019.08.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, B., Hwang, M., Song, J. S., Ryu, A. J., Joung, B., Shim, E. B., et al. (2017). Effectiveness of atrial fibrillation rotor ablation is dependent on conduction velocity: An in-silico 3-dimensional modeling study. PLoS One 12:e0190398. doi: 10.1371/journal.pone.0190398

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, B., Kim, J., Hwang, M., Song, J. S., Lee, J. K., Yu, H. T., et al. (2020a). In situ procedure for high-efficiency computational modeling of atrial fibrillation reflecting personal anatomy, fiber orientation, fibrosis, and electrophysiology. Sci. Rep. 10:2417. doi: 10.1038/s41598-020-59372-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, B., Park, J. W., Hwang, M., Ryu, A. J., Kim, I. S., Yu, H. T., et al. (2020b). Electrophysiological significance of the interatrial conduction including cavo-tricuspid isthmus during atrial fibrillation. J. Physiol. 598, 3597–3612. doi: 10.1113/JP279660

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, D., Santangeli, P., Zado, E. S., Bala, R., Hutchinson, M. D., Riley, M. P., et al. (2015). Electrophysiologic findings and long-term outcomes in patients undergoing third or more catheter ablation procedures for atrial fibrillation. J. Cardiovasc. Electrophysiol. 26, 371–377. doi: 10.1111/jce.12603

CrossRef Full Text | Google Scholar

Moreno, J. D., Zhu, Z. I., Yang, P. C., Bankston, J. R., Jeng, M. T., Kang, C., et al. (2011). A computational model to predict the effects of class I anti-arrhythmic drugs on ventricular rhythms. Sci. Transl. Med. 3:98ra83. doi: 10.1126/scitranslmed.3002588

CrossRef Full Text | Google Scholar

Nilsson, B., Chen, X., Pehrson, S., Kober, L., Hilden, J., and Svendsen, J. H. (2006). Recurrence of pulmonary vein conduction and atrial fibrillation after pulmonary vein isolation for atrial fibrillation: a randomized trial of the ostial versus the extraostial ablation strategy. Am. Heart J. 152, 537.e1–537.e8. doi: 10.1016/j.ahj.2006.05.029

CrossRef Full Text | Google Scholar

Oral, H., Pappone, C., Chugh, A., Good, E., Bogun, F., Pelosi, F. Jr., et al. (2006). Circumferential pulmonary-vein ablation for chronic atrial fibrillation. N. Engl. J. Med. 354, 934–941. doi: 10.1056/NEJMoa050955

CrossRef Full Text | Google Scholar

Park, J. W., Yang, P. S., Bae, H. J., Yang, S. Y., Yu, H. T., Kim, T. H., et al. (2019). Five-year change in the renal function After catheter ablation of atrial fibrillation. J. Am. Heart Assoc. 8:e013204. doi: 10.1161/JAHA.119.013204

CrossRef Full Text | Google Scholar

Pashakhanloo, F., Herzka, D. A., Ashikaga, H., Mori, S., Gai, N., Bluemke, D. A., et al. (2016). Myofiber architecture of the human atria as revealed by submillimeter diffusion tensor imaging. Circ. Arrhythm. Electrophysiol. 9:e004133. doi: 10.1161/CIRCEP.116.004133

CrossRef Full Text | Google Scholar

Po, S. S., Nakagawa, H., and Jackman, W. M. (2009). Localization of left atrial ganglionated plexi in patients with atrial fibrillation. J. Cardiovasc. Electrophysiol. 20, 1186–1189. doi: 10.1111/j.1540-8167.2009.01515.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Saguner, A. M., Maurer, T., Wissner, E., Santoro, F., Lemes, C., Mathew, S., et al. (2018). Catheter ablation of atrial fibrillation in very young adults: a 5-year follow-up study. Europace 20, 58–64. doi: 10.1093/europace/euw378

PubMed Abstract | CrossRef Full Text | Google Scholar

Shim, J., Hwang, M., Song, J. S., Lim, B., Kim, T. H., Joung, B., et al. (2017). Virtual in-Silico modeling guided catheter ablation predicts effective linear ablation lesion set for longstanding persistent atrial fibrillation: multicenter prospective randomized study. Front. Physiol. 8:792. doi: 10.3389/fphys.2017.00792

PubMed Abstract | CrossRef Full Text | Google Scholar

Sossalla, S., Kallmeyer, B., Wagner, S., Mazur, M., Maurer, U., Toischer, K., et al. (2010). Altered Na(+) currents in atrial fibrillation effects of ranolazine on arrhythmias and contractility in human atrial myocardium. J. Am. Coll. Cardiol. 55, 2330–2342. doi: 10.1016/j.jacc.2009.12.055

PubMed Abstract | CrossRef Full Text | Google Scholar

Sutanto, H., Laudy, L., Clerx, M., Dobrev, D., Crijns, H., and Heijman, J. (2019). Maastricht antiarrhythmic drug evaluator (MANTA): a computational tool for better understanding of antiarrhythmic drugs. Pharmacol. Res. 148:104444. doi: 10.1016/j.phrs.2019.104444

PubMed Abstract | CrossRef Full Text | Google Scholar

Turagam, M. K., Garg, J., Whang, W., Sartori, S., Koruth, J. S., Miller, M. A., et al. (2019). Catheter ablation of atrial fibrillation in patients With heart failure: a meta-analysis of randomized controlled trials. Ann. Intern. Med. 170, 41–50. doi: 10.7326/M18-0992

CrossRef Full Text | Google Scholar

Ugarte, J. P., Tobon, C., Orozco-Duque, A., Becerra, M. A., and Bustamante, J. (2015). Effect of the electrograms density in detecting and ablating the tip of the rotor during chronic atrial fibrillation: an in silico study. Europace 17(Suppl. 2), ii97–ii104. doi: 10.1093/europace/euv244

PubMed Abstract | CrossRef Full Text | Google Scholar

Velagapudi, P., Turagam, M. K., Leal, M. A., and Kocheril, A. G. (2013). Atrial fibrosis: a risk stratifier for atrial fibrillation. Expert. Rev. Cardiovasc. Ther. 11, 155–160. doi: 10.1586/erc.12.174

PubMed Abstract | CrossRef Full Text | Google Scholar

Verma, A., Jiang, C. Y., Betts, T. R., Chen, J., Deisenhofer, I., Mantovan, R., et al. (2015). Approaches to catheter ablation for persistent atrial fibrillation. N. Engl. J. Med. 372, 1812–1822. doi: 10.1056/NEJMoa1408288

CrossRef Full Text | Google Scholar

Wang, Z., Fermini, B., and Nattel, S. (1993). Mechanism of flecainide’s rate-dependent actions on action potential duration in canine atrial tissue. J. Pharmacol. Exp. Ther. 267, 575–581.

PubMed Abstract | Google Scholar

Zahid, S., Cochet, H., Boyle, P. M., Schwarz, E. L., Whyte, K. N., Vigmond, E. J., et al. (2016). Patient-derived models link re-entrant driver localization in atrial fibrillation to fibrosis spatial pattern. Cardiovasc. Res. 110, 443–454. doi: 10.1093/cvr/cvw073

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: atrial fibrillation, computational modeling, pulmonary vein, gap, dominant frequency

Citation: Jin Z, Hwang I, Lim B, Kwon O-S, Park J-W, Yu H-T, Kim T-H, Joung B, Lee M-H and Pak H-N (2022) Anti-atrial Fibrillation Effects of Pulmonary Vein Isolation With or Without Ablation Gaps: A Computational Modeling Study. Front. Physiol. 13:846620. doi: 10.3389/fphys.2022.846620

Received: 31 December 2021; Accepted: 24 February 2022;
Published: 17 March 2022.

Edited by:

Martin Bishop, King’s College London, United Kingdom

Reviewed by:

Mohit Turagam, Mount Sinai Hospital, United States
Jordi Heijman, Maastricht University, Netherlands

Copyright © 2022 Jin, Hwang, Lim, Kwon, Park, Yu, Kim, Joung, Lee and Pak. 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: Hui-Nam Pak, hnpak@yuhs.ac

Disclaimer: 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.