Evaluating intramural virtual electrodes in the myocardial wedge preparation: simulations of experimental conditions

Electrical Stimulation for Low-Energy Termination of Cardiac Arrhythmias: a Review

Cardiac arrhythmias are a leading cause of morbidity and mortality in the developed world, estimated to be responsible for hundreds of thousands of deaths annually. Our understanding of the electrophysiological mechanisms of such arrhythmias has grown since they were formally characterized in the late nineteenth century, and this has led to the development of numerous devices and therapies that have markedly improved outcomes for patients affected by such conditions. Despite these advancements, the application of a single large shock remains the clinical standard for treating deadly tachyarrhythmias. Such defibrillating shocks are undoubtedly effective in terminating such arrhythmias; however, they are applied without forewarning, contributing to the patient's stress and anxiety; they can be intensely painful; and they can have adverse psychological and physiological effects on patients. In recent years, there has been interest in developing defibrillation protocols that can terminate arrhythmias without crossing the human pain threshold for energy delivery, generally estimated to be between 0.1 and 1 J. In this article, we review existing literature on the development of such low-energy defibrillation methods and their underlying mechanisms, in an attempt to broadly describe the current landscape of these technologies.

Computational rabbit models to investigate the initiation, perpetuation, and termination of ventricular arrhythmia

Current understanding of cardiac electrophysiology has been greatly aided by computational work performed using rabbit ventricular models. This article reviews the contributions of multiscale models of rabbit ventricles in understanding cardiac arrhythmia mechanisms. This review will provide an overview of multiscale modeling of the rabbit ventricles. It will then highlight works that provide insights into the role of the conduction system, complex geometric structures, and heterogeneous cellular electrophysiology in diseased and healthy rabbit hearts to the initiation and maintenance of ventricular arrhythmia. Finally, it will provide an overview on the contributions of rabbit ventricular modeling on understanding the mechanisms underlying shock-induced defibrillation.

Optical Imaging of Cardiac Action Potential

This chapter reviews the major milestones and scientific achievements facilitated by optical imaging of the action potential in the heart over more than four decades since its introduction. We discuss the limitations of this technique, which sometimes are not fully recognized; the unresolved issues, such as motion artifacts, and the newest developments and future directions.

Three-dimensional histology: tools and application to quantitative assessment of cell-type distribution in rabbit heart

Aims

Cardiac histo-anatomical organization is a major determinant of function. Changes in tissue structure are a relevant factor in normal and disease development, and form targets of therapeutic interventions. The purpose of this study was to test tools aimed to allow quantitative assessment of cell-type distribution from large histology and magnetic resonance imaging- (MRI) based datasets.

Methods and results

Rabbit heart fixation during cardioplegic arrest and MRI were followed by serial sectioning of the whole heart and light-microscopic imaging of trichrome-stained tissue. Segmentation techniques developed specifically for this project were applied to segment myocardial tissue in the MRI and histology datasets. In addition, histology slices were segmented into myocytes, connective tissue, and undefined. A bounding surface, containing the whole heart, was established for both MRI and histology. Volumes contained in the bounding surface (called ‘anatomical volume’), as well as that identified as containing any of the above tissue categories (called ‘morphological volume’), were calculated. The anatomical volume was 7.8 cm³ in MRI, and this reduced to 4.9 cm³ after histological processing, representing an ‘anatomical’ shrinkage by 37.2%. The morphological volume decreased by 48% between MRI and histology, highlighting the presence of additional tissue-level shrinkage (e.g. an increase in interstitial cleft space). The ratio of pixels classified as containing myocytes to pixels identified as non-myocytes was roughly 6:1 (61.6 vs. 9.8%; the remaining fraction of 28.6% was ‘undefined’).

Conclusion

Qualitative and quantitative differentiation between myocytes and connective tissue, using state-of-the-art high-resolution serial histology techniques, allows identification of cell-type distribution in whole-heart datasets. Comparison with MRI illustrates a pronounced reduction in anatomical and morphological volumes during histology processing.

