arxiv: v1 [q-bio.to] 29 Oct 2018

Transcription

1 Unified Mechanism of Atrial Fibrillation in a Simple Model Max Falkenberg 1,2,3, Andrew J. Ford 1, Anthony C. Li 1, Alberto Ciacci 1,2,3, Nicholas S. Peters 3, and Kim Christensen 1,2,3 1 Blackett Laboratory, Imperial College London, London SW7 2AZ, United Kingdom 2 Centre for Complexity Science, Imperial College London, London SW7 2AZ, United Kingdom and 3 ElectroCardioMaths Programme, Imperial Centre for Cardiac Engineering, Imperial College London, London W12 0NN, United Kingdom (Dated: October 30, 2018) arxiv: v1 [q-bio.to] 29 Oct 2018 The mechanism of atrial fibrillation (AF) is poorly understood, resulting in disappointing success rates of ablative treatment. Different mechanisms defined largely by different atrial activation patterns have been proposed and, arguably, this dispute has slowed the progress of AF research. Recent clinical evidence suggests a unifying mechanism based on sustained re-entrant circuits in the complex atrial architecture. Here, we present a simple computational model showing spontaneous emergence of AF that strongly supports, and gives a theoretical explanation for, the clinically observed diversity of activation. We show that the difference in surface activation patterns is a direct consequence of the thickness of the discrete network of heart muscle cells through which electrical signals percolate to reach the imaged surface. The model naturally follows the clinical spectrum of AF spanning sinus rhythm, paroxysmal and persistent AF as the decoupling of myocardial cells results in the lattice approaching the percolation threshold. This allows the model to make additional predictions beyond the current clinical understanding, showing that for paroxysmal AF re-entrant circuits emerge near the endocardium, but in persistent AF they emerge deeper in the bulk of the atrial wall where endocardial ablation is less effective. If clinically confirmed, this may explain the lower success rate of ablation in long-lasting persistent AF. Atrial fibrillation (AF) is the most common cardiac arrhythmia and is projected to affect 1% of the world s population by 2030 [1]. The key disagreements in the mechanistic understanding of AF concern whether the sustaining mechanism of AF is due to localised drivers of fibrillatory activity or the random motion and interaction of multiple wavelets across the atria [2 6]. As a result, treatment by ablation still has disappointing success rates [7 9]. There is growing evidence in favour of local drivers as a sustaining mechanism of AF [4, 9 15], with a recent meta-analysis [16] demonstrating that ablating local drivers has an additional benefit over pulmonary vein isolation, the most common ablation strategy [17]. Ad- Corresponding author: k.christensen@imperial.ac.uk ditionally, it has been shown that persistent AF is only terminated when local drivers are targeted [18]. However, there is debate about the mechanistic origin of local drivers: Some studies have primarily identified the local drivers as rotors using phase analysis [11, 12], whereas another study has identified the drivers of AF as focal or breakthrough points [10]. Recently, a unifying mechanism has been proposed, suggesting that all local drivers may be explained by the presence of small re-entrant circuits [13, 19, 20]. These re-entrant circuits have a spatial extent of approximately 15mm [13]. Using high-resolution simultaneous optical mapping of the endocardium (inner layer of heart wall) and epicardium (outer layer of heart wall) in explanted, diseased human hearts, the study shows that stable transmural re-entrant circuits may project differently onto the endocardium and the epicardium. Projections may appear as rotational activity, or as stable or unstable breakthrough points. Hence, the apparently incompatible 2d projections onto the endocardium and epicardium are consistent with a single underlying mechanism of localised 3d re-entrant circuits in the transmural region [13 15]. In this paper, we demonstrate that these seemingly incompatible surface activation patterns are two sides of the same coin. We introduce a model in which activation wavefronts propagate in a 3d medium mimicking the complex discrete fibre structure of the atria. AF emerges spontaneously in the model through the formation of spatially stable but temporarily intermittent re-entrant circuits. The model predicts that these re-entrant circuits should have a spatial extent just above a lower bound of 12.5mm which is directly compatible with Ref. [13]. A strength of this model is that despite being very simple, it successfully demonstrates the range of activation patterns observed on the endocardium and the epicardium during AF, closely matching the key clinical observations in Ref. [13] as well as several other observations in the recent AF literature [10, 18, 21]. The observed activation patterns fundamentally change as a function of the depth from the imaged surface to the driving re-entrant circuit. If this depth is large, we observe focal or breakthrough activity. If this depth is small, we observe re-entrant or rotational activity. Hence, our model in combination with the clinical results in Ref. [13] provides strong evidence for this elegant unifying mechanism of AF. Furthermore, the model predicts that the drivers of paroxysmal AF are predominantly located near the endocardium whereas the

