Cardiac excitation mechanisms, wavefront dynamics and strength–interval curves predicted by 3D orthotropic bidomain simulations

Bidomain modeling of electrical and mechanical properties of cardiac tissue

Throughout the history of cardiac research, there has been a clear need to establish mathematical models to complement experimental studies. In an effort to create a more complete picture of cardiac phenomena, the bidomain model was established in the late 1970s to better understand pacing and defibrillation in the heart. This mathematical model has seen ongoing use in cardiac research, offering mechanistic insight that could not be obtained from experimental pursuits. Introduced from a historical perspective, the origins of the bidomain model are reviewed to provide a foundation for researchers new to the field and those conducting interdisciplinary research. The interplay of theory and experiment with the bidomain model is explored, and the contributions of this model to cardiac biophysics are critically evaluated. Also discussed is the mechanical bidomain model, which is employed to describe mechanotransduction. Current challenges and outstanding questions in the use of the bidomain model are addressed to give a forward-facing perspective of the model in future studies.

A clinical-in silico study on the effectiveness of multipoint bicathodic and cathodic-anodal pacing in cardiac resynchronization therapy

Up to one-third of patients undergoing cardiac resynchronization therapy (CRT) are nonresponders. Multipoint bicathodic and cathodic-anodal left ventricle (LV) stimulations could overcome this clinical challenge, but their effectiveness remains controversial. Here we evaluate the performance of such stimulations through both in vivo and in silico experiments, the latter based on computer electromechanical modeling. Seven patients, all candidates for CRT, received a quadripolar LV lead. Four stimulations were tested: right ventricular (RVS); conventional single point biventricular (S-BS); multipoint biventricular bicathodic (CC-BS) and multipoint biventricular cathodic-anodal (CA-BS). The following parameters were processed: QRS duration; maximal time derivative of arterial pressure (dPdt_max); systolic arterial pressure (P_sys); and stroke volume (SV). Echocardiographic data of each patient were then obtained to create an LV geometric model. Numerical simulations were based on a strongly coupled Bidomain electromechanical coupling model. Considering the in vivo parameters, when comparing S-BS to RVS, there was no significant decrease in SV (from 45 ± 11 to 44 ± 20 ml) and 6% and 4% increases of dPdt_max and P_sys, respectively. Focusing on in silico parameters, with respect to RVS, S-BS exhibited a significant increase of SV, dPdt_max and P_sys. Neither the in vivo nor in silico results showed any significant hemodynamic and electrical difference among S-BS, CC-BS and CA-BS configurations. These results show that CC-BS and CA-BS yield a comparable CRT performance, but they do not always yield improvement in terms of hemodynamic parameters with respect to S-BS. The computational results confirmed the in vivo observations, thus providing theoretical support to the clinical experiments.

Incorporating inductances in tissue-scale models of cardiac electrophysiology

In standard models of cardiac electrophysiology, including the bidomain and monodomain models, local perturbations can propagate at infinite speed. We address this unrealistic property by developing a hyperbolic bidomain model that is based on a generalization of Ohm's law with a Cattaneo-type model for the fluxes. Further, we obtain a hyperbolic monodomain model in the case that the intracellular and extracellular conductivity tensors have the same anisotropy ratio. In one spatial dimension, the hyperbolic monodomain model is equivalent to a cable model that includes axial inductances, and the relaxation times of the Cattaneo fluxes are strictly related to these inductances. A purely linear analysis shows that the inductances are negligible, but models of cardiac electrophysiology are highly nonlinear, and linear predictions may not capture the fully nonlinear dynamics. In fact, contrary to the linear analysis, we show that for simple nonlinear ionic models, an increase in conduction velocity is obtained for small and moderate values of the relaxation time. A similar behavior is also demonstrated with biophysically detailed ionic models. Using the Fenton-Karma model along with a low-order finite element spatial discretization, we numerically analyze differences between the standard monodomain model and the hyperbolic monodomain model. In a simple benchmark test, we show that the propagation of the action potential is strongly influenced by the alignment of the fibers with respect to the mesh in both the parabolic and hyperbolic models when using relatively coarse spatial discretizations. Accurate predictions of the conduction velocity require computational mesh spacings on the order of a single cardiac cell. We also compare the two formulations in the case of spiral break up and atrial fibrillation in an anatomically detailed model of the left atrium, and we examine the effect of intracellular and extracellular inductances on the virtual electrode phenomenon.

Effect of anisotropy on ventricular vulnerability to unidirectional block and reentry by single premature stimulation during normal sinus rhythm in rat heart

