# From Millimeters to Micrometers; Re-introducing Myocytes in Models of Cardiac Electrophysiology

^{1}Simula Research Laboratory, Lysaker, Norway^{2}Department of Physiology and Pharmacology, Faculty of Medicine, University of Calgary, Calgary, AB, Canada

Computational modeling has contributed significantly to present understanding of cardiac electrophysiology including cardiac conduction, excitation-contraction coupling, and the effects and side-effects of drugs. However, the accuracy of *in silico* analysis of electrochemical wave dynamics in cardiac tissue is limited by the homogenization procedure (spatial averaging) intrinsic to standard continuum models of conduction. Averaged models cannot resolve the intricate dynamics in the vicinity of individual cardiomyocytes simply because the myocytes are not present in these models. Here we demonstrate how recently developed mathematical models based on representing every myocyte can significantly increase the accuracy, and thus the utility of modeling electrophysiological function and dysfunction in collections of coupled cardiomyocytes. The present gold standard of numerical simulation for cardiac electrophysiology is based on the bidomain model. In the bidomain model, the extracellular (E) space, the cell membrane (M) and the intracellular (I) space are all assumed to be present everywhere in the tissue. Consequently, it is impossible to study biophysical processes taking place close to individual myocytes. The bidomain model represents the tissue by averaging over several hundred myocytes and this inherently limits the accuracy of the model. In our alternative approach both E, M, and I are represented in the model which is therefore referred to as the EMI model. The EMI model approach allows for detailed analysis of the biophysical processes going on in functionally important spaces very close to individual myocytes, although at the cost of significantly increased CPU-requirements.

## 1. Introduction

In this brief Perspective we present our approach for cell-based modeling of the cardiac syncytium. We point out the limitations of the standard, averaged models, and then illustrate how our approach can be used to reveal fundamental dynamics of a range of important, but incompletely understood mechanisms of arrhythmia. The main advantage of the cell-based approach is significantly increased accuracy, and the main disadvantage is significantly increased cost of the computations, and increased cost of the associated software development.

## 2. A Quantitative Understanding of Most Arrhythmias Remains Elusive

Optimal functioning of the heart relies on distinct electrochemical waves repetitively traversing the entire four chambered structure. Even transient perturbations to this wave, referred to as arrhythmias, can be life-threating and have therefore been subject to intense research efforts for at least a century. In general, both atrial and ventricular arrhythmias are well-defined clinical phenomena, but for many, precise explanations of their origin and/or maintenance—how the coordinated activity of myocytes becomes pathologically discoordinated—is lacking. This significantly compromises clinical efforts to restore normal cardiac conduction (Spach et al., 2007; Nattel and Dobrev, 2017; Heijman et al., 2021).

## 3. The Limit of Accuracy of the Standard Models of Cardiac Electrophysiology

The adult human heart contains several billion myocytes (Tirziu et al., 2010). Spatial averaging (homogenization) has been an invaluable tool for creating computationally tractable mathematical models of electrochemical wave dynamics including the action potential (AP) and details of change in intracellular calcium (Ca^{2+}). This has led to broad adoption of two standard models of cardiac electrophysiology; the monodomain and bidomain models. Since the late 1970's, the bidomain model (see e.g., Tung, 1978; Neu and Krassowska, 1993; Franzone et al., 2014) has been the gold standard for simulating electrochemical conduction in cardiac tissue, and has often been approximated by the slightly simpler monodomain model (Sundnes et al., 2007; Clayton and Panfilov, 2008; Vigmond et al., 2008). Coupled with models of myocyte membrane ion transport, these models have been very successful in accounting key aspects of the electrical activity in cardiac tissue. Furthermore, owing to dramatically increased computing power and substantially improved numerical methods, it is now possible to solve the bidomain model to convergence. A spatial resolution of Δ*x* ≈ 0.25 mm is generally considered to be sufficient to compute the solution of the model (Xie et al., 2004; Clayton and Panfilov, 2008; Niederer S. A. et al., 2011; Niederer S. et al., 2011). This mesh resolution represents tissue blocks containing ~980 myocytes each (assuming myocyte volume of 16 pL, including associated extracellular space, see numerical example below; Nygren et al., 1998). Since the converged solution can be computed, further mesh refinement does not provide greater insight to the physiological processes or dynamic characteristics of normal conduction. However, electrochemical conduction disturbances that underlie arrhythmia are not well-approximated by normal planar wave conduction. Rotor wave dynamics provide one important and ubiquitous example in which key details of the wave dynamics in smaller tissue regions are likely to be important. In this paradigm, the trajectory of a rotating wave is determined by a phase singularity and an adjacent region of maximum wavefront curvature. The singularity “meanders” around a central core and the maximum curvature dictates that meander in a manner that determines measurable clinical properties of the arrhythmia (e.g., dominant frequency of fibrillation). Importantly, even in large mammal (sheep) myocardium, rotor wave cores measured to be as small as 3 mm^{2} (Mandapati et al., 2000) are sufficient to support arrhythmic activity. The rotor wave core is one case in which understanding sub-millimeter-scale dynamics is likely to be important for defining and managing macroscopic outcomes. In addition, other dynamical conditions (e.g., ectopic focus), and heterogeneous innervation or non-uniform drug effects also likely rely on dynamics that are not adequately captured by models constructed to replicate only planar wave dynamics.

