Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 11 October 2017
Sec. Computational Physiology and Medicine
This article is part of the Research Topic Clinical Applications of Physiome Models View all 11 articles

Virtual In-Silico Modeling Guided Catheter Ablation Predicts Effective Linear Ablation Lesion Set for Longstanding Persistent Atrial Fibrillation: Multicenter Prospective Randomized Study

\r\nJaemin Shim&#x;Jaemin Shim1Minki Hwang&#x;Minki Hwang2Jun-Seop SongJun-Seop Song2Byounghyun LimByounghyun Lim2Tae-Hoon KimTae-Hoon Kim2Boyoung JoungBoyoung Joung2Sung-Hwan KimSung-Hwan Kim3Yong-Seog OhYong-Seog Oh3Gi-Byung NamGi-Byung Nam4Young Keun OnYoung Keun On5Seil OhSeil Oh6Young-Hoon KimYoung-Hoon Kim1Hui-Nam Pak*Hui-Nam Pak2*
  • 1Cardiovascular Center, Korea University, Seoul, South Korea
  • 2Division of Cardiology, Yonsei University Health System, Seoul, South Korea
  • 3Division of Cardiology, Catholic University of Korea, Seoul, South Korea
  • 4Asan Medical Center, University of Ulsan, Seoul, South Korea
  • 5Samsung Medical Center, Sungkyunkwan University, Seoul, South Korea
  • 6Division of Cardiology, Seoul National University, Seoul, South Korea

Objective: Radiofrequency catheter ablation for persistent atrial fibrillation (PeAF) still has a substantial recurrence rate. This study aims to investigate whether an AF ablation lesion set chosen using in-silico ablation (V-ABL) is clinically feasible and more effective than an empirically chosen ablation lesion set (Em-ABL) in patients with PeAF.

Methods: We prospectively included 108 patients with antiarrhythmic drug-resistant PeAF (77.8% men, age 60.8 ± 9.9 years), and randomly assigned them to the V-ABL (n = 53) and Em-ABL (n = 55) groups. Five different in-silico ablation lesion sets [1 pulmonary vein isolation (PVI), 3 linear ablations, and 1 electrogram-guided ablation] were compared using heart-CT integrated AF modeling. We evaluated the feasibility, safety, and efficacy of V-ABL compared with that of Em-ABL.

Results: The pre-procedural computing time for five different ablation strategies was 166 ± 11 min. In the Em-ABL group, the earliest terminating blinded in-silico lesion set matched with the Em-ABL lesion set in 21.8%. V-ABL was not inferior to Em-ABL in terms of procedure time (p = 0.403), ablation time (p = 0.510), and major complication rate (p = 0.900). During 12.6 ± 3.8 months of follow-up, the clinical recurrence rate was 14.0% in the V-ABL group and 18.9% in the Em-ABL group (p = 0.538). In Em-ABL group, clinical recurrence rate was significantly lower after PVI+posterior box+anterior linear ablation, which showed the most frequent termination during in-silico ablation (log-rank p = 0.027).

Conclusions: V-ABL was feasible in clinical practice, not inferior to Em-ABL, and predicts the most effective ablation lesion set in patients who underwent PeAF ablation.

Introduction

The prevalence of AF increases with age, and ~20% of all ischemic strokes are associated with AF (Goldstein et al., 2011). It is predicted that the prevalence of AF would more than double by the year 2050 (Ball et al., 2013). Although radiofrequency catheter ablation (RFCA) has been established as an effective rhythm control strategy, it still has a substantial long-term recurrence rate, especially in patients with persistent AF (Dewire and Calkins, 2013). AF catheter ablation introduces the risk of collateral damages, and a long duration of ablation procedure is associated with poor clinical outcomes (Shim et al., 2013). Moreover, the ablation strategy is selected largely based on the physician's experience, which bears the possibility of outcome differences depending on the physician. On the other hand, computer simulation has been widely used in the study of the mechanisms of AF (Hwang et al., 2014; Trayanova, 2014). Recently, a number of attempts have been made to use computer simulation for studying the effects of catheter ablation (Zhao et al., 2015). We previously conducted a retrospective study and reported that the earliest virtual AF terminating lesion set chosen by using a computer was commonly identical to that selected by the physician who performed the ablations in patients with persistent AF (Hwang et al., 2014). Because the computing speed for complicated in-silico modeling of AF is one of the limitations for the application of computer simulation in clinical practice, we have shortened the computing time for AF modeling by about 80 times compared with conventional modeling (Hwang et al., 2015; Zhao et al., 2015). This study was a prospective randomized trial to examine the applicability, safety, and efficacy of computer simulation for catheter ablation of AF. The ablation strategy was selected on the basis of computer simulation for one group of patients, and according to the physician's experience for the other group. Then, we evaluated the feasibility of in-silico ablation (computing time), and compared procedure time, procedure related complication rate, long-term rhythm outcome of both groups, and the most effective lesion set in in-silico ablation and empirical ablation (Em-ABL).

