Gap junctions are key mediators of the intercellular communication in cardiac tissue, and their function is vital to sustain normal cardiac electrical activity. Conduction through gap junctions strongly depends on the hemichannel arrangement and transjunctional voltage, rendering the intercellular conductance highly non-Ohmic.Despite this marked non-linear behavior, current tissue-level models of cardiac conduction are rooted on the assumption that gap-junctions conductance is constant (Ohmic), which results in inaccurate predictions of electrical propagation, particularly in the low junctional-coupling regime observed under pathological conditions. In this work, we present a novel non-Ohmic multiscale (NOM) model of cardiac conduction that is suitable for tissue-level simulations. Using non-linear homogenization theory, we develop June 28, 2019 1/17 a conductivity model that seamlessly upscales the voltage-dependent conductance of gap junctions, without the need of explicitly modeling gap junctions. The NOM model allows for the simulation of electrical propagation in tissue-level cardiac domains that accurately resemble that of cell-based microscopic models for a wide range of junctional coupling scenarios, recovering key conduction features at a fraction of the computational complexity. A unique feature of the NOM model is the possibility of upscaling the response of non-symmetric gap-junction conductance distributions, which result in conduction velocities that strongly depend on the direction of propagation, thus allowing to model the normal and retrograde conduction observed in certain regions of the heart. We envision that the NOM model will enable organ-level simulations that are informed by sub-and inter-cellular mechanisms, delivering an accurate and predictive in-silico tool for understanding the heart function.
Author summaryThe heart relies on the propagation of electrical impulses that are mediated gap junctions, whose conduction properties vary depending on the transjunctional voltage. Despite this non-linear feature, current mathematical models assume that cardiac tissue behaves like an Ohmic (linear) material, thus delivering inaccurate results when simulated in a computer. Here we present a novel mathematical multiscale model that explicitly includes the non-Ohmic response of gap junctions in its predictions. Our results show that the proposed model recovers important conduction features modulated by gap junctions at a fraction of the computational complexity. This contribution represents an important step towards constructing computer models of a whole heart that can predict organ-level behavior in reasonable computing times. 33 cable (monodomain) model of cardiac electrophysiology [9].34 Using two-scale asymptotic homogenization techniques, analytic expressions have 35 been obtained for the effective conductivity tensor, which is then used to model the 36 electrical current in an average macroscopic sense [10-12]. To this end, periodicity at 37 June 28, 2019 3/17 65 compared with those obtained...