As mentioned above, the bidomain model has been very successful in simulating propagation of electrochemical waves in cardiac tissue on the macroscale when the transmembrane conductances are assumed to vary only on that scale. However, key properties of gap junctions can vary individually. Therefore, averaging can give misleading results. In order to illustrate this, let us consider a case where a long, one-dimensional strand in which myocytes are connected via gap junctions. With normal gap junctions and a strong electrical stimulation at one end of the strand, an excitation wave will be initiated and move along all myocytes both for EMI and for the bidomain models. Also, if the parameters of the models are properly adjusted, the conduction velocity of this wave will be similar. Suppose next that one gap junction is significantly disturbed leading to almost infinite resistance from one myocyte to the next. At that point the wave will stop when using the EMI model since the impaired junction will cause a complete conduction block. For the bidomain model however, change of one gap junction connection will mean little when the average is computed and the wave will only move a little slower, but not stop. Although this is a specific single example and the bidomain model can probably be adjusted (e.g., by mesh alignment) to compute the correct solution, it does illustrate the fundamental difficulties of averaging.

## 4. Cell-Based Models of Electrophysiology

A second major assumption of the monodomain and bidomain models is that the extracellular space, the myocyte membrane and the intracellular space are all present everywhere in the cardiac tissue. This represents a very significant simplification since individual myocytes can be ignored, but it also imposes strict limits on the range of problems than can be realistically studied using these models. For example it excludes the ability to understand how spatially localized ion channel expression may impact macroscopic conduction, and this is known to be fundamental to the contribution of sodium channel activation to cardiac conduction (Rivaud et al., 2017; Jæger et al., 2019). In short, classical cardiac models of electrochemical conduction are relevant for phenomena at scales from several millimeters to centimeters, and under conditions approximating planar conduction; whereas dynamics occurring among groups of myocytes at or below the millimeter scale are likely to be key to many mechanisms of arrhythmia. Here, we further develop our approach described in Tveito et al. (2017a), Jæger et al. (2019), and Jæger and Tveito (2021) and show that it is feasible to devise mathematical models at the micrometer level and thereby include every individual myocyte in the *in silico* tissue. We refer to this as the EMI model since it explicitly represents the extracellular (E) space, the cell membrane (M), and the intracellular (I) space. Figure 1B illustrates the different components of the spatial EMI model domain. Our results demonstrate the utility of the EMI model by addressing key aspects of how arrhythmias can arise in cardiac tissue. Specifically, we show that re-entry can arise in very small, partially de-coupled, collections of myocytes.

**Figure 1**. Illustration of the EMI model domain. **(A)** Shows the approximate cylindrical geometry and associated mesh of a single myocyte of length 120 μm and diameter ranging from 13 to 14 μm (Nygren et al., 1998). **(B)** Shows an illustration of the different components of the EMI model domain for an example collection of four connected myocytes. The domain consists of a number of myocytes surrounded by an extracellular space. The cell membrane is defined at the interface between the intracellular and extracellular spaces and intercalated discs with gap junctions are defined at the interface between adjacent myocytes. All computations presented here are in 3D.