Methods

Study Design

The study protocol adhered to the Declaration of Helsinki and was approved by the institutional review board of the Yonsei University Health System (clinicaltrials.gov; NCT 02171364). The enrolled patients were randomly assigned to either the virtual ablation (V-ABL) group or the Em-ABL group (Figure 1). For all patients in both groups, virtual ablations were performed by applying the following 5 strategies: (i) circumferential pulmonary vein isolation (CPVI), (ii) CPVI + posterior box lesion, (iii) CPVI + posterior box lesion + anterior line, (iv) CPVI + roof line + left lateral isthmus line, and (v) CPVI + complex fractionated atrial electrogram (CFAE)-guided ablation (Hwang et al., 2014; Song et al., 2016). For patients in the V-ABL group, the best virtual ablation strategy was applied in real ablation. The virtual ablation strategy that resulted in the earliest termination of AF after in-silico ablation was determined as the best in-silico ablation lesion set among the five ablation strategies. For patients in the Em-ABL group, the ablation strategy selected by the physician based on experience was applied in real ablation. The ablation outcomes were compared between the V-ABL and Em-ABL groups.

FIGURE 1
www.frontiersin.org

Figure 1. Study flow chart. The enrolled patients were randomly assigned to either the virtual ablation (V-ABL) group or the empirical ablation (Em-ABL) group. Virtual ablations were performed for all patients in both groups, by applying five strategies. For patients in the V-ABL group, the best virtual ablation strategy was applied in real ablation. For patients in the Em-ABL group, the ablation strategy selected by the physician based on experience was applied in real ablation. The ablation outcomes were compared between the V-ABL and Em-ABL groups.

Study Population

In this study, we enrolled 118 patients with persistent AF (PeAF) who were eligible for RFCA. All patients provided written informed consent. The inclusion criteria were (i) age 19 years or older, (ii) antiarrhythmic drug resistance, and (iii) availability of three-dimensional (3D) computed tomography (CT) or magnetic resonance imaging (MRI) scans of the left atrium (LA). The exclusion criteria were (i) presence of significant structural heart disease, (ii) valvular AF, (iii) previous ablation or cardiac surgery, and (iv) presence of a cardiac implantable electronic device.

After receiving the patients' consent, DICOM files of heart CT were sent to the core laboratory a night before or in the early morning of the real AF ablation procedure. Patient randomization, in-silico AF modeling, and virtual AF ablation were done at the core laboratory during work hours (9.a.m.–6 p.m.). Among 118 patients, four were excluded owing to incorrect patient selection or patient request. Additional six patients were excluded because of failed in-silico AF ablation before the real ablation procedure, due to DICOM file error (n = 2) and communication error (n = 4). Finally, 108 patients (53 in the V-ABL group and 55 in the Em-ABL group) were enrolled in this study (Figure 1). The patient characteristics are shown in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Patient characteristics.

In-Silico Simulation of AF

Atrial geometries were reconstructed from the 3D CT merged NavX data (NavX; St. Jude Medical Inc., Minnetonka, MN, USA) defining the surface of the atrium. A triangular mesh was generated on the surface of the 3D atrial geometry. The final number of grid elements was between 400,000 and 500,000. The LA appendage and myocardial sleeves were included in the mesh. For the simulation of cardiac wave propagation in the atrial wall, the following reaction-diffusion equation was solved numerically (Zozor et al., 2003):

Vmt=1βCm{·DVm-β(Iion+Is)},    (1)

where Vm is the membrane potential (Unit: volt); β is the membrane surface-to-volume ratio (Unit: meter−1); Cm is the membrane capacitance per unit area (Unit: farad/meter2); D is conductivity tensor (Unit: siemens/meter); and Iion and Is are the ion current and stimulation current, respectively (Unit: ampere/meter2).

For the calculation of ionic currents, a mathematical model of human atrial action potential proposed by Courtemanche et al. (1998) was used. Electrical stimulation was applied at the location of Bachmann's bundle, and reentry was initiated by rapid pacing: a total of 24 pacings with pacing cycle lengths of 200, 190, and 180 ms. To induce self-sustained fibrillation, Ito, ICaL, and IKur were reduced by 80, 40, and 50%, respectively, and IK1 was increased by 50% (Li et al., 2016). We chose a conduction velocity of 0.4 m/s based on real human patient data (Yonsei AF ablation cohort data; n = 1,980; mean CV = 0.43 ± 0.24 m/s) (Park et al., 2014) and previous simulation studies (Hwang et al., 2014; Li et al., 2016).