2 2 (a) (b) (c) (d) y Sub-endo x Sub-epi FIG. 1. The emergence of AF in the inhomogeneous model with simultaneous endocardial (z = 0) and epicardial (z = 24) imaging. Gaussian smoothing has been applied for clarity. Red: excited nodes, blue: resting (excitable) nodes, other: refractory (unexcitable) nodes. (a) Planar wavefronts propagating during sinus rhythm. Endo view: Dotted box indicates gap in wavefront formed by conduction block. (b) Endo view: Arrow indicates excitation re-entering gap in wavefront. Epi view: Re-entry is not observed. (c) Emergence of fibrillatory activity. (d) Maintenance of fibrillatory activity. Epi view: Activity emerges on surface as point source located at the star. See videos A in supplementary information. drivers of persistent AF are uniformly distributed within the bulk of the atrial wall. Consistent with clinical experience, the former is more easily ablated than the latter. The simplicity of our model, as opposed to highly complex models of AF [22], makes it ideal for testing mechanistic hypotheses, proposed treatment and preventive strategies, and targeting and delivering of ablation. Therefore, in the future, the model may be used as a test bed for unlocking new promising intervention procedures. Model. Electrical signals in the heart are mediated by discrete heart cells arranged in long, intertwined fibres. This motivates the construction of a 3d cellular automata model of AF where each node acts as a single cell (or block of cells), and nodes are connected stochastically to their neighbours to mimic the branching structure of discrete heart muscle fibres, generalising a 2d model [23]. Discrete models of cardiac tissue have gained popularity recently [23 27]. In these models, the onset of re-entry has been associated with the approach from above of the bond occupation probability to the percolation threshold. Some models have been extended to 3d to highlight the importance of the thickness of the atrial wall to the probability of re-entry [26]. Additionally, leading clinicians have highlighted the importance of bond percolation to the emergence of AF, noting that the appearance of complex fractionated atrial electrograms can be explained by percolation [25]. In our model, we consider a simplified L x L y L z pipe topology of the atria with open boundaries in x and z and periodic boundaries in y. Nodes are connected longitudinally to their neighbours in the x direction with frequency ν and transversely in the y and z directions with frequency ν. Each node can take one of three states: resting, where the node can be excited by an active neighbour, excited, or refractory, where for τ time steps after excitation a node cannot be re-excited. We define the sinus node (pacemaker) as nodes at the boundary, x = 0, which excites every T time steps. Model parameters L x = L y = 200, L z = 25, T = 220 and τ = 50 are informed by clinically observed values, see supplementary information. We define a small fraction of nodes, δ, as representing locations where cells are susceptible to conduction block. With small probability ɛ, these nodes fail to activate when their neighbours excite. When the sinus node excites in normal conduction, wavefronts are initiated at x = 0 and propagate smoothly in the +x direction and terminate at x = L x. However, re-entrant circuits can form when nodes are sufficiently decoupled, through fibrosis [28] or otherwise, such that the shortest closed loop from a node back to itself is partially isolated from the remaining tissue, and the path length, in units of propagation time, exceeds the refractory period τ. This gives a minimum spatial extent of τ/2 which, when converted back into clinical units, corresponds to 12.5mm. The formation of these re-entrant circuits, wrapped around a line of conduction block, initiate and maintain (drive) fibrillatory activity where the electrical activity is disorganised within the atria [23, 29], see Fig. 1. A key observation in Ref. [13] is that re-entry activity is typically observed on the endocardial surface, whereas focal or breakthrough points are often found on the epicardial surface. Additionally, there is much stronger longitudinal coupling on the endocardial surface than on the epicardial surface. To test whether fibre orientation can account for the distribution of activation patterns observed clinically, we consider a homogeneous and inhomogeneous model. In the homogeneous model all nodes are connected to their neighbours with the same combination of frequencies ν and ν. In the inhomogeneous model, we allow the variation in fibre direction to change with depth. Here we fix the average decoupling of nodes,