## 5. Cell-Based Models Can Resolve Micro-Reentry in a Simulated Pulmonary Vein Sleeve

In atrial fibrillation, the arrhythmia often begins in the “sleeve” of the pulmonary veins of the left atrium (see e.g., Haissaguerre et al., 1998). In particular, the pulmonary vein/left atrial “junction” is assumed to be a driver of atrial fibrillation partially due to the highly heterogeneous intercellular coupling in this region (Koura et al., 2002; Pfenniger, 2020), and regions of structural conduction discontinuities (Hamabe et al., 2003). Heterogeneities in small collections of myocytes represented by cell-to-cell variations in membrane ion currents and intercellular coupling can initiate arrhythmias. Arrhythmias caused by heterogeneities can clearly not be accurately represented by models where the heterogeneities are accounted for only by averages. In Figure 2, we illustrate a collection of 25 × 25 myocytes covering an area of 1 mm^{2} with associated volume of 18 nL. We simulate the electrochemical wave (AP and Ca^{2+}-transient) in this heterogeneous syncytium using the EMI model, combined with a myocyte membrane model representing canine left pulmonary vein myocytes. To represent phenotypic heterogeneity in the tissue, the maximum conductance of each ionic current are in fact given as a random combination of these values in models of the left atrial free wall myocyte and the pure pulmonary vein myocyte. These differences are derived from detailed studies of the two myocyte populations (Ehrlich et al., 2003; Melnyk et al., 2005). Furthermore, to replicate one known property of diseased atrial tissue the intercellular resistance is increased (reducing gap-junction conductance) again chosen randomly for each myocyte-to-myocyte coupling in the right panel of Figure 2, similarly to Cherry et al. (2007). The simulations reveal that stable re-entrant waves can develop and propagate in very small collections of unhealthy myocytes under conditions of markedly reduced intercellular coupling. Further details of the model and method of solution are provided in Supplementary Material.

**Figure 2**. Re-entry in a cell-based model. The upper panels **(A,B)** show the results (intracellular potential) of S1+S2 stimulations (see, e.g., Spach et al., 2007) of a collection of 25×25 myocytes. The simulations are started by stimulation of the four leftmost rows of cells (marked by S1) and a resulting depolarization wave propagates through the tissue from left to right. After the center cell is repolarized to −60 mV, the lower left quarter of cells (marked by S2) are stimulated, potentially initiating a second depolarization wave through the tissue. In the healthy case **(A)**, the S2 stimulation does not produce a second wave, whereas in the unhealthy case **(B)**, the solution evolves into to a stable re-entrant wave that continues indefinitely. In the lower panels (**C,D)**, we show the solution of the single myocyte with the largest observed variation in intracellular potential during the simulation at the time point when the largest variation is observed. The maximum single cell variation is much larger in the healthy case **(C)** than in the unhealthy case **(D)**.

## 6. EMI Represents a Grand Challenge in High Performance Computing

In the finite element mesh used in the simulations reported here, the average edge length of a computational element is about 10 μm. This means that the average mesh block has a volume of 1 pL. This should be compared to a mesh block for the bidomain model which is (0.25mm)^{3} = 15, 600 pL, and to the volume of a myocyte which in the present model is 16 pL. Because of the significantly increased resolution, the EMI model is only applicable to relatively small collections of myocytes in a syncytium. Using the finite element code described in the Supplementary Material, we find that the computing times needed for simulating one time step of size 0.001 ms with 100, 625, or 2,500 myocytes are 19, 125, and 568 ms, respectively. Accordingly, the computing time requirements per myocyte are 0.19, 0.20, and 0.23 ms, respectively. These computing efforts fulfills the requirements given by Feynman (2018) who stated that the cost of a computation should be proportional to the space-time volume of the problem under consideration. For the EMI model, this means that the computational efforts should be proportional to the number of myocytes included in our simulations. This optimality criterion is met by the present finite element code and also by a finite difference method applied earlier (Tveito et al., 2017a; Jæger et al., 2021; Kuchta et al., 2021). However, these computing times remain prohibitive for many interesting and important applications, and further work on optimizing computing methods and application of larger HPC-facilities is definitely needed. Important steps toward improved solution technology is underway in the EU-funded *Microcard project (www.microcard.eu)*.