Cardiac response to low-energy field pacing challenges the standard theory of defibrillation

BACKGROUND

The electric response of myocardial tissue to periodic field stimuli has attracted significant attention as the basis for low-energy antifibrillation pacing, potentially more effective than traditional single high-energy shocks. In conventional models, an electric field produces a highly nonuniform response of the myocardial wall, with discrete excitations, or hot spots (HS), occurring at cathodal tissue surfaces or large coronary vessels. We test this prediction using novel 3-dimensional tomographic optical imaging.

METHODS AND RESULTS

Experiments were performed in isolated coronary perfused pig ventricular wall preparations stained with near-infrared voltage-sensitive fluorescent dye DI-4-ANBDQBS. The 3-dimensional coordinates of HS were determined using alternating transillumination. To relate HS formation with myocardial structures, we used ultradeep confocal imaging (interrogation depths, >4 mm). The peak HS distribution is located deep inside the heart wall, and the depth is not significantly affected by field polarity. We did not observe the strong colocalization of HS with major coronary vessels anticipated from theory. Yet, we observed considerable lateral displacement of HS with field polarity reversal. Models that de-emphasized lateral intracellular coupling and accounted for resistive heterogeneity in the extracellular space showed similar HS distributions to the experimental observations.

CONCLUSIONS

The HS distributions within the myocardial wall and the significant lateral displacements with field polarity reversal are inconsistent with standard theories of defibrillation. Extended theories based on enhanced descriptions of cellular scale electric mechanisms may be necessary. The considerable lateral displacement of HS with field polarity reversal supports the hypothesis of biphasic stimuli in low-energy antifibrillation pacing being advantageous.

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.

New insights into defibrillation of the heart from realistic simulation studies

Cardiac defibrillation, as accomplished nowadays by automatic, implantable devices, constitutes the most important means of combating sudden cardiac death. Advancing our understanding towards a full appreciation of the mechanisms by which a shock interacts with the heart, particularly under diseased conditions, is a promising approach to achieve an optimal therapy. The aim of this article is to assess the current state-of-the-art in whole-heart defibrillation modelling, focusing on major insights that have been obtained using defibrillation models, primarily those of realistic heart geometry and disease remodelling. The article showcases the contributions that modelling and simulation have made to our understanding of the defibrillation process. The review thus provides an example of biophysically based computational modelling of the heart (i.e. cardiac defibrillation) that has advanced the understanding of cardiac electrophysiological interaction at the organ level, and has the potential to contribute to the betterment of the clinical practice of defibrillation.

On boundary stimulation and optimal boundary control of the bidomain equations

The bidomain equations with Neumann boundary stimulation and optimal control of these stimuli are investigated. First an analytical framework for boundary control is provided. Then a parallel finite element based algorithm is devised and its efficiency is demonstrated not only for the direct problem but also for the optimal control problem. The computations realize a model configuration corresponding to optimal boundary defibrillation of a reentry phenomenon by applying current density stimuli.

Computational modeling of the human atrial anatomy and electrophysiology

This review article gives a comprehensive survey of the progress made in computational modeling of the human atria during the last 10 years. Modeling the anatomy has emerged from simple "peanut"-like structures to very detailed models including atrial wall and fiber direction. Electrophysiological models started with just two cellular models in 1998. Today, five models exist considering e.g. details of intracellular compartments and atrial heterogeneity. On the pathological side, modeling atrial remodeling and fibrotic tissue are the other important aspects. The bridge to data that are measured in the catheter laboratory and on the body surface (ECG) is under construction. Every measurement can be used either for model personalization or for validation. Potential clinical applications are briefly outlined and future research perspectives are suggested.

Representing cardiac bidomain bath-loading effects by an augmented monodomain approach: application to complex ventricular models