3 3 (a) (b) 1.0 Time in AF Depth z Frequency Depth z ν Frequency 0.2 ν = ν = FIG. 2. The average time spent in AF for a model with (a) homogeneous coupling parameters, and (b) inhomogeneous coupling parameters with a fixed average coupling, ν. Both models show a transition from sinus rhythm at large coupling (0% AF time) through paroxysmal AF (0% < AF time < 100%) to persistent AF (100% AF time). The grey parameter regions, indicate coupling parameters where there is no spanning network of muscle cells from the sinus node at x = 0 to x = L x, and are thus irrelevant to clinical AF. (a) Homogeneous coupling parameters. Blue: sinus rhythm, red: persistent AF. The white transition region corresponds to paroxysmal AF. Guidelines of positive (negative) gradient indicate constant θ ( ν). (b) Inhomogeneous coupling parameter risk curve (red graph). There is no AF for ν 0.4. Decreasing ν associated with decoupling cells, there is a transition from sinus rhythm through paroxysmal AF to persistent AF. Inset: Histograms showing the distribution of re-entrant circuits driving AF for paroxysmal AF (top) and persistent AF (bottom) as a function of depth z from the endocardium. Drivers cluster in the sub-endocardial region (z = 0) for paroxysmal AF but are uniformly distributed across the bulk for persistent AF. ν = (2ν +4ν )/6 and vary linearly the average fibre angle, θ = tan 1 (ν /ν ), in each layer from θ ENDO = 24 at the endocardium (z = 0) to θ EPI = 42 at the epicardium (z = 24), informed by [15]. Results. We have studied the emergence of AF in the model as a function of ν and ν for the homogeneous model and as a function of ν for the inhomogeneous model. The phase space for the risk of entering AF for homogeneous model is shown in Fig. 2a and for the inhomogeneous model in Fig. 2b. For large values of the coupling parameters (no fibrosis) the model exhibits sinus rhythm indefinitely. As the coupling reduces, due to fibrosis or otherwise, a transition takes place where a small number of re-entrant circuits can form. Here we observe paroxysmal AF with intermittent episodes of irregular activity. As the coupling is even further reduced, the model enters persistent AF where, once AF has been initiated, the model will never return to sinus rhythm without external intervention. These results are consistent with recent evidence showing that local drivers anchor at or near fibrotic lesions [12], with fibrosis increasing the number of re-entrant regions in the atria, and increasing the duration spent in AF [30]. In the model, we can quantify where re-entrant circuits emerge in the substrate and observe the associated activation patterns on the endocardial and epicardial surfaces. For the homogeneous model variant there is no distinction between the endocardium and the epicardium, and we find that in paroxysmal AF, the majority of drivers form equivalently on either surface and not in the bulk, whereas for persistent AF, the majority of drivers form in the bulk and not on the surfaces, see Fig. 11 in supplementary information. However, the inhomogeneous model breaks the symmetry between the epicardium and the endocardium, consistent with clinical observations [15]. Here we find that in paroxysmal AF, drivers form preferentially near the endocardial surface, with very few drivers forming on the epicardial surface and almost none in the bulk of the atria, see Fig. 2b, inset 1. However, in persistent AF, drivers are uniformly distributed throughout the atrial wall, see Fig. 2b, inset 2. Hence, as AF becomes more persistent, the average position of drivers moves away from the endocardium and into the bulk of the atrial wall. The surface activity driven from the atrial bulk will typically appear as a breakthrough point, as opposed to re-entry activity observed when drivers are close to the imaged surface. Indeed, in our model we find that ablating the active reentrant circuits in a tissue terminates AF, consistent with [13], see videos F&G in supplementary information. However, if the ablation does not penetrate far enough into the tissue, fibrillatory activity will continue, see Fig. 14 in supplementary information. In the 3d model, re-entrant circuits can form anywhere in the bulk of the atrial wall. If we visualise the activation patterns in the inhomogeneous model on the surfaces at z = 0 (endocardium) and z = 24 (epicardium), we observe a diversity of activation patterns consistent with