Virtual AF Ablation

Virtual ablation was conducted for all patients in both the V-ABL and Em-ABL groups. We developed a graphical user interface software with which the user can perform virtual ablation by mouse-clicking on the atrial geometry (CUVIA, Model: SH01, ver. 1.0; Laonmed Inc., Seoul, Korea; Figure 2). At the point where the mouse is clicked, a 2-mm-diameter circular region is virtually ablated, mimicking real ablation with a radiofrequency catheter. In the ablated region, no electrical flux condition was applied. The pattern of ablation was completed by serial mouse-clickings, similarly to real catheter ablation. For CFAE-guided ablations (Nademanee et al., 2004), the areas with CFAE cycle length < 120 ms were ablated as long as the total ablated area is < 5% of the whole atrial area. CFAE was observed in particular region in structurally homogenous simulation due to the three-dimensional geometrical variation. The area where the mitral annulus was located was also considered a non-conductive area. Virtual ablation was applied at 4 s after the end of pacings, and the simulation was run until 25 s from the start of the simulation. When all points in the atrial wall were repolarized by > 90%, the fibrillation was considered terminated. The duration until the termination of fibrillation was recorded for each ablation strategy for each patient. NavX potential map was not used in the virtual ablation. The physicians, however, utilized the endocardial potential map data in their decision makings. Fiber orientation was not considered in the simulation model, and conductivity was isotropic even though three-dimensional geometrical shape may have affected conductivity.

FIGURE 2
www.frontiersin.org

Figure 2. Virtual ablation procedure. Computed tomography images are segmented to construct three-dimensional atrial geometry. Virtual atrial fibrillation (AF) is induced. The maps of cardiac wave dynamics parameters such as complex fractionated atrial electrogram (CFAE) are generated. Virtual ablation is performed on the atrial geometry. AF simulation is resumed. The cardiac wave propagation pattern is observed.

Real AF Ablation

The electrophysiological mapping method and the RFCA technique were as described previously (Mun et al., 2012). Briefly, we used an open irrigated-tip catheter (Celsius [Johnson & Johnson Inc., Diamond Bar, CA, USA] or Coolflex [St. Jude Medical]; 30–35 W, 47°C) to deliver radiofrequency energy for ablation under 3D electroanatomical mapping (NavX, St. Jude Medical; CARTO3, Johnson & Johnson) merged with 3D spiral CT. The procedure was completed when there was no immediate recurrence of AF after cardioversion with isoproterenol infusion (5 μg/min). All RFCA procedures were performed by an operator with > 10 years of experience.

Post-ablation Management and Follow-Up

As per RFCA procedure protocols, antiarrhythmic drugs were discontinued in all patients. The patients visited the outpatient clinic at 1, 3, 6, and 12 months after RFCA, or whenever symptoms occurred. All patients underwent electrocardiography at each visit, and 24-h Holter recording was done at 3, 6, and 12 months, according to the 2012 Heart Rhythm Society/European Heart Rhythm Association/European Cardiac Arrhythmia Society Expert Consensus Statement (Calkins et al., 2012). However, whenever patients reported palpitations, Holter monitor or event monitor recordings were obtained and evaluated to check for arrhythmia recurrence. We defined recurrence of AF as any episode of AF or atrial tachycardia lasting for at least 30 s. Any electrocardiography documentation of AF recurrence after 3 months of the blanking period was diagnosed as clinical recurrence.

Statistical Analysis

Because the trial was conducted as a pilot study, the sample size was driven by the computation time and feasibility of recruiting eligible patients, and enrollment was open for 12 months. Results are expressed as mean values ± standard deviation for continuous variables, and absolute number and percentages for categorical variables. Continuous variables were compared by using Student's t-test, and categorical variables were compared by using either the chi-square test or Fisher's exact test as appropriate. The primary end point of the study was freedom from any atrial arrhythmias during follow-up after a 3-month blanking period. The time to recurrence and arrhythmia-free survival were assessed by using Kaplan-Meier analysis, and differences were calculated by using the log-rank test. A p-value of < 0.05 was considered statistically significant. All statistical analyses were performed by using SPSS version 18.0 software (SPSS Inc., Chicago, IL, USA).

Results

Patient Characteristics

