Do Intramural Virtual Electrodes Facilitate Successful Defibrillation? Model-Based Analysis of Experimental Evidence

An Introductory Overview of Image-Based Computational Modeling in Personalized Cardiovascular Medicine

Cardiovascular diseases account for the number one cause of deaths in the world. Part of the reason for such grim statistics is our limited understanding of the underlying mechanisms causing these devastating pathologies, which is made difficult by the invasiveness of the procedures associated with their diagnosis (e.g., inserting catheters into the coronal artery to measure blood flow to the heart). Likewise, it is also difficult to design and test assistive devices without implanting them in vivo. However, with the recent advancements made in biomedical scanning technologies and computer simulations, image-based modeling (IBM) has arisen as the next logical step in the evolution of non-invasive patient-specific cardiovascular medicine. Yet, due to its novelty, it is still relatively unknown outside of the niche field. Therefore, the goal of this manuscript is to review the current state-of-the-art and the limitations of the methods used in this area of research, as well as their applications to personalized cardiovascular investigations and treatments. Specifically, the modeling of three different physics – electrophysiology, biomechanics and hemodynamics – used in the cardiovascular IBM is discussed in the context of the physiology that each one of them describes and the mechanisms of the underlying cardiac diseases that they can provide insight into. Only the “bare-bones” of the modeling approaches are discussed in order to make this introductory material more accessible to an outside observer. Additionally, the imaging methods, the aspects of the unique cardiac anatomy derived from them, and their relation to the modeling algorithms are reviewed. Finally, conclusions are drawn about the future evolution of these methods and their potential toward revolutionizing the non-invasive diagnosis, virtual design of treatments/assistive devices, and increasing our understanding of these lethal cardiovascular diseases.

Conceptual Intra-Cardiac Electrode Configurations That Facilitate Directional Cardiac Stimulation for Optimal Electrotherapy

OBJECTIVE

Electrotherapy remains the most effective direct therapy against lethal cardiac arrhythmias. When an arrhythmic event is sensed, either strong electric shocks or controlled rapid pacing is automatically applied directly to the heart via an implanted cardioverter defibrillator (ICDs). Despite their success, ICDs remain a highly non-optimal therapy: the strong shocks required for defibrillation cause significant extra-cardiac stimulation, resulting in pain and long-term tissue damage, and can also limit battery life. When used in anti-tachycardia pacing mode, ICDs are also often ineffective, as the pacing electrode can be far away from the centre of the arrhythmia, making it hard for the paced wave to interrupt and terminate it.

METHODS

In this paper, we present two conceptual intra-cardiac directional electrode configurations in silico based on novel arrangements of pairs of positive-negative electrodes. Both configurations have the potential to cause preferential excitation on specific regions of the heart.

RESULTS

We demonstrate how the properties of the induced field varies spatially around the electrodes and how it depends upon the specific arrangements of dipole electrode pairs. The results show that when tested within anatomically-realistic rabbit ventricular models, both electrode configurations produce strong virtual electrodes on the targeted endocardial surfaces, with weaker virtual electrodes produced elsewhere.

CONCLUSIONS

The proposed electrode configurations may facilitate targeted far-field anti-tachycardia pacing and/or defibrillation, which may be useful in cases where conventional anti-tachycardia pacing fails. In addition, the conceptual electrode designs intrinsically confine the electric field to the immediate vicinity of the electrodes, and may, thus, minimize pain due to unnecessary extra-cardiac stimulation.

An efficient finite element approach for modeling fibrotic clefts in the heart

Advanced medical imaging technologies provide a wealth of information on cardiac anatomy and structure at a paracellular resolution, allowing to identify microstructural discontinuities which disrupt the intracellular matrix. Current state-of-the-art computer models built upon such datasets account for increasingly finer anatomical details, however, structural discontinuities at the paracellular level are typically discarded in the model generation process, owing to the significant costs which incur when using high resolutions for explicit representation. In this study, a novel discontinuous finite element (dFE) approach for discretizing the bidomain equations is presented, which accounts for fine-scale structures in a computer model without the need to increase spatial resolution. In the dFE method, this is achieved by imposing infinitely thin lines of electrical insulation along edges of finite elements which approximate the geometry of discontinuities in the intracellular matrix. Simulation results demonstrate that the dFE approach accounts for effects induced by microscopic size scale discontinuities, such as the formation of microscopic virtual electrodes, with vast computational savings as compared to high resolution continuous finite element models. Moreover, the method can be implemented in any standard continuous finite element code with minor effort.