Since the computing efforts needed to solve the EMI model is proportional to the number of myocytes in the simulations, it is obvious that there are well defined restrictions on the classes of application can be studied by this model. Restricted collections of myocytes like the special nodes of the cardiac conduction system (SA- and AV-nodes, Purkinje fibers etc.), areas like the outlet of the pulmonary veins, border zones associated with ischemic regions can be analyzed. In contrast, whole organ simulations are probably only accessible for very small model animals, such as the zebrafish and perhaps atrial or right ventricular tissue from adult mice.

## 7. Still Averaging

Finally, it should be pointed out that even with explicit representation of all myocytes in a simulation, the EMI model is still based on averaging of the processes taking place at a length scale that remains unreachable for tissue scale simulations. Specifically, Ca^{2+}-driven arrhythmia is known to involve criticality in the ensemble behavior of subcellular (micron scale) Ca^{2+}-handling structures. Detailed models of sub-cellular Ca^{2+}-dynamics have been developed (see, e.g., Colman et al., 2020) but application of these models in syncytium of myocytes remains challenging.

## 8. Discussion

At present, almost all tissue simulations of cardiac electrophysiology are based on models like the bidomain and monodomain models which entail spatial averaging and related limitations in spatial resolution. However, these simulations used to be performed using models where individual myocytes were coupled in cables and the electrical conduction was assumed to be one dimensional (see, e.g., Lieberman et al., 1973; Joyner, 1982; Rudy and Quan, 1987; Quan and Rudy, 1990; Shaw and Rudy, 1997). These models based on representation of individual myocytes have also been extended to two and three space dimensions (Roberts et al., 2008; Stinstra et al., 2010; Hubbard and Henriquez, 2012) while retaining the myocytes in the model. However, these models are complex from an implementational point of view, and also quite demanding in terms of CPU requirements. Our EMI model is based on formulations suggested by Krassowska and Neu (1994), Henríquez et al. (2013), and Agudelo-Toro and Neef (2013) and we first applied it to evaluate the accuracy of classical models of the transmembrane potential of neurons (see Tveito et al., 2017b).

As explained above, a major disadvantage of the EMI approach compared to the bidomain or monodoman model is the CPU requirements needed to perform simulations. It is, however, worthwhile to recall that this also used to be a disadvantage of the “averaged” models. In 1984, it was estimated (see Barr and Plonsey, 1984) that simulating 10 ms using 10^{6} mesh points, applying a brute force method, would take 3,000 years! But already in 2006, a 26 ×10^{6} nodes simulation (600 ms) was performed in only 2 days (Potse et al., 2006), and a few years later even more challenging computations were performed in 5 min (see Niederer S. A. et al., 2011). Today, such simulations are performed routinely. A comprehensive and interesting discussion of the history of tissue models is provided in Henriquez (2014).

We suggest that a series of important physiological questions can be addressed using the EMI model. In particular, we believe the EMI model can yield important new insights into mechanism of action of cardiac drugs, how early after depolarizations (EADs) and delayed after depolarizations (DADs) are generated in collections of myocytes, how deterioration of gap-junction coupling can modulate cardiac conduction, how reentry can be generated in very small collections of myocytes, and whether ephaptic coupling can maintain conduction in diseased tissue. Furthermore, we believe that a more complete understanding of the role of T-tubules in ventricular myocytes can be achieved if the EMI model is applied along with improved spatial modeling of ionic concentration changes in the associated restricted extracellular spaces (see Ellingsrud et al., 2020; Setterberg et al., 2021).

## Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

## Author Contributions

KJ and AT did the mathematical modeling, performed the numerical simulations, and wrote the paper. AE and WG reviewed and edited the paper. All authors contributed to the article and approved the submitted version.

## Conflict of Interest

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.

## Publisher's Note

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.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.763584/full#supplementary-material

