General Commentary ARTICLE
Commentary: Atrial Rotor Dynamics Under Complex Fractional Order Diffusion
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
A Commentary on
Atrial Rotor Dynamics Under Complex Fractional Order Diffusion
by Ugarte, J. P., Tobón, C., Lopes, A. M., and Tenreiro Machado, J. A. (2018). Front. Physiol. 9:975. doi: 10.3389/fphys.2018.00975
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. Impact of fractional diffusion operators on cardiac conduction for decreasing order α. Activation maps for central domain stimulation are shown (5 ms separation isochrones). (A) Fractional Laplacian (Bueno-Orovio et al., 2014b). (B) Fractional Riesz operator (Ugarte et al., 2018). Diffusion coefficients were optimized to match standard diffusion conduction velocity at the center of 5 cm fiber strands. Ionic term: Fenton–Karma (modified Beeler–Reuter) dynamics; domain size: 5 × 5 cm; space discretisation: 512 × 512 points; time resolution: 0.025 ms.
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.
The author confirms being the sole contributor of this work and has approved it for publication.
AB-O holds a British Heart Foundation Intermediate Basic Science Research Fellowship (FS/17/22/32644).
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Bueno-Orovio, A., Kay, D., Grau, V., Rodriguez, B., and Burrage, K. (2014b). Fractional diffusion models of cardiac electrical propagation: role of structural heterogeneity in dispersion of repolarization. J. R. Soc. Interface 11:20140352. doi: 10.1098/rsif.2014.0352
Bueno-Orovio, A., Teh, I., Schneider, J. E., Burrage, K., and Grau, V. (2016). Anomalous diffusion in cardiac tissue as an index of myocardial microstructure. IEEE Trans. Med. Imaging 35, 2200–2207. doi: 10.1109/TMI.2016.2548503
Comtois, P., Kneller, J., and Nattel, S. (2005). Of circles and spirals: bridging the gap between the leading circle and spiral wave concepts of cardiac reentry. Europace 7, S10–S20. doi: 10.1016/j.eupc.2005.05.011
Cusimano, N., Bueno-Orovio, A., Turner, I., and Burrage, K. (2015). On the order of the fractional Laplacian in determining the spatio-temporal evolution of a space-fractional model of cardiac electrophysiology. PLoS ONE 10:e0143938. doi: 10.1371/journal.pone.0143938
Cusimano, N., and Gerardo-Giorda, L. (2018). A space-fractional monodomain model for cardiac electrophysiology combining anisotropy and heterogeneity on realistic geometries. J. Comput. Phys. 362, 409–424. doi: 10.1016/j.jcp.2018.02.034
de Diego, C., Chen, F., Xie, Y., Pai, R. K., Slavin, L., Parker, J., et al. (2011). Anisotropic conduction block and reentry in neonatal rat ventricular myocyte monolayers. Am. J. Physiol. Heart Circ. Physiol. 300, H271–H278. doi: 10.1152/ajpheart.00758.2009
Koura, T., Hara, M., Takeuchi, S., Ota, K., Okada, Y., Miyoshi, S., et al. (2002). Anisotropic conduction properties in canine atria analyzed by high-resolution optical mapping: preferential direction of conduction block changes from longitudinal to transverse with increasing age. Circulation 105, 2092–2098. doi: 10.1161/01.CIR.0000015506.36371.0D
Liu, F., Zhuang, P., Turner, I., Anh, V., and Burrage, K. (2015). A semi-alternating direction method for a 2-D fractional FitzHugh–Nagumo monodomain model on an approximate irregular domain. J. Comput. Phys. 293, 252–263. doi: 10.1016/j.jcp.2014.06.001
Loppini, A., Gizzi, A., Cherubini, C., Cherry, E. M., Fenton, F. H., and Filippi, S. (2018). Spatiotemporal correlation uncovers fractional scaling in cardiac tissue. arXiv [preprint] arXiv:1806.04507.
Molitoris, J. M., Paliwal, S., Sekar, R. B., Blake, R., Park, J., Trayanova, N. A., et al. (2016). Precisely parameterized experimental and computational models of tissue organization. Integr. Biol. 8, 230–242. doi: 10.1039/C5IB00270B
Spach, M. S., Miller III, W. T., Geselowitz, D. B., Barr, R. C., Kootsey, J. M., and Johnson, E. A. (1981). The discontinuous nature of propagation in normal canine cardiac muscle. Evidence for recurrent discontinuities of intracellular resistance that affect the membrane currents. Circ. Res. 48, 39–54. doi: 10.1161/01.RES.48.1.39
Wilhelms, M., Hettmann, H., Maleckar, M. M., Koivumäki, J. T., Dössel, O., and Seemann, G. (2013). Benchmarking electrophysiological models of human atrial myocytes. Front. Physiol. 3:487. doi: 10.3389/fphys.2012.00487
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.
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
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, firstname.lastname@example.org