Simulating a dual-array electrode configuration to investigate the influence of skeletal muscle fatigue following functional electrical stimulation

A novel, anatomically-accurate model of a tibialis anterior muscle is used to investigate the electro-physiological properties of denervated muscles following functional electrical stimulation. The model includes a state-of-the-art description of cell electro-physiology. The main objective of this work is to develop a computational framework capable of predicting the effects of different stimulation trains and electrode configurations on the excitability and fatigue of skeletal muscle tissue. Utilizing a reduced but computationally amenable model, the effects of different electrode sizes and inter-electrode distances on the number of activated muscle fibers are investigated and qualitatively compared to existing literature. To analyze muscle fatigue, the sodium current, specifically the K+ ion concentrations within the t-tubule and the calcium release from the sarcoplasmic reticulum, is used to quantify membrane and metabolic fatigue. The simulations demonstrate that lower stimulation frequencies and biphasic pulse waveforms cause less fatigue than higher stimulation frequencies and monophasic pulses. A comparison between single and dual electrode configurations (with the same overall stimulation surface) is presented to locally investigate the differences in muscle fatigue. The dual electrode configuration causes the muscle tissue to fatigue quicker.

Probing field-induced tissue polarization using transillumination fluorescent imaging

Despite major successes of biophysical theories in predicting the effects of electrical shocks within the heart, recent optical mapping studies have revealed two major discrepancies between theory and experiment: 1), the presence of negative bulk polarization recorded during strong shocks; and 2), the unexpectedly small surface polarization under shock electrodes. There is little consensus as to whether these differences result from deficiencies of experimental techniques, artifacts of tissue damage, or deficiencies of existing theories. Here, we take advantage of recently developed near-infrared voltage-sensitive dyes and transillumination optical imaging to perform, for the first time that we know of, noninvasive probing of field effects deep inside the intact ventricular wall. This technique removes some of the limitations encountered in previous experimental studies. We explicitly demonstrate that deep inside intact myocardial tissue preparations, strong electrical shocks do produce considerable negative bulk polarization previously inferred from surface recordings. We also demonstrate that near-threshold diastolic field stimulation produces activation of deep myocardial layers 2-6 mm away from the cathodal surface, contrary to theory. Using bidomain simulations we explore factors that may improve the agreement between theory and experiment. We show that the inclusion of negative asymmetric current can qualitatively explain negative bulk polarization in a discontinuous bidomain model.

Shock induced electrical activation in structurally detailed models of pig left-ventricular tissue

Detailed models of sample specific structures in pig left-ventricular tissue have been constructed. These models include epicardial and endocardial surfaces, fiber and sheet orientations, vessels and cleavage planes with significant dimensions. This work shows that it is possible to extract from 3D tissue images reduced dimension descriptions of cleavage planes in the heart wall. These descriptions are used to analyze the response of tissue to electrical shocks of varying strengths. The presence of explicit discontinuities in the heart significantly reduces the time required for transmural activation and provides a basis for understanding successful defibrillation.

Development of an anatomically detailed MRI-derived rabbit ventricular model and assessment of its impact on simulations of electrophysiological function

Recent advances in magnetic resonance (MR) imaging technology have unveiled a wealth of information regarding cardiac histoanatomical complexity. However, methods to faithfully translate this level of fine-scale structural detail into computational whole ventricular models are still in their infancy, and, thus, the relevance of this additional complexity for simulations of cardiac function has yet to be elucidated. Here, we describe the development of a highly detailed finite-element computational model (resolution: approximately 125 microm) of rabbit ventricles constructed from high-resolution MR data (raw data resolution: 43 x 43 x 36 microm), including the processes of segmentation (using a combination of level-set approaches), identification of relevant anatomical features, mesh generation, and myocyte orientation representation (using a rule-based approach). Full access is provided to the completed model and MR data. Simulation results were compared with those from a simplified model built from the same images but excluding finer anatomical features (vessels/endocardial structures). Initial simulations showed that the presence of trabeculations can provide shortcut paths for excitation, causing regional differences in activation after pacing between models. Endocardial structures gave rise to small-scale virtual electrodes upon the application of external field stimulation, which appeared to protect parts of the endocardium in the complex model from strong polarizations, whereas intramural virtual electrodes caused by blood vessels and extracellular cleft spaces appeared to reduce polarization of the epicardium. Postshock, these differences resulted in the genesis of new excitation wavefronts that were not observed in more simplified models. Furthermore, global differences in the stimulus recovery rates of apex/base regions were observed, causing differences in the ensuing arrhythmogenic episodes. In conclusion, structurally simplified models are well suited for a large range of cardiac modeling applications. However, important differences are seen when behavior at microscales is relevant, particularly when examining the effects of external electrical stimulation on tissue electrophysiology and arrhythmia induction. This highlights the utility of histoanatomically detailed models for investigations of cardiac function, in particular for future patient-specific modeling.