The patients' characteristics are summarized in Table 1. The mean age was 60.8 ± 9.6 years, and 76.9% were male patients. Among 108 patients with PeAF, 77.8% had longstanding PeAF (lasting for >1 year). More than 28% of patients had a history of stroke or transient ischemic attack; the mean left atrial dimension was 45.1 ± 4.4 mm and the mean ejection fraction was 59.3 ± 9.7%. Randomization was done well, and there was no statistically significant difference between the V-ABL and Em-ABL groups.

In-Silico Ablation Results

In-silico AF ablations were done in all included patients; however, the Em-ABL group was blinded to the in-silico ablation results. The pre-procedural computing time for the five virtual ablation strategies described in the Methods section was 166 ± 11 min. Heart CT segmentation with a semi-automatic method and five different in-silico ablation lesion settings with the manual method took about 1 h. Another 1 h was taken for computing the five different in-silico ablation protocols simultaneously, and additional 30 min was required for the calculation of the CFAE area.

Table 2 summarizes the outcome of in-silico AF ablation. After virtual ablation, we waited for 25 s, measuring the termination time. Among the five in-silico ablation strategies, [CPVI+Posterior box+Anterior linear ablation] showed the highest AF termination rate (81.5%, 88 of 108) and shortest time to AF termination (16.792 ± 5.672 ms). The AF termination rates within 25 s were 73.1% after [CPVI+Roof line+Left lateral isthmus] in-silico ablation, 28.7% after [CPVI+Posterior box] after in-silico ablation, 11.1% after [CPVI only], and 8.3% after [CPVI+CFAE] virtual ablation. Compared with virtual [CPVI only], virtual [CPVI+Posterior box+Anterior line] (p < 0.001) and virtual [CPVI+Roof line+Left lateral isthmus line] (p < 0.001) showed significantly higher virtual AF termination rates in both the V-ABL and Em-ABL groups. There was no significant difference in the in-silico ablation outcome between the V-ABL group and the Em-ABL group.

TABLE 2
www.frontiersin.org

Table 2. Virtual ablation outcome.

Comparison of V-ABL and Em-ABL

In the real ablation procedure, CPVI only was more commonly chosen in the Em-ABL group (30.9%) than in the V-ABL group (1.9%, p < 0.001; Table 3). Among patients in the V-ABL group, [CPVI+Roof+Left lateral isthmus] ablation was most commonly chosen (43.4%), although in-silico AF termination rate was highest in [CPVI+Posterior box+Anterior line] (Table 2). This was because, when AF was terminated in both lesion sets (32 cases), we chose the earliest terminating virtual lesion set (16 in [CPVI+Roof+Left lateral isthmus] and 11 in [CPVI+Posterior box+Anterior line]). In the Em-ABL group, [CPVI+Posterior box+Anterior line] was selected most often (36.4%), followed by [CPVI only] (30.9%; p = 0.687) (Table 3). Among patients in the Em-ABL group, the earliest terminating in-silico lesion set, which was blinded to the operator, was identical to the empirically chosen real ablation lesion set in 21.8% (12 of 55).

TABLE 3
www.frontiersin.org

Table 3. Clinical outcome.

The procedure-related acute outcomes are summarized in Table 3. There was no significant difference in procedure time (p = 0.289), ablation time (p = 0.988), complication rates (p = 0.359), and early recurrence rate within 3 months (p = 1.000) between the V-ABL group and the Em-ABL group.

Comparison of Rhythm Outcomes

During 12.6 ± 3.8 months of follow-up, the clinical recurrence rates were 14.0% in the V-ABL group and 20.9% in the Em-ABL group (p = 0.411). The maintenance rates of antiarrhythmic drugs were 49.1% in the V-ABL group and 39.5% in the Em-ABL group (p = 0.840). Kaplan-Meier analyses showed consistent findings in the overall patients (log-rank p = 0.732; Figure 3A) and after excluding patients who maintained antiarrhythmic drugs (log-rank p = 0.751; Figure 3C).

FIGURE 3
www.frontiersin.org

Figure 3. Kaplan-Meier curves. (A) Overall patients. (B) Patients maintaining antiarrhythmic drug use were excluded.

However, when we compared the clinical recurrence rate according to the ablation lesion set, the clinical recurrence rates were significantly lower in patients in the Em-ABL group who underwent [CPVI+Posterior box+Anterior line] ablation, which showed the most frequent termination and shortest AF maintenance duration in virtual ablation (log-rank p = 0.027; Figure 4).

FIGURE 4
www.frontiersin.org

Figure 4. Kaplan-Meier curves comparing the clinical recurrence rate of patients who underwent circumferential pulmonary vein isolation (CPVI) + posterior box lesion + anterior line with that of the other patients. (A) Overall patients. (B) V-ABL group. (C) Em-ABL group.