## References

Agudelo-Toro, A., and Neef, A. (2013). Computationally efficient simulation of electrical activity at cell membranes interacting with self-generated and externally imposed electric fields. *J. Neural Eng*. 10:026019. doi: 10.1088/1741-2560/10/2/026019

Barr, R. C., and Plonsey, R. (1984). Propagation of excitation in idealized anisotropic two-dimensional tissue. *Biophys. J*. 45, 1191–1202. doi: 10.1016/S0006-3495(84)84268-X

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. *Heart Rhythm* 4, 1553–1562. doi: 10.1016/j.hrthm.2007.08.017

Clayton, R., and Panfilov, A. (2008). A guide to modelling cardiac electrical activity in anatomically detailed ventricles. *Prog. Biophys. Mol. Biol*. 96, 19–43. doi: 10.1016/j.pbiomolbio.2007.07.004

Colman, M. A., Holmes, M., Whittaker, D. G., Jayasinghe, I., and Benson, A. P. (2020). Multi-scale approaches for the simulation of cardiac electrophysiology: I-sub-cellular and stochastic calcium dynamics from cell to organ. *Methods* 185, 49–59. doi: 10.1016/j.ymeth.2020.02.011

Ehrlich, J. R., Cha, T.-J., Zhang, L., Chartier, D., Melnyk, P., Hohnloser, S. H., et al. (2003). Cellular electrophysiology of canine pulmonary vein cardiomyocytes: action potential and ionic current properties. *J. Physiol*. 551, 801–813. doi: 10.1113/jphysiol.2003.046417

Ellingsrud, A. J., Solbrå, A., Einevoll, G. T., Halnes, G., and Rognes, M. E. (2020). Finite element simulation of ionic electrodiffusion in cellular geometries. *Front. Neuroinformatics* 14:11. doi: 10.3389/fninf.2020.00011

Feynman, R. P. (2018). “Simulating physics with computers,” in *Feynman and Computation*, ed A. J. G. Hey (Boca Raton, FL: CRC Press), 133–153. doi: 10.1201/9780429500459-11

Franzone, P. C., Pavarino, L. F., and Scacchi, S. (2014). *Mathematical Cardiac Electrophysiology, Vol. 13*. Cham: Springer.

Haissaguerre, 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

Hamabe, A., Okuyama, Y., Miyauchi, Y., Zhou, S., Pak, H.-N., Karagueuzian, H. S., et al. (2003). Correlation between anatomy and electrical activation in canine pulmonary veins. *Circulation* 107, 1550–1555. doi: 10.1161/01.CIR.0000056765.97013.5E

Heijman, J., Sutanto, H., Crijns, H. J., Nattel, S., and Trayanova, N. A. (2021). Computational models of atrial fibrillation: achievements, challenges, and perspectives for improving clinical care. *Cardiovasc. Res*. 117, 1682–1699. doi: 10.1093/cvr/cvab138

Henriquez, C. S. (2014). A brief history of tissue models for cardiac electrophysiology. *IEEE Trans. Biomed. Eng*. 61, 1457–1465. doi: 10.1109/TBME.2014.2310515

Henríquez, F., Jerez-Hanckes, C., and Altermatt, M. F. R. (2013). “Dynamic finite-element model of axon extracellular stimulation,” in *2013 6th International IEEE/EMBS Conference on Neural Engineering (NER)* (San Diego, CA), 589–592. doi: 10.1109/NER.2013.6696003

Hubbard, M. L., and Henriquez, C. S. (2012). Microscopic variations in interstitial and intracellular structure modulate the distribution of conduction delays and block in cardiac tissue with source-load mismatch. *Europace* 14(Suppl. 5), v3–v9. doi: 10.1093/europace/eus267

Jæger, K. H., Edwards, A. G., McCulloch, A., and Tveito, A. (2019). Properties of cardiac conduction in a cell-based computational model. *PLoS Comput. Biol*. 15:e1007042. doi: 10.1371/journal.pcbi.1007042