Although the cardiac bidomain model has been widely used in the simulation of electrical activation, its relatively computationally expensive nature means that monodomain approaches are generally required for long-duration simulations (for example, investigations of arrhythmia mechanisms). However, the presence of a conducting bath surrounding the tissue is known to induce wavefront curvature (surface leading bulk), a phenomena absent in standard monodomain approaches. Here, we investigate the biophysical origin of the bidomain bath-loading induced wavefront curvature and present a novel augmented monodomain-equivalent bidomain approach faithfully replicating all aspects of bidomain wavefront morphology and conduction velocity, but with a fraction of the computational cost. Bath-loading effects are shown to be highly dependent upon specific conductivity parameters, but less dependent upon the thickness or conductivity of the surrounding bath, with even relatively thin surrounding fluid layers (~ 0.1 mm) producing significant wavefront curvature in bidomain simulations. We demonstrate that our augmented monodomain approach can be easily adapted for different conductivity sets and applied to anatomically complex models, thus facilitating fast and accurate simulation of cardiac wavefront dynamics during long-duration simulations, further aiding the faithful comparison of simulations with experiments.

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.

Modeling the role of the coronary vasculature during external field stimulation

The exact mechanisms by which defibrillation shocks excite cardiac tissue far from both the electrodes and heart surfaces require elucidation. Bidomain theory explains this phenomena through the existence of intramural virtual electrodes (VEs), caused by discontinuities in myocardial tissue structure. In this study, we assess the modeling components essential in constructing a finite-element cardiac tissue model including blood vessels from high-resolution magnetic resonance data and investigate the specific role played by coronary vasculature in VE formation, which currently remains largely unknown. We use a novel method for assigning histologically based fiber architecture around intramural structures and include an experimentally derived vessel lumen wall conductance within the model. Shock-tissue interaction in the presence of vessels is assessed through comparison with a simplified model lacking intramural structures. Results indicate that VEs form around blood vessels for shocks > 8 V/cm. The magnitude of induced polarizations is attenuated by realistic representation of fiber negotiation around vessel cavities, as well as the insulating effects of the vessel lumen wall. Furthermore, VEs formed around large subepicardial vessels reduce epicardial polarization levels. In conclusion, we have found that coronary vasculature acts as an important substrate for VE formation, which may help interpretation of optical mapping data.

High-throughput cardiac science on the Grid

Cardiac electrophysiology is a mature discipline, with the first model of a cardiac cell action potential having been developed in 1962. Current models range from single ion channels, through very complex models of individual cardiac cells, to geometrically and anatomically detailed models of the electrical activity in whole ventricles. A critical issue for model developers is how to choose parameters that allow the model to faithfully reproduce observed physiological effects without over-fitting. In this paper, we discuss the use of a parametric modelling toolkit, called Nimrod, that makes it possible both to explore model behaviour as parameters are changed and also to tune parameters by optimizing model output. Importantly, Nimrod leverages computers on the Grid, accelerating experiments by using available high-performance platforms. We illustrate the use of Nimrod with two case studies, one at the cardiac tissue level and one at the cellular level.

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.

Termination of atrial fibrillation using pulsed low-energy far-field stimulation

BACKGROUND

Electrically based therapies for terminating atrial fibrillation (AF) currently fall into 2 categories: antitachycardia pacing and cardioversion. Antitachycardia pacing uses low-intensity pacing stimuli delivered via a single electrode and is effective for terminating slower tachycardias but is less effective for treating AF. In contrast, cardioversion uses a single high-voltage shock to terminate AF reliably, but the voltages required produce undesirable side effects, including tissue damage and pain. We propose a new method to terminate AF called far-field antifibrillation pacing, which delivers a short train of low-intensity electric pulses at the frequency of antitachycardia pacing but from field electrodes. Prior theoretical work has suggested that this approach can create a large number of activation sites ("virtual" electrodes) that emit propagating waves within the tissue without implanting physical electrodes and thereby may be more effective than point-source stimulation.

METHODS AND RESULTS