4 4 (a) (b) (c) Sub-endo Removed from arxiv version Complete due to copyright. See Fig. 7A in [Hansen et al., re-entry European Heart Journal, Volume 36, Issue 35, 14 September Incomplete 2015, Pages ]. re-entry Link in caption. Unstable breakthrough Sub-Endo Sub-endo 1 4 Sub-endo 2 Stable breakthrough 3 Sub-epi Sub-epi Sub-Epi 1 Sub-epi Removed from arxiv version due to copyright. See Fig. 5D in [Hansen et al., European Heart Journal, Volume 36, Issue 35, 14 September 2015, Pages ]. Link in caption FIG. 3. (a) The distribution of activation patterns observed from simultaneous mapping of the endocardium and the epicardium. LRA# refer to the patient number in the study [13]. The figure shows that re-entry activity is predominantly observed on the endocardial surface, whereas focal sources are only seen on the epicardium. (b) An activation map from the clinical study [13] showing re-entry activity from sub-endo imaging and a focal source on the sub-epi side with four marked reference points. (c) The equivalent as (b) for our model. Gaussian smoothing has been applied to images for clarity. Red: excited nodes, blue: resting (excitable) nodes, other: refractory nodes. (a) and (b) have been removed from arxiv version. See Figs. 7A and 5D respectively in [13]. For links to figures click HERE and HERE. the clinical study [13] which is summarised in Fig. 3a. Figure 3b shows a specific example activation map from the clinical study. Figure 3c shows the equivalent from our model. In Fig. 3, the left (right) panel shows the activity when viewed from the endocardium (epicardium). For the activation maps in (b) and (c), electrical activity spreads across the tissue from point (1) through point (2), but is blocked from reaching the isolated fibre at point (4). The activity loops around the isolated fibre and re-enters the fibre from the far end at point (3). This back-propagating excitation now passes through point (4) before re-exciting point (1) which sustains the re-entrant circuit. Viewing the same region from the epicardium does not show the same re-entry activity. Instead, the excitation emerges as a point source indicated by a star. In incomplete re-entry, only part of the cycle shown in the left panel is visible from endocardial imaging. For unstable breakthrough, the point at which excitations emerge in the epicardial view vary from beat to beat. For example activation maps and videos, see supplementary information. Discussion. These model observations are directly compatible with recent clinical findings that identified rotational activity from endocardial mapping in patients with paroxysmal AF [11, 21], and another clinical observation identifying focal activity in patients with persistent AF from epicardial mapping [10]. Moreover, our model offers a natural explanation why ablation is more successful for paroxysmal than persistent AF. Clinical studies have shown that the radius of ablation lesions become smaller as they penetrate further into the tissue from the endocardial surface, and that ablation struggles to penetrate more than 2mm into the wall of the atria [31] which can be up to 7mm thick [13]. Hence, because the average driver position moves deeper into the atrial wall as AF becomes more persistent, ablation lesions need to be more accurately positioned and they must penetrate further into the tissue. This restricts the efficacy of AF ablation. The model may also explain the prevalence of paroxysmal AF in the pulmonary veins (PVs) [17, 20]. Sleeves of cardiac tissue extend into the PVs from the atria. The sleeves get thinner, further into the PVs and show many abrupt changes in fibre orientation [32 34]. These physiological observations all reduce local cell to cell coupling which the model has shown is the key requirement for forming the re-entrant circuits that initiate and sustain AF. Conclusions. The overarching aim was to create the simplest possible bio-inspired model that would recreate the variety of clinically observed activation patterns. The model that integrates a coarse grained dynamics of an electrical signal percolating through a complex discrete fibre network mimicking the real 3d structure of the human atria is very simple in its formulation. Nevertheless, we show that a variety of known key characteristics of clinical AF emerge naturally without much of the details thought to be essential in modern AF modelling. In addition to explaining the diversity of activation patterns, the model is consistent with the typical evolution of AF from sinus rhythm through paroxysmal AF to persistent AF. Combined with the clinical observations in Ref. [13], our model gives substantial evidence for the proposed unifying mechanism of AF, suggesting a potential resolution to long-standing debates. Complementing this finding, we predict that the inhomogeneity in the fibre orientation causes the average depth of drivers to move from the endocardium and into the bulk of the atrial wall as AF becomes more persistent. We hypothesise that this may explain why ablation is less successful for persistent AF

