Abstract
When modelling tissue-level cardiac electrophysiology, a continuum approximation to the discrete cell-level equations, known as the bidomain equations, is often used to maintain computational tractability. Analysing the derivation of the bidomain equations allows us to investigate how microstructure, in particular gap junctions that electrically connect cells, affect tissue-level conductivity properties. Using a one-dimensional cable model, we derive a modified form of the bidomain equations that take gap junctions into account, and compare results of simulations using both the discrete and continuum models, finding that the underlying conduction velocity of the action potential ceases to match up between models when gap junctions are introduced at physiologically realistic coupling levels. We show that this effect is magnified by: (i) modelling gap junctions with reduced conductivity; (ii) increasing the conductance of the fast sodium channel; and (iii) an increase in myocyte length. From this, we conclude that the conduction velocity arising from the bidomain equations may not be an accurate representation of the underlying discrete system. In particular, the bidomain equations are less likely to be valid when modelling certain diseased states whose symptoms include a reduction in gap junction coupling or an increase in myocyte length.











Similar content being viewed by others
References
Peracchia, C. (1999). Gap junctions: molecular basis of cell communication in health and disease. Current topics in membranes: Vol. 49 (1st ed.). San Diego: Academic Press.
Beeler, G. W., & Reuter, H. (1977). Reconstruction of the action potential of ventricular myocardial fibres. J. Physiol., 268(1), 177–210.
Bensoussan, A., Lions, J.-L., Papanicolaou, G., & Caughey, T. K. (1979). Asymptotic analysis of periodic structures. J. Appl. Mech., 46(2), 477.
Clerc, L. (1976). Directional differences of impulse spread in trabecular muscle from mammalian heart. J. Physiol., 255(2), 335–346.
Dhein, S. (2004). Cardiac gap junctions: physiology, regulation, pathophysiology and pharmacology (1st ed.). Basel: Karger.
Diaz, P. J., Rudy, Y., & Plonsey, R. (1983). Intercalated discs as a cause for discontinuous propagation in cardiac muscle: a theoretical simulation. Ann. Biomed. Eng., 11(3), 177–189.
Gutstein, D. E., Morley, G. E., Vaidya, D., Liu, F., Chen, F. L., Stuhlmann, H., & Fishman, G. I. (2001). Heterogeneous expression of gap junction channels in the heart leads to conduction defects and ventricular dysfunction. Circulation, 104(10), 1194–1199.
Hand, P. E., & Griffith, B. E. (2010). Adaptive multiscale model for simulating cardiac conduction. Proc. Natl. Acad. Sci., 107(33), 14603–14608.
Hand, P. E., & Griffith, B. E. (2011). Empirical study of an adaptive multiscale model for simulating cardiac conduction. Bull. Math. Biol., 73(12), 3071–3089.
Harris, A. L., & Locke, D. (2009). Connexins. Clifton: Humana Press.
Hoyt, R. H., Cohen, M. L., & Saffitz, J. E. (1989). Distribution and three-dimensional structure of intercellular junctions in canine myocardium. Circ. Res., 64(3), 563–574.
Hubbard, M. L. L., Ying, W., & Henriquez, C. S. (2007). Effect of gap junction distribution on impulse propagation in a monolayer of myocytes: a model study. Europace 9(6).
Jansen, J. A., van Veen, T. A. B., de Bakker, J. M. T., & van Rijen, H. V. M. (2010). Cardiac connexins and impulse propagation. J. Mol. Cell. Cardiol., 48(1), 76–82.
Jongsma, H. J., & Wilders, R. (2000). Gap junctions in cardiovascular disease. Circ. Res., 86(12), 1193–1197.
Keener, J., & Panfilov, A. (1996). A biophysical model for defibrillation of cardiac tissue. Biophys. J., 71(3), 1335–1345.
Keener, J., & Sneyd, J. (2001). Mathematical physiology. Berlin: Springer. Corrected edition.
Kléber, A. G., & Rudy, Y. (2004). Basic mechanisms of cardiac impulse propagation and associated arrhythmias. Physiol. Rev., 84(2), 431–488.
Levy, D., Garrison, R. J., Savage, D. D., Kannel, W. B., & Castelli, W. P. (1990). Prognostic implications of echocardiographically determined left ventricular mass in the framingham heart study. N. Engl. J. Med., 322(22), 1561–1566.
Neu, J. C., & Krassowska, W. (1993). Homogenization of syncytial tissues. Crit. Rev. Biomed. Eng., 21(2), 137–199.
Noble, D., & Rudy, Y. (2001). Models of cardiac ventricular action potentials: iterative interaction between experiment and simulation. Math. Phys. Eng. Sci., 359(1783), 1127–1142.
Pathmanathan, P., Bernabeu, M. O., Bordas, R., Cooper, J., Garny, A., Pitt-Francis, J. M., Whiteley, J. P., & Gavaghan, D. J. (2010). A numerical guide to the solution of the bidomain equations of cardiac electrophysiology. Prog. Biophys. Mol. Biol., 102(2–3), 136–155.
Pathmanathan, P., Mirams, G. R., Southern, J., & Whiteley, J. P. (2011). The significant effect of the choice of ionic current integration method in cardiac electro-physiological simulations. Int. J. Numer. Methods Biomed. Eng., 27(11), 1751–1770.
Reddy, J. N. (1993). Introduction to the finite element method (2nd ed.). New York: McGraw-Hill.
Richardson, G., & Chapman, S. J. (2011). Derivation of the bidomain equations for a beating heart with a general microstructure. SIAM J. Appl. Math., 71(3), 657–675.
Rohr, S. (2004). Role of gap junctions in the propagation of the cardiac action potential. Cardiovasc. Res., 62(2), 309–322.
Satoh, H., Delbridge, L. M., Blatter, L. A., & Bers, D. M. (1996). Surface:volume relationship in cardiac myocytes studied with confocal microscopy and membrane capacitance measurements: species-dependence and developmental effects. Biophys. J., 70, 1494–1504.
Seidel, T., Salameh, A., & Dhein, S. (2010). A simulation study of cellular hypertrophy and connexin lateralization in cardiac tissue. Biophys. J., 99(9), 2821–2830.
Sommer, J. R. (1983). Implications of structure and geometry on cardiac electrical activity. Ann. Biomed. Eng., 11(3–4), 149–157.
Yao, J. A., Gutstein, D. E., Liu, F., Fishman, G. I., & Wit, A. L. (2003). Cell coupling between ventricular myocyte pairs fromftable Connexin43-deficient murine hearts. Circ. Res., 93(8), 736–743.
Acknowledgements
Doug Bruce is supported by an EPSRC grant to the Life Sciences Interface Doctoral Training Centre.
Author information
Authors and Affiliations
Corresponding author
Appendix
Appendix
1.1 Proof of the Symmetry of the Homogenised Conductivity Tensors
We will show that the conductivity tensor given in Eq. (8) is symmetric by demonstrating that \(\mathbf{S_{i}}^{(1,2)}\) = \(\mathbf{S_{i}}^{(2,1)}\), with the symmetry of the other components following in a similar way. The extracellular conductivity tensor given in Eq. (12) can then be shown to be symmetric using the same logic. Throughout, we will write W for W (i). For reference, the intracellular tensor is given by
As the identity matrix is symmetric, it remains to show that \(\int \sigma_{i} \frac{\partial W}{\partial\mathbf{z}} \, \mathrm{d} V_{\mathbf {z}}\) is symmetric. Using integration by parts, we have
where σ i =σ(z) and, for discontinuous σ, the derivative is taken to be the weak derivative. Using Eq. (9), this becomes
which upon rearrangement can be expressed as
Applying the divergence theorem then gives
and by grouping terms, we see that this may be written
We see from Eq. (10) that this quantity is zero, so that \(\iiint\sigma \frac{ \partial W_{1}}{\partial z_{2}} \, \mathrm{d} V_{\mathbf{z}} = \iiint\sigma\frac{\partial W_{2}}{\partial z_{1}} \, \mathrm{d} V_{\mathbf{z,}} \) and hence \(\mathbf{S_{i}}^{(1,2)} = \mathbf{S_{i}}^{(2,1)}\) as required.
1.2 Derivation of the Intracellular Conductivity Tensor in the Case Where Gap Junctions are Present
We now give a derivation of Eq. (23). We consider the governing equations for the weight functions in the cell and the gap junction separately, and again write W for W (i). As each compartment has homogeneous conductivity, W 1 will satisfy Laplace’s equation in each, and so a solution that satisfies the governing PDE and boundary conditions has a derivative
where B 1 and B 2 are constants. To find B 1 and B 2, we can integrate the governing equation for \(W^{i}_{1}\), given in (9), across the boundary between the cell and the gap junction at x=δ to give
and, therefore, obtain the relation between B 1 and B 2:
We then use the fact that W 1 has zero mean in x to give us
Solving for B 1 and B 2 gives
and substituting the above into (8) gives
Rights and permissions
About this article
Cite this article
Bruce, D., Pathmanathan, P. & Whiteley, J.P. Modelling the Effect of Gap Junctions on Tissue-Level Cardiac Electrophysiology. Bull Math Biol 76, 431–454 (2014). https://doi.org/10.1007/s11538-013-9927-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11538-013-9927-1