Current injection into a two-dimensional anisotropic bidomain.

A two-dimensional sheet of anisotropic cardiac tissue is represented with the bidomain model, and the finite element method is used to solve the bidomain equations. When the anisotropy ratios of the intracellular and extracellular spaces are not equal, the injection of current into the tissue induces a transmembrane potential that has a complicated spatial dependence, including adjacent regions of depolarized and hyperpolarized tissue. This behavior may have important implications for the electrical stimulation of cardiac tissue and for defibrillation.     

Virtual electrodes in cardiac tissue: a common mechanism for anodal and cathodal stimulation.

Traditional cable analyses cannot explain complex patterns of excitation in cardiac tissue with unipolar, extracellular anodal, or cathodal stimuli. Epifluorescence imaging of the transmembrane potential during and after stimulation of both refractory and excitable tissue shows distinctive regions of simultaneous depolarization and hyperpolarization during stimulation that act as virtual cathodes and anodes. The results confirm bidomain model predictions that the onset (make) of a stimulus induces propagation from the virtual cathode, whereas stimulus termination (break) induces it from the virtual anode. In make stimulation, the virtual anode can delay activation of the underlying tissue, whereas in break stimulation this occurs under the virtual cathode. Thus make and break stimulations in cardiac tissue have a common mechanism that is the result of differences in the electrical anisotropy of the intracellular and extracellular spaces and provides clear proof of the validity of the bidomain model.     

Effect of nonuniform interstitial space properties on impulse propagation: a discrete multidomain model

This work presents a discrete multidomain model that describes ionic diffusion pathways between connected cells and within the interstitium. Unlike classical models of impulse propagation, the intracellular and extracellular spaces are represented as spatially distinct volumes with dynamic/static boundary conditions that electrically couple neighboring spaces. The model is used to investigate the impact of nonuniform geometrical and electrical properties of the interstitial space surrounding a fiber on conduction velocity and action potential waveshape. Comparison of the multidomain and bidomain models shows that although the conduction velocity is relatively insensitive to cases that confine 50% of the membrane surface by narrow extracellular depths (≥2 nm), the action potential morphology varies greatly around the fiber perimeter, resulting in changes in the magnitude of extracellular potential in the tight spaces. Results also show that when the conductivity of the tight spaces is sufficiently reduced, the membrane adjacent to the tight space is eliminated from participating in propagation, and the conduction velocity increases. Owing to its ability to describe the spatial discontinuity of cardiac microstructure, the discrete multidomain can be used to determine appropriate tissue properties for use in classical macroscopic models such as the bidomain during normal and pathophysiological conditions.     

High resolution magnetic images of planar wave fronts reveal bidomain properties of cardiac tissue

We magnetically imaged the magnetic action field and optically imaged the transmembrane potentials generated by planar wavefronts on the surface of the left ventricular wall of Langendorff-perfused isolated rabbit hearts. The magnetic action field images were used to produce a time series of two-dimensional action current maps. Overlaying epifluorescent images allowed us to identify a net current along the wavefront and perpendicular to gradients in the transmembrane potential. This is in contrast to a traditional uniform double-layer model where the net current flows along the gradient in the transmembrane potential. Our findings are supported by numerical simulations that treat cardiac tissue as a bidomain with unequal anisotropies in the intra- and extracellular spaces. Our measurements reveal the anisotropic bidomain nature of cardiac tissue during plane wave propagation. These bidomain effects play an important role in the generation of the whole-heart magnetocardiogram and cannot be ignored.     

Potential and current distributions in a cylindrical bundle of cardiac tissue.

The intracellular and interstitial potentials associated with each cell or fiber in multicellular preparations carrying a uniformly propagating wave are important for characterizing the electrophysiological behavior of the preparation and in particular, for evaluating the source contributed by each fiber. The aforementioned potentials depend on a number of factors including the conductivities characterizing the intracellular, interstitial, and extracellular domains, the thickness of the tissue, and the distance (depth) of the field point from the surface of the tissue. A model study is presented describing the extracellular and interstitial potential distribution and current flow in a cylindrical bundle of cardiac muscle arising from a planar wavefront. For simplicity, the bundle is considered as a bidomain. Using typical values of conductivity, the results show that the intracellular and interstitial potential of fibers near the center of a very large bundle (greater than 10 mm) may be approximated by the potentials of a single fiber surrounded by a limited extracellular space (a fiber in oil), hence justifying a core-conductor model. For smaller bundles, the peak interstitial potential is less than that predicted by the core-conductor model but still large enough to affect the overall source strength. The magnitude of the source strength is greatest for fibers lying near the center of the bundle and diminishes sharply for fibers within 50 microns of the surface.     