Construction and validation of a plunge electrode array for three-dimensional determination of conductivity in the heart

The heart's response to electrical shock, electrical propagation in sinus rhythm, and the spatiotemporal dynamics of ventricular fibrillation all depend critically on the electrical anisotropy of cardiac tissue. Analysis of the microstructure of the heart predicts that three unique intracellular electrical conductances can be defined at any point in the ventricular wall; however, to date, there has been no experimental confirmation of this concept. We report the design, fabrication, and validation of a novel plunge electrode array capable of addressing this issue. A new technique involving nylon coating of 24G hypodermic needles is performed to achieve nonconductive electrodes that can be combined to give moderate-density multisite intramural measurement of extracellular potential in the heart. Each needle houses 13 silver wires within a total diameter of 0.7 mm, and the combined electrode array gives 137 sites of recording. The ability of the electrode array to accurately assess conductances is validated by mapping the potential field induced by a point current source within baths of saline of varying concentration. A bidomain model of current injection in the heart is then used to test an approximate relationship between the monodomain conductivities measured by the array, and the full set of bidomain conductivities that describe cardiac tissue.

Multilevel homogenization applied to the cardiac bidomain equations

Accurate cardiac tissue-based modeling using the bidomain equations requires the incorporation of fine-scale structures observed at the 50-100 micron level. By including such features we can more easily observe how defibrillation shocks lead to total depolarization of the heart. Several modeling studies that have investigated the effect of fine scale structures on defibrillation success have been completed. Results have shown that such structures aid, through the creation of virtual electrodes, in total depolarization. An obstacle that occurs with this modeling style is the massive amount of data that must be incorporated into detailed tissue models for even a cubic millimeter sample of cardiac tissue. In this paper, we discuss our approach to generating upscaled, or homogenized, versions of these models that can be used to perform simulations at a more reasonable modeling scale. They have the advantage of incorporating fine scale structure into the model at a reduced modeling cost. We introduce and briefly explore the advantages of this upscaling method.

Laminar Arrangement of Ventricular Myocytes Influences Electrical Behavior of the Heart

The response of the heart to electrical shock, electrical propagation in sinus rhythm, and the spatiotemporal dynamics of ventricular fibrillation all depend critically on the electrical anisotropy of cardiac tissue. A long-held view of cardiac electrical anisotropy is that electrical conductivity is greatest along the myocyte axis allowing most rapid propagation of electrical activation in this direction, and that conductivity is isotropic transverse to the myocyte axis supporting a slower uniform spread of activation in this plane. In this context, knowledge of conductivity in two directions, parallel and transverse to the myofiber axis, is sufficient to characterize the electrical action of the heart. Here we present new experimental data that challenge this view. We have used a novel combination of intramural electrical mapping, and experiment-specific computer modeling, to demonstrate that left ventricular myocardium has unique bulk conductivities associated with three microstructurally-defined axes. We show that voltage fields induced by intramural current injection are influenced by not only myofiber direction, but also the transmural arrangement of muscle layers or myolaminae. Computer models of these experiments, in which measured 3D tissue structure was reconstructed in-silico, best matched recorded voltages with conductivities in the myofiber direction, and parallel and normal to myolaminae, set in the ratio 4:2:1, respectively. These findings redefine cardiac tissue as an electrically orthotropic substrate and enhance our understanding of how external shocks may act to successfully reset the fibrillating heart into a uniform electrical state. More generally, the mechanisms governing the destabilization of coordinated electrical propagation into ventricular arrhythmia need to be evaluated in the light of this discovery.

Myocardial segment-specific model generation for simulating the electrical action of the heart

BACKGROUND

Computer models of the electrical and mechanical actions of the heart, solved on geometrically realistic domains, are becoming an increasingly useful scientific tool. Construction of these models requires detailed measurement of the microstructural features which impact on the function of the heart. Currently a few generic cardiac models are in use for a wide range of simulation problems, and contributions to publicly accessible databases of cardiac structures, on which models can be solved, remain rare. This paper presents to-date the largest database of porcine left ventricular segment microstructural architecture, for use in both electrical and mechanical simulation.

METHODS