Discussion

Main Findings

In this prospective randomized trial, we investigated the feasibility of using a virtual ablation system in clinical practice and compared the clinical outcome of simulation-guided ablation (V-ABL) with empirical ablation (Em-ABL) in patients with PeAF. Patient-specific atrial geometry was used for finite element modeling, and five different virtual ablation patterns were tested for each patient. This is the first randomized study to validate the clinical outcomes of simulation-guided ablation, demonstrating the potential of in-silico computer models for planning individualized ablation strategies in clinical practice. We demonstrated that the most effective in-silico AF ablation lesion set showed the best clinical outcome in the Em-ABL group. The overall safety and efficacy of simulation-guided intervention was not inferior to empirical ablation of PeAF.

Unmet Clinical Needs in PeAF Ablation

Although pulmonary vein isolation is the cornerstone of catheter ablation for AF and effective in patients with paroxysmal AF, this was found to be insufficient to treat longstanding persistent AF. Even in the STAR AF II trial (Verma et al., 2015), which advocated CPVI alone for PeAF, a 41% recurrence rate within 18 months of CPVI procedure was not good enough for clinical practice. However, there is no consensus on further ablation after CPVI in PeAF (Verma et al., 2015). Linear ablation modifies the atrial substrates and compartmentalizes the atrium into smaller regions to reduce critical mass and AF maintenance (Pak et al., 2011). The question on what is the best linear lesion set for individual patients remains unanswered. CFAE ablation has been reported to be effective for substrate modification in PeAF by some investigators (Nademanee et al., 2004; Verma et al., 2010); however, this positive outcome was not reproduced in other investigations (Oral et al., 2008). Other various ablation strategies targeting rotor (Narayan et al., 2012), driver (Haissaguerre et al., 2014), dominant frequency (Atienza et al., 2014), or Shannon entropy (Ganesan et al., 2013) remain to be validated in patients with PeAF. There are potential explanations for the lack of a “one-fits-all ablation protocol.” First, PeAF is a progressive disease with various stages, and its nature is affected by the anatomy, histology, electrophysiology, and hemodynamic loading of individual patients. Second, the AF ablation procedure is operator dependent and the outcome is affected by the ablation lesion set; catheter stability; or ablation power, duration, and contact force. Third, it is not possible to test three different ablation designs in a patient because radiofrequency ablation causes irreversible tissue change. Therefore, in-silico AF ablation might be useful to test several different ablation lesion sets in a personalized modeling approach reflecting anatomy, histology, electrophysiology, and even genetic characteristics (Hwang et al., 2014).

Roles and Limitations of In-Silico AF Ablation

The advantage of computer simulation modeling is the potential for direct comparison of different ablation strategies for individual pre-procedural planning. As a result, modeling is emerging as a complementary approach to animal experiments and clinical trials in investigating more effective treatment (Winslow et al., 2012). Previously, Ruchat et al. (2007) compared the success rates in converting AF to sinus rhythm for 19 different linear ablation patterns by using a single identical human atrial model. There have been six AF simulation studies of catheter ablation thus far (Jacquemet, 2016), and we previously reported on multiple virtual ablation protocols in 20 different personalized human atrial models (Hwang et al., 2014). In this study, we conducted in-silico ablations prospectively in 108 personalized LA models reflecting individual anatomy; however, it was also a single-layer homogeneous modeling approach without integrating local histology. The main reason for our use of over-simplified simulation modeling was to reduce the computing time. Unless AF simulation and virtual intervention can generate the result within a few hours, it is not acceptable for clinical application during intervention. Although we already made the computing time 80 times faster by using a graphics processing unit system, further innovation is still needed for sophisticated realistic AF modeling reflecting detailed histology and electrophysiology.

Another discrepancy of in-silico AF ablation and real catheter ablation is completeness and maintenance of conduction block after linear ablation. Although [CPVI+Posterior box+Anterior linear ablation] showed the highest antifibrillatory effect in both virtual ablation and Em-ABL, it was difficult to achieve complete and permanent transmural conduction block by using endocardial radiofrequency energy delivery in the clinical setting with a beating heart and respiratory motion. A non-transmural lesion or ablation gap provides one of the main reasons for AF recurrence after an ablation procedure, discordant to in-silico AF ablation.

Study Limitations