Jæger, K. H., Hustad, K. G., Cai, X., and Tveito, A. (2021). Efficient numerical solution of the emi model representing the extracellular space (e), cell membrane (m) and intracellular space (i) of a collection of cardiac cells. *Front. Phys*. 8:579461. doi: 10.3389/fphy.2020.579461

Jæger, K. H., and Tveito, A. (2021). “Derivation of a cell-based mathematical model of excitable cells,” in *Modeling Excitable Tissue*, eds A. Tveito, K-A. Mardal, and M. E. Rogne (Cham: Springer), 1–13. doi: 10.1007/978-3-030-61157-6_1

Joyner, R. W. (1982). Effects of the discrete pattern of electrical coupling on propagation through an electrical syncytium. *Circul. Res*. 50, 192–200. doi: 10.1161/01.RES.50.2.192

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

Krassowska, W., and Neu, J. C. (1994). Response of a single cell to an external electric field. *Biophys. J*. 66, 1768–1776. doi: 10.1016/S0006-3495(94)80971-3

Kuchta, M., Mardal, K.-A., and Rognes, M. E. (2021). “Solving the EMI equations using finite element methods,” in *Modeling Excitable Tissue*, eds A. Tveito, K-A. Mardal, and M. E. Rognes (Cham: Springer), 56–69. doi: 10.1007/978-3-030-61157-6_5

Lieberman, M., Kootsey, J. M., Johnson, E. A., and Sawanobori, T. (1973). Slow conduction in cardiac muscle: a biophysical model. *Biophys. J*. 13, 37–55. doi: 10.1016/S0006-3495(73)85968-5

Mandapati, R., Skanes, A., Chen, J., Berenfeld, O., and Jalife, J. (2000). Stable microreentrant sources as a mechanism of atrial fibrillation in the isolated sheep heart. *Circulation* 101, 194–199. doi: 10.1161/01.CIR.101.2.194

Melnyk, P., Ehrlich, J. R., Pourrier, M., Villeneuve, L., Cha, T.-J., and Nattel, S. (2005). Comparison of ion channel distribution and expression in cardiomyocytes of canine pulmonary veins versus left atrium. *Cardiovasc. Res*. 65, 104–116. doi: 10.1016/j.cardiores.2004.08.014

Nattel, S., and Dobrev, D. (2017). Controversies about atrial fibrillation mechanisms: aiming for order in chaos and whether it matters. *Circul. Res*. 120, 1396–1398. doi: 10.1161/CIRCRESAHA.116.310489

Neu, J., and Krassowska, W. (1993). Homogenization of syncytial tissues. *Crit. Rev. Biomed. Eng*. 21, 137–199.

Niederer, S., Mitchell, L., Smith, N., and Plank, G. (2011). Simulating human cardiac electrophysiology on clinical time-scales. *Front. Physiol*. 2:14. doi: 10.3389/fphys.2011.00014

Niederer, S. A., Kerfoot, E., Benson, A. P., Bernabeu, M. O., Bernus, O., Bradley, C., et al. (2011). Verification of cardiac tissue electrophysiology simulators using an n-version benchmark. *Philos. Trans. R. Soc. A* 369, 4331–4351. doi: 10.1098/rsta.2011.0139

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. *Circul. Res*. 82, 63–81. doi: 10.1161/01.RES.82.1.63

Pfenniger, A. (2020). Irregularly irregular-nonuniform cardiac anisotropy as a substrate for atrial fibrillation. *J. Cardiovasc. Electrophysiol*. 31, 2210–2212. doi: 10.1111/jce.14596

Potse, M., Dubé, B., Richer, J., Vinet, A., and Gulrajani, R. M. (2006). A comparison of monodomain and bidomain reaction-diffusion models for action potential propagation in the human heart. *IEEE Trans. Biomed. Eng*. 53, 2425–2435. doi: 10.1109/TBME.2006.880875

Quan, W., and Rudy, Y. (1990). Unidirectional block and reentry of cardiac excitation: a model study. *Circul. Res*. 66, 367–382. doi: 10.1161/01.RES.66.2.367

Rivaud, M. R., Agullo-Pascual, E., Lin, X., Leo-Macias, A., Zhang, M., Rothenberg, E., et al. (2017). Sodium channel remodeling in subcellular microdomains of murine failing cardiomyocytes. *J. Am. Heart Assoc*. 6:e007622. doi: 10.1161/JAHA.117.007622