A numerical solution of the mechanical bidomain model

Introduction: The mechanical bidomain model predicts forces on integrin proteins in the membrane. It has been solved analytically for idealized examples, but a numerical algorithm is needed to address realistic problems. Methods: The bidomain equations are approximated using finite differences. An ischemic region is modeled as a circular area having no active tension, surrounded by normal tissue. Results: The membrane force is large in the ischemic border zone, but is small elsewhere. Strain is distributed widely throughout the ischemic region and surrounding tissue. Conclusion: This calculation provides a testable prediction for the mechanism of mechanotransduction and remodeling in cardiac tissue.     

Modification of a cylindrical bidomain model for cardiac tissue.

Previous models based on a cylindrical bidomain assumed either that the ratio of intracellular and interstitial conductivities in the principal directions were the same or that there was no radial variation in potential (i.e., a planar front, delta Vm/delta rho = 0). This paper presents a formulation and the expressions for the intracellular, interstitial, extracellular, and transmembrane potentials arising from nonplanar propagation along a cylindrical bundle of cardiac tissue represented as a bidomain with arbitrary anisotropy. For unequal anisotropy, the transmembrane current depends not only on the local change of the transmembrane potential but also on the nature of the transmembrane potential throughout the volume.     

How the anisotropy of the intracellular and extracellular conductivities influences stimulation of cardiac muscle

The bidomain model, which describes the behavior of many electrically active tissues, is equivalent to a multi-dimensional cable model and can be represented by a network of resistors and capacitors. For a two-dimensional sheet of tissue, the intracellular and extracellular conductivity tensors can be visualized as two ellipses. For any pair of conductivity tensors, a coordinate transformation can be found that reduces the extracellular ellipse to a circle and aligns the intracellular ellipse with the coordinate axes. The eccentricity of the intracellular ellipse in this new coordinate system is an important parameter. It can have two special values: zero (in which case the tissue has equal anisotropy ratios) or one (in which case the tissue is comprised of one-dimensional fibers coupled through the two-dimensional extracellular space). Thus the bidomain model provides a unifying framework within which the electrical behavior of a wide variety of nerve and muscle tissues can be studied.When the anisotropy ratios in the intracellular and extracellular domains are not equal, stimulation with an anode always causes depolarization of some region of tissue. An analogous effect occurs in models that describe one-dimensional fibers, in which an activating function determines the site of stimulation. Experiments indicate that cardiac muscle does not have equal anisotropy ratios. Therefore, models developed to describe stimulation of axons may also help in understanding stimulation of two- or three-dimensional cardiac tissue, and may explain the concept of anodal stimulation of cardiac tissue through a virtual cathode.     

Genesis of the monophasic action potential: role of interstitial resistance and boundary gradients

The extracellular potential at the site of a mechanical deformation has been shown to resemble the underlying transmembrane action potential, providing a minimally invasive way to access membrane dynamics. The biophysical factors underlying the genesis of this signal, however, are still poorly understood. With the use of data from a recent experimental study in a murine heart, a three-dimensional anisotropic bidomain model of the mouse ventricular free wall was developed to study the currents and potentials resulting from the application of a point mechanical load on cardiac tissue. The applied pressure is assumed to open nonspecific pressure-sensitive channels depolarizing the membrane, leading to monophasic currents at the electrode edge that give rise to the monophasic action potential (MAP). The results show that the magnitude and the time course of the MAP are reproduced only for certain combinations of local or global intracellular and interstitial resistances that form a resting tissue length constant that, if applied over the entire domain, is smaller than that required to match the wave speed. The results suggest that the application of pressure not only causes local depolarization but also changes local tissue properties, both of which appear to play a critical role in the genesis of the MAP.     

The effect of simplifying assumptions in the bidomain model of cardiac tissue: application to ST segment shifts during partial ischaemia

In this study various electrical conductivity approximations used in bidomain models of cardiac tissue are considered. Comparisons are based on epicardial surface potential distributions arising from regions of subendocardial ischaemia situated within the cardiac tissue. Approximations studied are a single conductivity bidomain model, an isotropic bidomain model and equal and reciprocal anisotropy ratios both with and without fibre rotation. It is demonstrated both analytically and numerically that the approximations involving a single conductivity bidomain, an isotropic bidomain or equal anisotropy ratios (ignoring fibre rotation) results in identical epicardial potential distributions for all degrees of subendocardial ischaemia. This result is contrary to experimental observations. It is further shown that by assuming reciprocal anisotropy ratios, epicardial potential distributions vary with the degree of subendocardial ischaemia. However, it is concluded that unequal anisotropy ratios must be used to obtain the true character of experimental observations.     

