Comparison between the role of discontinuities in cardiac conduction and in a one-dimensional hardware model

Reentrant excitation in an analog-digital hybrid circuit model of cardiac tissue

We propose an analog-digital hybrid circuit model of one-dimensional cardiac tissue with hardware implementation that allows us to perform real-time simulations of spatially conducting cardiac action potentials. Each active nodal compartment of the tissue model is designed using analog circuits and a dsPIC microcontroller, by which the time-dependent and time-independent nonlinear current-voltage relationships of six types of ion channel currents employed in the Luo-Rudy phase I (LR-I) model for a single mammalian cardiac ventricular cell can be reproduced quantitatively. Here, we perform real-time simulations of reentrant excitation conduction in a ring-shaped tissue model that includes eighty nodal compartments. In particular, we show that the hybrid tissue model can exhibit real-time dynamics for initiation of reentries induced by uni-directional block, as well as those for phase resetting that leads to annihilation of the reentry in response to impulsive current stimulations at appropriate nodes and timings. The dynamics of the hybrid model are comparable to those of a spatially distributed tissue model with LR-I compartments. Thus, it is conceivable that the hybrid model might be a useful tool for large scale simulations of cardiac tissue dynamics, as an alternative to numerical simulations, leading toward further understanding of the reentrant mechanisms.

Continuous description of lattice discreteness effects in front propagation

Models describing microscopic or mesoscopic phenomena in physics are inherently discrete, where the lattice spacing between fundamental components, such as in the case of atomic sites, is a fundamental physical parameter. The effect of spatial discreteness over front propagation phenomenon in an overdamped one-dimensional periodic lattice is studied. We show here that the study of front propagation leads in a discrete description to different conclusions that in the case of its, respectively, continuous description, and also that the results of the discrete model, can be inferred by effective continuous equations with a supplementary spatially periodic term that we have denominated Peierls-Nabarro drift, which describes the bifurcation diagram of the front speed, the appearance of particle-type solutions and their snaking bifurcation diagram. Numerical simulations of the discrete equation show quite good agreement with the phenomenological description.

Transverse propagation in an expanded PSpice model for cardiac muscle with gap-junction ion channels

Transverse propagation was previously found to occur in a two-dimensional model of cardiac muscle using the PSpice software program for electronic circuit design and analysis. Longitudinal propagation within each chain, and transverse propagation between parallel chains, occurred even when there were no gap-junction (g-j) channels inserted between the simulated myocardial cells either longitudinally or transversely. In those studies, there were pronounced edge (boundary) effects and end-effects even within single chains. Transverse velocity increased with increase in model size. The present study was performed to examine boundary effects on transverse propagation velocity when the length of the chains was held constant at 10 cells and the number of parallel chains was varied from 3 to 5, to 7, to 10, and to 20. The number of g-j channels was either zero, both longitudinally and transversely (0/0), or 100/100. Some experiments were also made at 100/0, 1/1, and 10/10. Transverse velocity and overall velocity (both longitudinal and transverse components) was calculated from the measured total propagation time (TPT), i.e., the elapsed time between when the first action potential (AP) and the last AP crossed the zero potential level. The transverse g-j channels were placed only at the ends of each chain, such that propagation would occur in a zigzag pattern. Electrical stimulation was applied intracellularly between cells A1 and A2. It was found that, with no g-j channels (0/0), overall velocity increased almost linearly when more and more chains were placed in parallel. In contrast, with many g-j channels (100/100), there was a much flatter relationship between overall velocity and number of parallel chains. The difference in velocities with 0/0 channels and 100/100 channels was reduced as the number of chains was increased. In conclusion, edges have important effects on propagation velocity (overall and transverse) in cardiac muscle simulations.

Theoretical and experimental study of two discrete coupled Nagumo chains

We analyze front wave (kink and antikink) propagation and pattern formation in a system composed of two coupled discrete Nagumo chains using analytical and numerical methods. In the case of homogeneous interaction among the chains, we show the possibility of the effective control on wave propagation. In addition, physical experiments on electrical chains confirm all theoretical behaviors.

Propagation failure in discrete bistable reaction-diffusion systems: theory and experiments

Wave front propagation failure is investigated in discrete bistable reaction-diffusion systems. We present a theoretical approach including dissipative effects and leading to an analytical expression of the critical coupling beyond which front propagation can occur as a function of the nonlinearity threshold parameter. Our theoretical predictions are confirmed by numerical simulations and experimental results on an equivalent electrical diffusive lattice.

Universal scaling of wave propagation failure in arrays of coupled nonlinear cells

We study the onset of the propagation failure of wave fronts in systems of coupled cells. We introduce a new method to analyze the scaling of the critical external field at which fronts cease to propagate, as a function of intercellular coupling. We find the universal scaling of the field throughout the range of couplings and show that the field becomes exponentially small for large couplings. Our method is generic and applicable to a wide class of cellular dynamics in chemical, biological, and engineering systems. We confirm our results by direct numerical simulations.