Roberts, S. F., Stinstra, J. G., and Henriquez, C. S. (2008). Effect of nonuniform interstitial space properties on impulse propagation: a discrete multidomain model. *Biophys. J*. 95, 3724–3737. doi: 10.1529/biophysj.108.137349

Rudy, Y., and Quan, W. (1987). A model study of the effects of the discrete cellular structure on electrical propagation in cardiac tissue. *Circul. Res*. 61, 815–823. doi: 10.1161/01.RES.61.6.815

Setterberg, I. E., Le, C., Frisk, M., Li, J., and Louch, W. E. (2021). The physiology and pathophysiology of t-tubules in the heart. *Front. Physiol*. 12:718404. doi: 10.3389/fphys.2021.718404

Shaw, R. M., and Rudy, Y. (1997). Electrophysiologic effects of acute myocardial ischemia: a mechanistic investigation of action potential conduction and conduction failure. *Circul. Res*. 80, 124–138. doi: 10.1161/01.RES.80.1.124

Spach, M. S., Heidlage, J. F., Dolber, P. C., and Barr, R. C. (2007). Mechanism of origin of conduction disturbances in aging human atrial bundles: experimental and model study. *Heart Rhythm* 4, 175–185. doi: 10.1016/j.hrthm.2006.10.023

Stinstra, J., MacLeod, R., and Henriquez, C. (2010). Incorporating histology into a 3d microscopic computer model of myocardium to study propagation at a cellular level. *Ann. Biomed. Eng*. 38, 1399–1414. doi: 10.1007/s10439-009-9883-y

Sundnes, J., Lines, G. T., Cai, X., Nielsen, B. F., Mardal, K.-A., and Tveito, A. (2007). *Computing the Electrical Activity in the Heart, Vol. 1*. Heidelberg: Springer Science & Business Media.

Tirziu, D., Giordano, F. J., and Simons, M. (2010). Cell communications in the heart. *Circulation* 122, 928–937. doi: 10.1161/CIRCULATIONAHA.108.847731

Tung, L. (1978). *A BI-domain model for describing ischemic myocardial DC potentials* (Ph.D. thesis). Massachusetts Institute of Technology, Cambridge, MA, United States.

Tveito, A., Jæger, K. H., Kuchta, M., Mardal, K.-A., and Rognes, M. E. (2017a). A cell-based framework for numerical modeling of electrical conduction in cardiac tissue. *Front. Phys*. 5:48. doi: 10.3389/fphy.2017.00048

Tveito, A., Jæger, K. H., Lines, G. T., Paszkowski, Ł., Sundnes, J., Edwards, A. G., et al. (2017b). An evaluation of the accuracy of classical models for computing the membrane potential and extracellular potential for neurons. *Front. Comput. Neurosci*. 11:27. doi: 10.3389/fncom.2017.00027

Vigmond, E., Dos Santos, R. W., Prassl, A., Deo, M., and Plank, G. (2008). Solvers for the cardiac bidomain equations. *Prog. Biophys. Mol. Biol*. 96, 3–18. doi: 10.1016/j.pbiomolbio.2007.07.012

Keywords: computational modeling, computational electrophysiology, cardiac modeling, action potential propagation, cell based model, EMI model, cardiac disease, conduction abnormalities

Citation: Jæger KH, Edwards AG, Giles WR and Tveito A (2021) From Millimeters to Micrometers; Re-introducing Myocytes in Models of Cardiac Electrophysiology. *Front. Physiol.* 12:763584. doi: 10.3389/fphys.2021.763584

Received: 24 August 2021; Accepted: 30 September 2021;

Published: 27 October 2021.

Edited by:

Eun Bo Shim, Kangwon National University, South KoreaReviewed by:

Steven Alexander Niederer, King's College London, United KingdomEdward Joseph Vigmond, Université de Bordeaux, France

Copyright © 2021 Jæger, Edwards, Giles and Tveito. 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: Karoline Horgmo Jæger, karolihj@simula.no