- 1Lab of Computing, Medical Informatics and Biomedical Imaging Technologies, School of Medicine, Aristotle University of Thessaloniki, Thessaloniki, Greece
- 23rd Cardiology Department, Aristotle University of Thessaloniki, Thessaloniki, Greece
- 3Department of Industrial Engineering and Management Sciences, Northwestern University, Evanston, IL, United States
The remarkable advances in high-performance computing and the resulting increase of the computational power have the potential to leverage computational cardiology toward improving our understanding of the pathophysiological mechanisms of arrhythmias, such as Atrial Fibrillation (AF). In AF, a complex interaction between various triggers and the atrial substrate is considered to be the leading cause of AF initiation and perpetuation. In electrocardiography (ECG), P-wave is supposed to reflect atrial depolarization. It has been found that even during sinus rhythm (SR), multiple P-wave morphologies are present in AF patients with a history of AF, suggesting a higher dispersion of the conduction route in this population. In this scoping review, we focused on the mechanisms which modify the electrical substrate of the atria in AF patients, while investigating the existence of computational models that simulate the propagation of the electrical signal through different routes. The adopted review methodology is based on a structured analytical framework which includes the extraction of the keywords based on an initial limited bibliographic search, the extensive literature search and finally the identification of relevant articles based on the reference list of the studies. The leading mechanisms identified were classified according to their scale, spanning from mechanisms in the cell, tissue or organ level, and the produced outputs. The computational modeling approaches for each of the factors that influence the initiation and the perpetuation of AF are presented here to provide a clear overview of the existing literature. Several levels of categorization were adopted while the studies which aim to translate their findings to ECG phenotyping are highlighted. The results denote the availability of multiple models, which are appropriate under specific conditions. However, the consideration of complex scenarios taking into account multiple spatiotemporal scales, personalization of electrophysiological and anatomical models and the reproducibility in terms of ECG phenotyping has only partially been tackled so far.
Introduction
Over a century ago, Atrial Fibrillation (AF) had been recognized as the most common arrhythmia in adults (Nishida and Nattel, 2014). It is characterized by chaotic atrial activation leading to impaired atrial myocardial function. It is easily recognized on a surface electrocardiogram by lack of atrial depolarization represented by P-wave, quivering isoelectric line and irregular ventricular activation represented by QRS-complexes, which lead to contractile dysfunction (Pellman and Sheikh, 2015). In the case that AF converts to normal sinus rhythm (SR) within 7 days, it is classified as “paroxysmal,” while in the case of it lasting more than 7 days or more than a year it is classified as “persistent” or “long-standing persistent,” respectively (Camm et al., 2012a; Kirchhof et al., 2016).
Atrial fibrillation is associated with increased morbidity and mortality as it is widely recognized as a risk factor for embolic stroke and heart failure (HF) exacerbation (Camm et al., 2012b). On the other hand, numerous risk factors have been correlated with the development of AF. Aging doubles AF risk per decade, whereas gender influences the incidence of AF with males having 1.5 fold risk to develop AF (Andrade et al., 2014). Most cases of persistent and permanent AF are related to hypertensive, valvular, ischaemic or other types of structural heart disease, while lone AF represents only 15% of the cases (Markides and Schilling, 2003).
Initiation and maintenance of the abnormal rhythm in both lone AF and AF secondary to structural heart disease are supposed to require pathophysiological remodeling of the atria (Pellman and Sheikh, 2015). Remodeling can be grouped into three categories that include:
(i) electrical remodeling, including modulation of L-type Ca2+ current, various K+ currents, and gap junction function,
(ii) structural remodeling, including changes in tissues properties, size, and ultrastructure, and
(iii) autonomic remodeling, including altered sympathovagal activity, and hyperinnervation (Pellman and Sheikh, 2015).
Electrical, structural, and autonomic remodeling all contribute to creating an AF-prone substrate which can produce AF-associated electrical phenomena including a rapidly firing focus, complex multiple reentrant circuits, or rotors (Nattel, 2002). Since the 1960s, the most popular theory has held that AF consists of multiple wavelets of functional re-entry (Moe et al., 1964). Conditions that increase the atrial size or decrease the wavelength (by decreasing the conduction velocity or refractory period or both) permit multiple wavelets and promote AF (Rensma et al., 1988).
Multiple factors have been found to relate and contribute to AF initiation or perpetuation mechanisms. Ionic currents and gap junction function are considered to be fundamental parts of electrical remodeling occurring in atrial cardiomyocytes, which leads to reduced action potential (AP) duration and refractoriness and affects conduction velocity and wavelength, both known determinants of AF initiation and perpetuation (Gaspo et al., 1997; Pellman and Sheikh, 2015). Atrial enlargement has been associated with AF many decades ago (Henry et al., 1976). Moreover, scarring, fibrosis, and increased atrial size observed in HF patients, are related to conduction delay and refractory period prolongation (Sanders et al., 2003). The autonomic nervous system (ANS) exerts significant control over cardiac electrophysiology, while it has been proposed that patterns of baseline autonomic nerve activity are essential in the development of pacing-induced sustained AF (Shen et al., 2011). There is a suggestion that the increased sympathetic activity leads to heterogeneous changes in atrial refractoriness, which in turn favors reentrant waves that generally contribute to the maintenance of AF (Olgin et al., 1998). On the other hand, ablation of the ganglionated plexi can also improve long-term AF symptoms (Katritsis et al., 2013). Additionally, the muscular sleeves of the pulmonary veins (PVs) have been identified as a source of tachyarrhythmias and atrial premature beats that could trigger paroxysms of AF (Haïssaguerre et al., 1998). These findings lead to a unifying theory that focal tachycardias (mostly originating from the PVs) promote atrial remodeling and they are required to trigger and maintain a substrate capable of multiple wavelet reentry (Veenhuyzen et al., 2004). AF causes electrophysiological changes in the atrial myocardium itself, which might explain the progressive nature of the arrhythmia. In a landmark study, the perpetuation of AF was accompanied by a shortening of the atrial refractory period inducing electrophysiologic changes that promote further AF. These include electrical, contractile and structural modifications to the atria that have collectively become known as atrial remodeling and the authors concluded that “AF begets AF” (Wijffels et al., 1995). However, most of the findings are based on electrophysiological studies making the reproducibility of the experiments a challenging case.
The interplay of different factors contributing to AF and the inherent complexity of the biological systems, necessitates the development of computational models across different levels, from the cell, tissue, and patient level. Computational modeling can help study and shed light on the mechanisms of AF, and thus new therapeutic approaches can be developed and tested before applying them to patients (Jacquemet, 2016). The main advantage of computational modeling is that researchers have access to all levels of interest and thus they can perform experiments in controlled and repeatable conditions, which is a limitation of experimental and clinical research. The increased evidence on the importance and the effectiveness of computational modeling in several fields of science, including cardiac electrophysiology, in association with the advances in the computational speed (Virag et al., 2012) and the availability of such infrastructure, were the primary catalysts on the increased interest in the field of computational modeling during the recent years. Overall, computational modeling can link the phenomenological organ level findings, such as findings in the electrocardiography (ECG), and the actual mechanisms at smaller scales.
As regards the AF patients, the analysis of P-wave during SR is of clinical value, since P-wave morphology is affected by the site of origin of the SR beat and the atrial conduction routes (Platonov, 2012). The analysis of P-waves can reveal information related to the prediction of AF initiation (Martínez et al., 2014a) or the success of the pulmonary vein isolation (Huo et al., 2015). Also, the analysis of the P-wave morphology variability, on a beat-to-beat basis, also revealed differences between healthy subjects and AF patients (Censi et al., 2016; Filos et al., 2017). However, the direct correlation of those ECG findings with the underlying substrate modifications is not well understood. The motivation behind this research was to find what mechanisms may explain the difference in the percentage of the primary/secondary P-wave morphologies observed in ECG signals of paroxysmal AF vs. normal subjects, and more importantly whether these mechanisms were studied, quantified and reproduced via simulation studies. To serve this need, we designed a scoping review, (a) to map the existing literature in the field of atrial models and their multifaceted properties/components related to AF, and (b) to synthesize this knowledge toward explaining the link between observable AF-related P-wave morphologies and computational atrial models of AF. Finally, this work attempts to identify the research gaps and to make recommendations for future research.
Methodology
We conducted a scoping review adopting the methodology proposed by Arksey and O’Malley (2005) and Peters et al. (2015). The aim of scoping reviews is to map existing literature in a particular field as well as to synthesize the knowledge. The definition proposed by Colquhoun et al. (2014) states that “A scoping review or scoping study is a form of knowledge synthesis that addresses an exploratory research question aimed at mapping key concepts, types of evidence, and gaps in research related to a defined area or field by systematically searching, selecting, and synthesizing existing knowledge.”
Objective
We focused on the identification of the mechanisms, either on cell, tissue or organ level, which can describe the variation of the P-wave morphologies in AF patients during SR while we tried to investigate the existence of computational models which can reproduce those mechanisms. The following research questions were posed:
(1) How was AF modeled?
(2) Which computational models were proposed for the description of patients prone to AF during SR?
(3) How was the initiation of AF episodes modeled?
(4) Were those models able to reproduce the multiple P-wave morphologies that appear in AF patients during SR?
The primary outcomes that are of interest in this scoping review are the characteristics of the AF substrate which describe atrial remodeling occurring as a result of AF. The comprehension of such mechanisms may lead to get more insights on the prediction of the AF onset and thus facilitate the management of patients in terms of the prevention and adaptation of the therapeutic approach.
Inclusion Criteria
We applied the PICO (Population, Intervention, Comparator, and Outcome) eligibility criteria. Regarding the population, we considered all AF studies without distinction between human and animal models. All types of AF were considered, including Asymptomatic, Paroxysmal, Persistent, or Permanent. As for the intervention, any level of computational models that enabled the simulation of the atrial activation was eligible for inclusion. Models focusing on the cell or tissue level, as well as on the organ level were examined.
Studies comparing simulated and real ECG phenotyping or comparison of different computational models were included in the literature review as well as studies without a comparator, such as theoretical models. Finally, any publication, such as articles, conference proceedings, editorials and chapters in textbooks, published after 1998 until 1 of March 2018, written in English was included in the review, whereas all review articles were excluded.
Search Strategy
A two steps search strategy was adopted to map the existing research on this topic. Two databases, which are extensively used in biomedical sciences, were searched (PubMed and Scopus). Initially, a limited search in both databases was undertaken followed by an analysis of the text words contained in the title and the abstract as well as on the indexed terms of the article. The second search step included the extensive literature search based on the keywords extracted in the first step and the manual inspection of the studies.
Regarding the text analysis made during the first step, review articles, published during the last 10 years, in which the words “cardiac” and “simulation” were included in their title or abstract were selected, as they are considered to cover the field of the current review. As a result, 115 and 163 articles were found in PubMed and Scopus, respectively. After removing the duplicates, 216 review papers were screened, whereas only 99 review articles were found to be relevant to our study. Subsequently, a text analysis of the titles and the abstract was conducted. A table of the most frequent words and phrases was produced, while equivalent words were identified (e.g., modeling and modeling). The word cloud depicted in Figure 1 presents in a visual way the most frequent words. The most common phrases were found to be “computation model*,” “mathematical model*,” or “computer model*,” which describe the same process. On the other hand, the word “model*” and “simulat*” appeared in the same abstract, but not necessarily accompanied, to describe the same approach as the phrases above. The asterisk (*) following a word describes all of the words starting with the phrase before it [e.g., model(s), modeling are described as model*]. Based on these findings we formulated the first two blocks of Figure 2, which are connected to an OR logical gate. The terms “cardiac,” “atrial,” and “heart” were extracted as the most frequent words, that are related to the description of the model or the tissue of interest. Finally, the term “Atrial Fibrillation” was not extracted by the text analysis but it was used since each article was focused on AF.
Figure 1. Word cloud of the most frequent words, where the highest the font the higher the frequency. For the purposes of the visualization the terms model, models, modeling, and modeling were grouped under the term “model.” As depicted, the most used words are model, cardiac, computational, simulation, heart, and with frequencies 16.6, 11.16, 7.2, 6.5, and 5.7%, respectively.
Figure 2. Visual representation of the final query using logic gates. The final query is described as: {[(“computational model*”[Title/Abstract] OR “mathematical model*”[Title/Abstract] OR “computer model*”[Title/Abstract]) OR (simulat*[Title/Abstract] AND model*[Title/Abstract])] AND (“cardiac”[Title/Abstract] OR “atrial”[Title/Abstract] OR “heart”[Title/Abstract])} AND “atrial fibrillation”[Text Word].
These findings led to the creation of the search query that was used in the second step of the literature review (Figure 2). In this step, the words and phrases identified were searched in the abstract, the title or the indexed keywords. The whole literature search, its evaluation, and categorization were conducted independently by two reviewers (DF and DT). After the completion of the evaluation, in case of disagreement, the two reviewers worked together to reach a consensus.
Results
A total number of 721 articles (excluding duplicates) was raised. Titles and abstracts were screened to check whether they are relevant to the research question or not. Approximately half of them (377) were eligible for full-text review. Seventy-six articles were excluded, resulting in a total of 301 articles which were finally included in the review (Figure 3). The observation of Figure 4 reveals a continuous increase in the number of articles published over time.
Figure 4. Evolution of the number of publications in the field of AF modeling from 1998 till 1st March 2018. A parabolic line was found to fit better the data depicted in the bars ().
Modeling Components
The interplay of multiple factors can lead to the initiation and perpetuation of the arrhythmia (Kirchhof et al., 2016). In this section, the existing literature in the field of atrial models is presented while a categorization of the articles was made based on the factors that are addressed in each simulation study.
Model Experimental Context
Animal studies
While in the first years, the ratio between the studies focusing on humans or animals was approximately 1:1, during the last decade, an increased interest in human studies is observed (Figure 5). The animal studies can improve our understanding of the pathophysiological properties of AF while their findings can be useful in human studies.
Figure 5. Comparison of the studies based on human or animal models for AF for each year. For 2018, until 1st of March 2018, the studies targeting humans and animals were 3 and 0, respectively. To avoid any false perception those data have not been visualized.
A detailed model of the canine atrial cell is presented by Ramirez et al. (2000) which is applicable in both normal and electrically remodeled cellular substrate. The following years, additional animal models were proposed able to describe the sheep (Goodman et al., 2005; Butters et al., 2013), rat (Majumder et al., 2016) or other animals’ atrial myocardium (Aslanidi et al., 2009b). A novel computational model of the canine left atrium (LA) and PVs is proposed by Cherry et al. (2007) and Aslanidi et al. (2011a) through which an attempt was made to understand the arrhythmogenicity of PV and how re-entries and their geometry can promote their excitation. The differences between the right atrium (RA) and LA in terms of the electrophysiological properties in the animals are described in Aslanidi et al. (2009a) and in Xia et al. (2010).
Aslanidi et al. (2008) studied the effects of Ca2+ concentration in the cell of the RA in terms of AF initiation. According to the simulation performed in a cell and tissue level, no crucial effect of intracellular Ca2+ concentration on reentries has been observed as its results are reversible. Furthermore, Kneller et al. (2005) studied the influence of antiarrhythmic drugs in the canine atrium, and found that the inhibition of pure Na current during AF can terminate the arrhythmia.
Other mechanisms were also addressed in animal models before being incorporated in the human ones, such as the effect of ANS (Kneller et al., 2002) and of the antiarrhythmic drugs (Comtois et al., 2008; Almquist et al., 2010; Aguilar-Shardonofsky et al., 2012; Colman et al., 2014; Aguilar et al., 2015; Varela et al., 2016).
Vigmond et al. (2001) presented a morphologically realistic atrial model where all the major anatomical structures of the atrium, including fiber orientation, muscle structures such as crista terminalis and pectinate muscles and the orifices of the veins and valves, were considered and scaled to a canine atrium level. The resulting model did not consist of an exact representation of the atria. However, it reflected the interrelations between atrial anatomical structures. Zhao et al. (2012) proposed a detailed structural model of the sheep atria which was based on ex vivo analysis of the atrial chambers. The details of the fiber orientations and myofiber architecture were incorporated into an anatomical model. The simulation performed in normal and electrically remodeled conditions, confirm the unique roles of crista terminalis, pectinate muscles and Bachman’s bundle on the activation time as well as the differences on the electrical propagation through the posterior LA.
In addition to animal models, many studies are taking into account findings from experiments performed on animals, to modify the human model. The goal of those studies extends from studying the effects of drug therapy (Gomez et al., 2005; Ehrlich et al., 2008) to impact of gene mutations (Hancox et al., 2014; Syeda et al., 2016).
Evaluation of model output
Computational modeling is by nature an approximation of reality and several approaches have been adopted to evaluate their output. However, most of the studies employed in silico experiments, rather than a direct evaluation of the simulation outputs with real recordings. On the other hand, some of the articles included a comparison between the in silico and the in vivo findings, based on electrophysiological findings (Courtemanche et al., 1998, 1999; Ramirez et al., 2000; Wettwer et al., 2004; Kneller et al., 2005; Kuijpers et al., 2007; Comtois et al., 2008; Lombardo et al., 2016). Furthermore, the reproducibility of the simulated signals, electrograms or ECG, with the real recordings is discussed in the context of solving the forward (Ogawa et al., 2007; Burdumy et al., 2012) or inverse problem (Pedron-Torrecilla et al., 2016). Finally, the need for personalized anatomical and electrical models is highlighted toward the highest reproducibility of the results (Krueger et al., 2013c; Roney et al., 2016).
Model Electrophysiological Components at Different Scales
Electrophysiological cell models in humans
The core of each study is the electrophysiological model which describes the cellular ionic properties and transmembrane currents. A plethora of different approaches has been proposed to simulate electrical activity in human atria.
Initially, a common approach was the use of models describing the ventricular activity, such as those of Beeler and Reuter (1977) or Luo and Rudy (1991), to simulate atrial excitation (Virag et al., 2002). On the other hand, membrane kinetics describing the atrial myocardium in animals, such as in canines, have been applied in the human atrium (Ruchat et al., 2007b, 2007c), hypothesizing similar ionic properties.
The CRN (Courtemanche et al., 1998) and Nygren et al. (1998) models are the first attempts to describe the ionic mechanism of the atrial myocytes in humans. The CRN model was based on adaptations made on a model describing ventricular myocytes. The Maleckar et al. (2009) model refined the description of Nygren model’s K current, while the Koivumäki et al. (2011) extended the Nygren model in order to improve the description of the Ca2+ dynamics. The Grandi et al. (2011) model focuses on patients with chronic AF (cAF) and in SR. Figure 6 depicts the number of articles focusing on the human atria that where based on each of the most frequently used electrophysiological models. The most frequently used model is the one proposed by Courtemanche et al. (1998). Compared to the first decade (1998–2007), the number of articles for AF which adopt ventricular models [such as that of Luo and Rudy (1991)] decreased to half in the second decade (2008–2018). On the other side, the newly proposed models by Koivumäki et al. (2011), Grandi et al. (2011), and Colman et al. (2016) seem to gain the research community’s acceptance.
Cherry and Evans (2008), as well as Krueger et al. (2013a) focused on the comparison of electrophysiological models. In particular, Krueger et al. (2013a) performed a comparison between CRN, Nygren, and Maleckar electrophysiological models adjusted to describe 13 different atrial regions heterogeneities, including electrical remodeling. A detailed and personalized MRI-based anatomical model, including fiber orientation detail and thorax model, was estimated for each of the eight patients to reproduce P and Ta-wave accurately. According to the results, the CRN model simulates better action potential duration (APD) in each of the atrial regions compared with the remaining two models. Furthermore, the simulated P-waves and Ta-waves morphology calculated using the CRN models fitted well with the experimental recordings in both physiological and AF-remodeled conditions. Wilhelms et al. (2013) compared the five most frequently used electrophysiological models of human atria and presented the main strengths and weaknesses of each of the models. Lombardo et al. (2016) performed a comparison between two electrophysiological models, a detailed and a simpler one. The results showed that both models reproduce the clinical recordings and the spiral waves dynamics were similar and thus, in the case that computationa cost is crucial, simple models could be used to model arrhythmias in spatially extended domains.
On the other hand, some articles (Ciaccio et al., 2017; Lin et al., 2017) adopted the use of a cellular automaton model (CA) like the one proposed by Moe et al. (1964), which was the first CA model for AF. In general, CA models describe the depolarization and repolarization of a cell based on rules of its present state as well as on its neighbors while more detailed CA models have been proposed recently (Manani et al., 2016). The main property of a CA model is that it is suitable for observation of macroscopic properties of the electrical propagation since it is a simplified model of cardiac electrophysiology and the required computational load is low.
Genes and mutations
Several studies reveal the occurrence of AF on a family basis which implies a genetic cause of this arrhythmia. Genetic mutations have been found to lead to modification of ionic channel functionality. The first study identified is the one of Seemann et al. (2004) where a 2D model of the RA was used to investigate the effect of a mutation in the KCNQ1 gene on the initiation of AF. In the following years, the effect of additional gene mutations or gene expressions on the initiation of AF was studied (Table 1). In all these studies, the parameters of the atrial models have been modified accordingly, to reproduce the findings from the electrophysiological studies. Most of the studies addressed the effect on the slow delayed outward rectifier potassium channel () (Seemann et al., 2004; Ehrlich et al., 2005; Kharche et al., 2012a) as well as on other potassium (Seemann et al., 2009) or sodium channels (Ziyadeh-Isleem et al., 2014). While most articles focus on the effects of gene expression on a cell or 2D tissue model, Whittaker et al. (2017) studied the stability of the re-entrant waves using a realistic atrial model. Finally, Ni et al. (2017a) examined six observed mutations of KCNA5 regarding their role in the electromechanical function of the atrium. Half of those mutations lead to gain-of-function resulting in a worsening of the contractile function of the atrium while the rest lead to loss-of-function of that mediated the positive inotropic effects.
Modeling drug effect on atrial cells and AF
Several types of antiarrhythmic drugs have been analyzed in silico regarding their performance on the termination of AF (Ehrlich et al., 2008; Sanchez et al., 2017), under SR and AF (Aslanidi et al., 2012a). Table 2 summarizes the studies focusing on each anti-arrhythmic drug category (class I to v).
The first study identified in this context, Namba et al. (1999) investigated the effects of a specific sodium channel blocker drug on the stability versus meandering of spiral waves in a canine atrial model.
One of the major challenges of antiarrhythmic drug therapy is the minimization of side effects, like QT prolongation and ventricular pro-arrhythmic conditions (Aguilar-Shardonofsky et al., 2012; Aguilar et al., 2015). A key concept to address these challenges is the design of drugs with atrial selectivity. Since ultrarapid delayed rectifying potassium current () is present in atrial and absent from ventricular myocardium, it is a justified target for AF antiarrhythmic therapy (Almquist et al., 2010). However, the effect of such a () blocker may be pro-arrhythmic during sinus rhythm and anti-arrhythmic during AF (Law et al., 2010). Moreover, the effect of selective (fast and slow onset) blockers of depends on the level of electrical remodeling (Wettwer et al., 2004; Tsujimae et al., 2008). Although is diminished in cAF, the blockers show prominent significance in cAF treatment, as investigated in silico in a number of theoretical drugs in both SR and cAF conditions (Ellinwood et al., 2017a, 2017b).
Along this line, Ni et al. (2017b) studied in silico the effects of the multi-channel blockers on atrial and ventricular activation, and found that the synergistic use of a specific and potentially selective potassium channel blocker (acacetin) and a blocker seems to have no significant impact on QT prolongation or on ventricular activation. As suggested, the synergetic antiarrhythmic treatment leads to more effective termination of atrial reentries (Aguilar et al., 2015; Ni et al., 2016), while ventricular side effects in the ventricles are limited. Furthermore, the relationship between drug concentration and AF termination is also considered in other studies (Duarte et al., 2013; Tobón et al., 2014a, 2017).
Colman et al. (2014) and Varela et al. (2016) studied the effectiveness of the antiarrhythmic therapy on termination of the wave breaks near the PVs. The latter study discusses the effect of the multi-channel class III antiarrhythmic drugs in combination with atrial heterogeneities on the termination of AF, using detailed 3D anatomical models of human atria. Additionally, the role of the antiarrhythmic drugs in the prevention of AF initiation due to bradycardia (Cacciani and Zaniboni, 2015) or how isoproterenol can induce AF (Krummen et al., 2012) were also studied.
A fundamental approach to AF treatment is β-blockers usage to control ventricular response. The chronic use of such treatment can lead to drug-induced remodeling. Marshall et al. (2012) studied these pharmacological remodeling effects of β-blockers on specific potassium currents in the human atrium. Kharche et al. (2014b) discussed this pharmaceutical remodeling and observed that chronic use of β-blockers suppresses AF by APD and effective refractory period (ERP) prolongation, with potential anti-arrhythmic consequences.
Autonomic nervous system and the role of ganglia
The effect of vagal tone on atrial arrhythmogenesis has been reported in several studies (Vigmond et al., 2004). Modifying the ion channel conductance, Ashihara et al. (2002) simulated the two branches of the ANS and used the Luo-Rudy-1 model to examine why vagally mediated AF lasts more than sympathetic mediated AF. Interestingly, it was found that sympathetic tone promotes spiral waves and restrains their breakup. In contrast, vagal tone promotes spiral wave breakup. A significant limitation of this study is that the modifications on the ion channels were speculative as there was insufficient information from clinical data.
A significant number of studies incorporates, in the electrophysiological models, the effect of Acetylcholine to model the vagal activation and study its effect on AF (Comtois and Nattel, 2011; Voigt et al., 2013). Acetylcholine (Ach) is a neurotransmitter used in vagal action studying, while its concentration affects the APD (Vigmond et al., 2004). Experimental findings reveal increased vagal activity during hemodialysis (HD) sessions suggesting that enhanced acetylcholine concentration can lead to ERP shortening and depolarization prolongation (Vincenti et al., 2014). Kneller et al. (2002) used the RNC canine atrial model, and the Ach effect was incorporated based on a novel description of the Ach-induced current () that is included in all the studies which consider the vagal influence on AF. Hwang et al. (2016) considered three different approaches for the distribution of the ACh concentration on the atrium: (1) random distribution, (2) ACh concentration within the four ganglionated plexi (GP) areas, and (3) using the octopus hypothesis where eight nerves with gradient ACh concentration originate from the GPs and spread through the atrium. The octopus hypothesis was adopted in order to examine the effect of both branches of ANS (sympathetic and parasympathetic nervous system) (Hwang et al., 2017). The parasympathetic system was modeled as an increase of the in the ganglionated plexi and their nerves in the LA 3D model while the sympathetic system was modeled as an increase of the L-type calcium current. ANS stimulation can induce triggering activity and PV automaticity, leading to local re-entries in the LA-PV junction.
AF Remodeling
Atrial fibrillation is a progressive disease, where at the first stages the episodes are short and infrequent (paroxysmal AF), and later on, they are more frequent and more prolonged (persistent or permanent AF). As AF progresses to more stable types, different parameters of the electrophysiological models are modified to reflect the AF-induced electrical remodeling. These changes on ion currents and conductance are based on findings from experimental studies (Courtemanche et al., 1999; Zhang et al., 2005). The majority of the 301 studies included in this review focused on more stable types of AF while only about 10% of them directly referred to Paroxysmal AF (Ashihara et al., 2002; Gong et al., 2007; Severi et al., 2010; Duarte et al., 2013; Calvo et al., 2014; Vincenti et al., 2014; Voigt et al., 2014; Berenfeld, 2016).
The electrically remodeled substrate favors the initiation and the perpetuation of AF. APD shortening and decreased conduction velocity lead to reduced wavelength. However, the mechanisms for AF initiation and perpetuation may differ (Chang and Trayanova, 2016) which is the subject of several studies (Krogh-Madsen et al., 2012; Colman et al., 2013; Koivumäki et al., 2014; Lee et al., 2016). Furthermore, studies included the extent of the electrical remodeling in each atrium (Luca et al., 2015), concerning its influence on the duration of AF episodes. Among the articles included in this review, 44 articles were found to study AF initiation and 25 AF termination, while the majority of the articles examine the AF perpetuation mechanisms. Finally, 37 articles investigate atrial activation under SR (Weber et al., 2009a, 2009b; Krueger et al., 2013c, 2013d; Voigt et al., 2013).
Apart from electrical remodeling, structural remodeling has also been found to influence the susceptibility to AF. It can be the result of structural heart disease or other conditions such as hypertension, while AF itself can also modify the atrial substrate. As AF progresses, apart from the electrical remodeling, structural remodeling is also present (Wijffels et al., 1995) and interestingly larger atrial size favors reentrant circuits, possibly due to more area available for rotor formation (Zou et al., 2005). The structural remodeling is modeled by considering fibrosis (including collagenous septa, remodeled gap junctions and proliferation of myofibroblast) (McDowell et al., 2013), scars as a result of previous ablation procedures (Gonzales et al., 2014), endo/epicardium dissociation (Gharaviri et al., 2017) and absence of t-tubules (Li et al., 2012). The role of structural remodeling in AF initiation and perpetuation has been thoroughly examined. Zhao et al. (2013b) confirmed that structural remodeling facilitates reentries and multiple wavelets, increasing AF susceptibility. Zahid et al. (2016) described the effect of atrial fibrosis on the appearance of reentrant drivers in patients with persistent AF, and by the use of patient-specific models of the atria, demonstrated that the reentrant activity in fibrotic zones perpetuates AF. The role of the fibroblasts, which can serve both as triggers and substrate to cardiac arrhythmias is described by Koivumaki et al. (2014).
Furthermore, Aguilar et al. (2014) investigated fibroblasts involvement in the maintenance of AF, in HF patients, where they are usually activated. Based on electrophysiological studies in dogs, it was found that a novel fibroblast K+ current must be taken into account in the mathematical models to describe the effects of such type of structural remodeling. Finally, Jacquemet and Henriquez (2009) described the impact of microfibrosis progression on electrograms’ fractionation as structural remodeling increases. However, fibrosis distribution can remarkably alter the activation pathways, so there is a great need for accurate atrial substrate mapping using advanced mapping techniques, such as Late Gadolinium (Zahid et al., 2016) and contrast-enhanced mapping (Zhao et al., 2017).
Integrating With Structure, Geometry, and Anatomy
Detailed anatomical structures
Atrial tissue presents a complex morphology both in terms of structural and electrical heterogeneities. Fiber orientation constitutes a key aspect of atrial anatomy while electrical heterogeneities (such as the fast conduction systems of Bachman bundle, Pectinate muscles, and Crista Terminalis) alter active potential propagation and the way this can influence the initiation and maintenance of AF (Jacquemet and Henriquez, 2009). Vigmond et al. (2001) presented a simple geometrical representation of atrial anatomy toward assessing of the role of the various anatomical structures on the initiation of reentries. In that study, all major anatomical structures of the atrium were considered, including fiber orientation, muscle structures such as crista terminalis and pectinate muscles as well as the orifices of the veins and valves. The outcome model did not constitute an exact representation of the atria, although it reflected the interrelations between the anatomical structures of the atria.
In most studies, fiber orientation was linear whereas, in another study, a sigmoid representation of fiber orientation was adopted (Ashihara et al., 2002). However, in 3D atrial models, fiber orientation is extracted based on imaging data and ex vivo tissue analysis (Zhao et al., 2012). Based on successive images from sheep atria, the myofiber orientation was determined by the purpose-development structure tensor analysis whereas the role of the fiber architecture on AF initiation has been studied (Zhao et al., 2013a). Kharche et al. (2014a) constructed a detailed model of rabbit atria from micro-CT images and found that fiber orientation anisotropy serves as a regulator of atrial activation while any modification of myofiber architecture may be pro-arrhythmic.
To develop realistic anatomical models of the atria, Zhao et al. (2008) suggested a methodology, based on ex vivo samples from pigs. During the following years, several models were designed to describe a realistic human or animal atrial anatomy (Deng and Xia, 2010; Aslanidi et al., 2011b; Tobón et al., 2013; Li et al., 2014; Ferrer et al., 2015; Zhao et al., 2017), whereas Zhao et al. (2017) proposed a novel 3D computational and structural model of human atria. This model is the output of a systematic analysis. It includes wall thickness, fiber orientation, and transmural fibrosis, and it is the most complete model of the human atrial structure.
Gharaviri et al. (2012a) introduced the novel approach of a dual-layer model for the study of AF. This model has been identified in 11 articles in the current review. Its basic hypothesis is that structural remodeling does not only lead to dissociation within the epicardial layer but also dissociation between the epicardial and endocardial layer promoting AF. AF induces this type of structural remodeling, based on experimental findings from goats (Eckstein et al., 2011), while it was revealed that this dissociation increases AF stability (Gharaviri et al., 2012b, 2017; Verheule et al., 2013). A methodology for the creation of a bilayer model of human atria, using CT images, is proposed by Vigmond et al. (2013), where the difference between RA and LA thickness, as well as the main anatomical structures and fiber orientation, is taken into account. The dual layer approach has also been adopted to examine the effectiveness of different radiofrequency ablation strategies (Bayer et al., 2016).
Simulation levels
Electrical activation of each cell results in modification of the electrophysiological properties spanning from the cell to the tissue, organ and body level. According to the scope of each study different observational outcomes were considered, such as AP, activation map, electrogram or ECG. In the cell level simulations, a single cell is considered. In the tissue level studies, the atrial excitation is simulated in a 2D level including several anatomical structures (Seemann et al., 2005; Carrillo et al., 2008). In the organ level studies, 3D models of the atria are considered. 3D structures, like icosahedron (Ellis et al., 2000), or simple atrial representations (Vigmond et al., 2001) have been proposed, while realistic geometrical models of human atria were constructed using MRI (Virag et al., 2002; Aslanidi et al., 2011b; Colman et al., 2013) or CT images (Burdumy et al., 2012). Furthermore, advanced imaging techniques, like late gadolinium-enhanced magnetic resonance imaging (LGE-MRI) allow the construction of patient-specific atrial substrate with fibrosis (McDowell et al., 2012). In body level articles, a realistic model of a human torso is adopted to simulate body surface potentials and the produced ECG. One of the first works in this field is the one of Jacquemet et al. (2006) where an ECG signal representing atrial activation is created, using a geometrical model atrial and torso model, based on MR images. Finally, Krueger et al. (2013a) proposed torso models to simulate surface ECG recordings.
Studying the electrophysiological changes on AP and ion kinetics (Jones et al., 2012; Sanchez et al., 2012; Voigt et al., 2013, 2014) is the main subject of cell level simulations. 2D level studies target at identifying AF drivers, such as rotors and spiral waves (Becerra et al., 2014; Ganesan et al., 2016; Salmin et al., 2016; Duarte-Salazar et al., 2017; Aronis and Ashikaga, 2018) or the effect of spatial resolution on AF drivers detection (Rappel and Narayan, 2013; Roney et al., 2017). The use of patient-specific 3D models can give more insight into the mechanism of AF (Zahid et al., 2016). Those simulations are used for the analysis of the mechanisms for AF initiation (Virag et al., 2002; Chang and Trayanova, 2016; Alcaine et al., 2017), AF stability (Rottmann et al., 2015) or the rotor drifting toward PVs (Calvo et al., 2014) and their anchoring in obstacles, such as fibrotic regions (Morgan et al., 2016; Nguyen et al., 2016; Deng et al., 2017). The detection of rotors is also studied using 3D simulations of human atria (Ugarte et al., 2014; Berenfeld, 2016; Martinez et al., 2016) to investigate the effectiveness of an ablation procedure (Dang et al., 2005; Hwang et al., 2014).
Apart from activation patterns analysis, signals can also be simulated (such as electrograms and ECG recordings). In the 2D or 3D case, the simulated electrograms can serve as a tool for AF drivers identification (Benson et al., 2014; Hummel et al., 2017; Sahli Costabal et al., 2018), atrial activation analysis and complex fractionated atrial electrogram (CFAE) regions detection (Jacquemet et al., 2003; Jacquemet and Henriquez, 2009; Ugarte et al., 2013, 2015; Roney et al., 2016). Jacquemet et al. (2005, 2009) constructed synthetic ECG signals to evaluate ventricular activity cancelation techniques. 3D simulations have been conducted to analyze the ECG characteristics, such as P-wave morphology (van Oosterom and Jacquemet, 2005; Ogawa et al., 2007; Krueger et al., 2013a, 2013c). On the other hand, the inverse problem solutions and body surface recordings were employed for the detection of AF drivers (Pedrón-Torrecilla et al., 2014; Figuera et al., 2016b; Guillem et al., 2016; Pedron-Torrecilla et al., 2016). Finally, an adaptation of the standard ECG systems for the optimization of the atrial signal recording was proposed by Ihara et al. (2007) and van Oosterom et al. (2007).
As shown in Figure 7, only a small minority of the studies adopt the body level simulations, whereas, more and more studies adopt a 3D realistic model of human atria to study the electrical activation during AF (organ level). Furthermore, the interest in the cell level analysis is persistent, as this constitutes the basis for any simulation protocol.
Figure 7. Evolution of the number of studies focusing each of the simulation levels (bar chart) and the distribution of the studies for the last 20 years (pie chart). For 2018 1 study targeted cell level and 2 studies the tissue level (until 1st of March 2018).
Electromechanics focus on the estimation of the electrical activation of the atria under the consideration of their mechanical function. Eight articles were found to include the electromechanical function in the simulations. In general, atrial contractility is associated with tissue depolarization as well as with the structural remodeling co-existence. Kuijpers et al. (2007) published the first model that integrates cardiac electromechanics with physiological details. According to this model, stretching of an atrial fiber influences the stretch-activated current (). The model complies with experimental observations of conduction slowing and blocking as a result of an acutely dilated atria. This model was integrated into a 3D geometrical model of the atria in order to investigate the role of stretch-activated current on the onset and termination of AF episodes (Kuijpers et al., 2009). The role of acutely dilated atria on AF initiation was later reported (Kuijpers et al., 2011). Ni et al. (2014) used an electromechanical model of the ventricles (Adeniran et al., 2013) to study the effect of KCNA5 mutation on the electromechanical function of the human atrial cell while Adeniran et al. (2015) studied the effect of AF-induced electrical remodeling on the electromechanical function of the heart. Finally, a more recent study focused on the analysis of spiral waves when the atrial mechanoelectrical function is considered (Brocklehurst et al., 2017). In this study, one electrophysiological model of human atria (Colman et al., 2013) was modified to include the stretch-activated current, while combined with a mechanical model (Rice et al., 2008) which describes the active force generated in cellular level in response to the electrical signal. This study highlights the effect of mechanoelectrical feedback on AF.
Modeling Disease Context and Purpose
Rhythm status (SR or AF)
In a healthy heart, the electrical pulse originates from the sinus node and spreads through the atrium (Platonov, 2012), while during an AF episode chaotic activation of atrial myocardium occurs (Yamazaki and Jalife, 2012). In the majority of the papers included in this review, AF was induced and the analysis focused on the identification of AF drivers (Ashihara et al., 2004; Martínez et al., 2014b; Tobón et al., 2014b; Ugarte et al., 2014, 2015; Felix et al., 2015; Martinez et al., 2016; Salmin et al., 2016; Roney et al., 2017; Duarte-Salazar et al., 2017) or on the study of the electrical activation during AF in general. However, analysis of the electrical activation of the remodeled atria during SR was the subject of 52 studies, as this can also improve our understanding of the pathology.
The main property of Sinoatrial Node (SAN) models is their auto-oscillatory behavior (Fink et al., 2011). In most cases, they are used under normal conditions instead of a remodeled substrate. The effect of SAN dysfunction on AF initiation is addressed in several studies. Kharche et al. (2017) developed a 3D electro-anatomical model of the SAN, based on the Fenton-Karma model. The model was based on histological studies, and the role of the SAN exit pathways was examined in terms of arrhythmias initiation. Glynn et al. (2014) adapted a detailed murine SAN model to confirm the restitution hypothesis, which states that the APD is associated with the duration of the previous diastolic intervals. They concluded that cycle length restitution analysis can enhance our understanding on SAN dysfunction and the initiation of arrhythmias. The simulations performed denote that the restitution curve can be useful for the analysis of cell dynamics and dysfunction in SAN. Lin et al. (2017) used a simplified spherical model of LA to investigate the mechanism of SR to AF transition. It was found that during SR, in healthy subjects, SAN to AV node conduction occurs through the interatrial septum via the fast pathway, while during AF, SAN long refractory period protects it from overdrive and the AV node filters the APs which would reach the ventricles (Li et al., 2014).
A common reason of SR usage is to test the reproducibility of the computational models, in both AF and normal SR (Goodman et al., 2005; van Oosterom and Jacquemet, 2005; Cherry and Evans, 2008; Grandi et al., 2011; Burdumy et al., 2012; Krueger et al., 2013a, 2013d; Colman et al., 2016) or to examine differences in activation patterns with or without electrical remodeling (Zhao et al., 2012; Voigt et al., 2014). Simulation under SR, using patient-specific models, including fibrosis, can accurately reproduce local activation (Krueger et al., 2014). The reconstruction of the ECG signals by solving the forward problem of electrophysiology, to understand the genesis of the P-waves, was also found to be very accurate when the patient is in SR (Ferrer et al., 2015; Jacquemet, 2015; Figuera et al., 2016b). Finally, Zhao et al. (2012, 2013a) highlighted that the structurally remodeled substrate and the anatomical heterogeneities can increase the vulnerability to rhythm disturbances.
In addition, considering patients during SR, the effect of antiarrhythmic drugs, such as potassium or sodium channel blockers, on the APD or ionic current characteristics, has been studied (Wettwer et al., 2004; Ehrlich et al., 2008; Law et al., 2010; Aslanidi et al., 2012a; Marshall et al., 2012; Voigt et al., 2013; Ni et al., 2016; Skibsbye et al., 2016; Ellinwood et al., 2017b; Tobón et al., 2017). Ablation success, applied during SR, has been presented (Tobón et al., 2015; Bayer et al., 2016) and the changes occurring in the P-waves after the PV isolation procedure have also been examined (Saha et al., 2016).
AF initiation
The conditions, under which AF can be induced, were extensively studied and several approaches have been proposed. In general, a triggering factor, in the presence of an arrhythmogenic substrate, can favor the initiation of PAF. PVs have been reported to present pacemaking automaticity and triggering activity, and thus their role on arrhythmogenicity has been examined (Reumann et al., 2006; Cherry et al., 2007), whereas mathematical models have also been implemented to describe PVs in both animals (Lo et al., 2006; Seol et al., 2008) and humans. Jones et al. (2012) modeled the PVs under the experimental evidence that two types of cells exist in PVs, pacemaking, and non-pacemaking.
Abramovich-Sivan and Akselrod (1999) initiated a new approach, where a model of the atrial cells, able to describe the transition from normal rhythm to AF, is proposed. According to this study, the atrial cells are considered as “pacemaker-like” cells with different cycle length and different type of intercellular connection, where, the effect of the junctions between the different types of cell are taken into account to investigate how the initiation of atrial rhythm disturbances occurs and maintains.
The role of the coupling between anatomical heterogeneous regions, such as CT (Ellis et al., 2000; Reumann et al., 2007; Dorn et al., 2012; Butters et al., 2013; Zhao et al., 2013a), or between the PV and the LA (Aslanidi et al., 2012b, 2013) on the initiation of AF is discussed, whereas atrial contraction and electromechanics are proposed as a mechanism of arrhythmogenesis (Kuijpers et al., 2009). Fibrosis (McDowell et al., 2012; Zhao et al., 2013b; O’Connell et al., 2015) or other anatomical obstacles (Ciaccio et al., 2017), as well as the imperfect ablation lines (Dang et al., 2005) have also been proposed to facilitate AF initiation.
The effect of specific ion currents and electrolyte concentration on arrhythmogeneity was studied (Kuijpers et al., 2006; Aslanidi et al., 2008; Law et al., 2009; Varela et al., 2014; Voigt et al., 2014; Onal et al., 2017) whereas electrical remodeling, induced by AF, and its role on episodes induction is discussed by Colman et al. (2013) and Lee et al. (2016). Furthermore, Morgan et al. (2014) reported that co-existence of structural and electrical remodeling can promote AF initiation.
AF initiation by specific diseases, such as ankyrin-B syndrome (Wolf et al., 2013) and end-stage renal disease, was also examined. In particular, HD effect on AF is highlighted in several studies (Severi et al., 2008, 2010; Krueger et al., 2011) suggesting that HD sessions may be a triggering factor for the onset of the arrhythmia (Vincenti et al., 2014).
Atrial ectopies and heart rate (HR) variations have reported to induce AF. Virag et al. (2002) used three different stimulation protocols. The first one is the S1-S2 protocol, where the SA node stimulation is followed by an ectopic beat located in the RA. The second one is the S1-S2-S3 protocol where an additional ectopic beat is applied in the same areas as the previous one, while the third protocol is a burst pacing protocol of 20 Hz applied in the SA node for several seconds. These protocols, applied in both normal and electrically remodeled substrate, simulate both self-terminated AF episodes as well as more sustained types of AF. It is worth mentioning that most of the studies, analyzing AF drivers, use the S1–S2 protocol for AF initiation (Ramirez et al., 2000; Liu et al., 2007; Tobón et al., 2014b; Ugarte et al., 2015). Shusterman et al. (2003) mentioned that abrupt short-term variations in cycle length, such as atrial ectopic activity and HR variations, lead to spatial irregularities of the wavefront propagation and thus initiate reentries in a simple two-dimensional tissue model. The role of alternans as a proarrhythmic state has also been discussed (Chang et al., 2014; Chang and Trayanova, 2016). Gong et al. (2007) proposed that ectopies originating from the PVs facilitate reentry induction by prolonging the vulnerability window. Hwang et al. (2016) and Hwang et al. (2017) examined the effect of the ANS on arrhythmogenesis where the role of ganglia on the initiation of AF is highlighted. The concept of early warning that can forecast a change of system’s stability is applied in AF electrophysiology (Garcia-Gudino et al., 2017). In this study, a cellular automaton is used to describe atrial conduction on the tissue. Applying different pacing approaches, along with real RR interval data derived from Fantasia database (Iyengar et al., 1996), it was shown that even when complex pacing is introduced, this approach can serve as a predictor of transition from normal rhythm to AF.
Finally, genetic background and gene mutations represent an additional field of research on AF initiation (Carrillo et al., 2008; Seemann et al., 2009; Kharche et al., 2012a; Whittaker et al., 2017).
Modeling the interventions for AF termination
Atrial fibrillation termination approaches include pharmacological treatment, ablation procedure, and electrical cardioversion. The first study identified on defibrillation usage for AF termination is the one by Boegli et al. (2000) where different defibrillation techniques were tested for their effectiveness. Kharche et al. (2012b), using a 3D realistic model of human atria, observed that low energy cardioversion efficacy depends on atrial anatomy. It has also been mentioned that in patients with no electrical remodeling lower cardioversion energy and shorter duration must be applied (Kharche et al., 2013). Rusu et al. (2014) proposed that pacing therapy must be patient-specific, due to its influence by the presence of heterogeneities in vagal activation and repolarization. On the other hand, it has been suggested patient-specific management of defibrillation (Kharche et al., 2015; Luca et al., 2016).
Dang et al. (2005) evaluated the outcome of the radiofrequency ablation and studied the effect of imperfect lines. In this study, the significance of the biophysical models’ application on the development of new ablation approaches is highlighted. The simulation results were confirmed with clinical data. The identification of the optimal ablation line is the subject of several studies which point out the need for tailored ablation lines for each patient (Rotter et al., 2007; Ruchat et al., 2007a, 2007b, 2007c; Reumann et al., 2008; Kwon et al., 2013; Hwang et al., 2014; Zhu et al., 2014; Rappel et al., 2015; Bayer et al., 2016). During the ablation procedure, in most cases, the PVs are isolated. However, fibrotic region ablation is also proposed (Ashihara et al., 2012) whereas other studies suggest other approaches which can guide the ablation procedure (Benson et al., 2014; Tobón et al., 2014b; Carrick et al., 2016; Kuklik et al., 2016; Li et al., 2016). In a 3D model of LA, the ablation of ganglia seems to terminate AF (Hwang et al., 2017). Interestingly, the effect of the poor contact of the electrodes with the tissue on AF drivers identification is discussed, and a methodology for signals reconstruction is proposed by Green et al. (2017).
Apart from the approaches mentioned above, drug treatment is also followed for AF termination. Most of the antiarrhythmic drugs, acting on atrial myocardium, may also have an effect on ventricular ion channels. State-dependent -channel blockers (class I antiarrhythmic drugs) maximize their effect on the atrium while at the same time the influence on the ventricles is minimal. AF-selectivity of such drugs is the topic of Aguilar-Shardonofsky et al. (2012). According to this study, increased ventricular proarrhythmia correlates with drug effectiveness close to 100%, while in case of minimum ventricular proarrhythmia, AF termination rates were low. Furthermore, several potassium channel blockers were also analyzed, in terms of computer models, in order to examine their effect on AF termination (Scholz et al., 2013; Tobón et al., 2014a; Aguilar et al., 2017; Ellinwood et al., 2017a). Ni et al. (2017b) compared the effects on atrial and ventricular activation of these two types of blockers combination to single channel blockers usage. The human heart simulation results showed that synergetic antiarrhythmic treatment leads to more effective termination of reentries, compared to single channels block, while the side effects in the ventricles are limited. Moreover, analyzing inter-subject variability on AP, distinguishes drug responders from non-responders (Liberos et al., 2016).
The self-terminating nature of the arrhythmia has also been studied (Uldry et al., 2012), however, its transient nature is the main barrier while an analysis framework has been proposed (Lin et al., 2017). It was found that in low complexity AF, the spontaneous termination process lasted 3.2 s and was initially started in the LA, while the last fibrillatory signal was detected in the RA. On the other hand, in complex AF, the termination process is less predictable. Additionally, the role of electro-mechanical feedback on AF self-termination is discussed by Kuijpers et al. (2009).
Evolution of Models
Computational modeling evolution can be studied by network theory. The components, involved in AF modeling, are considered as nodes of a network, whereas, the connections between them represent the studies they are reported in. Two distinct periods (1998–2007 and 2007–2018) were considered, and two metrics, describing the network complexity were extracted. The first metric is the average degree (AD) of the graph, which is defined as the average connections between the nodes, while the second one is the network density, defined as , where E is the number of edges and N the number of nodes.
According to the analysis AD was 16.96 during the first decade its value increased to 21.52 during the second decade, demonstrating an increase in network complexity as more factors were involved in the computational modeling. D was also increased during the second decade, from 0.707 to 0.897. These results suggest that the simulation experiments are getting more and more complex, taking into account additional factors which influence AF.
As seen in Figure 8, the study of subjects such as the dual layer, the electromechanics and the involvement of ANS, gained the attention of the research community during the second decade. Additionally, 3D model simulations also doubled their frequency which can be attributed to the improvement of imaging techniques and the creation of detailed anatomical models of the atria. Finally, it is observed that more studies focus on atrial activation analysis during SR and not only while on AF.
Figure 8. The network chart depicted the evolution of AF modeling complexity during the last 20 years (1998–2008: left, 2009–2018: right). The nodes of the network represent the basic components of the AF modeling concept, and the edges represent their connections. The stronger the prevalence of the nodes and edges, the higher their thickness, and the darker their color.
On the other hand, the proportion of studies, conducting an AP model analysis seems to be the same between the two decades, implying a stable trend in the analysis of APs.
AF Models: What Do We Learn About ECG and P-Wave Morphology
One of the great challenges in the field of computational modeling is to translate the output of the simulation studies into clinical practice (Trayanova and Chang, 2016). ECG recordings constitute an accurate non-invasive tool for heart activity assessment which is extensively useful in daily practice. ECG signals extraction from computational studies and vice versa would improve AF pathophysiology understanding and thus optimize clinical decision and treatment strategy. In the present review, several approaches for ECG signal reconstruction have been identified.
An attempt to extract a pseudoelectrocardiogram (pECG) was made by Namba et al. (1999) where a point in a realistic distance of the excitable tissue was considered and the pECG was calculated as the integration of all the membrane potentials, assuming a homogeneous and zero resistance area. In this work, it was found that the spiral wave results on a low Atrial Fibrillation Interval (FF) while under sodium channel blockers the FF intervals increase. Aslanidi et al. (2011b) described the development of the 3D virtual human atria which is a multi-scale detailed computational model. The anatomical model of the atria was based on MR images where the main anatomical structures, with different electrical characteristics, and fiber orientation were included. The 3D virtual atria model was adjusted in a torso model, enabling several activation patterns simulation and ECG extraction. A detailed anatomical model of human atria and torso, based on histological data and MR images, is also described by Ferrer et al. (2015). The concept of this work was to create a multi-scale framework used to explain, from an electrophysiological point of view, findings that are detectable in P-waves and body surface potential maps. P-wave morphology and duration were observed to be highly influenced by the conductivity properties of the torso and the fiber orientation in the atria.
Kuijpers et al. (2011) tested the effect of mechanical function on AF initiation and its effect on the extracted ECG was estimated. In another work from Krueger et al. (2013d) the personalized anatomical model in terms of the location of the interatrial bridges is studied. MRI images were used to create a personalized anatomical model while local activation time (LAT) maps were recorded using electroanatomical mapping systems [Ensite NavX (St. Jude Medical, St. Paul, MN, United States) or Carto (Biosense Webster, Haifa, Israel)] before the ablation procedure. Those LAT maps were used to localize the interatrial bridges in each patient. The simulation results, during SR, denote personalized interatrial bridges consideration can improve simulated and real body surface potentials concordance.
In an ECG signal atrial activity, especially atrial repolarization is hidden by the much higher amplitude of the ventricular electrical activity. In this respect, the cancelation of the ventricular wave is of great importance. Validation of ventricular activity cancelation algorithms is feasible employing computational modeling and ECG signal reconstruction. Jacquemet et al. (2009) created a realistic 3D model of human atria and torso, based on MR images, and the signal, reflecting atrial activity, was reconstructed in the points that a standard ECG lead is measured. That signal has later been summed with QRST complexes recorded during SR, adapting the RR intervals accordingly. This work does not focus on the analysis of the F-waves, but it provides the methodological framework for the evaluation of QRST cancelation techniques.
The standard ECG-lead systems, such as the 12-lead ECG, were designed in such a way to capture the ventricular activity effectively. Thus, these lead systems may be sub-optimal in terms of atrial activity detection. A few studies have been conducted (Ihara et al., 2005, 2007) proposing the modification of these systems toward the effective capturing of atrial activity. Additionally, van Oosterom et al. (2007) suggested a modification of electrodes’ placing vectorcardiographic (VCG) signals acquisition in AF patients, as, according to the authors, the Frank lead system seems suboptimal for atrial activity monitoring during AF. These studies are based on 3D models of atria and torso, where the ECG on the body surface is calculated using simulations.
Electrocardiography and vectorcardiographic recordings have been analyzed to identify AF drivers as well as to characterize atrial activation complexity. An interesting visual representation of the atrial activation, using the time course of the VCG dipole components, is based on a sphere where the main vessel and vein orifices were included. The trajectory of the atrial activation demonstrates the differences between SR and different types of AF (Jacquemet et al., 2006) while the extracted spatial characteristics can localize stable AF drivers in the atria (Lemay et al., 2007b, 2009). Lemay et al. (2007a) proposed that the existence of multiple AF drivers in the atrium can be revealed through the analysis of the full spectrum of the ECG signal instead of the analysis of the dominant frequency only. Duchêne et al. (2009) used a multiple frequency tracking algorithm to examine whether ECG signals analysis can provide information related to AF drivers complexity. Furthermore, rotor activity detection using spherical atrial and torso model and surface phase maps extraction is proposed in Rodrigo et al. (2014a, b). Finally, Ni et al. (2017b) examined the effect of antiarrhythmic drugs on the power spectrum of the simulated ECG.
In the surface ECG, atrial depolarization is mainly represented by P-wave. While AF is characterized by the absence of P-waves in the ECG signals, during SR, this complex can provide information related to the underline substrate. It is known that the first part reflects the activation of the RA and the second part the activation of the LA (Michelucci et al., 2002). According to Platonov (2012), the origin of the sinus rhythm, the shape and the size of the atria can affect the morphological characteristics of the P-waves. Solving of the forward or the inverse problem of the electrophysiology helps to translate the P-wave characteristics on tissue level phenomena. However, it is worth mentioning that both the inverse and forward approaches require a detailed representation of the torso as this affects the precision of the results (Bear et al., 2015).
Forward Problem – P-Wave Estimation
The estimation of the body surface potentials based on the epicardial ones, define the forward problem of electrophysiology. The key work from van Oosterom and Jacquemet (2005) used the Equivalent Double Layer (EDL) source model to simulate the atrial activity of the human body. The proposed approach was applied both during SR and AF. During SR, the application of the EDL model presented high correlation with the actual recorded signal, primarily related to the P-wave morphology. Furthermore, during AF, simulated signals were also highly correlated with the signals that are typically recorded in AF patients.
Krueger et al. (2013c) presented a framework for the generation and validation of multiscale cardiac modeling. This study describes how a priori knowledge, such as fiber orientation and electrophysiological models, can be enhanced by information derived from ECG characteristics like P-wave duration or HR, electrolytes concentrations and information from cardiac and thoracic MRI. The availability of LATs and ECG recordings can be used for models evaluation. As highlighted, HR must be considered during the simulation process as it has a significant impact on atrial repolarization. Additionally, P-wave morphology was properly reproduced only in specific leads, while P-wave duration is accurately computed in any lead. Krueger et al. (2013a) simulated P-wave as well as T-wave, based on the patient-specific detailed anatomical model. The simulated P-waves and T-waves morphology fitted well with the experimental recordings in both physiological and AF-remodeled conditions.
A typical example of macroscopic phenomena translation is to predict AF recurrences by analyzing the P-wave shortening observed in patients undergoing PVs isolation. P-wave duration reflects the interatrial conduction time while in AF patients the heterogeneous atrial propagation results in prolonged P-waves. Ogawa et al. (2007) used signal-averaged P-wave duration, P-wave dispersion and root-mean-square voltage of the terminal 20 ms (RMS20) of the averaged P-wave, to predict AF recurrences after successful PV ablation. Comparing real data from patient underwent cardiac ablation with P-wave generated by mathematical simulation of realistic human atria, it was found that the elimination of PV sleeve results on shortening of the P-wave duration and modification of its morphology. Jacquemet (2015) discussed LA and RA substrate contribution on P-wave morphology confirming that P-wave reflects any structural modification of the atrial substrate.
An effort to interpret the atrial substrate mechanisms, which lead to changes in P-wave duration after the ablation procedure, was performed by Saha et al. (2016). Using mathematical models of human atria and torso, VCG and 16-lead ECG were reconstructed before and after PVs isolation. According to the results, P-wave shortening is detected in leads V3 and V4 whereas these changes were more apparent when larger PV areas were isolated. Furthermore, the P-wave area was increased in V1, while decreased in V6. Additionally, possible reconnections in the PVs could be identified, and thus the 16-lead ECG can be used as a tool for PVs reconnections detection in patients undergone ablation.
P-wave terminal force in lead V1 (PTF-V1) has been proposed as an AF risk marker. PTF-V1 is the product of the duration and the amplitude of the terminal negative part of the P-wave (Morris et al., 1964). Loewe et al. (2016) studied the effect of the Early Activation Sites (EAS), as well as, their proximity to inter-atrial connections (IACs) on the PTF-V1. The term EAS is defined as the area where the sinus node excitation is captured by the myocardium. An inter-personal variability, as well as, a variation over time on the same individual has been observed. Based on experimental findings from dogs, this variation over time is suggested to reflect the vagal activation. On the other hand, the location and the conductive properties of the IACs are known to vary among subjects. Several anatomical models and multiple EAS consideration have shown that these two factors affect PTF-V1.
Alday et al. (2014, 2015) examined the effect of ectopic foci on P-wave morphology. In the latter study, the P-wave morphology is examined in terms of polarity (positive, negative, and biphasic). Under the assumption that a negative deflection on an ECG depicts a wavefront traveling away from the location where the electrodes are placed, the torso and the atria were separated into four quadrants. The analysis of the P-wave polarity in each of the quadrant can reveal the area from which the ectopic foci originates, while the implemented algorithm demonstrated high success rates on the localization of the ectopic source.
As mentioned, alternations on the ionic currents during HD sessions can induce PAF episodes. Observational studies have shown P-waves prolongation during HD sessions implying intra-atrial conduction slowing. A computer model describing the ionic concentrations of a patient during HD has been used to explain this macroscopic phenomenon (Severi et al., 2010). alternations (hypokalemia) result in conduction velocity reduction and ERP shortening. Krueger et al. (2011) confirmed these findings, where simulated ECGs were extracted, using realistic models of human atria and torso.
Finally, Zhu et al. (2014) examined P-wave morphology to check the effectiveness of 4 ablation protocols, on cAF treatment. Without providing any details on the evaluation of P-wave morphology, it is mentioned that simulated P-waves were similar after the application of different ablation protocols. The authors concluded that the Maze III approach is effective, and the transmural ablation has the same effectiveness as the non-transmural one. Table 3 summarizes the features mentioned above for the analysis of the P-wave using computational modeling approaches.
Table 3. Features used toward the analysis of the P-wave, during SR, using computational modeling approach.
Inverse Problem
Estimation of bioelectric source potentials using remote electric field measurements is known as the inverse problem of electrophysiology, and it has been applied in both SR and AF scenarios (Pedrón-Torrecilla et al., 2014). The accuracy of the extracted results is a crucial topic and several studies focused on the implementation of algorithms able to reconstruct activation patterns matching with the observed ones (Pedron-Torrecilla et al., 2016). Zemzemi et al. (2012) proposed one such approach based in machine learning techniques where the simulated ECGs, computed by the consideration of atria and torso model, were used for the training of the algorithms while the results were examined in terms of reproducibility. One of the application areas of this field is the identification of AF drivers in the atria (Guillem et al., 2016). On the other hand, different methods of estimating epicardial potentials, dominant frequency, phase maps, and singularity points locations, were examined for their accuracy (Figuera et al., 2016a, 2016b).
Discussion
Computational modeling is an emerging area of research useful to understand the mechanisms involved with pathophysiology, diagnosis, and treatment of human diseases (Bers and Grandi, 2011). This article presents the results of systematic scoping review related to multi-scale computational modeling of AF. In the context of this review, two databases were searched and 301 articles, published during the last 20 years, were identified as relevant for our purposes.
Atrial fibrillation is a complex arrhythmia, and the electrophysiological mechanisms for its initiation and maintenance are not clearly understood. In general, a triggering factor, as well as an arrhythmogenic substrate, must coexist for AF initiation. The complicated nature of AF is reflected in the significant number of studies included in the current review and the numerous factors to be considered on simulation protocols designing. The advances in imaging technologies led to the creation of detailed anatomical models, depicting the structural and electrical heterogeneities of the atria whereas, several mechanisms involved in arrhythmogenesis have been modeled related to electrical and structural remodeling. The computational approaches move toward more detailed models, incorporating multiple factors, such ANS effect, atrial electromechanical, the endo and epicardium layer dissociation, fibrosis and the genetic substrate. This progress will lead to patient-specific models creation (Krueger et al., 2013b) which can be useful in real-life cases analysis.
The observation of the number of articles published over time reveals an increased interest in the field of computational modeling of AF while this interest is projected to continue. Furthermore, the number of papers focusing on in silico AF models is increasing, and multiple scales are being analyzed, ranging from cell to tissue and organ level. In time to come, the advances of mapping and imaging techniques as well as high-performance computing engineering, will force computational models to adopt new elements to improve the repeatability and consistency with clinical observations (Trayanova and Chang, 2016). The main challenge is to transform the experimental knowledge into clinical practice in order to serve as a tool for physicians to perform for clinical interventions (Trayanova and McDowell, 2018). The importance of Physiome project (Hunter and Borg, 2003) and the Virtual Physiological Human (VPH) initiative (Viceconti and Hunter, 2016) and the positive effect they had on the increased interest on computational modeling must be highlighted.
Apart from atrial substrate analysis during AF, some articles focus on atrial excitation during SR. Computer simulations can help interpret the findings from ECG and P-wave analysis. As Dössel et al. (2017) mentioned, merging of biosignal analysis with computer modeling can enhance our understanding of AF. Atrial geometry (Ihara et al., 2007) as well as the modification of anatomical IACs have also shown the ability to change the morphology of the P-waves (Loewe et al., 2016) while the effect of structural remodeling on the P-wave (Jacquemet, 2015) such as its shortening after the PV isolation procedure (Saha et al., 2016), has also been reported. Furthermore, P-wave is supposed to be reformed by ectopic foci (Alday et al., 2015), suggesting that any deflection from the SA node excitation can alter the conduction routes. P-wave morphology depends on (1) the origin of the SR as well as (2) the conduction route (Platonov, 2012). However, many studies have revealed the existence of a distinct secondary P-wave morphology in PAF patients (Filos et al., 2017), while the percentage of beats matching the main or the secondary P-wave morphology predict the outcome of PV isolation (Huo et al., 2015) or the chance for AF initiation (Martínez et al., 2012). Monfredi et al. (2010) investigated whether the variability in the P-wave morphology can be attributed either on a transient switch from one group of SAN pacemaker cells to another or on multiple activation routes. Studies in dogs revealed that ANS, and in particular the different branches shifted the EAS in different regions (Platonov, 2012). Apart from the effect of the ANS, additional factors which alter on different time scales, such as electromechanics and subcellular phenomena like metabolism, must also be included in the analysis of the P-wave morphology variability. However, there is a need to consider longer simulations times and analysis of more beats to capture the transient modification of the P-wave morphology.
However, atria must not be regarded as being isolated from the rest of the heart. The effect of AV node on AF, as well as the role of accessory pathways in conditions, such as Wolff–Parkinson–White Syndrome, must also be considered. Recently, the predictive value of corrected QT interval (QTc) prolongation, a metric of ventricular repolarization, on the AF incidence was highlighted (Zhang et al., 2018). However, this metric can be considered as a surrogate marker. Future simulation protocols should also incorporate the role of other anatomical structures of the heart to investigate the distributed conduction in the atrial substrate which is implied by the variable P-wave morphologies.
Any effort made in the field of computational modeling must seriously consider how the outcome can be used in clinical practice (Grandi and Maleckar, 2016). Guided ablation is one of the main topics influenced by the advances in cardiac simulations (Zhao et al., 2015; Jacquemet, 2016). Computer models can identify AF drivers effectively when the effect of multiple factors on AF, such as structural heterogeneities and electrical remodeling, is considered. In the future, the analysis will be based on personalized atrial models since the evaluation of theoretical ablation procedures in such models showed their effectiveness on AF treatment. Upon application of these techniques in real patients, computational modeling is expected to optimize ablation procedures (Boyle et al., 2017). Furthermore, the development of more effective drugs with substantial atrial selectivity is foreseen (Grandi and Maleckar, 2016).
Limitations
Studies were limited to the AF while other pro-arrhythmic states were not included in the review. Furthermore, since one of the inclusion criteria was that the term “AF” needed to be indexed as a keyword, we recognize that some pioneering works on the field of computational modeling may be excluded. A typical example is the work from Vigmond et al. (2004), which was excluded because it was not indexed with the term “AF” as a keyword neither in PubMed nor in Scopus libraries, despite the fact its focus is on atrial arrhythmogenesis, and AF can be considered as a part of it. However, the fact that all the remaining inclusion criteria are met denotes the effectiveness of the query for the description of the field of computational modeling. Finally, the modeling of the AV node, which serves as a filter for the atrial AP to prevent the arrhythmia to be diffused into the ventricles, was not considered here, and thus its effect on the initiation and maintenance of AF was not studied.
Conclusion
Computer models become more and more complex as they include an increasing number of parameters. The evolution of the complexity, depicted in Figure 8, clearly reflects this pattern. The consideration of ANS branches, in terms of ganglia modeling, and the development of detailed representations of the anatomical heterogeneities, using MRI techniques, improved the reproducibility of the simulations compared to clinical observations. The advances on genomics is also an aspect which is under consideration. The advances on high-performance computing can overcome the computation load and the barriers introduced and thus it can allow the observation of more macroscopic events, such as the transition of the conduction routes which are reflected in the existence of multiple P-wave morphologies in AF patients during SR. Eventually, electrophysiological models of the heart will be used to optimize AF treatment and improve the quality of life in this group of patients.
Author Contributions
The study was initiated by the previous works of all the authors. DF and IC defined the inclusion/exclusion rules. DF, IC, and DT initiated the design, development, and synthesis of the studies and received feedback from NM and VV. DF and IC coordinated the writing of all drafts of the manuscript and all other authors contributed to the submitted versions of the manuscript. All authors have read and agreed to the manuscript being submitted as it is.
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
Abramovich-Sivan, S., and Akselrod, S. (1999). Simulation of atrial activity by a phase response curve based model of a two-dimensional pacemaker cells array: the transition from a normal activation pattern to atrial fibrillation. Biol. Cybern. 80, 141–153. doi: 10.1007/s004220050512
Adeniran, I., Hancox, J. C., and Zhang, H. (2013). In silico investigation of the short QT syndrome, using human ventricle models incorporating electromechanical coupling. Front. Physiol. 4:166. doi: 10.3389/fphys.2013.00166
Adeniran, I., Maciver, D. H., Garratt, C. J., Ye, J., Hancox, J. C., and Zhang, H. (2015). Effects of persistent atrial fibrillation- induced electrical remodeling on atrial electro-mechanics - insights from a 3D model of the human atria. PLoS One 10:e0142397. doi: 10.1371/journal.pone.0142397
Aguilar, M., Feng, J., Vigmond, E., Comtois, P., and Nattel, S. (2017). Rate-dependent role of ikur in human atrial repolarization and atrial fibrillation maintenance. Biophys. J. 112, 1997–2010. doi: 10.1016/j.bpj.2017.03.022
Aguilar, M., Qi, X. Y., Huang, H., and Nattel, S. (2014). Fibroblast electrical remodeling in heart failure and potential effects on atrial fibrillation. Biophys. J. 107, 2444–2455. doi: 10.1016/j.bpj.2014.10.014
Aguilar, M., Xiong, F., Qi, X. Y., Comtois, P., and Nattel, S. (2015). Potassium channel blockade enhances atrial fibrillation-selective antiarrhythmic effects of optimized state-dependent sodium channel blockade. Circulation 132, 2203–2211. doi: 10.1161/CIRCULATIONAHA.115.018016
Aguilar-Shardonofsky, M., Vigmond, E. J., Nattel, S., and Comtois, P. (2012). In silico optimization of atrial fibrillation-selective sodium channel blocker pharmacodynamics. Biophys. J. 102, 951–960. doi: 10.1016/j.bpj.2012.01.032
Alcaine, A., Masè, M., Cristoforetti, A., Ravelli, F., Nollo, G., Laguna, P., et al. (2017). A multi-variate predictability framework to assess invasive cardiac activity and interactions during atrial fibrillation. IEEE Trans. Biomed. Eng. 64, 1157–1168. doi: 10.1109/TBME.2016.2592953
Alday, E. A. P., Colman, M. A., Langley, P., Butters, T. D., Higham, J., Workman, A. J., et al. (2015). A new algorithm to diagnose atrial ectopic origin from multi lead ecg systems - insights from 3d virtual human atria and torso. PLoS Comput. Biol. 11:e1004026. doi: 10.1371/journal.pcbi.1004026
Alday, E. A. P., Colman, M. A., Langley, P., and Zhang, H. (2014). “Spatial refinement of a new algorithm to identify focus of atrial ectopic activity from 64-lead ECGs,” in Computing in Cardiology, (Washington, D.C: IEEE Computer Society), 501–504.
Almquist, J., Wallman, M., Jacobson, I., and Jirstrand, M. (2010). Modeling the effect of Kv1.5 block on the canine action potential. Biophys. J. 99, 2726–2736. doi: 10.1016/j.bpj.2010.08.062
Andrade, J., Khairy, P., Dobrev, D., and Nattel, S. (2014). The clinical profile and pathophysiology of atrial fibrillation: relationships among clinical features, epidemiology, and mechanisms. Circ. Res. 114, 1453–1468. doi: 10.1161/CIRCRESAHA.114.303211
Arksey, H., and O’Malley, L. (2005). Scoping studies: towards a methodological framework. Int. J. Soc. Res. Methodol. Theory Pract. 8, 19–32. doi: 10.1080/1364557032000119616
Aronis, K. N., and Ashikaga, H. (2018). Impact of number of co-existing rotors and inter-electrode distance on accuracy of rotor localization. J. Electrocardiol. 51, 82–91. doi: 10.1016/j.jelectrocard.2017.08.032
Ashihara, T., Haraguchi, R., Nakazawa, K., Namba, T., Ikeda, T., Nakazawa, Y., et al. (2012). The role of fibroblasts in complex fractionated electrograms during persistent/permanent atrial fibrillation: implications for electrogram-based catheter ablation. Circ. Res. 110, 275–284. doi: 10.1161/CIRCRESAHA.111.255026
Ashihara, T., Namba, T., Ikeda, T., Ito, M., Nakazawa, K., and Trayanova, N. (2004). Mechanisms of Myocardial Capture and Temporal Excitable Gap during Spiral Wave Reentry in a Bidomain Model. Circulation 109, 920–925. doi: 10.1161/01.CIR.0000118331.13524.75
Ashihara, T., Yao, T., Namba, T., Kawase, A., Ikeda, T., Nakazawa, K., et al. (2002). Differences in sympathetic and vagal effects on paroxysmal atrial fibrillation: a simulation study. Biomed. Pharmacother. 56, 359s–363s.
Aslanidi, O. V., Al-Owais, M., Benson, A. P., Colman, M., Garratt, C. J., Gilbert, S. H., et al. (2012a). Virtual tissue engineering of the human atrium: modelling pharmacological actions on atrial arrhythmogenesis. Eur. J. Pharm. Sci. 46, 209–221. doi: 10.1016/j.ejps.2011.08.014
Aslanidi, O. V., Colman, M. A., Zhao, J., Smaill, B. H., Gilbert, S. H., Hancox, J. C., et al. (2012b). Arrhythmogenic substrate for atrial fibrillation: insights from an integrative computational model of pulmonary veins. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2012, 203–206. doi: 10.1109/EMBC.2012.6345906
Aslanidi, O. V., Boyett, M. R., and Zhang, H. (2008). Effects of the intracellular Ca2+ dynamics on restitution properties and stability of reentry in rabbit atrial tissue model. Comput. Cardiol. 35, 295–297. doi: 10.1109/CIC.2008.4749036
Aslanidi, O. V., Boyett, M. R., and Zhang, H. (2009a). Left to right atrial electrophysiological differences: substrate for a dominant reentrant source during atrial fibrillation. Lect. Notes Comput. Sci. 5528, 154–161. doi: 10.1007/978-3-642-01932-6_17
Aslanidi, O. V., Robinson, R., Cheverton, D., Boyett, M. R., and Zhang, H. (2009b). Electrophysiological substrate for a dominant reentrant source during atrial fibrillation. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2009, 2819–2822. doi: 10.1109/IEMBS.2009.5333573
Aslanidi, O. V., Butters, T. D., Ren, C. X., Ryecroft, G., and Zhang, H. (2011a). Electrophysiological models for the heterogeneous canine atria: computational platform for studying rapid atrial arrhythmias. Conf. IEEE Eng. Med. Biol. Soc. 2011, 1693–1696.
Aslanidi, O. V., Colman, M. A., Stott, J., Dobrzynski, H., Boyett, M. R., Holden, A. V., et al. (2011b). 3D virtual human atria: a computational platform for studying clinical atrial fibrillation. Prog. Biophys. Mol. Biol. 107, 156–168. doi: 10.1016/j.pbiomolbio.2011.06.011
Aslanidi, O. V., Colman, M. A., Varela, M., Zhao, J., Smaill, B. H., Hancox, J. C., et al. (2013). Heterogeneous and anisotropic integrative model of pulmonary veins: computational study of arrhythmogenic substrate for atrial fibrillation. Inter. Focus 3:20120069. doi: 10.1098/rsfs.2012.0069
Bayer, J. D., Roney, C. H., Pashaei, A., Jaïs, P., and Vigmond, E. J. (2016). Novel radiofrequency ablation strategies for terminating atrial fibrillation in the left atrium: a simulation study. Front. Physiol. 7:108. doi: 10.3389/fphys.2016.00108
Bear, L. R., Cheng, L. K., LeGrice, I. J., Sands, G. B., Lever, N. A., Paterson, D. J., et al. (2015). Forward problem of electrocardiography: is it solved? Circ. Arrhythmia Electrophysiol. 8, 677–684. doi: 10.1161/CIRCEP.114.001573
Becerra, M. A., Murillo-Escobar, J., Palacio, L. C., and Zuluaga, C. T. (2014). “Noise and spatial-resolution effect of electrode array on rotor tip location during atrial fibrillation: A simulation study,” in Computing in Cardiology, (Washington, DC: IEEE Computer Society), 797–800.
Beeler, G. W., and Reuter, H. (1977). Reconstruction of the action potential of ventricular myocardial fibres. J. Physiol. 268, 177–210. doi: 10.1113/jphysiol.1977.sp011853
Benson, B. E., Carrick, R., Habel, N., Bates, O., Bates, J. H. T., Bielau, P., et al. (2014). Mapping multi-wavelet reentry without isochrones: an electrogram-guided approach to define substrate distribution. Europace 16, iv102–iv109. doi: 10.1093/europace/euu254
Berenfeld, O. (2016). The major role of IK1 in mechanisms of rotor drift in the atria: a computational study. Clin. Med. Insights Cardiol. 10, (Suppl. 1), 71–79. doi: 10.4137/CMC.S39773
Bers, D. M., and Grandi, E. (2011). Human atrial fibrillation: insights from computational electrophysiological models. Trends Cardiovasc. Med. 21, 145–150. doi: 10.1016/j.tcm.2012.04.004
Boegli, M., Blanc, O., Virag, N., Vesin, J.-M., and Kappenberger, L. (2000). Study of the defibrillation process in a computer model of human atria. Annu. Int. Conf. IEEE Eng. Med. Biol. Proc. 3, 1848–1849. doi: 10.1109/IEMBS.2000.900446
Boyle, P. M., Zahid, S., and Trayanova, N. A. (2017). Using personalized computer models to custom-tailor ablation procedures for atrial fibrillation patients: are we there yet? Expert Rev. Cardiovasc. Ther. 15, 339–341. doi: 10.1080/14779072.2017.1317593
Brocklehurst, P., Ni, H., Zhang, H., and Ye, J. (2017). Electro-mechanical dynamics of spiral waves in a discrete 2D model of human atrial tissue. PLoS One 12:e0176607. doi: 10.1371/journal.pone.0176607
Burdumy, M., Luik, A., Neher, P., Hanna, R., Krueger, M. W., Schilling, C., et al. (2012). Comparing measured and simulated wave directions in the left atrium a workflow for model personalization and validation. Biomed. Tech. 57, 79–87. doi: 10.1515/bmt-2011-0059
Butters, T. D., Aslanidi, O., Zhao, J., Smaill, B., and Zhang, H. (2013). A novel computational sheep atria model for the study of atrial fibrillation. Inter. focus 3:20120067. doi: 10.1098/rsfs.2012.0067
Cacciani, F., and Zaniboni, M. (2015). Chronotropic modulation of the source-sink relationship of sinoatrial-atrial impulse conduction and its significance to initiation of AF: a one-dimensional model study. Biomed. Res. Int. 2015:18. doi: 10.1155/2015/496418
Calvo, C. J., Deo, M., Zlochiver, S., Millet, J., and Berenfeld, O. (2014). Attraction of rotors to the pulmonary veins in paroxysmal atrial fibrillation: a modeling study. Biophys. J. 106, 1811–1821. doi: 10.1016/j.bpj.2014.02.030
Camm, A. J., Al-Khatib, S. M., Calkins, H., Halperin, J. L., Kirchhof, P., Lip, G. Y. H., et al. (2012a). A proposal for new clinical concepts in the management of atrial fibrillation. Am. Heart J. 164, 292.e–302.e. doi: 10.1016/j.ahj.2012.05.017
Camm, A. J., Lip, G. Y. H., De Caterina, R., Savelieva, I., Atar, D., Hohnloser, S. H., et al. (2012b). 2012 focused update of the ESC Guidelines for the management of atrial fibrillation: An update of the 2010 ESC Guidelines for the management of atrial fibrillation * Developed with the special contribution of the European Heart Rhythm Association. Europace 14, 1385–1413. doi: 10.1093/europace/eus305
Carrick, R. T., Benson, B. E., Bates, J. H. T., and Spector, P. S. (2016). Prospective, tissue-specific optimization of ablation for multiwavelet reentry: predicting the required amount, location, and configuration of lesions. Circ. Arrhythmia Electrophysiol. 9:e003555. doi: 10.1161/CIRCEP.115.003555
Carrillo, P., Seemann, G., Scholz, E., Weiss, D. L., and Dössel, O. (2008). Impact of the hERG channel mutation N588K on the electrical properties of the human atrium. IFMBE 22, 2583–2586. doi: 10.1007/978-3-540-89208-3_620
Censi, F., Corazza, I., Reggiani, E., Calcagnini, G., Mattei, E., Triventi, M., et al. (2016). P-wave variability and atrial fibrillation. Sci. Rep. 6:26799. doi: 10.1038/srep26799
Chang, K. C., Bayer, J. D., and Trayanova, N. A. (2014). Disrupted calcium release as a mechanism for atrial alternans associated with human atrial fibrillation. PLoS Comput. Biol. 10:e1004011. doi: 10.1371/journal.pcbi.1004011
Chang, K. C., and Trayanova, N. A. (2016). Mechanisms of arrhythmogenesis related to calcium-driven alternans in a model of human atrial fibrillation. Sci. Rep. 6:36395. doi: 10.1038/srep36395
Cherry, E. M., Ehrlich, J. R., Nattel, S., and Fenton, F. H. (2007). Pulmonary vein reentry-properties and size matter: insights from a computational analysis. Hear. Rhythm 4, 1553–1562. doi: 10.1016/j.hrthm.2007.08.017
Cherry, E. M., and Evans, S. J. (2008). Properties of two human atrial cell models in tissue: Restitution, memory, propagation, and reentry. J. Theor. Biol. 254, 674–690. doi: 10.1016/j.jtbi.2008.06.030
Ciaccio, E. J., Biviano, A. B., Wan, E. Y., Peters, N. S., and Garan, H. (2017). Development of an automaton model of rotational activity driving atrial fibrillation. Comput. Biol. Med. 83, 166–181. doi: 10.1016/j.compbiomed.2017.02.008
Colman, M. A., Aslanidi, O. V., Kharche, S., Boyett, M. R., Garratt, C., Hancox, J. C., et al. (2013). Pro-arrhythmogenic effects of atrial fibrillation-induced electrical remodelling: insights from the three-dimensional virtual human atria. J. Physiol. 591, 4249–4272. doi: 10.1113/jphysiol.2013.254987
Colman, M. A., Sarathy, P. P., Macquiaide, N., and Workman, A. J. (2016). “A new model of the human atrial myocyte with variable T-tubule organization for the study of atrial fibrillation,” in Computing in Cardiology, (Washington, DC: IEEE Computer Society), 221–224.
Colman, M. A., Varela, M., Hancox, J. C., Zhang, H., and Aslanidi, O. V. (2014). Evolution and pharmacological modulation of the arrhythmogenic wave dynamics in canine pulmonary vein model. Europace 16, 416–423. doi: 10.1093/europace/eut349
Colquhoun, H. L., Levac, D., O’Brien, K. K., Straus, S., Tricco, A. C., Perrier, L., et al. (2014). Scoping reviews: Time for clarity in definition, methods, and reporting. J. Clin. Epidemiol. 67, 1291–1294. doi: 10.1016/j.jclinepi.2014.03.013
Comtois, P., and Nattel, S. (2011). Impact of tissue geometry on simulated cholinergic atrial fibrillation: a modeling study. Chaos 21:13108. doi: 10.1063/1.3544470
Comtois, P., Sakabe, M., Vigmond, E. J., Munoz, M., Texier, A., Shiroshita-Takeshita, A., et al. (2008). Mechanisms of atrial fibrillation termination by rapidly unbinding Na + channel blockers: Insights from mathematical models and experimental correlates. Am. J. Physiol. Hear. Circ. Physiol. 295, H1489–H1504. doi: 10.1152/ajpheart.01054.2007
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, H301–H321. doi: 10.1152/ajpheart.1998.275.1.H301
Courtemanche, M., Ramirez, R. J., and Nattel, S. (1999). Ionic targets for drug therapy and atrial fibrillation-induced electrical remodeling: Insights from a mathematical model. Cardiovasc. Res. 42, 477–489. doi: 10.1016/S0008-6363(99)00034-36
Dang, L., Virag, N., Ihara, Z., Jacquemet, V., Vesin, J.-M., Schlaepfer, J., et al. (2005). Evaluation of ablation patterns using a biophysical model of atrial fibrillation. Ann. Biomed. Eng. 33, 465–474. doi: 10.1007/s10439-005-2502-7
Deng, D., Murphy, M. J., Hakim, J. B., Franceschi, W. H., Zahid, S., Pashakhanloo, F., et al. (2017). Sensitivity of reentrant driver localization to electrophysiological parameter variability in image-based computational models of persistent atrial fibrillation sustained by a fibrotic substrate. Chaos 27:93932. doi: 10.1063/1.5003340
Deng, D., and Xia, L. (2010). “Study the effect of tissue heterogeneity and anisotropy in atrial fibrillation based on a human atrial model,” in Proceedings of the Computing in Cardiology, (Belfast: IEEE), 433–436.
Dorn, A., Krueger, M. W., Seemann, G., and Doessel, O. (2012). Modelling of heterogeneous human atrial electrophysiology. Biomed. Tech. 57, 350–353. doi: 10.1515/bmt-2012-4215
Dössel, O., Lenis, G., Loewe, A., Pollnow, S., Rottmann, M., Verma, B., et al. (2017). “Atrial Signals – Modeling Meets Biosignal Analysis,” in IFMBE Proceedings, eds H. Eskola and O. Vaisanen (Tampere: Spinger), 723–726. doi: 10.1007/978-981-10-5122-7_181
Duarte, M., Restrepo, A. L., Tobon, C., Cardona, K., and Saiz, J. (2013). “Chloroquine effect on human atrial action potential under normal conditions and during paroxysmal and chronic atrial fibrillation. A simulation study,” in Proceedings of the Pan American Health Care Exchanges, PAHCE, (Medellin: IEEE).
Duarte-Salazar, C. A., Orozco-Duque, A., Tobón, C., Peluffo-Ordóñez, D. H., Luna, J. A. G., and Becerra, M. A. (2017). “Comparison Between Unipolar And Bipolar Electrograms For Detecting Rotor Tip From 2D Fibrillation Model Using Image Fusion. A Simulation Study,” in 2016 IEEE Latin American Conference on Computational Intelligence, LA-CCI 2016 - Proceedings, ed. G. J. B. Rodriguez (New Jersey: Institute of Electrical and Electronics Engineers Inc.).
Duchêne, C., Lemay, M., Vesin, J. M., and Van Oosterom, A. (2009). Adaptive multiple frequency tracking algorithm: detection of stable atrial fibrillation sources from standard 12-lead ECG. Comput. Cardiol. 36, 505–508.
Eckstein, J., Maesen, B., Linz, D., Zeemering, S., Van Hunnik, A., Verheule, S., et al. (2011). Time course and mechanisms of endo-epicardial electrical dissociation during atrial fibrillation in the goat. Cardiovasc. Res. 89, 816–824. doi: 10.1093/cvr/cvq336
Ehrlich, J. R., Ocholla, H., Ziemek, D., Rütten, H., Hohnloser, S. H., and Gögelein, H. (2008). Characterization of human cardiac Kv1.5 inhibition by the novel atrial-selective antiarrhythmic compound AVE1231. J. Cardiovasc. Pharmacol. 51, 380–387. doi: 10.1097/FJC.0b013e3181669030
Ehrlich, J. R., Zicha, S., Coutu, P., Hébert, T. E., and Nattel, S. (2005). Atrial fibrillation-associated minK38G/S polymorphism modulates delayed rectifier current and membrane localization. Cardiovasc. Res. 67, 520–528. doi: 10.1016/j.cardiores.2005.03.007
Ellinwood, N., Dobrev, D., Morotti, S., and Grandi, E. (2017a). In silico assessment of efficacy and safety of Ikur inhibitors in chronic atrial fibrillation: role of kinetics and state-dependence of drug binding. Front. Pharmacol. 8:799. doi: 10.3389/fphar.2017.00799
Ellinwood, N., Dobrev, D., Morotti, S., and Grandi, E. (2017b). Revealing kinetics and state-dependent binding properties of IKur-targeting drugs that maximize atrial fibrillation selectivity. Chaos 27:93918. doi: 10.1063/1.5000226
Ellis, W. S., Sippensgroenewegen, A., Auslander, D. M., and Lesh, M. D. (2000). The role of the crista terminalis in atrial flutter and fibrillation: a computer modeling study. Ann. Biomed. Eng. 28, 742–754. doi: 10.1114/1.1289456
Felix, J., Jacquemet, V., Alcaraz, R., and Rieta, J. J. (2015). Study on the trustability of phase mapping methods to represent atrial potentials in atrial fibrillation. Comput. Cardiol. 42, 49–52. doi: 10.1109/CIC.2015.7408583
Ferrer, A., Sebastián, R., Sánchez-Quintana, D., Rodríguez, J. F., Godoy, E. J., Martínez, L., et al. (2015). Detailed anatomical and electrophysiological models of human atria and torso for the simulation of atrial activation. PLoS One 10:e0141573. doi: 10.1371/journal.pone.0141573
Figuera, C., Suarez-Gutierrez, V., Barquero-Perez, O., Goya-Esteban, R., Rodrigo, M., Hernandez, I., et al. (2016a). “Performance of inverse problem regularization methods for driver location during atrial fibrillation,” in Computing in Cardiology, (Vancouver, BC: IEEE Computer Society), 693–696.
Figuera, C., Suárez-Gutiérrez, V., Hernández-Romero, I., Rodrigo, M., Liberos, A., Atienza, F., et al. (2016b). Regularization techniques for ECG imaging during atrial fibrillation: A computational study. Front. Physiol. 7:466. doi: 10.3389/fphys.2016.00466
Filos, D., Chouvarda, I., Tachmatzidis, D., Vassilikos, V., and Maglaveras, N. (2017). Beat-to-beat P-wave morphology as a predictor of paroxysmal atrial fibrillation. Comput. Methods Programs Biomed. 151, 111–121. doi: 10.1016/j.cmpb.2017.08.016
Fink, M., Niederer, S. A., Cherry, E. M., Fenton, F. H., Koivumäki, J. T., Seemann, G., et al. (2011). Cardiac cell modelling: observations from the heart of the cardiac physiome project. Prog. Biophys. Mol. Biol. 104, 2–21. doi: 10.1016/j.pbiomolbio.2010.03.002
Ganesan, P., Salmin, A., Cherry, E. M., and Ghoraani, B. (2016). “Development of a novel probabilistic algorithm for localization of rotors during atrial fibrillation,” in Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, EMBS, (New Jersey: Institute of Electrical and Electronics Engineers Inc.), 493–496.
Garcia-Gudino, D., Landa, E., Mendoza-Temis, J., Albarado-Ibanez, A., Toledo-Roy, J. C., Morales, I. O., et al. (2017). Enhancement of early warning properties in the Kuramoto model and in an atrial fibrillation model due to an external perturbation of the system. PLoS One 12:e0181953. doi: 10.1371/journal.pone.0181953
Gaspo, R., Bosch, R. F., Talajic, M., and Nattel, S. (1997). Functional mechanisms underlying tachycardia-induced sustained atrial fibrillation in a chronic dog model. Circulation 96, 4027–4035. doi: 10.1161/01.CIR.96.11.4027
Gharaviri, A., Verheule, S., Eckstein, J., Potse, M., Kuijpers, N. H. L., and Schotten, U. (2012a). A computer model of endo-epicardial electrical dissociation and transmural conduction during atrial fibrillation. Europace 14, 10–16. doi: 10.1093/europace/eus270
Gharaviri, A., Verheule, S., Kuijpers, N., and Schotten, U. (2012b). “Mutual influence between dyssynchrony and transmural conduction maintains atrial fibrillation,” in Proceedings of the Computing in Cardiology, (Krakow: IEEE), 897–900.
Gharaviri, A., Verheule, S., Eckstein, J., Potse, M., Kuklik, P., Kuijpers, N. H. L., et al. (2017). How disruption of endo-epicardial electrical connections enhances endo-epicardial conduction during atrial fibrillation. Europace 19, 308–318. doi: 10.1093/europace/euv445
Glynn, P., Onal, B., and Hund, T. J. (2014). Cycle length restitution in sinoatrial node cells: a theory for understanding spontaneous action potential dynamics. PLoS One 9:e89049. doi: 10.1371/journal.pone.0089049
Gomez, R., Nunez, L., Caballero, R., Vaquero, M., Tamargo, J., and Delpon, E. (2005). Spironolactone and its main metabolite canrenoic acid block hKv1.5. Kv4.3 and Kv7.1 + minK channels. Br. J. Pharmacol. 146, 146–161. doi: 10.1038/sj.bjp.0706302
Gong, Y., Xie, F., Stein, K. M., Garfinkel, A., Culianu, C. A., Lerman, B. B., et al. (2007). Mechanism underlying initiation of paroxysmal atrial flutter/atrial fibrillation by ectopic foci: A simulation study. Circulation 115, 2094–2102. doi: 10.1161/CIRCULATIONAHA.106.656504
Gonzales, M. J., Vincent, K. P., Rappel, W. J., Narayan, S. M., and McCulloch, A. D. (2014). Structural contributions to fibrillatory rotors in a patient-derived computational model of the atria. Europace 16, (Suppl. 4), iv3–iv10. doi: 10.1093/europace/euu251
Goodman, A. M., Oliver, R. A., Henriquez, C. S., and Wolf, P. D. (2005). A membrane model of electrically remodelled atrial myocardium derived from in vivo measurements. Europace 7, S135–S145. doi: 10.1016/j.eupc.2005.04.010
Grandi, E., and Maleckar, M. M. (2016). Anti-arrhythmic strategies for atrial fibrillation: The role of computational modeling in discovery, development, and optimization. Pharmacol. Ther. 168, 126–142. doi: 10.1016/j.pharmthera.2016.09.012
Grandi, E., Pandit, S. V., Voigt, N., Workman, A. J., Dobrev, D., Jalife, J., et al. (2011). Human atrial action potential and Ca 2+ model: Sinus rhythm and chronic atrial fibrillation. Circ. Res. 109, 1055–1066. doi: 10.1161/CIRCRESAHA.111.253955
Green, H. D., Thomas, G., and Terry, J. R. (2017). Signal reconstruction of pulmonary vein recordings using a phenomenological mathematical model: application to pulmonary vein isolation therapy. Front. Physiol. 8:496. doi: 10.3389/fphys.2017.00496
Guillem, M. S., Climent, A. M., Rodrigo, M., Hernandez-Romero, I., Liberos, A., Fernandez-Aviles, F., et al. (2016). “Noninvasive identification of atrial fibrillation drivers: Simulation and patient data evaluation,” in Computing in Cardiology, (Vancouver, BC: IEEE Computer Society), 121–124.
Haïssaguerre, M., Jaïs, P., Shah, D. C., Takahashi, A., Hocini, M., Quiniou, G., et al. (1998). Spontaneous initiation of atrial fibrillation by ectopic beats originating in the pulmonary veins. N. Engl. J. Med. 339, 659–666. doi: 10.1056/NEJM199809033391003
Hancox, J. C., Kharche, S., El Harchi, A., Stott, J., Law, P., and Zhang, H. (2014). In silico investigation of a KCNQ1 mutation associated with familial atrial fibrillation. J. Electrocardiol. 47, 158–165. doi: 10.1016/j.jelectrocard.2013.12.004
Henry, W. L., Morganroth, J., Pearlman, A. S., Clark, C. E., Redwood, D. R., Itscoitz, S. B., et al. (1976). Relation between echocardiographically determined left atrial size and atrial fibrillation. Circulation 53, 273–279. doi: 10.1161/01.CIR.53.2.273
Hong, K., Piper, D. R., Diaz-Valdecantos, A., Brugada, J., Oliva, A., Burashnikov, E., et al. (2005). De novo KCNQ1 mutation responsible for atrial fibrillation and short QT syndrome in utero. Cardiovasc. Res. 68, 433–440. doi: 10.1016/j.cardiores.2005.06.023
Hummel, J. P., Baher, A., Buck, B., Fanarjian, M., Webber, C. L. J., and Akar, J. G. (2017). A method for quantifying recurrent patterns of local wavefront direction during atrial fibrillation. Comput. Biol. Med. 89, 497–504. doi: 10.1016/j.compbiomed.2017.08.027
Hunter, P., and Borg, T. (2003). Integration from protiens to organs: the Physiome Project. Nature 4, 237–243. doi: 10.1038/nrm1054
Huo, Y., Holmqvist, F., Carlson, J., Gaspar, T., Hindricks, G., Piorkowski, C., et al. (2015). Variability of P-wave morphology predicts the outcome of circumferential pulmonary vein isolation in patients with recurrent atrial fibrillation. J. Electrocardiol. 48, 218–225. doi: 10.1016/j.jelectrocard.2014.11.011
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
Hwang, M., Lee, H.-S., Pak, H.-N., and Shim, E. B. (2016). Inducibility of human atrial fibrillation in an in silico model reflecting local acetylcholine distribution and concentration. Korean J. Physiol. Pharmacol. 20, 111–117. doi: 10.4196/kjpp.2016.20.1.111
Hwang, M., Lim, B., Song, J.-S. J. S., Yu, H. T., Ryu, A.-J. A. J., Lee, Y. S. Y.-S., et al. (2017). Ganglionated plexi stimulation induces pulmonary vein triggers and promotes atrial arrhythmogenecity: in silico modeling study. PLoS One 12:e0172931. doi: 10.1371/journal.pone.0172931
Ihara, Z., Jacquemet, V., Vesin, J.-M., and Van Oosterom, A. (2005). “Adaptation of the standard 12-lead EGG system focusing on atrial electrical activity,” in Proceedings of the Computers in Cardiology, (Lyon: IEEE), 203–205.
Ihara, Z., van Oosterom, A., Jacquemet, V., and Hoekema, R. (2007). Adaptation of the standard 12-lead electrocardiogram system dedicated to the analysis of atrial fibrillation. J. Electrocardiol. 40, .e1–.e68. doi: 10.1016/j.jelectrocard.2006.04.006
Iyengar, N., Peng, C. K., Morin, R., Goldberger, A. L., and Lipsitz, L. A. (1996). Age-related alterations in the fractal scaling of cardiac interbeat interval dynamics. Am. J. Physiol. 271, R1078–R1084. doi: 10.1152/ajpregu.1996.271.4.R1078
Jacquemet, V. (2015). Modeling left and right atrial contributions to the ECG: A dipole-current source approach. Comput. Biol. Med. 65, 192–199. doi: 10.1016/j.compbiomed.2015.06.007
Jacquemet, V. (2016). Lessons from computer simulations of ablation of atrial fibrillation. J. Physiol. 594, 2417–2430. doi: 10.1113/JP271660
Jacquemet, V., and Henriquez, C. S. (2009). Genesis of complex fractionated atrial electrograms in zones of slow conduction: a computer model of microfibrosis. Hear. Rhythm 6, 803–810. doi: 10.1016/j.hrthm.2009.02.026
Jacquemet, V., Lemay, M., Uldry, L., Duchêne, C., Van Oosterom, A., Kappenberger, L., et al. (2009). The role of atrial modeling in the development of ECG processing tools. IFMBE 25, 429–432. doi: 10.1007/978-3-642-03882-2-114
Jacquemet, V., Lemay, M., Van Oosterom, A., and Kappenberger, L. (2006). The equivalent dipole used to characterize atrial fibrillation. Comput. Cardiol. 33, 149–152.
Jacquemet, V., Lemay, M., Vesin, J. M., Van Oosterom, A., and Kappenberger, L. (2005). A biophysical model of ECG signals during atrial fibrillation used to evaluate the performance of QRST cancellation algorithms. Comput. Cardiol. 2005, 343–346. doi: 10.1109/CIC.2005.1588107
Jacquemet, V., Virag, N., Ihara, Z., Dang, L., Blanc, O., Zozor, S., et al. (2003). Study of Unipolar Electrogram Morphology in a Computer Model of Atrial Fibrillation. J. Cardiovasc. Electrophysiol. 14, S172–S179.
Jones, G., Spencer, B. D., Adeniran, I., and Zhang, H. (2012). Development of biophysically detailed electrophysiological models for pacemaking and non-pacemaking human pulmonary vein cardiomyocytes. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2012, 199–202. doi: 10.1109/EMBC.2012.6345905
Katritsis, D. G., Pokushalov, E., Romanov, A., Giazitzoglou, E., Siontis, G. C. M., Po, S. S., et al. (2013). Autonomic denervation added to pulmonary vein isolation for paroxysmal atrial fibrillation: A randomized clinical trial. J. Am. Coll. Cardiol. 62, 2318–2325. doi: 10.1016/j.jacc.2013.06.053
Kharche, S., Adeniran, I., Stott, J., Law, P., Boyett, M. R., Hancox, J. C., et al. (2012a). Pro-arrhythmogenic effects of the S140G KCNQ1 mutation in human atrial fibrillation - insights from modelling. J. Physiol. 590, 4501–4514. doi: 10.1113/jphysiol.2012.229146
Kharche, S., Biktasheva, I., Seemann, G., Zhang, H. G., and Biktashev, V. (2012b). Cardioversion using feedback stimuli in human atria. Comput. Cardiol. 39, 133–136.
Kharche, S., Castro, S., Thomas, B., Colman, M., Jarvis, J., Smaill, B., et al. (2014a). “Role of fiber orientation in atrial arrythmogenesis,” in Computing in Cardiology, (Washington, D.C: IEEE Computer Society), 1041–1044.
Kharche, S. R., Stary, T., Colman, M. A., Biktasheva, I. V., Workman, A. J., Rankin, A. C., et al. (2014b). Effects of human atrial ionic remodelling by β-blocker therapy onmechanisms of atrial fibrillation: A computer simulation. Europace 16, 1524–1533. doi: 10.1093/europace/euu084
Kharche, S., Garratt, C. J., Boyett, M. R., Inada, S., Holden, A. V., Hancox, J. C., et al. (2008). Atrial proarrhythmia due to increased inward rectifier current (IK1) arising from KCNJ2 mutation - A simulation study. Prog. Biophys. Mol. Biol. 98, 186–197. doi: 10.1016/j.pbiomolbio.2008.10.010
Kharche, S. R., Biktasheva, I. V., Seemann, G., Zhang, H., Zhao, J., and Biktashev, V. N. (2015). Computational modelling of low voltage resonant drift of scroll waves in the realistic human atria. Lect. Notes Comput. Sci. 9126, 421–429. doi: 10.1007/978-3-319-20309-6_48
Kharche, S. R., Biktasheva, I. V., Zhang, H., and Biktashev, V. N. (2013). Simulating the role of anisotropy in human atrial cardioversion. Conf. IEEE Eng. Med. Biol. Soc. 2013, 6838–6841. doi: 10.1109/EMBC.2013.6611128
Kharche, S. R., Vigmond, E., Efimov, I. R., and Dobrzynski, H. (2017). Computational assessment of the functional role of sinoatrial node exit pathways in the human heart. PLoS One 12:e0183727. doi: 10.1371/journal.pone.0183727
Ki, C.-S., Jung, C. L., Kim, H.-J., Baek, K.-H., Park, S. J., On, Y. K., et al. (2014). A KCNQ1 mutation causes age-dependant bradycardia and persistent atrial fibrillation. Pflugers Arch. Eur. J. Physiol. 466, 529–540. doi: 10.1007/s00424-013-1337-1336
Kirchhof, P., Benussi, S., Kotecha, D., Ahlsson, A., Atar, D., Casadei, B., et al. (2016). 2016 ESC Guidelines for the management of atrial fibrillation developed in collaboration with EACTS. Eur. Heart J. 37, 2893–2962. doi: 10.1093/eurheartj/ehw210
Kneller, J., Kalifa, J., Zou, R., Zaitsev, A. V., Warren, M., Berenfeld, O., et al. (2005). Mechanisms of atrial fibrillation termination by pure sodium channel blockade in an ionically-realistic mathematical model. Circ. Res. 96, e35–e47.
Kneller, J., Zou, R., Vigmond, E. J., Wang, Z., Leon, L. J., and Nattel, S. (2002). Cholinergic atrial fibrillation in a computer model of a two-dimensional sheet of canine atrial cells with realistic ionic properties. Circ. Res. 90, E73–E87.
Koivumaki, J. T., Clark, R. B., Belke, D., Kondo, C., Fedak, P. W. M., Maleckar, M. M. C., et al. (2014). Na(+) current expression in human atrial myofibroblasts: identity and functional roles. Front. Physiol. 5:275. doi: 10.3389/fphys.2014.00275
Koivumäki, J. T., Korhonen, T., and Tavi, P. (2011). Impact of sarcoplasmic reticulum calcium release on calcium dynamics and action potential morphology in human atrial myocytes: A computational study. PLoS Comput. Biol. 7:e1001067. doi: 10.1371/journal.pcbi.1001067
Koivumäki, J. T., Seemann, G., Maleckar, M. M., and Tavi, P. (2014). In silico screening of the key cellular remodeling targets in chronic atrial fibrillation. PLoS Comput. Biol. 10:e1003620. doi: 10.1371/journal.pcbi.1003620
Krogh-Madsen, T., Abbott, G. W., and Christini, D. J. (2012). Effects of electrical and structural remodeling on atrial fibrillation maintenance: a simulation study. PLoS Comput. Biol. 8:e1002390. doi: 10.1371/journal.pcbi.1002390
Krueger, M. W., Dorn, A., Keller, D. U. J., Holmqvist, F., Carlson, J., Platonov, P. G., et al. (2013a). In-silico modeling of atrial repolarization in normal and atrial fibrillation remodeled state. Med. Biol. Eng. Comput. 51, 1105–1119. doi: 10.1007/s11517-013-1090-1
Krueger, M. W., Schulze, W. H. W., Rhode, K. S., Razavi, R., Seemann, G., and Dössel, O. (2013b). Towards personalized clinical in-silico modeling of atrial anatomy and electrophysiology. Med. Biol. Eng. Comput. 51, 1251–1260. doi: 10.1007/s11517-012-0970-0
Krueger, M. W., Seemann, G., Rhode, K., Keller, D. U. J., Schilling, C., Arujuna, A., et al. (2013c). Personalization of atrial anatomy and electrophysiology as a basis for clinical modeling of radio-frequency ablation of atrial fibrillation. IEEE Trans. Med. Imaging 32, 73–84. doi: 10.1109/TMI.2012.2201948
Krueger, M. W., Seemann, G., Rhode, K. S., Weber, F. M., Linton, N., Williams, S., et al. (2013d). Fusion of local activation time maps and image data to personalize anatomical atrial models. Lect. Notes Comput. Sci. 7945, 1–10. doi: 10.1007/978-3-642-38899-6_1
Krueger, M. W., Rhode, K. S., O’Neill, M. D., Rinaldi, C. A., Gill, J., Razavi, R., et al. (2014). Patient-specific modeling of atrial fibrosis increases the accuracy of sinus rhythm simulations and may explain maintenance of atrial fibrillation. J. Electrocardiol. 47, 324–328. doi: 10.1016/j.jelectrocard.2013.11.003
Krueger, M. W., Severi, S., Rhode, K., Genovesi, S., Weber, F. M., Vincenti, A., et al. (2011). Alterations of atrial electrophysiology related to hemodialysis session: Insights from a multiscale computer model. J. Electrocardiol. 44, 176–183. doi: 10.1016/j.jelectrocard.2010.11.016
Krummen, D. E., Bayer, J. D., Ho, J., Ho, G., Smetak, M. R., Clopton, P., et al. (2012). Mechanisms of human atrial fibrillation initiation clinical and computational studies of repolarization restitution and activation latency. Circ. Arrhythmia Electrophysiol. 5, 1149–1159. doi: 10.1161/CIRCEP.111.969022
Kuijpers, N., Ten Eikelder, H., and Verheule, S. (2009). Atrial anatomy influences onset and termination of atrial fibrillation: a computer model study. Lect. Notes Comput. Sci. 5528, 285–294. doi: 10.1007/978-3-642-01932-6_31
Kuijpers, N. H. L., Keldermann, R. H., Ten Eikelder, H. M. M., Arts, T., and Hilbers, P. A. J. (2006). The role of the hyperpolarization-activated inward current If in arrhythmogenesis: A computer model study. IEEE Trans. Biomed. Eng. 53, 1499–1511. doi: 10.1109/TBME.2006.877801
Kuijpers, N. H. L., Potse, M., Van Dam, P. M., Ten Eikelder, H. M. M., Verheule, S., Prinzen, F. W., et al. (2011). Mechanoelectrical coupling enhances initiation and affects perpetuation of atrial fibrillation during acute atrial dilation. Hear. Rhythm 8, 429–436. doi: 10.1016/j.hrthm.2010.11.020
Kuijpers, N. H. L., Ten Eikelder, H. M. M., Bovendeerd, P. H. M., Verheule, S., Arts, T., and Hilbers, P. A. J. (2007). Mechanoelectric feedback leads to conduction slowing and block in acutely dilated atria: A modeling study of cardiac electromechanics. Am. J. Physiol. - Hear. Circ. Physiol. 292, H2832–H2853. doi: 10.1152/ajpheart.00923.2006
Kuklik, P., Schäffer, B., Hoffmann, B. A., Ganesan, A. N., Schreiber, D., Moser, J. M., et al. (2016). Local electrical dyssynchrony during atrial fibrillation: theoretical considerations and initial catheter ablation results. PLoS One 11:e0164236. doi: 10.1371/journal.pone.0164236
Kwon, S.-S., Yun, Y. H., Hong, S.-B., Pak, H.-N., and Shim, E. B. (2013). A patient-specific model of virtual ablation for atrial fibrillation. Conf. IEEE Eng. Med. Biol. Soc. 2013, 1522–1525. doi: 10.1109/EMBC.2013.6609802
Law, P., Kharche, S., Higham, J., and Zhang, H. (2010). Anti-arrhythmic effects of atrial specific ikur block: A simulation study. Comput. Cardiol. 37, 429–432.
Law, P., Kharche, S., Stott, J., and Zhang, H. (2009). Effects of elevated homocysteine hormone on electrical activity in the human atrium: A simulation study. in proceedings of the 31st annual international conference of the IEEE engineering in medicine and biology society: engineering the future of biomedicine. EMBC 2009, 3936–3939. doi: 10.1109/IEMBS.2009.5333530
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
Lemay, M., Jacquemet, V., Duchêne, C., Van Oosterom, A., Abächerli, R., and Vesin, J. M. (2009). “Activity level of an atrial ectopic focus observed through the atrial vectorcardiogram: A biophysical model,” in Proceedings of the 36th Annual Computers in Cardiology Conference (CinC), (Park City, UT: IEEE), 221–224.
Lemay, M., Jacquemet, V., Jousset, F., Vesin, J. M., and Van Oosterom, A. (2007a). “The mean firing rate of atrial fibrillation as estimated from the ECG evaluation using a biophysical model,” in Proceedings of the Computers in Cardiology, (Durham, NC: IEEE), 37–40.
Lemay, M., Vesin, J. M., Jacquemet, V., Forclaz, A., Kappenberger, L., and van Oosterom, A. (2007b). Spatial dynamics of atrial activity assessed by the vectorcardiogram: from sinus rhythm to atrial fibrillation. Europace 9, (Suppl. 6), vi109–vi118. doi: 10.1093/europace/eum215
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
Li, J., Inada, S., Schneider, J. E., Zhang, H., Dobrzynski, H., and Boyett, M. R. (2014). Three-dimensional computer model of the right atrium including the sinoatrial and atrioventricular nodes predicts classical nodal behaviours. PLoS One 9:e112547. doi: 10.1371/journal.pone.0112547
Li, Q., O’Neill, S. C., Tao, T., Li, Y., Eisner, D., and Zhang, H. (2012). Mechanisms by which cytoplasmic calcium wave propagation and alternans are generated in cardiac atrial myocytes lacking T-tubules-insights from a simulation study. Biophys. J. 102, 1471–1482. doi: 10.1016/j.bpj.2012.03.007
Liberos, A., Bueno-Orovio, A., Rodrigo, M., Ravens, U., Hernandez-Romero, I., Fernandez-Aviles, F., et al. (2016). Balance between sodium and calcium currents underlying chronic atrial fibrillation termination: an in silico intersubject variability study. Hear. Rhythm 13, 2358–2365. doi: 10.1016/j.hrthm.2016.08.028
Lin, Y. T., Chang, E. T., Eatock, J., Galla, T., and Clayton, R. H. (2017). Mechanisms of stochastic onset and termination of atrial fibrillation studied with a cellular automaton model. J. R. Soc. Inter. 14:20160968. doi: 10.1098/rsif.2016.0968
Liu, J., Holden, A. V., and Zhang, H. (2007). Actions of an external electrical shock on human atrial excitation - a computer model study. Lect. Notes Comput. Sci. 4561, 659–667. doi: 10.1007/978-3-540-73321-8_76
Lo, C.-P., Horng, T.-L., Luk, H.-N., Tien, H.-C., Yuan, J.-M., and Lee, D. (2006). A computational model of the canine pulmonary veins sleeves. J. Med. Biol. Eng. 26, 125–130. doi: 10.1109/EMBC.2012.6345906
Loewe, A., Krueger, M. W., Holmqvist, F., Dossel, O., Seemann, G., and Platonov, P. G. (2016). Influence of the earliest right atrial activation site and its proximity to interatrial connections on P-wave morphology. Europace 18, iv35–iv43. doi: 10.1093/europace/euw349
Loewe, A., Wilhelms, M., Fischer, F., Scholz, E. P., Dössel, O., and Seemann, G. (2014). Arrhythmic potency of human ether-à-go-go-related gene mutations L532P and N588K in a computational model of human atrial myocytes. Europace 16, 435–443. doi: 10.1093/europace/eut375
Lombardo, D. M., Fenton, F. H., Narayan, S. M., and Rappel, W. J. (2016). Comparison of detailed and simplified models of human atrial myocytes to recapitulate patient specific properties. PLoS Comput. Biol. 12:e1005060. doi: 10.1371/journal.pcbi.1005060
Luca, A., Jacquemet, V., Virag, N., and Vesin, J.-M. (2015). “Influence of right and left atrial tissue heterogeneity on atrial fibrillation perpetuation,” in Computing in Cardiology, (Washington, DC: IEEE Computer Society), 449–452.
Luca, A., Kallmyer, T., and Virag, N. (2016). Atrial fibrillation septal pacing: translation of modelling results. Europace 18, iv53–iv59. doi: 10.1093/europace/euw360
Luo, C. H., and Rudy, Y. (1991). A model of the ventricular cardiac action potential. Circ. Res. 68, 1501–1526. doi: 10.1161/01.RES.68.6.1501
Lutz, Y., Loewe, A., Wilhelms, M., Dössel, O., and Seemann, G. (2014). Specific antiarrhythmic therapy for familial atrial fibrillation in a numerical model of human atrial electrophysiology. Biomed. Tech. 59, S933–S936. doi: 10.1515/bmt-2014-5012
Majumder, R., Jangsangthong, W., Feola, I., Ypey, D. L., Pijnappels, D. A., and Panfilov, A. V. (2016). A mathematical model of neonatal rat atrial monolayers with constitutively active acetylcholine-mediated K+ current. PLoS Comput. Biol. 12:e1004946. doi: 10.1371/journal.pcbi.1004946
Maleckar, M. M., Greenstein, J. L., Giles, W. R., and Trayanova, N. A. (2009). K+ current changes account for the rate dependence of the action potential in the human atrial myocyte. AJP Hear. Circ. Physiol. 297, H1398–H1410. doi: 10.1152/ajpheart.00411.2009
Manani, K. A., Christensen, K., and Peters, N. S. (2016). Myocardial architecture and patient variability in clinical patterns of atrial fibrillation. Phys. Rev. E. 94, 1–7. doi: 10.1103/PhysRevE.94.042401
Mann, S. A., Otway, R., Guo, G., Soka, M., Karlsdotter, L., Trivedi, G., et al. (2012). Epistatic effects of potassium channel variation on cardiac repolarization and atrial fibrillation risk. J. Am. Coll. Cardiol. 59, 1017–1025. doi: 10.1016/j.jacc.2011.11.039
Markides, V., and Schilling, R. J. (2003). Atrial fibrillation: classification, pathophysiology, mechanisms and drug treatment. Heart 89, 939–943. doi: 10.1136/heart.89.8.939
Marshall, G. E., Russell, J. A., Tellez, J. O., Jhund, P. S., Currie, S., Dempster, J., et al. (2012). Remodelling of human atrial K + currents but not ion channel expression by chronic β-blockade. Pflugers Arch. Eur. J. Physiol. 463, 537–548. doi: 10.1007/s00424-011-1061-z
Martínez, A., Alcaraz, R., and Rieta, J. J. (2012). Study on the P-wave feature time course as early predictors of paroxysmal atrial fibrillation. Physiol. Meas. 33, 1959–1974. doi: 10.1088/0967-3334/33/12/1959
Martínez, A., Alcaraz, R., and Rieta, J. J. (2014a). Morphological variability of the P-wave for premature envision of paroxysmal atrial fibrillation events. Physiol. Meas. 35, 1–14. doi: 10.1088/0967-3334/35/1/1
Martínez, L., Romero, L., Tobón, C., Ferrero, J. M., Jalife, J., Berenfeld, O., et al. (2014b). “Accurate characterization of rotor activity during atrial fibrillation depends on the properties of the multi-electrode grid,” in Proceedings of the Computing in Cardiology, (Washington, DC: IEEE Computer Society), 757–760.
Martinez, L., Romero, L., Ferrer, A., Jalife, J., Berenfeld, O., and Saiz, J. (2016). “Effect of multi-electrode configurations on accuracy of rotor detection in the atria,” in Proceedings of the Computing in Cardiology, (Washington, DC: IEEE Computer Society), 1085–1088.
McDowell, K. S., Vadakkumpadan, F., Blake, R., Blauer, J., Plank, G., MacLeod, R. S., et al. (2012). Methodology for patient-specific modeling of atrial fibrosis as a substrate for atrial fibrillation. J. Electrocardiol. 45, 640–645. doi: 10.1016/j.jelectrocard.2012.08.005
McDowell, K. S., Vadakkumpadan, F., Blake, R., Blauer, J., Plank, G., Macleod, R. S., et al. (2013). Mechanistic inquiry into the role of tissue remodeling in fibrotic lesions in human atrial fibrillation. Biophys. J. 104, 2764–2773. doi: 10.1016/j.bpj.2013.05.025
Michelucci, A., Bagliani, G., Colella, A., Pieragnoli, P., Porciani, M. C., Gensini, G., et al. (2002). P wave assessment: State of the art update. Card. Electrophysiol. Rev. 6, 215–220. doi: 10.1023/A:1016368723033
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-90370
Monfredi, O., Dobrzynski, H., Mondal, T., Boyett, M. R., and Morris, G. M. (2010). The anatomy and physiology of the sinoatrial node-A contemporary review. Pacing Clin. Electrophysiol. 33, 1392–1406. doi: 10.1111/j.1540-8159.2010.02838.x
Morgan, R., Colman, M., Kruger, M., Seemann, G., Rhode, K., and Aslanidi, O. (2014). “Evaluating effects of fibrosis in atrial arrhythmogenesis using 3D computational modelling,” in Computing in Cardiology, (Washington, DC: IEEE Computer Society), 765–768.
Morgan, R., Colman, M. A., Chubb, H., Seemann, G., and Aslanidi, O. V. (2016). Slow conduction in the border zones of patchy fibrosis stabilizes the drivers for atrial fibrillation: insights from multi-scale human atrial modeling. Front. Physiol. 7:474. doi: 10.3389/fphys.2016.00474
Morris, J. J., Estes, E. J., Whalen, R., Thompson, H. J., and McIntosh, H. (1964). P-wave analysis in valvular heart disease. Circulation 29, 242–252. doi: 10.1161/01.CIR.29.2.242
Namba, T., Ashihara, T., Nakazawa, K., and Ohe, T. (1999). Effect of pilsicainide, a pure sodium channel blocker, on spiral waves during atrial fibrillation: Theoretical analysis by numerical simulation. J. Electrocardiol. 32, 321–334. doi: 10.1016/S0022-0736(99)90003-90001
Nattel, S. (2002). New ideas about atrial fibrillation 50 years on. Nature 415, 219–226. doi: 10.1038/415219a
Nguyen, H. Q., Bailey, A., Coleman, K. J., Desai, S., Fan, V. S., Gould, M. K., et al. (2016). Patient-centered physical activity coaching in COPD (Walk On!): A study protocol for a pragmatic randomized controlled trial. Contemp. Clin. Trials 46, 18–29. doi: 10.1016/j.cct.2015.10.010
Ni, H., Adeniran, I., and Zhang, H. (2017a). In-silico investigations of the functional impact of KCNA5 mutations on atrial mechanical dynamics. J. Mol. Cell. Cardiol. 111, 86–95. doi: 10.1016/j.yjmcc.2017.08.005
Ni, H., Whittaker, D. G., Wang, W., Giles, W. R., Narayan, S. M., and Zhang, H. (2017b). Synergistic anti-arrhythmic effects in human atria with combined use of sodium blockers and acacetin. Front. Physiol. 8:946. doi: 10.3389/fphys.2017.00946
Ni, H., Colman, M. A., and Zhang, H. (2014). “Modelling the functional impact of KCNA5 mutations on the electrical and mechanical activities of human atrial cells,” in Computing in Cardiology, (Washington, DC: IEEE Computer Society), 57–60.
Ni, H., Whittaker, D. G., Wang, W., and Zhang, H. (2016). “Synergistic anti-arrhythmic effects of combined blockade of sodium and ultra-rapid delayed rectifier potassium channels in human atria,” in Proceedings of the Computing in Cardiology, (Washington, DC: IEEE Computer Society), 837–840.
Nishida, K., and Nattel, S. (2014). Atrial fibrillation compendium: Historical context and detailed translational perspective on an important clinical problem. Circ. Res. 114, 1447–1452. doi: 10.1161/CIRCRESAHA.114.303466
Nygren, A., Fiset, C., Firek, L., Clark, J. W., Lindblad, D. S., Clark, R. B., et al. (1998). Mathematical model of an adult human atrial cell: the role of K+ currents in repolarization. Circ. Res. 82, 63–81. doi: 10.1161/01.RES.82.1.63
O’Connell, R. P., Musa, H., Gomez, M. S. M., Avula, U. M., Herron, T. J., Kalifa, J., et al. (2015). Free fatty acid effects on the atrial myocardium: membrane ionic currents are remodeled by the disruption of t-tubular architecture. PLoS One 10:e0133052. doi: 10.1371/journal.pone.0133052
Ogawa, M., Kumagai, K., Vakulenko, M., Yasuda, T., Siegerman, C., Garfinkel, A., et al. (2007). Reduction of P-wave duration and successful pulmonary vein isolation in patients with atrial fibrillation. J. Cardiovasc. Electrophysiol. 18, 931–938. doi: 10.1111/j.1540-8167.2007.00890.x
Olgin, J. E., Sih, H. J., Hanish, S., Jayachandran, J. V., Wu, J., Zheng, Q. H., et al. (1998). Heterogeneous atrial denervation creates substrate for sustained atrial fibrillation. Circulation 98, 2608–2614. doi: 10.1161/01.CIR.98.23.2608
Onal, B., Gratz, D., and Hund, T. J. (2017). Ca(2+)/calmodulin-dependent kinase II-dependent regulation of atrial myocyte late Na(+) current, Ca(2+) cycling, and excitability: a mathematical modeling study. Am. J. Physiol. Heart Circ. Physiol. 313, H1227–H1239. doi: 10.1152/ajpheart.00185.2017
Pedrón-Torrecilla, J., Climent, A. M., Liberos, A., Rodrigo, M., Pérez-David, E., Millet, J., et al. (2014). “Accuracy of inverse solution computation of dominant frequencies and phases during atrial fibrillation,” in Proceedings of the Computing in Cardiology, (Washington, DC: IEEE Computer Society), 537–540.
Pedron-Torrecilla, J., Rodrigo, M., Climent, A. M., Liberos, A., Perez-David, E., Bermejo, J., et al. (2016). Noninvasive estimation of epicardial dominant high-frequency regions during atrial fibrillation. J. Cardiovasc. Electrophysiol. 27, 435–442. doi: 10.1111/jce.12931
Pellman, J., and Sheikh, F. (2015). Atrial fibrillation: mechanisms, therapeutics, and future directions. Compr. Physiol. 5, 649–665. doi: 10.1002/cphy.c140047
Peters, M. D. J., Godfrey, C. M., Khalil, H., McInerney, P., Parker, D., and Soares, C. B. (2015). Guidance for conducting systematic scoping reviews. Int. J. Evid. Based. Healthc. 13, 141–146. doi: 10.1097/XEB.0000000000000050
Platonov, P. G. (2012). P-wave morphology: underlying mechanisms and clinical implications. Ann. Noninvasive Electrocardiol. 17, 161–169. doi: 10.1111/j.1542-474X.2012.00534.x
Ramirez, R. J., Nattel, S., and Courtemanche, M. (2000). Mathematical analysis of canine atrial action potentials: Rate, regional factors, and electrical remodeling. Am. J. Physiol. - Hear. Circ. Physiol. 279, H1767–H1785.
Rappel, W., and Narayan, S. M. (2013). Theoretical considerations for mapping activation in human cardiac fibrillation. Chaos 23:023113. doi: 10.1063/1.4807098
Rappel, W.-J., Zaman, J. A. B., and Narayan, S. M. (2015). Mechanisms for the termination of atrial fibrillation by localized ablation: computational and clinical studies. Circ. Arrhythmia Electrophysiol. 8, 1325–1333. doi: 10.1161/CIRCEP.115.002956
Rensma, P. L., Allessie, M. A., Lammers, W. J. E. P., Bonke, F. I. M., and Schalij, M. J. (1988). Length of excitation wave and susceptibility to reentrant atrial arrhythmias in normal conscious dogs. Circ. Res. 62, 395–410. doi: 10.1161/01.RES.62.2.395
Reumann, M., Bohnert, J., and Doessel, O. (2006). Simulating pulmonary vein activity leading to atrial fibrillation using a rule-based approach on realistic anatomical data. Conf. Proc. IEEE Eng. Med. Biol. Soc. 1, 3943–3946. doi: 10.1109/IEMBS.2006.260383
Reumann, M., Bohnert, J., Osswald, B., Hagl, S., and Doessel, O. (2007). Multiple wavelets, rotors, and snakes in atrial fibrillation-a computer simulation study. J. Electrocardiol. 40, 328–334. doi: 10.1016/j.jelectrocard.2006.12.016
Reumann, M., Bohnert, J., Seemann, G., Osswald, B., and Dössel, O. (2008). Preventive ablation strategies in a biophysical model of atrial fibrillation based on realistic anatomical data. IEEE Trans. Biomed. Eng. 55, 399–406. doi: 10.1109/TBME.2007.912672
Rice, J. J., Wang, F., Bers, D. M., and De Tombe, P. P. (2008). Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations. Biophys. J. 95, 2368–2390. doi: 10.1529/biophysj.107.119487
Rodrigo, M., Clirnent, A. M., Liberos, A., Pedrón-Torrecilla, J., Millet, J., Fernández-Avilés, F., et al. (2014a). “Non-invasive detection of reentrant drivers during atrial fibrillation: A clinical-computational study,” in Proceedings of the Computing in Cardiology, (Washington, DC: IEEE Computer Society), 9–12.
Rodrigo, M., Guillem, M. S., Climent, A. M., Pedrón-Torrecilla, J., Liberos, A., Millet, J., et al. (2014b). Body surface localization of left and right atrial high-frequency rotors in atrial fibrillation patients: a clinical-computational study. Hear. Rhythm 11, 1584–1591. doi: 10.1016/j.hrthm.2014.05.013
Roney, C. H., Bayer, J. D., Zahid, S., Meo, M., Boyle, P. M. J., Trayanova, N. A., et al. (2016). Modelling methodology of atrial fibrosis affects rotor dynamics and electrograms. Europace 18, iv146–iv155. doi: 10.1093/europace/euw365
Roney, C. H., Cantwell, C. D., Bayer, J. D., Qureshi, N. A., Lim, P. B., Tweedy, J. H., et al. (2017). Spatial resolution requirements for accurate identification of drivers of atrial fibrillation. Circ. Arrhythm. Electrophysiol. 10:e004899. doi: 10.1161/CIRCEP.116.004899
Rotter, M., Dang, L., Jacquemet, V., Virag, N., Kappenberger, L., and Haïssaguerre, M. (2007). Impact of varying ablation patterns in a simulation model of persistent atrial fibrillation. Pacing Clin. Electrophysiol. 30, 314–321. doi: 10.1111/j.1540-8159.2007.00671.x
Rottmann, M., Unger, L., Kaltenbacher, W., Seemann, G., Loewe, A., Krueger, M. W., et al. (2015). “Methods for analyzing signal characteristics of stable and unstable rotors in a realistic heart model,” in Proceedings of the Computing in Cardiology, (Washington, DC: IEEE Computer Society), 485–488.
Ruchat, P., Dang, L., Schlaepfer, J., Virag, N., von Segesser, L. K., and Kappenberger, L. (2007a). Use of a biophysical model of atrial fibrillation in the interpretation of the outcome of surgical ablation procedures. Eur. J. Cardio-thoracic Surg. 32, 90–95. doi: 10.1016/j.ejcts.2007.02.031
Ruchat, P., Dang, L., Virag, N., Schlaepfer, J., von Segesser, L. K., and Kappenberger, L. (2007b). A biophysical model of atrial fibrillation to define the appropriate ablation pattern in modified maze. Eur. J. Cardio-thoracic Surg. 31, 65–69. doi: 10.1016/j.ejcts.2006.10.015
Ruchat, P., Virag, N., Dang, L., Schlaepfer, J., Pruvot, E., and Kappenberger, L. (2007c). A biophysical model of atrial fibrillation ablation: what can a surgeon learn from a computer model? Europace 9, (Suppl. 6), vi71–vi76. doi: 10.1093/europace/eum209
Rusu, A., Jacquemet, V., Vesin, J.-M., and Virag, N. (2014). Influence of atrial substrate on local capture induced by rapid pacing of atrial fibrillation. Europace 16, 766–773. doi: 10.1093/europace/euu003
Saha, M., Conte, G., Caputo, M. L., Regoli, F., Krause, R., Auricchio, A., et al. (2016). Changes in P-wave morphology after pulmonary vein isolation: insights from computer simulations. Europace 18, iv23–iv34. doi: 10.1093/europace/euw348
Sahli Costabal, F., Zaman, J. A. B., Kuhl, E., and Narayan, S. M. (2018). Interpreting activation mapping of atrial fibrillation: a hybrid computational/physiological study. Ann. Biomed. Eng. 46, 257–269. doi: 10.1007/s10439-017-1969-3
Salmin, A. J., Ganesan, P., Shillieto, K. E., Cherry, E. M., Huang, D. T., Pertsov, A. M., et al. (2016). A novel catheter-guidance algorithm for localization of atrial fibrillation rotor and focal sources. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2016, 501–504. doi: 10.1109/EMBC.2016.7590749
Sampson, K. J., Terrenoire, C., Cervantes, D. O., Kaba, R. A., Peters, N. S., and Kass, R. S. (2008). Adrenergic regulation of a key cardiac potassium channel can contribute to atrial fibrillation: Evidence from an IKs transgenic mouse. J. Physiol. 586, 627–637. doi: 10.1113/jphysiol.2007.141333
Sanchez, C., Bueno-Orovio, A., Pueyo, E., and Rodriguez, B. (2017). Atrial fibrillation dynamics and ionic block effects in six heterogeneous human 3d virtual atria with distinct repolarization dynamics. Front. Bioeng. Biotechnol. 5:29. doi: 10.3389/fbioe.2017.00029
Sanchez, C., Corrias, A., Bueno-Orovio, A., Davies, M., Swinton, J., Jacobson, I., et al. (2012). The Na+/K+ pump is an important modulator of refractoriness and rotor dynamics in human atrial tissue. Am. J. Physiol. Heart Circ. Physiol. 302, H1146–H1159. doi: 10.1152/ajpheart.00668.2011
Sanders, P., Morton, J. B., Davidson, N. C., Spence, S. J., Vohra, J. K., Sparks, P. B., et al. (2003). Electrical remodeling of the atria in congestive heart failure: Electrophysiological and electroanatomic mapping in humans. Circulation 108, 1461–1468. doi: 10.1161/01.CIR.0000090688.49283.67
Scholz, E. P., Carrillo-Bustamante, P., Fischer, F., Wilhelms, M., Zitron, E., Dössel, O., et al. (2013). Rotor termination is critically dependent on kinetic properties of I Kur inhibitors in an In Silico model of chronic atrial fibrillation. PLoS One 8:e83179. doi: 10.1371/journal.pone.0083179
Seemann, G., Carillo, P., Weiss, D. L., Krueger, M. W., Dössel, O., and Scholz, E. P. (2009). Investigating arrhythmogenic effects of the herg mutation n588k in virtual human atria. Lect. Notes Comput. Sci. 5528, 144–153. doi: 10.1007/978-3-642-01932-6_16
Seemann, G., Weiß, D. L., Sachse, F. B., and Dössel, O. (2004). “Familial atrial fibrillation: Simulation of the mechanisms and effects of a slow rectifier potassium channel mutation in human atrial tissue,” in Proceedings of the Computers in Cardiology, (Chicago, IL: IEEE), 125–128.
Seemann, G., Ying, H., Weiss, D. L., Sachse, F. B., and Dössel, O. (2005). “Effects of electrophysiological remodeling in human right atrium: A simulation study,” in Proceedings of the Computers in Cardiology, (Lyon: IEEE), 69–72. doi: 10.1109/CIC.2005.1588036
Seol, C. A., Kim, J., Kim, W. T., Ha, J. M., Choe, H., Jang, Y. J., et al. (2008). Simulation of spontaneous action potentials of cardiomyocytes in pulmonary veins of rabbits. Prog. Biophys. Mol. Biol. 96, 132–151. doi: 10.1016/j.pbiomolbio.2007.07.013
Severi, S., Fantini, G., Corsi, C., Vincenti, A., and Genovesi, S. (2008). Role of hemodialysis in atrial fibrillation onset: preliminary results from a combined computational and experimental analysis. Comput. Cardiol. 35, 877–879. doi: 10.1109/CIC.2008.4749182
Severi, S., Pogliani, D., Fantini, G., Fabbrini, P., Viganò, M. R., Galbiati, E., et al. (2010). Alterations of atrial electrophysiology induced by electrolyte variations: Combined computational and P-wave analysis. Europace 12, 842–849. doi: 10.1093/europace/euq042
Shen, M. J., Choi, E. K., Tan, A. Y., Han, S., Shinohara, T., Maruyama, M., et al. (2011). Patterns of baseline autonomic nerve activity and the development of pacing-induced sustained atrial fibrillation. Hear. Rhythm 8, 583–589. doi: 10.1016/j.hrthm.2010.11.040
Shusterman, V., Aysin, B., Ermentrout, G. B., London, B., and Schwartzman, D. (2003). Detecting instabilities of cardiac rhythm. J. Electrocardiol. 36(Suppl.), 219–226. doi: 10.1016/j.jelectrocard.2003.09.063
Skibsbye, L., Jespersen, T., Christ, T., Maleckar, M. M., van den Brink, J., Tavi, P., et al. (2016). Refractoriness in human atria: time and voltage dependence of sodium channel availability. J. Mol. Cell. Cardiol. 101, 26–34. doi: 10.1016/j.yjmcc.2016.10.009
Syed, Z. F., Vigmond, E., and Leon, L. J. (2005). Suitability of genetic algorithm generated models to simulate atrial fibrillation and K+ channel blockades. Conf. Proc. IEEE Eng. Med. Biol. Soc. 7, 7087–7090.
Syeda, F., Holmes, A. P., Yu, T. Y., Tull, S., Kuhlmann, S. M., Pavlovic, D., et al. (2016). PITX2 modulates atrial membrane potential and the antiarrhythmic effects of sodium-channel blockers. J. Am. Coll. Cardiol. 68, 1881–1894. doi: 10.1016/j.jacc.2016.07.766
Tobón, C., Cardona, E. A., Palacio, L. C., Duque, J. E., Ugarte, J. P., Orozco-Duque, A., et al. (2015). “Ablation pattern guided by approximate entropy maps to prevent chronic atrial fibrillation: A simulation study,” in IFMBE Proceedings, ed. H. A. Braidot (Berlin: Springer Verlag), 560–563. doi: 10.1007/978-3-319-13117-7_143
Tobón, C., Duarte, M., Duque, J. E., Becerra, M. A., Arango, S. S., Cardona, K., et al. (2014a). “Chloroquine effect on rotor termination under paroxysmal and chronic atrial fibrillation. 2D simulation study,” in Proceedings of the 2014 IEEE Central America and Panama Convention, CONCAPAN 2014, (Panama: Institute of Electrical and Electronics Engineers Inc.).
Tobón, C., Palacio, L. C., Duque, J. E., Cardona, E. A., Ugarte, J. P., Orozco-Duque, A., et al. (2014b). “Simple ablation guided by ApEn mapping in a 2D model during permanent atrial fibrillation,” in Computing in Cardiology, (Washington, DC: IEEE Computer Society), 1029–1032.
Tobón, C., Pachajoa, D., Ugarte, J. P., Saiz, J., Pérez, S., Ugarte, J. P., et al. (2017). “Dofetilide effect on human atrial action potential under normal and atrial fibrillation conditions,” in Silico Study. in IFMBE Proceedings, eds T. I. Bustamante and J. Sierra (Berlin: Springer Verlag), 38–41. doi: 10.1007/978-981-10-4086-3_10
Tobón, C., Ruiz-Villa, C. A., Heidenreich, E., Romero, L., Hornero, F., and Saiz, J. (2013). A three-dimensional human atrial model with fiber orientation. Electrograms and arrhythmic activation patterns relationship. PLoS One 8:e50883. doi: 10.1371/journal.pone.0050883
Trayanova, N. A., and Chang, K. C. (2016). How computer simulations of the human heart can improve anti-arrhythmia therapy. J. Physiol. 594, 2483–2502. doi: 10.1113/JP270532
Trayanova, N. A., and McDowell, K. S. (2018). Personalized Computational Modeling of Atrial Fibrosis to Guide Catheter Ablation of Atrial Fibrillation. Baltimore: The Johns Hopkins University.
Tsujimae, K., Murakami, S., and Kurachi, Y. (2008). In silico study on the effects of IKur block kinetics on prolongation of human action potential after atrial fibrillation-induced electrical remodeling. Am. J. Physiol. Hear. Circ. Physiol. 294, H793–H800. doi: 10.1152/ajpheart.01229.2007
Ugarte, J. P., Orozco-Duque, A., Tobón, C., Kremen, V., Novak, D., Saiz, J., et al. (2014). Dynamic approximate entropy electroanatomic maps detect rotors in a simulated atrial fibrillation model. PLoS One 9:e114577. doi: 10.1371/journal.pone.0114577
Ugarte, J. P., Orozco-Duque, A., Tobon, C., Saiz, J., Morillo, C. A., and Bustamante, J. (2013). “Localization of complex fractionated atrial electrograms by approximate entropy in a 3D model of human atria,” in Proceedings of the Pan American Health Care Exchanges, PAHCE, (Medellin: IEEE).
Ugarte, J. P., Tobón, 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, ii97–ii104. doi: 10.1093/europace/euv244
Uldry, L., Jacquemet, V., Virag, N., Kappenberger, L., and Vesin, J.-M. (2012). Estimating the time scale and anatomical location of atrial fibrillation spontaneous termination in a biophysical model. Med. Biol. Eng. Comput. 50, 155–163. doi: 10.1007/s11517-011-0859-3
van Oosterom, A., Ihara, Z., Jacquemet, V., and Hoekema, R. (2007). Vectorcardiographic lead systems for the characterization of atrial fibrillation. J. Electrocardiol. 40, e1–e343. doi: 10.1016/j.jelectrocard.2006.08.002
van Oosterom, A., and Jacquemet, V. (2005). Genesis of the P wave: Atrial signals as generated by the equivalent double layer source model. Europace 7, S21–S29. doi: 10.1016/j.eupc.2005.05.001
Varela, M., Colman, M. A., Hancox, J. C., and Aslanidi, O. V. (2016). Atrial heterogeneity generates re-entrant substrate during atrial fibrillation and anti-arrhythmic drug action: mechanistic insights from canine atrial models. PLoS Comput. Biol. 12:e1005245. doi: 10.1371/journal.pcbi.1005245
Varela, M., Morgan, R., Ghavami, N., James, S., and Aslanidi, O. (2014). “Ionic mechanisms of triggered activity in atrial cell models,” in Computing in Cardiology, (Washington, DC: IEEE Computer Society), 1021–1024.
Veenhuyzen, G. D., Simpson, C. S., and Abdollah, H. (2004). Atrial fibrillation. CMAJ 171, 755–760. doi: 10.1503/cmaj.1031364
Verheule, S., Tuyls, E., Gharaviri, A., Hulsmans, S., Van Hunnik, A., Kuiper, M., et al. (2013). Loss of continuity in the thin epicardial layer because of endomysial fibrosis increases the complexity of atrial fibrillatory conduction. Circ. Arrhythmia Electrophysiol. 6, 202–211. doi: 10.1161/CIRCEP.112.975144
Viceconti, M., and Hunter, P. (2016). The virtual physiological human: ten years after. Annu. Rev. Biomed. Eng. 18, 103–123. doi: 10.1146/annurev-bioeng-110915-114742
Vigmond, E., Labarthe, S., Cochet, H., Coudiere, Y., Henry, J., and Jais, P. (2013). A bilayer representation of the human atria. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2013, 1530–1533. doi: 10.1109/EMBC.2013.6609804
Vigmond, E. J., Ruckdeschel, R., and Trayanova, N. (2001). Reentry in a morphologically realistic atrial model. J. Cardiovasc. Electrophysiol. 12, 1046–1054. doi: 10.1046/j.1540-8167.2001.01046.x
Vigmond, E. J., Tsoi, V., Kuo, S., Arevalo, H., Kneller, J., Nattel, S., et al. (2004). The effect of vagally induced dispersion of action potential duration on atrial arrhythmogenesis. Hear. Rhythm 1, 334–344. doi: 10.1016/j.hrthm.2004.03.077
Vincenti, A., Passini, E., Fabbrini, P., Luise, M. C., Severi, S., and Genovesi, S. (2014). Recurrent intradialytic paroxysmal atrial fibrillation: hypotheses on onset mechanisms based on clinical data and computational analysis. Europace 16, 396–404. doi: 10.1093/europace/eut346
Virag, N., Jacquemet, V., Henriquez, C. S., Zozor, S., Blanc, O., Vesin, J.-M., et al. (2002). Study of atrial arrhythmias in a computer model based on magnetic resonance images of human atria. Chaos 12, 754–763. doi: 10.1063/1.1483935
Virag, N., Jacquemet, V., and Kappenberger, L. (2012). Modeling of Atrial Fibrillation. New Jersey: Wiley-Blackwell.
Voigt, N., Heijman, J., Trausch, A., Mintert-Jancke, E., Pott, L., Ravens, U., et al. (2013). Impaired Na+-dependent regulation of acetylcholine-activated inward-rectifier K+ current modulates action potential rate dependence in patients with chronic atrial fibrillation. J. Mol. Cell. Cardiol. 61, 142–152. doi: 10.1016/j.yjmcc.2013.03.011
Voigt, N., Heijman, J., Wang, Q., Chiang, D. Y., Li, N., Karck, M., et al. (2014). Cellular and molecular mechanisms of atrial arrhythmogenesis in patients with paroxysmal atrial fibrillation. Circulation 129, 145–156. doi: 10.1161/CIRCULATIONAHA.113.006641
Weber, F. M., Schilling, C., Straub, D., Gurm, S., Seemann, G., Lorenz, C., et al. (2009a). Extracting clinically relevant circular mapping and coronary sinus catheter potentials from atrial simulations. Lect. Notes Comput. Sci. 5528, 30–38. doi: 10.1007/978-3-642-01932-6_4
Weber, F. M., Schilling, C., Straub, D., Seemann, G., Lorenz, C., and Dössel, O. (2009b). Localizing ectopic foci in the pulmonary veins from intracardiac ECGs - A simulation study. IFMBE 25, 645–648. doi: 10.1007/978-3-642-03882-2-172
Wettwer, E., Hála, O., Christ, T., Heubach, J. F., Dobrev, D., Knaut, M., et al. (2004). Role of IKur in controlling action potential shape and contractility in the human atrium: Influence of chronic atrial fibrillation. Circulation 110, 2299–2306. doi: 10.1161/01.CIR.0000145155.60288.71
Whittaker, D. G., Ni, H., Harchi, A. E., Hancox, J. C., and Zhang, H. (2017). Atrial arrhythmogenicity of KCNJ2 mutations in short QT syndrome: Insights from virtual human atria. PLoS Comput. Biol. 13:e1005593. doi: 10.1371/journal.pcbi.1005593
Wijffels, M. C., Kirchhof, C. J., Dorland, R., and Allessie, M. A. (1995). Atrial fibrillation begets atrial fibrillation?: a study in awake chronically instrumented goats. Circulation 92, 1954–1968. doi: 10.1161/01.CIR.92.7.1954
Wilhelms, M., Hettmann, H., Maleckar, M. M., Koivumaki, J. T., Dossel, O., Seemann, G., et al. (2013). Benchmarking electrophysiological models of human atrial myocytes. Front. Physiol. 3:487. doi: 10.3389/fphys.2012.00487
Wolf, R. M., Glynn, P., Hashemi, S., Zarei, K., Mitchell, C. C., Anderson, M. E., et al. (2013). Atrial fibrillation and sinus node dysfunction in human ankyrin-B syndrome: A computational analysis. Am. J. Physiol. Hear. Circ. Physiol. 304, H1253–H1266. doi: 10.1152/ajpheart.00734.2012
Xia, L., Gong, Y.-L., Zhu, X.-W., Zhang, Y., Sun, Q., and Zhang, H.-G. (2010). Mathematical models of canine right and left atria cardiomyocytes. J. Zhejiang Univ. Sci. B 11, 402–416. doi: 10.1631/jzus.B0900346
Yamazaki, M., and Jalife, J. (2012). Pathophysiology of atrial fibrillation: From initiation to maintenance. J. Arrhythmia 28, 129–139. doi: 10.1016/j.joa.2012.05.002
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
Zemzemi, N., Labarthe, S., Dubois, R., and Coudiere, Y. (2012). “From body surface potential to activation maps on the atria: A machine learning technique,” in Proceedings of the Computing in Cardiology, (Krakow: IEEE), 125–128.
Zhang, H., Garratt, C. J., Zhu, J., and Holden, A. V. (2005). Role of up-regulation of IK1 in action potential shortening associated with atrial fibrillation in humans. Cardiovasc. Res. 66, 493–502. doi: 10.1016/j.cardiores.2005.01.020
Zhang, N., Gong, M., Tse, G., Zhang, Z., Meng, L., Yan, B. P., et al. (2018). Prolonged corrected QT interval in predicting atrial fibrillation: a systematic review and meta-analysis. Pacing Clin. Electrophysiol. 41, 321–327. doi: 10.1111/pace.13292
Zhao, J., Amiri, A., Sands, G. B., Trew, M., LeGrice, I., Smaill, B. H., et al. (2008). Structure specific models of electrical function in the right atrial appendage. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2008, 145–148. doi: 10.1109/IEMBS.2008.4649111
Zhao, J., Butters, T. D., Zhang, H., Legrice, I. J., Sands, G. B., and Smaill, B. H. (2013a). Image-based model of atrial anatomy and electrical activation: A computational platform for investigating atrial arrhythmia. IEEE Trans. Med. Imaging 32, 18–27. doi: 10.1109/TMI.2012.2227776
Zhao, J., Stephenson, R. S., Sands, G. B., LeGrice, I. J., Zhang, H., Jarvis, J. C., et al. (2013b). Atrial fibrosis and atrial fibrillation: a computer simulation in the posterior left atrium. Lect. Notes Comput. Sci. 7945, 400–408. doi: 10.1007/978-3-642-38899-6_47
Zhao, J., Butters, T. D., Zhang, H., Pullan, A. J., LeGrice, I. J., Sands, G. B., et al. (2012). An image-based model of atrial muscular architecture effects of structural anisotropy on electrical activation. Circ. Arrhythmia Electrophysiol. 5, 361–370. doi: 10.1161/CIRCEP.111.967950
Zhao, J., Hansen, B. J., Wang, Y., Csepe, T. A., Sul, L. V., Tang, A., et al. (2017). Three-dimensional integrated functional, structural, and computational mapping to define the structural “fingerprints” of heart-specific atrial fibrillation drivers in human heart ex vivo. J. Am. Heart Assoc. 6, e005922. doi: 10.1161/JAHA.117.005922
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
Zhu, X., Yang, D., Lu, W., Chen, W., Wei, D., Fukuda, K., et al. (2014). “Computer simulation of cathode ablation for atrial fibrillation,” in Proceedings - 2014 IEEE International Conference on Computer and Information Technology, CIT 2014, (New Jersey: Institute of Electrical and Electronics Engineers Inc.), 721–725. doi: 10.1109/CIT.2014.37
Ziyadeh-Isleem, A., Clatot, J., Duchatelet, S., Gandjbakhch, E., Denjoy, I., Hidden-Lucet, F., et al. (2014). A truncating SCN5A mutation combined with genetic variability causes sick sinus syndrome and early atrial fibrillation. Hear. Rhythm 11, 1015–1023. doi: 10.1016/j.hrthm.2014.02.021
Keywords: atrial fibrillation, multiple P-wave morphologies, sinus rhythm, computational models, simulation, scoping review
Citation: Filos D, Tachmatzidis D, Maglaveras N, Vassilikos V and Chouvarda I (2019) Understanding the Beat-to-Beat Variations of P-Waves Morphologies in AF Patients During Sinus Rhythm: A Scoping Review of the Atrial Simulation Studies. Front. Physiol. 10:742. doi: 10.3389/fphys.2019.00742
Received: 15 May 2018; Accepted: 28 May 2019;
Published: 18 June 2019.
Edited by:
Edward Joseph Vigmond, Université de Bordeaux, FranceReviewed by:
Mark Potse, Inria Bordeaux – Sud-Ouest Research Centre, FranceMartin Bishop, King’s College London, United Kingdom
Copyright © 2019 Filos, Tachmatzidis, Maglaveras, Vassilikos and Chouvarda. 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: Dimitrios Filos, ZGltZmlsb3NAYXV0aC5ncg==