There are limitations in our simulation modeling. The mechanism of AF in this particular computer model is based exclusively on multiple wavelet theories (Moe et al., 1964), and other well-described mechanisms like rotors can be masked in this model. Although we used patient-specific atrial geometry for the AF ablation simulation, individual electrophysiological and structural characteristics such as regional differences in action potential morphology, and fiber orientation were not included. Clinical studies showed that atrial wall thickness affected AF wave dynamics (Song et al., 2017), but the wall thickness was not considered in the present model. We also did not model the right atrium, because AF drivers have been known to exist mostly in the LA in humans (Sanders et al., 2005). We assume that the predictability of AF intervention by using simulation will be improved remarkably by integrating high-resolution cardiac imaging modalities reflecting the remodeled scar area, pathological electro-anatomical information acquired during procedure, and more realistic multilayer bi-atrial modeling study. Although our modeling protocol for virtual AF ablation may not reflect the mechanism of AF ablation precisely (whether virtual ablation stops an ongoing AF or whether it hinders the onset of new episode of AF, ectopic centers, rotor, or micro-reentry, etc.), it exhibits the antiarrhythmic effects of appropriate critical mass reduction, which is the main antiarrhythmic mechanism of linear ablation of AF or a surgical maze procedure. While CV may also change AF wave dynamics, we used a fixed CV based on real human patient data (Yonsei AF ablation cohort) (Park et al., 2014) and previous simulation studies (Hwang et al., 2014; Li et al., 2016). Because we tested five different virtual ablation protocols in every patient's LA modeling at the same CV condition, it did not affect the mechanism of critical mass reduction. For the model to show superiority to empirical ablation, individual patient's structural heterogeneity such as scar/fibrosis, fiber orientation, and wall thickness variation should be included in the model.

Conclusion

This is the first randomized multicenter study to demonstrate that simulation-guided ablation by using personalized models of LA was feasible, and the clinical outcome was not inferior to that of empirical ablation. Although there are limitations in applying virtual AF ablation to a personalized ablation strategy, our results showed that the virtual AF ablation system is capable of identifying the most effective ablation lesion set for individual patients. Advances in modeling technology might provide useful clinical insights for planning therapeutic interventions in the near future.

Author Contributions

Conception and design of the study: HP. Acquisition of data: JS, TK, BJ, SK, YSO, GN, YKO, SO, YK, and HP. Analysis of data: JS, MH, JSS, and BL. Interpretation of data: JS and HP. Drafting the work: JS and MH. All authors contributed in revising the work, approved the final version to be published, and agreed to be accountable for all aspects of the work.

Funding

This research was supported by grants (A085136) from the Korea Health 21 R&D Project, the Ministry of Health and Welfare, [NRF-2017R1A2B4003983] from the Basic Science Research Program run by the National Research Foundation of Korea (NRF), and Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2017R1D1A1B03030495).

Conflict of Interest Statement

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.

References

Atienza, F., Almendral, J., Ormaetxe, J. M., Moya, A., Martinez-Alday, J. D., Hernandez-Madrid, A., et al. (2014). Comparison of radiofrequency catheter ablation of drivers and circumferential pulmonary vein isolation in atrial fibrillation: a noninferiority randomized multicenter RADAR-AF trial. J. Am. Coll. Cardiol. 64, 2455–2467. doi: 10.1016/j.jacc.2014.09.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Ball, J., Carrington, M. J., McMurray, J. J., and Stewart, S. (2013). Atrial fibrillation: profile and burden of an evolving epidemic in the 21st century. Int. J. Cardiol. 167, 1807–1824. doi: 10.1016/j.ijcard.2012.12.093

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: a report of the Heart Rhythm Society (HRS) Task Force on Catheter and Surgical Ablation of Atrial Fibrillation. Developed in partnership with the European Heart Rhythm Association (EHRA), a registered branch of the European Society of Cardiology (ESC) and the European Cardiac Arrhythmia Society (ECAS); and in collaboration with the American College of Cardiology (ACC), American Heart Association (AHA), the Asia Pacific Heart Rhythm Society (APHRS), and the Society of Thoracic Surgeons (STS). Endorsed by the governing bodies of the American College of Cardiology Foundation, the American Heart Association, the European Cardiac Arrhythmia Society, the European Heart Rhythm Association, the Society of Thoracic Surgeons, the Asia Pacific Heart Rhythm Society, and the Heart Rhythm Society. Heart Rhythm 9, 632 e621–696 e621. doi: 10.1016/j.hrthm.2011.12.016

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. Physiol. 275(1 Pt 2), H301–H321.

PubMed Abstract | Google Scholar

Dewire, J., and Calkins, H. (2013). Update on atrial fibrillation catheter ablation technologies and techniques. Nat. Rev. Cardiol. 10, 599–612. doi: 10.1038/nrcardio.2013.121