Electric and magnetic fields from two-dimensional anisotropic bisyncytia.

Cardiac tissue can be considered macroscopically as a bidomain, anisotropic conductor in which simple depolarization wavefronts produce complex current distributions. Since such distributions may be difficult to measure using electrical techniques, we have developed a mathematical model to determine the feasibility of magnetic localization of these currents. By applying the finite element method to an idealized two-dimensional bisyncytium with anisotropic conductivities, we have calculated the intracellular and extracellular potentials, the current distributions, and the magnetic fields for a circular depolarization wavefront. The calculated magnetic field 1 mm from the tissue is well within the sensitivity of a SQUID magnetometer. Our results show that complex bisyncytial current patterns can be studied magnetically, and these studies should provide valuable insight regarding the electrical anisotropy of cardiac tissue.     

Microdomain Effects on Transverse Cardiac Propagation

The effect of gap junctional coupling, sodium ion channel distribution, and extracellular conductivity on transverse conduction in cardiac tissue is explored using a microdomain model that incorporates aspects of the inhomogeneous cellular structure. The propagation velocities found in our model are compared to those in the classic bidomain model and indicate a strong ephaptic microdomain contribution to conduction depending on the parameter regime. We show that ephaptic effects can be quite significant in the junctional spaces between cells, and that the cell activation sequence is modified substantially by these effects. Further, we find that transverse propagation can be maintained by ephaptic effects, even in the absence of gap junctional coupling. The mechanism by which this occurs is found to be cablelike in that the junctional regions act like inverted cables. Our results provide insight into several recent experimental studies that indirectly indicate a mode of action potential propagation that does not rely exclusively on gap junctions.     

Microdomain Effects on Transverse Cardiac Propagation

The effect of gap junctional coupling, sodium ion channel distribution, and extracellular conductivity on transverse conduction in cardiac tissue is explored using a microdomain model that incorporates aspects of the inhomogeneous cellular structure. The propagation velocities found in our model are compared to those in the classic bidomain model and indicate a strong ephaptic microdomain contribution to conduction depending on the parameter regime. We show that ephaptic effects can be quite significant in the junctional spaces between cells, and that the cell activation sequence is modified substantially by these effects. Further, we find that transverse propagation can be maintained by ephaptic effects, even in the absence of gap junctional coupling. The mechanism by which this occurs is found to be cablelike in that the junctional regions act like inverted cables. Our results provide insight into several recent experimental studies that indirectly indicate a mode of action potential propagation that does not rely exclusively on gap junctions.     

Optimal control approach to termination of re-entry waves in cardiac electrophysiology

This work proposes an optimal control approach for the termination of re-entry waves in cardiac electrophysiology. The control enters as an extracellular current density into the bidomain equations which are well established model equations in the literature to describe the electrical behavior of the cardiac tissue. The optimal control formulation is inspired, in part, by the dynamical systems behavior of the underlying system of differential equations. Existence of optimal controls is established and the optimality system is derived formally. The numerical realization is described in detail and numerical experiments, which demonstrate the capability of influencing and terminating reentry phenomena, are presented.     

Quantal currents and potential in the three-dimensional anisotropic bidomain model of smooth muscle

The potential generated in the smooth muscle of the vas deferens on release of a quantum of transmitter from a varicosity was analyzed using a three-dimensional bidomain continuum model. Current was injected at the origin of the bidomain; this current had the temporal characteristics of the junctional current. The membrane potential, intracellular potential, and extracellular potential, as well as the extracellular current, were then calculated throughout the bidomain at different times. Calculations were performed to show the effect of changing the anisotropy ratios of the intracellular and extracellular conductivities on the spread of current and potential in each of the three dimensions. These results provide a theoretical framework for ascertaining the time course of transmitter interaction at a varicosity following the secretion of a quantum of transmitter.     

Mathematical modeling of the excitation process in myocardial tissue: influence of fiber rotation on wavefront propagation and potential field