5 5 than for paroxysmal AF. This insight is the first of its kind and may have key clinical implications. In the future, the model will allow us to test the hypothesis that adjusting a small number of parameters underpins the mechanistic origin of AF. Additionally, coupled with newly available high resolution atrial fibre maps [35], the model will enable rigorous testing and development of treatment and prevention strategies. Author Contributions. KC conceived the model. MFM proposed the connection to [13] and adapted the model for variable fibre orientation. AJF & ACL implemented the model and ran simulations, supervised by MFM. AC developed the theoretical aspects of the model. NSP provided the medical expertise for the project. KC supervised all aspects of project. All authors contributed to the writing of the manuscript. Acknowledgements. MFM thanks Timotej Kapus and Éamonn Murray for computational support. MFM and AC gratefully acknowledge support from the EP- SRC. NSP acknowledges funding from the British Heart Foundation (RG/16/3/32175 and Centre of Research Excellence), and the National Institute for Health Research (UK) Biomedical Research Centre. KC and NSP acknowledge funding from the Rosetrees Trust. [1] Patel, N. J., Atti, V., Mitrani, R. D., Viles-Gonzalez, J. F. & Goldberger, J. J. Global rising trends of atrial fibrillation: a major public health concern. Heart (2018). [2] Nattel, S. New ideas about atrial fibrillation 50 years on. Nature 415, (2002). [3] Schotten, U., Dobrev, D., Platonov, P. G., Kottkamp, H. & Hindricks, G. Current controversies in determining the main mechanisms of atrial fibrillation. J. Intern. Med. 279, (2016). [4] Fedorov, V. V. & Hansen, B. J. Now you see a rotor, now you dont. JACC: Clin. Electrophysiol. 4, (2018). [5] Nattel, S. & Dobrev, D. Controversies about atrial fibrillation mechanisms: aiming for order in chaos and whether it matters. Circ. Res. 120, (2017). [6] Mann, I., Sandler, B., Linton, N. & Kanagaratnam, P. Drivers of atrial fibrillation: Theoretical considerations and practical concerns. Arrhythm. Electrophysiol. Rev. 7, (2018). [7] Verma, A. et al. Approaches to catheter ablation for persistent atrial fibrillation. New Engl. J. Med. 372, (2015). [8] Macle, L. et al. Adenosine-guided pulmonary vein isolation for the treatment of paroxysmal atrial fibrillation: an international, multicentre, randomised superiority trial. Lancet 386, (2015). [9] Narayan, S. M. et al. The continuous challenge of af ablation: From foci to rotational activity. Rev. Port. Cardiol. 36, 9 17 (2017). [10] Lee, S. et al. Simultaneous Bi-Atrial High Density ( electrodes) Epicardial Mapping of Persistent and Long-Standing Persistent Atrial Fibrillation in Patients: New Insights into the Mechanism of Its Maintenance. Circulation 132, (2015). [11] Narayan, S. M., Krummen, D. E., Clopton, P., Shivkumar, K. & Miller, J. M. Direct or coincidental elimination of stable rotors or focal sources may explain successful atrial fibrillation ablation. J. Am. Coll. Cardiol. 62, (2013). [12] Haïssaguerre, M. et al. Driver domains in persistent atrial fibrillation. Circulation 130, (2014). [13] Hansen, B. J. et al. Atrial fibrillation driven by microanatomic intramural re-entry revealed by simultaneous sub-epicardial and sub-endocardial optical mapping in explanted human hearts. Eur. Heart J. 36, (2015). [14] Hansen, B. J. et al. Maintenance of Atrial Fibrillation: Are Reentrant Drivers with Spatial Stability the Key? Circ. Arrhythm. Electrophysiol. 9, e (2016). [15] Csepe, T. A., Hansen, B. J. & Fedorov, V. V. Atrial fibrillation driver mechanisms: Insight from the isolated human heart. Trends Cardiovas. Med. 27, 1 11 (2017). [16] Baykaner, T. et al. Clinical implications of ablation of drivers for atrial fibrillation: A systematic review and meta-analysis. Circ. Arrhythm. Electrophysiol. 11, e (2018). [17] Haïssaguerre, M. et al. Spontaneous initiation of atrial fibrillation by ectopic beats originating in the pulmonary veins. New Engl. J. Med. 339, (1998). [18] Zaman, J. A. B. et al. Identification and characterization of sites where persistent atrial fibrillation is terminated by localized ablation. Circ. Arrhythm. Electrophysiol. 11, e (2018). [19] Nattel, S., Xiong, F. & Aguilar, M. Demystifying rotors and their place in clinical translation of atrial fibrillation mechanisms. Nat. Rev. Cardiol. 14, (2017). [20] Calkins, H. et al HRS/EHRA/ECAS/APHRS/ SOLAECE expert consensus statement on catheter and surgical ablation of atrial fibrillation. Heart Rhythm 14, e275 e444 (2017). [21] Swarup, V. et al. Stability of rotors and focal sources for human atrial fibrillation: focal impulse and rotor mapping (firm) of af sources and fibrillatory conduction. J. Cardiovasc. Electr. 25, (2014). [22] Clayton, R. H. et al. Models of cardiac tissue electrophysiology: progress, challenges and open questions. Prog. Biophys. Mol. Bio. 104, (2011). [23] Christensen, K., Manani, K. A. & Peters, N. S. Simple model for identifying critical regions in atrial fibrillation. Phys. Rev. Lett. 114, (2015). [24] Alonso, S. & Bär, M. Reentry near the percolation threshold in a heterogeneous discrete model for cardiac tissue. Phys. Rev. Lett. 110, (2013). [25] Vigmond, E., Pashaei, A., Amraoui, S., Cochet, H. & Ha M. Percolation as a mechanism to explain atrial fractionated electrograms and reentry in a fibrosis model based on imaging data. Heart Rhythm 13, (2016). [26] Alonso, S., dos Santos, R. W. & Bär, M. Reentry and ectopic pacemakers emerge in a three-dimensional model for a slab of cardiac tissue with diffuse microfibrosis near the percolation threshold. PloS One 11, e (2016). [27] Lin, Y. T., Chang, E. T. Y., Eatock, J., Galla, T. & Clay-

