Commentary: Atrial Rotor Dynamics Under Complex Fractional Order Diffusion

Even at healthy states, cardiac tissue conforms one of the most representative cases of a highly heterogeneous and composite biological medium, whose spatial complexity has for long been known to modulate electrical conduction (Frank and Langer, 1974; Spach et al., 1981). How to capture its intricate structural heterogeneity at a tractable cost remains an open challenge in computational physiology and medicine, as traditional approaches such as the monodomain or bidomain equations inherently assume that the tissue behaves as an averaged syncytium with negligible contribution of its composite microstructure.

To overcome some of these limitations, we recently pioneered the use of fractional diffusion for the description of cardiac conduction (Bueno-Orovio et al., 2014b). Our proposed framework took the form

where −(−∇ · D∇V)^α/2 is the so-called fractional Laplacian of real order 1 < α ≤ 2. For α = 2, the model clearly recovers the standard monodomain equation, and could equally be extended to the bidomain setting. The well-founded potential theory around the fractional Laplacian allowed us to establish its biophysical interpretation, showing it represents the modulation of the electrical field of a homogeneous conductor by the secondary electrical sources associated with its inhomogeneities. The model further helped elucidating formerly unrelated effects of tissue microstructure on cardiac conduction, including widespread of the action potential foot during depolarization, action potential shortening along the activation pathway, and the modulated dispersion of repolarization. Experimentally, the model has been supported by diffusion spectrum imaging in ex-vivo hearts, indicating fractional diffusion metrics as indices of myocardial microstructure (Bueno-Orovio et al., 2016), as well as by high-resolution optical mapping on cardiac tissue preparations, demonstrating fractional scaling in the propagation of the cardiac wavelength (Loppini et al., 2018).

In the work under comment, Ugarte et al. build and expand on these ideas to present a two-dimensional isotropic fractional diffusion framework of complex order, of the form

with operators and involving pairs of complex-conjugate fractional derivatives defined by

where γ = α + jβ is the complex fractional order, and its complex conjugate. The authors' interpretation of their complex-order model newly builds on potential theory, which connects potential distributions over fractal domains and the complex-order fractional Laplacian. Indeed, the inclusion of the imaginary part β implies that cardiac tissue must satisfy a discrete-scale fractal structure (self-similarity at discrete scales). Whilst such a self-similarity could perhaps be arguable for the main volumetric constituents of cardiac tissue (cardiomyocytes), other components might very well exhibit a fractal structure (e.g., microvasculature). Importantly, such a complex-order fractional framework holds a great potential for consideration of structural remodeling, shall the associated structures (e.g. fibrotic clefts) are proved to have a fractal nature.

However, an important limitation of Ugarte et al. (2018) is that their proposed model is not consistent with the fractional Laplacian in which the authors base their analysis. Taking β = 0 for simplicity, Equations (2), (3) then reduce to

known as a fractional Riesz operator (fractional derivatives independently applied in each spatial coordinate). Conversely, under two-dimensional isotropic conditions, the fractional Laplacian model given by Equation (1) becomes

where for clarity the same notation κ_α has been used for the equivalent diffusion coefficient. Comparing (4) and (5), it becomes evident that the proposed fractional model is only equivalent to the fractional Laplacian under the standard diffusion case, given by α = 2.

The implications of these subtle but important discrepancies on cardiac conduction are exemplified in Figure 1. Simulations illustrate isotropic conduction for both models under decreasing fractional order α, with ion dynamics described for simplicity by Fenton and Karma (1998). Whereas the fractional Laplacian (Figure 1A) correctly replicates for all α the circular propagation patterns observed on isotropic cardiac monolayers as the simplest yet inhomogeneous in-vitro model of cardiac tissue (Badie and Bursac, 2009; Bian et al., 2014; Molitoris et al., 2016), the fractional Riesz operator (Figure 1B) induces increasingly larger curvature artifacts on wavefront conduction for decreasing α. Such curvature artifacts indeed translate into the results of Ugarte et al. (2018), as evidenced by their square-like spiral wavefronts and rotor trajectories. Given the well-known curvature-related modulation of conduction velocity and therefore wavefront-waveback interactions (Fast and Kléber, 1997; Comtois and Vinet, 1999; Comtois et al., 2005; Kadota et al., 2012), their results on vulnerability to re-entry and associated rotor biomarkers thus must be cautiously interpreted.

Figure 1

It is nevertheless relevant to note that more squared propagation patterns have been reported in both optical mapping (Koura et al., 2002; de Diego et al., 2011) and computational (He and Liu, 2010) studies. This was however under marked anisotropic conduction, not accounted in the isotropic model by Ugarte et al. (2018). In addition, fractional Riesz operators have been also used in modeling electrical propagation (Liu et al., 2013, 2015; Zeng et al., 2014). Such works, more centered in numerical analysis than in gaining physiological insights, might be additionally contributing to spreading the inconsistencies between these two types of fractional diffusion operators. Finally, a too coarse spatial resolution for atrial dynamics compared to previous studies (Wilhelms et al., 2013) could also contribute to partially unresolved re-entrant patterns. Although minimized by the high-order approach on which the authors base their numerical methods (Bueno-Orovio et al., 2014a), allowing considerably larger space steps than traditional stencils, this aspect certainly deserves further consideration.

As previously discussed, the ideas presented in Ugarte et al. (2018) hold a great potential for advancing the field of fractional diffusion applied to cardiac tissue, in order to promote our understanding of the role of tissue microstructure and structural remodeling in modulating wavefront propagation. However, this contribution raises awareness on the definition of suitable fractional diffusion models, exemplifying that simply recovering standard diffusion for a specific value of the considered tissue parameters is not a sufficient condition for realistic cardiac conduction. In this regard, frameworks that are consistent with the fractional Laplacian (Bueno-Orovio et al., 2014b; Cusimano et al., 2015; Cusimano and Gerardo-Giorda, 2018) seem a more suitable modeling approach to correctly capture the characteristic electrotonic loading of cardiac tissue.

Statements

Summary

Keywords

cardiac tissue, structural heterogeneity, electrical propagation, fractional diffusion, wavefront curvature

Citation

Bueno-Orovio A (2018) Commentary: Atrial Rotor Dynamics Under Complex Fractional Order Diffusion. Front. Physiol. 9:1386. doi: 10.3389/fphys.2018.01386

Received

30 August 2018

Accepted

11 September 2018

Published

04 October 2018

Volume

9 - 2018

Edited by

Jichao Zhao, University of Auckland, New Zealand

Reviewed by

Juan Pablo Ugarte, Universidad de san Buenaventura, Medellin, Colombia; Arun V. Holden, University of Leeds, United Kingdom

Updates

Check for updates

Copyright

© 2018 Bueno-Orovio.

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: Alfonso Bueno-Orovio alfonso.bueno@cs.ox.ac.uk

This article was submitted to Computational Physiology and Medicine, a section of the journal Frontiers in Physiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.