In our macroscopic model the heart tissue is represented as a bidomain coupling the intra- and extracellular media. Owing to the fiber structure of the myocardium, these media are anisotropic, and their conductivity tensors have a principal axis parallel to the local fiber direction. A reaction-diffusion system is derived that governs the distribution and evolution of the extracellular and transmembrane potentials during the depolarization phase of the heart beat. To investigate frontlike solutions, the system is rescaled and transformed into a system dependent on a small parameter. Subsequently a perturbation analysis is carried out that yields zero- and first-order approximations called eikonal equations. The effects of the transmural fiber rotation on wavefront propagation and the corresponding potential field, elicited by point stimulations, are investigated by means of numerical simulations.     

A sensitivity study of conductivity values in the passive bidomain equation

There is a complex interplay between the four conductivity values used in the bidomain equation and the resulting electric potential distribution in cardiac tissue arising from subendocardial ischaemia. Based on the three commonly used experimentally derived conductivity data sets, a non-dimensional formulation of the passive bidomain equation is derived, which gives rise naturally to several dimensionless conductivity ratios. The data sets are then used to define a parameter space of these ratios, which is studied by considering the correlation coefficients between different epicardial potential distributions.From this study, it is shown that the ratio of the intracellular longitudinal conductivity to the intracellular transverse conductivity is the key parameter in explaining the differences between the epicardial potential distributions observed with these three data sets.     

A solution method for the determination of cardiac potential distributions with an alternating current source

A recently presented solution method for the bidomain model (Johnston et al. 2006), which involves the application of direct current for studying electrical potential in a slab of cardiac tissue, is extended here to allow the use of an applied alternating current. The advantage of using AC current, in a four-electrode method for determining cardiac conductivities, is that instead of using 'close' and 'wide' electrode spacings to make potential measurements, increasing the frequency of the AC current redirects a fraction of the current from the extracellular space into the intracellular space. The model is based on the work of Le Guyader et al. (2001), but is able to include the effects of the fibre rotation between the epicardium and the endocardium on the potentials. Also, rather than using a full numerical technique, the solution method uses Fourier series and a simple one dimensional finite difference scheme, which has the advantage of allowing the potentials to be calculated only at points, such as the measuring electrodes, where they are required. The new alternating current model, which includes intracellular capacitance, is used with a particular four-electrode configuration, to show that the potential measured is affected by changes in fibre rotation. This is significant because it indicates that it is necessary to include fibre rotation in models, which are to be used in conjunction with measuring arrays that are more complex than those involving simply surface probes or a single vertical probe.     

An Analysis and Comparison of Convergence and Uniqueness of Time-Independent Bone Adaptation Models

A recently presented solution method for the bidomain model (Johnston et al. 2006), which involves the application of direct current for studying electrical potential in a slab of cardiac tissue, is extended here to allow the use of an applied alternating current. The advantage of using AC current, in a four-electrode method for determining cardiac conductivities, is that instead of using ‘close’ and ‘wide’ electrode spacings to make potential measurements, increasing the frequency of the AC current redirects a fraction of the current from the extracellular space into the intracellular space.

The model is based on the work of Le Guyader et al. (2001), but is able to include the effects of the fibre rotation between the epicardium and the endocardium on the potentials. Also, rather than using a full numerical technique, the solution method uses Fourier series and a simple one dimensional finite difference scheme, which has the advantage of allowing the potentials to be calculated only at points, such as the measuring electrodes, where they are required.

The new alternating current model, which includes intracellular capacitance, is used with a particular four-electrode configuration, to show that the potential measured is affected by changes in fibre rotation. This is significant because it indicates that it is necessary to include fibre rotation in models, which are to be used in conjunction with measuring arrays that are more complex than those involving simply surface probes or a single vertical probe.     

An approximate solution to the periodic bidomain equations in one dimension

An approximate, computationally tractable solution is proposed for the potentials in the bidomain model with periodic intracellular junctions (the periodic bidomain model). This new approach is based on the one-dimensional rigorous spectral method described previously by Trayanova and Pilkington (IEEE Trans. Biomed. Eng., May 1993). The total solution to the one-dimensional periodic bidomain problem is decomposed in the spectral domain into solutions to (1) the single-fiber classical bidomain problem in which the intracellular conductivity value incorporates the average contribution from cytoplasm and junction and (2) the “junctional” potential problem due to the presence of junctions at discrete locations alone. Solving for the junctional term rigorously requires most of the numerical effort in the solution for the periodic bidomain potentials. Here the junctional potential is found approximately with little numerical effort. A comparison between the rigorous and the approximate solutions serves as a justification for the proposed approximate solution procedure. The procedure outlined in this paper is applicable to higher spatial dimensions where both tissue anisotropy and junctional inhomogeneities play a role in establishing the transmembrane potential distribution.