PubMed Abstract | CrossRef Full Text | Google Scholar

Ganesan, A. N., Kuklik, P., Lau, D. H., Brooks, A. G., Baumert, M., Lim, W. W., et al. (2013). Bipolar electrogram shannon entropy at sites of rotational activation: implications for ablation of atrial fibrillation. Circ. Arrhythm. Electrophysiol. 6, 48–57. doi: 10.1161/CIRCEP.112.976654

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldstein, L. B., Bushnell, C. D., Adams, R. J., Appel, L. J., Braun, L. T., Chaturvedi, S., et al. (2011). Guidelines for the primary prevention of stroke: a guideline for healthcare professionals from the American Heart Association/American Stroke Association. Stroke 42, 517–584. doi: 10.1161/STR.0b013e3181fcb238

PubMed Abstract | CrossRef Full Text

Haissaguerre, M., Hocini, M., Denis, A., Shah, A. J., Komatsu, Y., Yamashita, S., et al. (2014). Driver domains in persistent atrial fibrillation. Circulation 130, 530–538. doi: 10.1161/CIRCULATIONAHA.113.005421

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, M., Kwon, S. S., Wi, J., Park, M., Lee, H. S., Park, J. S., et al. (2014). Virtual ablation for atrial fibrillation in personalized in-silico three-dimensional left atrial modeling: comparison with clinical catheter ablation. Prog. Biophys. Mol. Biol. 116, 40–47. doi: 10.1016/j.pbiomolbio.2014.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, M., Park, J., Lee, Y. S., Park, J. H., Choi, S. H., Shim, E. B., et al. (2015). Fibrillation number based on wavelength and critical mass in patients who underwent radiofrequency catheter ablation for atrial fibrillation. IEEE Trans. Biomed. Eng. 62, 673–679. doi: 10.1109/TBME.2014.2363669

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacquemet, V. (2016). Lessons from computer simulations of ablation of atrial fibrillation. J. Physiol. 594, 2417–2430. doi: 10.1113/JP271660

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Lim, B., Hwang, M., Song, J. S., Lee, Y. S., Joung, B., et al. (2016). The Spatiotemporal stability of dominant frequency sites in in-silico modeling of 3-dimensional left atrial mapping of atrial fibrillation. PLoS ONE 11:e0160017. doi: 10.1371/journal.pone.0160017

PubMed Abstract | CrossRef Full Text | Google Scholar

Moe, G. K., Rheinboldt, W. C., and Abildskov, J. A. (1964). A computer model of atrial fibrillation. Am. Heart J. 67, 200–220. doi: 10.1016/0002-8703(64)90371-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Mun, H. S., Joung, B., Shim, J., Hwang, H. J., Kim, J. Y., Lee, M. H., et al. (2012). Does additional linear ablation after circumferential pulmonary vein isolation improve clinical outcome in patients with paroxysmal atrial fibrillation? Prospective randomised study. Heart 98, 480–484. doi: 10.1136/heartjnl-2011-301107

PubMed Abstract | CrossRef Full Text | Google Scholar

Nademanee, K., McKenzie, J., Kosar, E., Schwab, M., Sunsaneewitayakul, B., Vasavakul, T., et al. (2004). A new approach for catheter ablation of atrial fibrillation: mapping of the electrophysiologic substrate. J. Am. Coll. Cardiol. 43, 2044–2053. doi: 10.1016/j.jacc.2003.12.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Narayan, S. M., Krummen, D. E., Shivkumar, K., Clopton, P., Rappel, W. J., and Miller, J. M. (2012). Treatment of atrial fibrillation by the ablation of localized sources: CONFIRM (Conventional Ablation for Atrial Fibrillation With or Without Focal Impulse and Rotor Modulation) trial. J. Am. Coll. Cardiol. 60, 628–636. doi: 10.1016/j.jacc.2012.05.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Oral, H., Chugh, A., Good, E., Crawford, T., Sarrazin, J. F., Kuhne, M., et al. (2008). Randomized evaluation of right atrial ablation after left atrial ablation of complex fractionated atrial electrograms for long-lasting persistent atrial fibrillation. Circ. Arrhythm. Electrophysiol. 1, 6–13. doi: 10.1161/CIRCEP.107.748780

PubMed Abstract | CrossRef Full Text | Google Scholar