Single high-intensity premature stimuli when applied to the ventricles during ventricular drive of an ectopic site, as in Winfree's "pinwheel experiment," usually induce reentry arrhythmias in the normal heart, while single low-intensity stimuli barely do. Yet ventricular arrhythmia vulnerability during normal sinus rhythm remains largely unexplored. With a view to define the role of anisotropy on ventricular vulnerability to unidirectional conduction block and reentry, we revisited the pinwheel experiment with reduced constraints in the in situ rat heart. New features included single premature stimulation during normal sinus rhythm, stimulation and unipolar potential mapping from the same high-resolution epicardial electrode array, and progressive increase in stimulation strength and prematurity from diastolic threshold until arrhythmia induction. Measurements were performed with 1-ms cathodal stimuli at multiple test sites (n = 26) in seven rats. Stimulus-induced virtual electrode polarization during sinus beat recovery phase influenced premature ventricular responses. Specifically, gradual increase in stimulus strength and prematurity progressively induced make, break, and graded-response stimulation mechanisms. Hence unidirectional conduction block occurred as follows: 1) along fiber direction, on right and left ventricular free walls (n = 23), initiating figure-eight reentry (n = 17) and tachycardia (n = 12), and 2) across fiber direction, on lower interventricular septum (n = 3), initiating spiral wave reentry (n = 2) and tachycardia (n = 1). Critical time window (55.1 ± 4.7 ms, 68.2 ± 6.0 ms) and stimulus strength lower limit (4.9 ± 0.6 mA) defined vulnerability to reentry. A novel finding of this study was that ventricular tachycardia evolves and is maintained by episodes of scroll-like wave and focal activation couplets. We also found that single low-intensity premature stimuli can induce repetitive ventricular response (n = 13) characterized by focal activations.NEW & NOTEWORTHY We performed ventricular cathodal point stimulation during sinus rhythm by progressively increasing stimulus strength and prematurity. Virtual electrode polarization and recovery gradient progressively induced make, break, and graded-response stimulation mechanisms. Unidirectional conduction block occurred along or across fiber direction, initiating figure-eight or spiral wave reentry, respectively, and tachycardia sustained by scroll wave and focal activations.

Effects of premature anodal stimulations on cardiac transmembrane potential and intracellular calcium distributions computed by anisotropic Bidomain models

AIMS

Cardiac unipolar electrode stimulations induce a particular structure of the transmembrane potential distribution (Vm), called virtual electrode polarization (VEP), which plays an important role in the mechanisms of cardiac excitation, reentry induction, and ventricular defibrillation. Recent experimental studies, based on the optical mapping techniques, have shown that premature stimulations also induce significant changes in the intracellular calcium (Cai) spatial distribution. The aim of this work is to investigate and compare by means of numerical simulations the morphology of the Vm and Cai patterns, generated by applying an S1-S2 stimulation protocol with a premature S2 anodal pulse.

METHODS AND RESULTS

We perform parallel finite element simulations of a three-dimensional orthotropic Bidomain model on a block of ventricular tissue by using four membrane models of two species (guinea pig and rabbit), that incorporate the phenomenological or more detailed mechanistic descriptions of the calcium dynamics. During the S2 anodal stimulus, the Cai spatial distribution, computed with all the considered models, presents a configuration similar to the typical VEP pattern of Vm, with a minimum inside the virtual anode and two maxima in the virtual cathodes. After the S2 stimulus turns off, the anode break excitation mechanism yields a Vm pattern exhibiting a clearly propagating wavefront. Differently, the Cai patterns do not show a clear separation between the resting and the activated regions, with the exception of one of the phenomenological models considered, but they show warped dog-bone shaped equi-level lines around an elevation in the virtual anode region.

CONCLUSION

The VEP pattern of the Cai spatial distribution during the S2 stimulus is in agreement with the previous experimental studies. Moreover, the Cai minimum in the virtual anode can be mainly attributable to the outflow of calcium ions produced by the sodium-calcium (NCX) exchanger, without a significant contribution of the ICaL current.

The strength-interval curve in cardiac tissue

The bidomain model describes the electrical properties of cardiac tissue and is often used to simulate the response of the heart to an electric shock. The strength-interval curve summarizes how refractory tissue is excited. This paper analyzes calculations of the strength-interval curve when a stimulus is applied through a unipolar electrode. In particular, the bidomain model is used to clarify why the cathodal and anodal strength-interval curves are different, and what the mechanism of the “dip” in the anodal strength-interval curve is.