Cryosectioning techniques were used to reconstruct the myofibre and myosheet orientations in tissue blocks of size ~15 x 15 x 15 mm, taken from the mid-anterior left ventricular freewall, of seven hearts. Tissue sections were gathered on orthogonal planes, and the angles of intersection of myofibres and myosheets with these planes determined automatically with a gradient intensity based algorithm. These angles were then combined to provide a description of myofibre and myosheet variation throughout the tissue, in a form able to be input to biophysically based computational models of the heart.

RESULTS

Several microstructural features were common across all hearts. Myofibres rotated through 141 +/- 18 degrees (mean +/- SD) from epicardium to endocardium, in near linear fashion. In the outer two-thirds of the wall sheet angles were predominantly negative, however, in the inner one-third an abrupt change in sheet angle, with reversal in sign, was seen in six of the seven hearts. Two distinct populations of sheets with orthogonal orientations often co-existed, usually with one population dominating. The utility of the tissue structures was demonstrated by simulating the passive and active electrical responses of two of the tissue blocks to current injection. Distinct patterns of electrical response were obtained in the two tissue blocks, illustrating the importance of testing model based predictions on a variety of tissue architectures.

CONCLUSION

This study significantly expands the set of geometries on which models of cardiac function can be solved.

Modeling cardiac electrical activity at the cell and tissue levels

Significant tissue structures exist in cardiac ventricular tissue, which are of supracellular dimension. It is hypothesized that these tissue structures contribute to the discontinuous spread of electrical activation, may contribute to arrhythmogenesis, and also provide a substrate for effective cardioversion. However, the influences of these mesoscale tissue structures in intact ventricular tissue are difficult to understand solely on the basis of experimental measurement. Current measurement technology is able to record at both the macroscale tissue level and the microscale cellular or subcellular level, but to date it has not been possible to obtain large volume, direct measurements at the mesoscales. To bridge this scale gap in experimental measurements, we use tissue-specific structure and mathematical modeling. Our models, which can incorporate ion channel models at the cell level into the reaction-diffusion equations at the tissue level, have enabled us to consider key hypotheses regarding discontinuous activation.

Modeling cardiac activation and the impact of a discontinuous myocardium

Methods have been developed for modeling cardiac activation which accounts for detailed myocardial geometric structures derived from specific tissue samples. This modeling allows both study and analysis of the effects of cleavage planes and other structural barriers to myocardial current flow. Specialized numerical and computational procedures have been developed to enable this modeling. The results that have been obtained clearly show the impact that discontinuities have in the formation of transmural virtual sources and also assist in better understanding experimental recordings. We are continuing to increase our capacity for modeling larger tissue samples, particularly those capable of sustaining reentry.

Solving the cardiac bidomain equations for discontinuous conductivities

Fast simulations of cardiac electrical phenomena demand fast matrix solvers for both the elliptic and parabolic parts of the bidomain equations. It is well known that fast matrix solvers for the elliptic part must address multiple physical scales in order to show robust behavior. Recent research on finding the proper solution method for the bidomain equations has addressed this issue whereby multigrid preconditioned conjugate gradients has been used as a solver. In this paper, a more robust multigrid method, called Black Box Multigrid, is presented as an alternative to conventional geometric multigrid, and the effect of discontinuities on solver performance for the elliptic and parabolic part is investigated. Test problems with discontinuities arising from inserted plunge electrodes and naturally occurring myocardial discontinuities are considered. For these problems, we explore the advantages to using a more advanced multigrid method like Black Box Multigrid over conventional geometric multigrid. Results will indicate that for certain discontinuous bidomain problems Black Box Multigrid provides 60% faster simulations than using conventional geometric multigrid. Also, for the problems examined, it will be shown that a direct usage of conventional multigrid leads to faster simulations than an indirect usage of conventional multigrid as a preconditioner unless there are sharp discontinuities.

Cardiac electrophysiology and tissue structure: bridging the scale gap with a joint measurement and modelling paradigm

Significant tissue structures exist in cardiac ventricular tissue that are of supracellular dimension. It is hypothesized that these tissue structures contribute to the discontinuous spread of electrical activation, may contribute to arrhymogenesis and also provide a substrate for effective cardioversion. However, the influences of these mesoscale tissue structures in intact ventricular tissue are difficult to understand solely on the basis of experimental measurement. Current measurement technology is able to record at both the macroscale tissue level and the microscale cellular or subcellular level, but to date it has not been possible to obtain large volume, direct measurements at the mesoscales. To bridge this scale gap in experimental measurements, we use tissue-specific structure and mathematical modelling. Our models have enabled us to consider key hypotheses regarding discontinuous activation. We also consider the future developments of our intact tissue experimental programme.