Pak, H. N., Oh, Y. S., Lim, H. E., Kim, Y. H., and Hwang, C. (2011). Comparison of voltage map-guided left atrial anterior wall ablation versus left lateral mitral isthmus ablation in patients with persistent atrial fibrillation. Heart Rhythm 8, 199–206. doi: 10.1016/j.hrthm.2010.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, J., Joung, B., Uhm, J. S., Young Shim, C., Hwang, C., Hyoung Lee, M., et al. (2014). High left atrial pressures are associated with advanced electroanatomical remodeling of left atrium and independent predictors for clinical recurrence of atrial fibrillation after catheter ablation. Heart Rhythm 11, 953–960. doi: 10.1016/j.hrthm.2014.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruchat, P., Virag, N., Dang, L., Schlaepfer, J., Pruvot, E., and Kappenberger, L. (2007). A biophysical model of atrial fibrillation ablation: what can a surgeon learn from a computer model? Europace (Suppl. 6), 6, vi71–vi76. doi: 10.1093/europace/eum209

CrossRef Full Text | Google Scholar

Sanders, P., Berenfeld, O., Hocini, M., Jais, P., Vaidyanathan, R., Hsu, L. F., et al. (2005). Spectral analysis identifies sites of high-frequency activity maintaining atrial fibrillation in humans. Circulation 112, 789–797. doi: 10.1161/CIRCULATIONAHA.104.517011

PubMed Abstract | CrossRef Full Text | Google Scholar

Shim, J., Joung, B., Park, J. H., Uhm, J. S., Lee, M. H., and Pak, H. N. (2013). Long duration of radiofrequency energy delivery is an independent predictor of clinical recurrence after catheter ablation of atrial fibrillation: over 500 cases experience. Int. J. Cardiol. 167, 2667–2672. doi: 10.1016/j.ijcard.2012.06.120

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, J. S., Lee, Y. S., Hwang, M., Lee, J. K., Li, C., Joung, B., et al. (2016). Spatial reproducibility of complex fractionated atrial electrogram depending on the direction and configuration of bipolar electrodes: an in-silico modeling study. Korean J. Physiol. Pharmacol. 20, 507–514. doi: 10.4196/kjpp.2016.20.5.507

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, J. S., Wi, J., Lee, H. J., Hwang, M., Lim, B., Kim, T. H., et al. (2017). Role of atrial wall thickness in wave-dynamics of atrial fibrillation. PLoS ONE 12:e0182174. doi: 10.1371/journal.pone.0182174

PubMed Abstract | CrossRef Full Text | Google Scholar

Trayanova, N. A. (2014). Mathematical approaches to understanding and imaging atrial fibrillation: significance for mechanisms and management. Circ. Res. 114, 1516–1531. doi: 10.1161/CIRCRESAHA.114.302240

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Verma, A., Mantovan, R., Macle, L., De Martino, G., Chen, J., Morillo, C. A., et al. (2010). Substrate and Trigger Ablation for Reduction of Atrial Fibrillation (STAR AF): a randomized, multicentre, international trial. Eur. Heart J. 31, 1344–1356. doi: 10.1093/eurheartj/ehq041

PubMed Abstract | CrossRef Full Text | Google Scholar

Winslow, R. L., Trayanova, N., Geman, D., and Miller, M. I. (2012). Computational medicine: translating models to clinical care. Sci. Transl. Med. 4, 158rv111. doi: 10.1126/scitranslmed.3003528

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Kharche, S. R., Hansen, B. J., Csepe, T. A., Wang, Y., Stiles, M. K., et al. (2015). Optimization of catheter ablation of atrial fibrillation: insights gained from clinically-derived computer models. Int. J. Mol. Sci. 16, 10834–10854. doi: 10.3390/ijms160510834

PubMed Abstract | CrossRef Full Text | Google Scholar

Zozor, S., Blanc, O., Jacquemet, V., Virag, N., Vesin, J. M., Pruvot, E., et al. (2003). A numerical scheme for modeling wavefront propagation on a monolayer of arbitrary geometry. IEEE Trans. Biomed. Eng. 50, 412–420. doi: 10.1109/TBME.2003.809505

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: atrial fibrillation, catheter ablation, virtual ablation, in-silico modeling, recurrence

Citation: Shim J, Hwang M, Song J-S, Lim B, Kim T-H, Joung B, Kim S-H, Oh Y-S, Nam G-B, On YK, Oh S, Kim Y-H and Pak H-N (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

Received: 25 May 2017; Accepted: 27 September 2017;
Published: 11 October 2017.

Edited by:

Eun Bo Shim, Kangwon National University, South Korea

Reviewed by:

Jae Boum Youm, Inje University, South Korea
Seiryo Sugiura, The University of Tokyo, Japan

Copyright © 2017 Shim, Hwang, Song, Lim, Kim, Joung, Kim, Oh, Nam, On, Oh, Kim 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) or licensor 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

These authors have contributed equally to this work.

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.