6 6 ton, R. H. Mechanisms of stochastic onset and termination of atrial fibrillation studied with a cellular automaton model. J. Royal Soc. Interface 14, (2017). [28] Spach, M. S. & Boineau, J. P. Microfibrosis produces electrical load variations due to loss of side-to-side cell connections; a major mechanism of structural heart disease arrhythmias. Pacing Clin. Electrophysio. 20, (2006). [29] McGillivray, M. F., Cheng, W., Peters, N. S. & Christensen, K. Machine learning methods for locating reentrant drivers from electrograms in a model of atrial fibrillation. Roy. Soc. Open Sci. 5, (2018). [30] Cochet, H. et al. Relationship Between Fibrosis Detected on Late Gadolinium-Enhanced Cardiac Magnetic Resonance and Re-Entrant Activity Assessed With Electrocardiographic Imaging in Human Persistent Atrial Fibrillation. JACC: Clin. Electrophysiol. 4, (2018). [31] Kumar, S. et al. Better lesion creation and assessment during catheter ablation. J. Atr. Fibrillation 8, 1189 (2015). [32] Yen Ho, S., Sanchez-Quintana, D., Cabrera, J. A. & Anderson, R. H. Anatomy of the left atrium. J. Cardiovasc. Electrophysiol. 10, (1999). [33] Weiss, C. et al. Impact of the distribution and structure of myocardium in the pulmonary veins for radiofrequency ablation of atrial fibrillation. Pacing Clin. Electrophysio. 25, (2002). [34] Weerasooriya, R. et al. Dissociated pulmonary vein arrhythmia: Incidence and characteristics. J. Cardiovasc. Electrophysiol. 14, (2003). [35] Pashakhanloo, F. et al. Myofiber architecture of the human atria as revealed by submillimeter diffusion tensor imaging. Circ. Arrhythm. Electrophysiol. 9, e (2016).