Using optical mapping in isolated perfused canine atrial preparations, we show that a series of pulses at low field strength (0.9 to 1.4 V/cm) is sufficient to entrain and subsequently extinguish AF with a success rate of 93% (69 of 74 trials in 8 preparations). We further demonstrate that the mechanism behind far-field antifibrillation pacing success is the generation of wave emission sites within the tissue by the applied electric field, which entrains the tissue as the field is pulsed.

CONCLUSIONS

AF in our model can be terminated by far-field antifibrillation pacing with only 13% of the energy required for cardioversion. Further studies are needed to determine whether this marked reduction in energy can increase the effectiveness and safety of terminating atrial tachyarrhythmias clinically.

Generation of histo-anatomically representative models of the individual heart: tools and application

This paper presents methods to build histo-anatomically detailed individualized cardiac models. The models are based on high-resolution three-dimensional anatomical and/or diffusion tensor magnetic resonance images, combined with serial histological sectioning data, and are used to investigate individualized cardiac function. The current state of the art is reviewed, and its limitations are discussed. We assess the challenges associated with the generation of histo-anatomically representative individualized in silico models of the heart. The entire processing pipeline including image acquisition, image processing, mesh generation, model set-up and execution of computer simulations, and the underlying methods are described. The multifaceted challenges associated with these goals are highlighted, suitable solutions are proposed, and an important application of developed high-resolution structure-function models in elucidating the effect of individual structural heterogeneity upon wavefront dynamics is demonstrated.

From mitochondrial ion channels to arrhythmias in the heart: computational techniques to bridge the spatio-temporal scales

Computer simulations of electrical behaviour in the whole ventricles have become commonplace during the last few years. The goals of this article are (i) to review the techniques that are currently employed to model cardiac electrical activity in the heart, discussing the strengths and weaknesses of the various approaches, and (ii) to implement a novel modelling approach, based on physiological reasoning, that lifts some of the restrictions imposed by current state-of-the-art ionic models. To illustrate the latter approach, the present study uses a recently developed ionic model of the ventricular myocyte that incorporates an excitation-contraction coupling and mitochondrial energetics model. A paradigm to bridge the vastly disparate spatial and temporal scales, from subcellular processes to the entire organ, and from sub-microseconds to minutes, is presented. Achieving sufficient computational efficiency is the key to success in the quest to develop multiscale realistic models that are expected to lead to better understanding of the mechanisms of arrhythmia induction following failure at the organelle level, and ultimately to the development of novel therapeutic applications.

Polarity reversal lowers activation time during diastolic field stimulation of the rabbit ventricles: insights into mechanisms

To fully characterize the mechanisms of defibrillation, it is necessary to understand the response, within the three-dimensional (3D) volume of the ventricles, to shocks given in diastole. Studies that have examined diastolic responses conducted measurements on the epicardium or on a transmural surface of the left ventricular (LV) wall only. The goal of this study was to use optical imaging experiments and 3D bidomain simulations, including a model of optical mapping, to ascertain the shock-induced virtual electrode and activation patterns throughout the rabbit ventricles following diastolic shocks. We tested the hypothesis that the locations of shock-induced regions of hyperpolarization govern the different diastolic activation patterns for shocks of reversed polarity. In model and experiment, uniform-field monophasic shocks of reversed polarities (cathode over the right ventricle is RV-, reverse polarity is LV-) were applied to the ventricles in diastole. Experiments and simulations revealed that RV- shocks resulted in longer activation times compared with LV- shocks of the same strength. 3D simulations demonstrated that RV- shocks induced a greater volume of hyperpolarization at shock end compared with LV- shocks; most of these hyperpolarized regions were located in the LV. The results of this study indicate that ventricular geometry plays an important role in both the location and size of the shock-induced virtual anodes that determine activation delay during the shock and subsequently affect shock-induced propagation. If regions of hyperpolarization that develop during the shock are sufficiently large, activation delay may persist until shock end.