The role of photon scattering in optical signal distortion during arrhythmia and defibrillation

Personalized Heart Digital Twins Detect Substrate Abnormalities in Scar-Dependent Ventricular Tachycardia

BACKGROUND

Current outcomes from catheter ablation for scar-dependent ventricular tachycardia (VT) are limited by high recurrence rates and long procedure durations. Personalized heart digital twin technology presents a noninvasive method of predicting critical substrate in VT, and its integration into clinical VT ablation offers a promising solution. The accuracy of the predictions of digital twins to detect invasive substrate abnormalities is unknown. We present the first prospective analysis of digital twin technology in predicting critical substrate abnormalities in VT.

METHODS

Heart digital twin models were created from 18 patients with scar-dependent VT undergoing catheter ablation. Contrast-enhanced cardiac magnetic resonance images were used to reconstruct finite-element meshes, onto which regional electrophysiological properties were applied. Rapid-pacing protocols were used to induce VTs and to define the VT circuits. Predicted optimum ablation sites to terminate all VTs in the models were identified. Invasive substrate mapping was performed, and the digital twins were merged with the electroanatomical map. Electrogram abnormalities and regions of conduction slowing were compared between digital twin-predicted sites and nonpredicted areas.

RESULTS

Electrogram abnormalities were significantly more frequent in digital twin-predicted sites compared with nonpredicted sites (468/1029 [45.5%] versus 519/1611 [32.2%]; P<0.001). Electrogram duration was longer at predicted sites compared with nonpredicted sites (82.0±25.9 milliseconds versus 69.7±22.3 milliseconds; P<0.001). Digital twins correctly identified 21 of 26 (80.8%) deceleration zones seen on isochronal late activation mapping.

CONCLUSIONS

Digital twin-predicted sites display a higher prevalence of abnormal and prolonged electrograms compared with nonpredicted sites and accurately identify regions of conduction slowing. Digital twin technology may help improve substrate-based VT ablation.

REGISTRATION

URL: https://www.clinicaltrials.gov; Unique identifier: NCT04632394.

POPDC1 Variants Cause Atrioventricular Node Dysfunction and Arrhythmogenic Changes in Cardiac Electrophysiology and Intracellular Calcium Handling in Zebrafish

Popeye domain-containing (POPDC) proteins selectively bind cAMP and mediate cellular responses to sympathetic nervous system (SNS) stimulation. The first discovered human genetic variant (POPDC1S201F) is associated with atrioventricular (AV) block, which is exacerbated by increased SNS activity. Zebrafish carrying the homologous mutation (popdc1S191F) display a similar phenotype to humans. To investigate the impact of POPDC1 dysfunction on cardiac electrophysiology and intracellular calcium handling, homozygous popdc1S191F and popdc1 knock-out (popdc1KO) zebrafish larvae and adult isolated popdc1S191F hearts were studied by functional fluorescent analysis. It was found that in popdc1S191F and popdc1KO larvae, heart rate (HR), AV delay, action potential (AP) and calcium transient (CaT) upstroke speed, and AP duration were less than in wild-type larvae, whereas CaT duration was greater. SNS stress by β-adrenergic receptor stimulation with isoproterenol increased HR, lengthened AV delay, slowed AP and CaT upstroke speed, and shortened AP and CaT duration, yet did not result in arrhythmias. In adult popdc1S191F zebrafish hearts, there was a higher incidence of AV block, slower AP upstroke speed, and longer AP duration compared to wild-type hearts, with no differences in CaT. SNS stress increased AV delay and led to further AV block in popdc1S191F hearts while decreasing AP and CaT duration. Overall, we have revealed that arrhythmogenic effects of POPDC1 dysfunction on cardiac electrophysiology and intracellular calcium handling in zebrafish are varied, but already present in early development, and that AV node dysfunction may underlie SNS-induced arrhythmogenesis associated with popdc1 mutation in adults.

Optogenetic Stimulation Using Anion Channelrhodopsin (GtACR1) Facilitates Termination of Reentrant Arrhythmias With Low Light Energy Requirements: A Computational Study

Optogenetic defibrillation of hearts expressing light-sensitive cation channels (e.g., ChR2) has been proposed as an alternative to conventional electrotherapy. Past modeling work has shown that ChR2 stimulation can depolarize enough myocardium to interrupt arrhythmia, but its efficacy is limited by light attenuation and high energy needs. These shortcomings may be mitigated by using new optogenetic proteins like Guillardia theta Anion Channelrhodopsin (GtACR1), which produces a repolarizing outward current upon illumination. Accordingly, we designed a study to assess the feasibility of GtACR1-based optogenetic arrhythmia termination in human hearts. We conducted electrophysiological simulations in MRI-based atrial or ventricular models (n = 3 each), with pathological remodeling from atrial fibrillation or ischemic cardiomyopathy, respectively. We simulated light sensitization via viral gene delivery of three different opsins (ChR2, red-shifted ChR2, GtACR1) and uniform endocardial illumination at the appropriate wavelengths (blue, red, or green light, respectively). To analyze consistency of arrhythmia termination, we varied pulse timing (three evenly spaced intervals spanning the reentrant cycle) and intensity (atrial: 0.001–1 mW/mm²; ventricular: 0.001–10 mW/mm²). In atrial models, GtACR1 stimulation with 0.005 mW/mm² green light consistently terminated reentry; this was 10–100x weaker than the threshold levels for ChR2-mediated defibrillation. In ventricular models, defibrillation was observed in 2/3 models for GtACR1 stimulation at 0.005 mW/mm² (100–200x weaker than ChR2 cases). In the third ventricular model, defibrillation failed in nearly all cases, suggesting that attenuation issues and patient-specific organ/scar geometry may thwart termination in some cases. Across all models, the mechanism of GtACR1-mediated defibrillation was voltage forcing of illuminated tissue toward the modeled channel reversal potential of −40 mV, which made propagation through affected regions impossible. Thus, our findings suggest GtACR1-based optogenetic defibrillation of the human heart may be feasible with ≈2–3 orders of magnitude less energy than ChR2.

OptoGap is an optogenetics-enabled assay for quantification of cell-cell coupling in multicellular cardiac tissue

Intercellular electrical coupling is an essential means of communication between cells. It is important to obtain quantitative knowledge of such coupling between cardiomyocytes and non-excitable cells when, for example, pathological electrical coupling between myofibroblasts and cardiomyocytes yields increased arrhythmia risk or during the integration of donor (e.g., cardiac progenitor) cells with native cardiomyocytes in cell-therapy approaches. Currently, there is no direct method for assessing heterocellular coupling within multicellular tissue. Here we demonstrate experimentally and computationally a new contactless assay for electrical coupling, OptoGap, based on selective illumination of inexcitable cells that express optogenetic actuators and optical sensing of the response of coupled excitable cells (e.g., cardiomyocytes) that are light-insensitive. Cell-cell coupling is quantified by the energy required to elicit an action potential via junctional current from the light-stimulated cell(s). The proposed technique is experimentally validated against the standard indirect approach, GapFRAP, using light-sensitive cardiac fibroblasts and non-transformed cardiomyocytes in a two-dimensional setting. Its potential applicability to the complex three-dimensional setting of the native heart is corroborated by computational modelling and proper calibration. Lastly, the sensitivity of OptoGap to intrinsic cell-scale excitability is robustly characterized via computational analysis.

An evaluation of some assumptions underpinning the bidomain equations of electrophysiology

Tissue level cardiac electrophysiology is usually modelled by the bidomain equations, or the monodomain simplification of the bidomain equations. One assumption made when deriving the bidomain equations is that both the intracellular and extracellular spaces are in electrical equilibrium. This assumption neglects the disturbance of this equilibrium in thin regions close to the cell membrane known as Debye layers. We first demonstrate that the governing equations at the cell, or microscale, level may be adapted to take account of these Debye layers with little additional complexity, provided the permittivity within the Debye layers satisfies certain conditions that are believed to be satisfied for biological cells. We then homogenize the microscale equations using a technique developed for an almost periodic microstructure. Cardiac tissue is usually modelled as sheets of cardiac fibres stacked on top of one another. A common assumption is that an orthogonal coordinate system can be defined at each point of cardiac tissue, where the first axis is in the fibre direction, the second axis is orthogonal to the first axis but lies in the sheet of cardiac fibres and the third axis is orthogonal to the cardiac sheet. It is assumed further that both the intracellular and extracellular conductivity tensors are diagonal with respect to these axes and that the diagonal entries of these tensors are constant across the whole tissue. Using the homogenization technique we find that this assumption is usually valid for cardiac tissue, but highlight situations where the assumption may not be valid.

Tissue-Specific Optical Mapping Models of Swine Atria Informed by Optical Coherence Tomography

Computational models and experimental optical mapping of cardiac electrophysiology serve as powerful tools to investigate the underlying mechanisms of arrhythmias. Modeling can also aid the interpretation of optical mapping signals, which may have different characteristics with respect to the underlying electrophysiological signals they represent. However, despite the prevalence of atrial arrhythmias such as atrial fibrillation, models of optical electrical mapping incorporating realistic structure of the atria are lacking. Therefore, we developed image-based models of atrial tissue using structural information extracted from optical coherence tomography (OCT), which can provide volumetric tissue characteristics in high resolution. OCT volumetric data of four swine atrial tissue samples were used to develop models incorporating tissue geometry, tissue-specific myofiber orientation, and ablation lesion regions. We demonstrated the use of these models through electrophysiology and photon scattering simulations. Changes in transmural electrical conduction were observed with the inclusion of OCT-derived, depth-resolved fiber orientation. Additionally, the amplitude of optical mapping signals were not found to correspond with lesion transmurality because of lesion geometry and electrical propagation occurring beyond excitation light penetration. This work established a framework for the development of tissue-specific models of atrial tissue derived from OCT imaging data, which can be useful in future investigations of electrophysiology and optical mapping signals with respect to realistic atrial tissue structure.

Virtual electrophysiological study as a tool for evaluating efficacy of MRI techniques in predicting adverse arrhythmic events in ischemic patients

Myocardial infarct (MI) related indices determined by late gadolinium enhancement (LGE) MRI have been widely investigated in determining patients suitable for implantable cardiovascular-defibrillator (ICD) therapy to complement left ventricular ejection fraction (LV EF). In comparison to LGE-MRI using inversion-recovery fast-gradient-echo (IR-FGRE), T1 mapping techniques, such as multi contrast late enhancement (MCLE), have been shown to provide more quantitative and reproducible estimates of infarct regions. The objective of this study is to use individualized heart computer models in determining the efficacy of IR-FGRE and MCLE techniques in predicting the occurrence of post-MI ventricular tachycardia (VT). Twenty-seven patients with MI underwent LGE-MRI using IR-FGRE and MCLE prior to ICD implantation and were followed up for 6-46 months. Individualized image-based computational models were built separately for each imaging technique; simulations of propensity to VT were conducted with each model. The imaging methods were evaluated by comparing simulated inducibility of VT to clinical outcome (appropriate ICD therapy) in patients. Twelve patients had at least one appropriate ICD therapy for VT at follow-up. For both MCLE and IR-FGRE, the outcomes of the simulations of VT were significantly associated with the events of appropriate ICD therapy. This indicates that, as compared to conventional measurements such as LV EF, the simulations of VT corresponding to both MCLE and IR-FGRE were more sensitive in predicting appropriate ICD therapy in post-MI patients.

Personalized virtual-heart technology for guiding the ablation of infarct-related ventricular tachycardia

Ventricular tachycardia (VT), which can lead to sudden cardiac death, occurs frequently in patients with myocardial infarction. Catheter-based radiofrequency ablation of cardiac tissue has achieved only modest efficacy, owing to the inaccurate identification of ablation targets by current electrical mapping techniques, which can lead to extensive lesions and to a prolonged, poorly tolerated procedure. Here we show that personalized virtual-heart technology based on cardiac imaging and computational modelling can identify optimal infarct-related VT ablation targets in retrospective animal (5 swine) and human studies (21 patients) and in a prospective feasibility study (5 patients). We first assessed in retrospective studies (one of which included a proportion of clinical images with artifacts) the capability of the technology to determine the minimum-size ablation targets for eradicating all VTs. In the prospective study, VT sites predicted by the technology were targeted directly, without relying on prior electrical mapping. The approach could improve infarct-related VT ablation guidance, where accurate identification of patient-specific optimal targets could be achieved on a personalized virtual heart prior to the clinical procedure.

Validation and Trustworthiness of Multiscale Models of Cardiac Electrophysiology

Computational models of cardiac electrophysiology have a long history in basic science applications and device design and evaluation, but have significant potential for clinical applications in all areas of cardiovascular medicine, including functional imaging and mapping, drug safety evaluation, disease diagnosis, patient selection, and therapy optimisation or personalisation. For all stakeholders to be confident in model-based clinical decisions, cardiac electrophysiological (CEP) models must be demonstrated to be trustworthy and reliable. Credibility, that is, the belief in the predictive capability, of a computational model is primarily established by performing validation, in which model predictions are compared to experimental or clinical data. However, there are numerous challenges to performing validation for highly complex multi-scale physiological models such as CEP models. As a result, credibility of CEP model predictions is usually founded upon a wide range of distinct factors, including various types of validation results, underlying theory, evidence supporting model assumptions, evidence from model calibration, all at a variety of scales from ion channel to cell to organ. Consequently, it is often unclear, or a matter for debate, the extent to which a CEP model can be trusted for a given application. The aim of this article is to clarify potential rationale for the trustworthiness of CEP models by reviewing evidence that has been (or could be) presented to support their credibility. We specifically address the complexity and multi-scale nature of CEP models which makes traditional model evaluation difficult. In addition, we make explicit some of the credibility justification that we believe is implicitly embedded in the CEP modeling literature. Overall, we provide a fresh perspective to CEP model credibility, and build a depiction and categorisation of the wide-ranging body of credibility evidence for CEP models. This paper also represents a step toward the extension of model evaluation methodologies that are currently being developed by the medical device community, to physiological models.

Termination of re-entrant atrial tachycardia via optogenetic stimulation with optimized spatial targeting: insights from computational models

KEY POINTS

Optogenetics has emerged as a potential alternative to electrotherapy for treating heart rhythm disorders, but its applicability for terminating atrial arrhythmias remains largely unexplored. We used computational models reconstructed from clinical MRI scans of fibrotic patient atria to explore the feasibility of optogenetic termination of atrial tachycardia (AT), comparing two different illumination strategies: distributed vs. targeted. We show that targeted optogenetic stimulation based on automated, non-invasive flow-network analysis of patient-specific re-entry morphology may be a reliable approach for identifying the optimal illumination target in each individual (i.e. the critical AT isthmus). The above-described approach yields very high success rates (up to 100%) and requires dramatically less input power than distributed illumination We conclude that simulations in patient-specific models show that targeted light pulses lasting longer than the AT cycle length can efficiently and reliably terminate AT if the human atria can be successfully light-sensitized via gene delivery of ChR2.

ABSTRACT

Optogenetics has emerged as a potential alternative to electrotherapy for treating arrhythmia, but feasibility studies have been limited to ventricular defibrillation via epicardial light application. Here, we assess the efficacy of optogenetic atrial tachycardia (AT) termination in human hearts using a strategy that targets for illumination specific regions identified in an automated manner. In three patient-specific models reconstructed from late gadolinium-enhanced MRI scans, we simulated channelrhodopsin-2 (ChR2) expression via gene delivery. In all three models, we attempted to terminate re-entrant AT (induced via rapid pacing) via optogenetic stimulation. We compared two strategies: (1) distributed illumination of the endocardium by multi-optrode grids (number of optrodes, N_opt  = 64, 128, 256) and (2) targeted illumination of the critical isthmus, which was identified via analysis of simulated activation patterns using an algorithm based on flow networks. The illuminated area and input power were smaller for the targeted approach (19-57.8 mm² ; 0.6-1.8 W) compared to the sparsest distributed arrays (N_opt  = 64; 124.9 ± 6.3 mm² ; 3.9 ± 0.2 W). AT termination rates for distributed illumination were low, ranging from <5% for short pulses (1/10 ms long) to ∼20% for longer stimuli (100/1000 ms). When we attempted to terminate the same AT episodes with targeted illumination, outcomes were similar for short pulses (1/10 ms long: 0% success) but improved for longer stimuli (100 ms: 54% success; 1000 ms: 90% success). We conclude that simulations in patient-specific models show that light pulses lasting longer than the AT cycle length can efficiently and reliably terminate AT in atria light-sensitized via gene delivery. We show that targeted optogenetic stimulation based on analysis of AT morphology may be a reliable approach for defibrillation and requires less power than distributed illumination.

Simultaneous Quantification of Spatially Discordant Alternans in Voltage and Intracellular Calcium in Langendorff-Perfused Rabbit Hearts and Inconsistencies with Models of Cardiac Action Potentials and Ca Transients

Rationale: Discordant alternans, a phenomenon in which the action potential duration (APDs) and/or intracellular calcium transient durations (CaDs) in different spatial regions of cardiac tissue are out of phase, present a dynamical instability for complex spatial dispersion that can be associated with long-QT syndrome (LQTS) and the initiation of reentrant arrhythmias. Because the use of numerical simulations to investigate arrhythmic effects, such as acquired LQTS by drugs is beginning to be studied by the FDA, it is crucial to validate mathematical models that may be used during this process. Objective: In this study, we characterized with high spatio-temporal resolution the development of discordant alternans patterns in transmembrane voltage (Vm) and intracellular calcium concentration ([Cai]+2) as a function of pacing period in rabbit hearts. Then we compared the dynamics to that of the latest state-of-the-art model for ventricular action potentials and calcium transients to better understand the underlying mechanisms of discordant alternans and compared the experimental data to the mathematical models representing Vm and [Cai]+2 dynamics. Methods and Results: We performed simultaneous dual optical mapping imaging of Vm and [Cai]+2 in Langendorff-perfused rabbit hearts with higher spatial resolutions compared with previous studies. The rabbit hearts developed discordant alternans through decreased pacing period protocols and we quantified the presence of multiple nodal points along the direction of wave propagation, both in APD and CaD, and compared these findings with results from theoretical models. In experiments, the nodal lines of CaD alternans have a steeper slope than those of APD alternans, but not as steep as predicted by numerical simulations in rabbit models. We further quantified several additional discrepancies between models and experiments. Conclusions: Alternans in CaD have nodal lines that are about an order of magnitude steeper compared to those of APD alternans. Current action potential models lack the necessary coupling between voltage and calcium compared to experiments and fail to reproduce some key dynamics such as, voltage amplitude alternans, smooth development of calcium alternans in time, conduction velocity and the steepness of the nodal lines of APD and CaD.

Effects of mechanical feedback on the stability of cardiac scroll waves: A bidomain electro-mechanical simulation study

In this work, we investigate the influence of cardiac tissue deformation on re-entrant wave dynamics. We have developed a 3D strongly coupled electro-mechanical Bidomain model posed on an ideal monoventricular geometry, including fiber direction anisotropy and stretch-activated currents (SACs). The cardiac mechanical deformation influences the bioelectrical activity with two main mechanical feedback: (a) the geometric feedback (GEF) due to the presence of the deformation gradient in the diffusion coefficients and in a convective term depending on the deformation rate and (b) the mechano-electric feedback (MEF) due to SACs. Here, we investigate the relative contribution of these two factors with respect to scroll wave stability. We extend the previous works [Keldermann et al., Am. J. Physiol. Heart Circ. Physiol. 299, H134-H143 (2010) and Hu et al., PLoS One 8(4), e60287 (2013)] that were based on the Monodomain model and a simple non-selective linear SAC, while here we consider the full Bidomain model and both selective and non-selective components of SACs. Our simulation results show that the stability of cardiac scroll waves is influenced by MEF, which in case of low reversal potential of non-selective SACs might be responsible for the onset of ventricular fibrillation; GEF increases the scroll wave meandering but does not determine the scroll wave stability.

Highly trabeculated structure of the human endocardium underlies asymmetrical response to low-energy monophasic shocks

Novel low-energy defibrillation therapies are thought to be driven by virtual-electrodes (VEs), due to the interaction of applied monophasic electric shocks with fine-scale anatomical structures within the heart. Significant inter-species differences in the cardiac (micro)-anatomy exist, however, particularly with respect to the degree of endocardial trabeculations, which may underlie important differences in response to low-energy defibrillation protocols. Understanding the interaction of monophasic electric fields with the specific human micro-anatomy is therefore imperative in facilitating the translation and optimisation of these promising experimental therapies to the clinic. In this study, we sought to investigate how electric fields from implanted devices interact with the highly trabeculated human endocardial surface to better understand shock success in order to help optimise future clinical protocols. A bi-ventricular human computational model was constructed from high resolution (350 μm) ex-vivo MR data, including anatomically accurate endocardial structures. Monophasic shocks were applied between a basal right ventricular catheter and an exterior ground. Shocks of varying strengths were applied with both anodal [positive right ventricle (RV) electrode] and cathodal (negative RV electrode) polarities at different states of tissue refractoriness and during induced arrhythmias. Anodal shocks induced isolated positive VEs at the distal side of "detached" trabeculations, which rapidly spread into hyperpolarised tissue on the surrounding endocardial surfaces following the shock. Anodal shocks thus depolarised more tissue 10 ms after the shock than cathodal shocks where the propagation of activation from VEs induced on the proximal side of "detached" trabeculations was prevented due to refractory endocardium. Anodal shocks increased arrhythmia complexity more than cathodal shocks during failed anti-arrhythmia shocks. In conclusion, multiple detached trabeculations in the human ventricle interact with anodal stimuli to induce multiple secondary sources from VEs, facilitating more rapid shock-induced ventricular excitation compared to cathodal shocks. Such a mechanism may help explain inter-species differences in response to shocks and help to develop novel defibrillation strategies.

Determining the light scattering and absorption parameters from forward-directed flux measurements in cardiac tissue

We describe a method to accurately measure the light scattering model parameters from forward-directed flux (FDF) measurements carried out with a fiber-optic probe (optrode). Improved determination of light scattering parameters will, in turn, permit better modeling and interpretation of optical mapping in the heart using voltage-sensitive dyes. Using our optrode-based system, we carried out high spatial resolution measurements of FDF in intact and homogenized cardiac tissue, as well as in intralipid-based tissue phantoms. The samples were illuminated with a broad collimated beam at 660 and 532 nm. Measurements were performed with a plunge fiber-optic probe (NA=0.22) at a spatial resolution of up to 10  μm. In the vicinity of the illuminated surface, the FDF consistently manifested a fast decaying exponent with a space constant comparable with the decay rate of ballistic photons. Using a Monte Carlo model, we obtained a simple empirical formula linking the rate of the fast exponent to the scattering coefficient, the anisotropy parameter g, and the numerical aperture of the probe. The estimates of scattering coefficient based on this formula were validated in tissue phantoms. Potential applications of optical fiber-based FDF measurements for the evaluation of optical parameters in turbid media are discussed.

Bidomain Predictions of Virtual Electrode-Induced Make and Break Excitations around Blood Vessels

Introduction and background

Virtual electrodes formed by field stimulation during defibrillation of cardiac tissue play an important role in eliciting activations. It has been suggested that the coronary vasculature is an important source of virtual electrodes, especially during low-energy defibrillation. This work aims to further the understanding of how virtual electrodes from the coronary vasculature influence defibrillation outcomes.

Methods

Using the bidomain model, we investigated how field stimulation elicited activations from virtual electrodes around idealized intramural blood vessels. Strength–interval curves, which quantify the stimulus strength required to elicit wavefront propagation from the vessels at different states of tissue refractoriness, were computed for each idealized geometry.

Results

Make excitations occurred at late diastolic intervals, originating from regions of depolarization around the vessel. Break excitations occurred at early diastolic intervals, whereby the vessels were able to excite surrounding refractory tissue due to the local restoration of excitability by virtual electrode-induced hyperpolarizations. Overall, strength–interval curves had similar morphologies and underlying excitation mechanisms compared with previous experimental and numerical unipolar stimulation studies of cardiac tissue. Including the presence of the vessel wall increased the field strength required for make excitations but decreased the field strength required for break excitations, and the field strength at which break excitations occurred was generally greater than 5 V/cm. Finally, in a more realistic ventricular slice geometry, the proximity of virtual electrodes around subepicardial vessels was seen to cause break excitations in the form of propagating unstable wavelets to the subepicardial layer.

Conclusion

Representing the blood vessel wall microstructure in computational bidomain models of defibrillation is recommended as it significantly alters the electrophysiological response of the vessel to field stimulation. Although vessels may facilitate excitation of relatively refractory tissue via break excitations, the field strength required for this is generally greater than those used in the literature on low-energy defibrillation. However, the high-intensity shocks used in standard defibrillation may elicit break excitation propagation from the coronary vasculature.

Virtual electrodes around anatomical structures and their roles in defibrillation

BACKGROUND

Virtual electrodes from structural/conductivity heterogeneities are known to elicit wavefront propagation, upon field-stimulation, and are thought to be important for defibrillation. In this work we investigate how the constitutive and geometrical parameters associated with such anatomical heterogeneities, represented by endo/epicardial surfaces and intramural surfaces in the form of blood-vessels, affect the virtual electrode patterns produced.

METHODS AND RESULTS

The steady-state bidomain model is used to obtain, using analytical and numerical methods, the virtual electrode patterns created around idealized endocardial trabeculations and blood-vessels. The virtual electrode pattern around blood-vessels is shown to be composed of two dominant effects; current traversing the vessel surface and conductivity heterogeneity from the fibre-architecture. The relative magnitudes of these two effects explain the swapping of the virtual electrode polarity observed, as a function of the vessel radius, and aid in the understanding of the virtual electrode patterns predicted by numerical bidomain modelling. The relatively high conductivity of blood, compared to myocardium, is shown to cause stronger depolarizations in the endocardial trabeculae grooves than the protrusions.

CONCLUSIONS

The results provide additional quantitative understanding of the virtual electrodes produced by small-scale ventricular anatomy, and highlight the importance of faithfully representing the physiology and the physics in the context of computational modelling of field stimulation.

Optogenetic defibrillation terminates ventricular arrhythmia in mouse hearts and human simulations

Ventricular arrhythmias are among the most severe complications of heart disease and can result in sudden cardiac death. Patients at risk currently receive implantable defibrillators that deliver electrical shocks to terminate arrhythmias on demand. However, strong electrical shocks can damage the heart and cause severe pain. Therefore, we have tested optogenetic defibrillation using expression of the light-sensitive channel channelrhodopsin-2 (ChR2) in cardiac tissue. Epicardial illumination effectively terminated ventricular arrhythmias in hearts from transgenic mice and from WT mice after adeno-associated virus-based gene transfer of ChR2. We also explored optogenetic defibrillation for human hearts, taking advantage of a recently developed, clinically validated in silico approach for simulating infarct-related ventricular tachycardia (VT). Our analysis revealed that illumination with red light effectively terminates VT in diseased, ChR2-expressing human hearts. Mechanistically, we determined that the observed VT termination is due to ChR2-mediated transmural depolarization of the myocardium, which causes a block of voltage-dependent Na+ channels throughout the myocardial wall and interrupts wavefront propagation into illuminated tissue. Thus, our results demonstrate that optogenetic defibrillation is highly effective in the mouse heart and could potentially be translated into humans to achieve nondamaging and pain-free termination of ventricular arrhythmia.

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.

Arrhythmia risk stratification of patients after myocardial infarction using personalized heart models

Sudden cardiac death (SCD) from arrhythmias is a leading cause of mortality. For patients at high SCD risk, prophylactic insertion of implantable cardioverter defibrillators (ICDs) reduces mortality. Current approaches to identify patients at risk for arrhythmia are, however, of low sensitivity and specificity, which results in a low rate of appropriate ICD therapy. Here, we develop a personalized approach to assess SCD risk in post-infarction patients based on cardiac imaging and computational modelling. We construct personalized three-dimensional computer models of post-infarction hearts from patients' clinical magnetic resonance imaging data and assess the propensity of each model to develop arrhythmia. In a proof-of-concept retrospective study, the virtual heart test significantly outperformed several existing clinical metrics in predicting future arrhythmic events. The robust and non-invasive personalized virtual heart risk assessment may have the potential to prevent SCD and avoid unnecessary ICD implantations.

Arrhythmia risk stratification of patients after myocardial infarction using personalized heart models

Sudden cardiac death (SCD) from arrhythmias is a leading cause of mortality. For patients at high SCD risk, prophylactic insertion of implantable cardioverter defibrillators (ICDs) reduces mortality. Current approaches to identify patients at risk for arrhythmia are, however, of low sensitivity and specificity, which results in a low rate of appropriate ICD therapy. Here, we develop a personalized approach to assess SCD risk in post-infarction patients based on cardiac imaging and computational modelling. We construct personalized three-dimensional computer models of post-infarction hearts from patients' clinical magnetic resonance imaging data and assess the propensity of each model to develop arrhythmia. In a proof-of-concept retrospective study, the virtual heart test significantly outperformed several existing clinical metrics in predicting future arrhythmic events. The robust and non-invasive personalized virtual heart risk assessment may have the potential to prevent SCD and avoid unnecessary ICD implantations.

Image-based reconstruction of three-dimensional myocardial infarct geometry for patient-specific modeling of cardiac electrophysiology

PURPOSE

Accurate three-dimensional (3D) reconstruction of myocardial infarct geometry is crucial to patient-specific modeling of the heart aimed at providing therapeutic guidance in ischemic cardiomyopathy. However, myocardial infarct imaging is clinically performed using two-dimensional (2D) late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) techniques, and a method to build accurate 3D infarct reconstructions from the 2D LGE-CMR images has been lacking. The purpose of this study was to address this need.

METHODS

The authors developed a novel methodology to reconstruct 3D infarct geometry from segmented low-resolution (Lo-res) clinical LGE-CMR images. Their methodology employed the so-called logarithm of odds (LogOdds) function to implicitly represent the shape of the infarct in segmented image slices as LogOdds maps. These 2D maps were then interpolated into a 3D image, and the result transformed via the inverse of LogOdds to a binary image representing the 3D infarct geometry. To assess the efficacy of this method, the authors utilized 39 high-resolution (Hi-res) LGE-CMR images, including 36 in vivo acquisitions of human subjects with prior myocardial infarction and 3 ex vivo scans of canine hearts following coronary ligation to induce infarction. The infarct was manually segmented by trained experts in each slice of the Hi-res images, and the segmented data were downsampled to typical clinical resolution. The proposed method was then used to reconstruct 3D infarct geometry from the downsampled images, and the resulting reconstructions were compared with the manually segmented data. The method was extensively evaluated using metrics based on geometry as well as results of electrophysiological simulations of cardiac sinus rhythm and ventricular tachycardia in individual hearts. Several alternative reconstruction techniques were also implemented and compared with the proposed method.

RESULTS

The accuracy of the LogOdds method in reconstructing 3D infarct geometry, as measured by the Dice similarity coefficient, was 82.10% ± 6.58%, a significantly higher value than those of the alternative reconstruction methods. Among outcomes of electrophysiological simulations with infarct reconstructions generated by various methods, the simulation results corresponding to the LogOdds method showed the smallest deviation from those corresponding to the manual reconstructions, as measured by metrics based on both activation maps and pseudo-ECGs.

CONCLUSIONS

The authors have developed a novel method for reconstructing 3D infarct geometry from segmented slices of Lo-res clinical 2D LGE-CMR images. This method outperformed alternative approaches in reproducing expert manual 3D reconstructions and in electrophysiological simulations.

Filament Dynamics during Simulated Ventricular Fibrillation in a High-Resolution Rabbit Heart

The mechanisms underlying ventricular fibrillation (VF) are not well understood. The electrical activity on the heart surface during VF has been recorded extensively in the experimental setting and in some cases clinically; however, corresponding transmural activation patterns are prohibitively difficult to measure. In this paper, we use a high-resolution biventricular heart model to study three-dimensional electrical activity during fibrillation, focusing on the driving sources of VF: “filaments,” the organising centres of unstable reentrant scroll waves. We show, for the first time, specific 3D filament dynamics during simulated VF in a whole heart geometry that includes fine-scale anatomical structures. Our results suggest that transmural activity is much more complex than what would be expected from surface observations alone. We present examples of complex intramural activity, including filament breakup and reattachment, anchoring to the thin right ventricular apex; rapid transitions among various filament shapes; and filament lengths much greater than wall thickness. We also present evidence for anatomy playing a major role in VF development and coronary vessels and trabeculae influencing filament dynamics. Overall, our results indicate that intramural activity during simulated VF is extraordinarily complex and suggest that further investigation of 3D filaments is necessary to fully comprehend recorded surface patterns.

Evaluation of excitation propagation in the rabbit heart: optical mapping and transmural microelectrode recordings

Background

Because of the optical features of heart tissue, optical and electrical action potentials are only moderately associated, especially when near-infrared dyes are used in optical mapping (OM) studies.

Objective

By simultaneously recording transmural electrical action potentials (APs) and optical action potentials (OAPs), we aimed to evaluate the contributions of both electrical and optical influences to the shape of the OAP upstroke.

Methods and Results

A standard glass microelectrode and OM, using an near-infrared fluorescent dye (di-4-ANBDQBS), were used to simultaneously record transmural APs and OAPs in a Langendorff-perfused rabbit heart during atrial, endocardial, and epicardial pacing. The actual profile of the transmural AP upstroke across the LV wall, together with the OAP upstroke, allowed for calculations of the probing-depth constant (k ~2.1 mm, n = 24) of the fluorescence measurements. In addition, the transmural AP recordings aided the quantitative evaluation of the influences of depth-weighted and lateral-scattering components on the OAP upstroke. These components correspond to the components of the propagating electrical wave that are transmural and parallel to the epicardium. The calculated mean values for the depth-weighted and lateral-scattering components, whose sum comprises the OAP upstroke, were (in ms) 10.18 ± 0.62 and 0.0 ± 0.56 for atrial stimulation, 9.37 ± 1.12 and 3.01 ± 1.30 for endocardial stimulation, and 6.09 ± 0.79 and 8.16 ± 0.98 for epicardial stimulation; (n = 8 for each). For this dye, 90% of the collected fluorescence originated up to 4.83 ± 0.18 mm (n = 24) from the epicardium.

Conclusions

The co-registration of OM and transmural microelectrode APs enabled the probing depth of fluorescence measurements to be calculated and the OAP upstroke to be divided into two components (depth-weighted and lateral-scattering), and it also allowed the relative strengths of their effects on the shape of the OAP upstroke to be evaluated.

Simulating photon scattering effects in structurally detailed ventricular models using a Monte Carlo approach

Light scattering during optical imaging of electrical activation within the heart is known to significantly distort the optically-recorded action potential (AP) upstroke, as well as affecting the magnitude of the measured response of ventricular tissue to strong electric shocks. Modeling approaches based on the photon diffusion equation have recently been instrumental in quantifying and helping to understand the origin of the resulting distortion. However, they are unable to faithfully represent regions of non-scattering media, such as small cavities within the myocardium which are filled with perfusate during experiments. Stochastic Monte Carlo (MC) approaches allow simulation and tracking of individual photon “packets” as they propagate through tissue with differing scattering properties. Here, we present a novel application of the MC method of photon scattering simulation, applied for the first time to the simulation of cardiac optical mapping signals within unstructured, tetrahedral, finite element computational ventricular models. The method faithfully allows simulation of optical signals over highly-detailed, anatomically-complex MR-based models, including representations of fine-scale anatomy and intramural cavities. We show that optical action potential upstroke is prolonged close to large subepicardial vessels than further away from vessels, at times having a distinct “humped” morphology. Furthermore, we uncover a novel mechanism by which photon scattering effects around vessels cavities interact with “virtual-electrode” regions of strong de-/hyper-polarized tissue surrounding cavities during shocks, significantly reducing the apparent optically-measured epicardial polarization. We therefore demonstrate the importance of this novel optical mapping simulation approach along with highly anatomically-detailed models to fully investigate electrophysiological phenomena driven by fine-scale structural heterogeneity.

Methodology for image-based reconstruction of ventricular geometry for patient-specific modeling of cardiac electrophysiology

Patient-specific modeling of ventricular electrophysiology requires an interpolated reconstruction of the 3-dimensional (3D) geometry of the patient ventricles from the low-resolution (Lo-res) clinical images. The goal of this study was to implement a processing pipeline for obtaining the interpolated reconstruction, and thoroughly evaluate the efficacy of this pipeline in comparison with alternative methods. The pipeline implemented here involves contouring the epi- and endocardial boundaries in Lo-res images, interpolating the contours using the variational implicit functions method, and merging the interpolation results to obtain the ventricular reconstruction. Five alternative interpolation methods, namely linear, cubic spline, spherical harmonics, cylindrical harmonics, and shape-based interpolation were implemented for comparison. In the thorough evaluation of the processing pipeline, Hi-res magnetic resonance (MR), computed tomography (CT), and diffusion tensor (DT) MR images from numerous hearts were used. Reconstructions obtained from the Hi-res images were compared with the reconstructions computed by each of the interpolation methods from a sparse sample of the Hi-res contours, which mimicked Lo-res clinical images. Qualitative and quantitative comparison of these ventricular geometry reconstructions showed that the variational implicit functions approach performed better than others. Additionally, the outcomes of electrophysiological simulations (sinus rhythm activation maps and pseudo-ECGs) conducted using models based on the various reconstructions were compared. These electrophysiological simulations demonstrated that our implementation of the variational implicit functions-based method had the best accuracy.

Images as drivers of progress in cardiac computational modelling

Computational models have become a fundamental tool in cardiac research. Models are evolving to cover multiple scales and physical mechanisms. They are moving towards mechanistic descriptions of personalised structure and function, including effects of natural variability. These developments are underpinned to a large extent by advances in imaging technologies. This article reviews how novel imaging technologies, or the innovative use and extension of established ones, integrate with computational models and drive novel insights into cardiac biophysics. In terms of structural characterization, we discuss how imaging is allowing a wide range of scales to be considered, from cellular levels to whole organs. We analyse how the evolution from structural to functional imaging is opening new avenues for computational models, and in this respect we review methods for measurement of electrical activity, mechanics and flow. Finally, we consider ways in which combined imaging and modelling research is likely to continue advancing cardiac research, and identify some of the main challenges that remain to be solved.

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.

Advances in modeling ventricular arrhythmias: from mechanisms to the clinic

Modern cardiovascular research has increasingly recognized that heart models and simulation can help interpret an array of experimental data and dissect important mechanisms and interrelationships, with developments rooted in the iterative interaction between modeling and experimentation. This article reviews the progress made in simulating cardiac electrical behavior at the level of the organ and, specifically, in the development of models of ventricular arrhythmias and fibrillation, as well as their termination (defibrillation). The ability to construct multiscale models of ventricular arrhythmias, representing integrative behavior from the molecule to the entire organ, has enabled mechanistic inquiry into the dynamics of ventricular arrhythmias in the diseased myocardium, in understanding drug-induced proarrhythmia, and in the development of new modalities for defibrillation, to name a few. In this article, we also review the initial use of ventricular models of arrhythmia in personalized diagnosis, treatment planning, and prevention of sudden cardiac death. Implementing individualized cardiac simulations at the patient bedside is poised to become one of the most thrilling examples of computational science and engineering approaches in translational medicine.

Diastolic field stimulation: the role of shock duration in epicardial activation and propagation

Detailed knowledge of tissue response to both systolic and diastolic shock is critical for understanding defibrillation. Diastolic field stimulation has been much less studied than systolic stimulation, particularly regarding transient virtual anodes. Here we investigated high-voltage-induced polarization and activation patterns in response to strong diastolic shocks of various durations and of both polarities, and tested the hypothesis that the activation versus shock duration curve contains a local minimum for moderate shock durations, and it grows for short and long durations. We found that 0.1-0.2-ms shocks produced slow and heterogeneous activation. During 0.8-1 ms shocks, the activation was very fast and homogeneous. Further shock extension to 8 ms delayed activation from 1.55 ± 0.27 ms and 1.63 ± 0.21 ms at 0.8 ms shock to 2.32 ± 0.41 ms and 2.37 ± 0.3 ms (N = 7) for normal and opposite polarities, respectively. The traces from hyperpolarized regions during 3-8 ms shocks exhibited four different phases: beginning negative polarization, fast depolarization, slow depolarization, and after-shock increase in upstroke velocity. Thus, the shocks of >3 ms in duration created strong hyperpolarization associated with significant delay (P < 0.05) in activation compared with moderate shocks of 0.8 and 1 ms. This effect appears as a dip in the activation-versus-shock-duration curve.

Defibrillation success with high frequency electric fields is related to degree and location of conduction block

BACKGROUND

We recently demonstrated that high frequency alternating current (HFAC) electric fields can reversibly block propagation in the heart by inducing an oscillating, elevated transmembrane potential (Vm) that maintains myocytes in a refractory state for the field duration and can terminate arrhythmias, including ventricular fibrillation (VF).

OBJECTIVES

To quantify and characterize conduction block (CB) induced by HFAC fields and to determine whether the degree of CB can be used to predict defibrillation success.

METHODS

Optical mapping was performed in adult guinea pig hearts (n = 14), and simulations were performed in an anatomically accurate rabbit ventricular model. HFAC fields (50-500 Hz) were applied to the ventricles. A novel power spectrum metric of CB-the loss of spectral power in the 1-30 Hz range, termed loss of conduction power (LCP)-was assessed during the HFAC field and compared with defibrillation success and VF vulnerability.

RESULTS

LCP increased with field strength and decreased with frequency. Optical mapping experiments conducted on the epicardial surface showed that LCP and the size of CB regions were significantly correlated with VF initiation and termination. In simulations, subsurface myocardial LCP and CB sizes were more closely correlated with VF termination than surface values. Multilinear regression analysis of simulation results revealed that while CB on both the surface and the subsurface myocardium was predictive, subsurface myocardial CB was the better predictor of defibrillation success.

CONCLUSIONS

HFAC fields induce a field-dependent state of CB, and defibrillation success is related to the degree and location of the CB.

Modeling defibrillation of the heart: approaches and insights

Cardiac defibrillation, as accomplished nowadays by automatic, implantable devices (ICDs), constitutes the most important means of combating sudden cardiac death. While ICD therapy has proved to be efficient and reliable, defibrillation is a traumatic experience. Thus, research on defibrillation mechanisms, particularly aimed at lowering defibrillation voltage, remains an important topic. Advancing our understanding towards a full appreciation of the mechanisms by which a shock interacts with the heart is the most promising approach to achieve this goal. The aim of this paper is to assess the current state-of-the-art in ventricular defibrillation modeling, focusing on both numerical modeling approaches and major insights that have been obtained using defibrillation models, primarily those of realistic ventricular geometry. The paper showcases the contributions that modeling and simulation have made to our understanding of the defibrillation process. The review thus provides an example of biophysically based computational modeling 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.

Mapping cardiac surface mechanics with structured light imaging

Cardiovascular disease often manifests as a combination of pathological electrical and structural heart remodeling. The relationship between mechanics and electrophysiology is crucial to our understanding of mechanisms of cardiac arrhythmias and the treatment of cardiac disease. While several technologies exist for describing whole heart electrophysiology, studies of cardiac mechanics are often limited to rhythmic patterns or small sections of tissue. Here, we present a comprehensive system based on ultrafast three-dimensional (3-D) structured light imaging to map surface dynamics of whole heart cardiac motion. Additionally, we introduce a novel nonrigid motion-tracking algorithm based on an isometry-maximizing optimization framework that forms correspondences between consecutive 3-D frames without the use of any fiducial markers. By combining our 3-D imaging system with nonrigid surface registration, we are able to measure cardiac surface mechanics at unprecedented spatial and temporal resolution. In conclusion, we demonstrate accurate cardiac deformation at over 200,000 surface points of a rabbit heart recorded at 200 frames/s and validate our results on highly contrasting heart motions during normal sinus rhythm, ventricular pacing, and ventricular fibrillation.

Three-dimensional mechanisms of increased vulnerability to electric shocks in myocardial infarction: altered virtual electrode polarizations and conduction delay in the peri-infarct zone

Defibrillation efficacy is decreased in infarcted hearts, but the mechanisms by which infarcted hearts are more vulnerable to electric shocks than healthy hearts remain poorly understood. The goal of this study was to provide insight into the 3D mechanisms for the increased vulnerability to electric shocks in infarcted hearts. We hypothesized that changes in virtual electrode polarizations (VEPs) and propagation delay through the peri-infarct zone (PZ) were responsible. We developed a micro anatomically detailed rabbit ventricular model with chronic myocardial infarction from magnetic resonance imaging and enriched the model with data from optical mapping experiments. We further developed a control model without the infarct. The simulation protocol involved apical pacing followed by biphasic shocks. Simulation results from both models were compared.The upper limit of vulnerability(ULV) was 8 V cm(-1) in the infarction model and 4 V cm(-1) in the control model. VEPs were less pronounced in the infarction model, providing a larger excitable area for postshock propagation but smaller transmembrane potential gradients to initiate new wavefronts. Initial post-shock transmural activation occurred at a later time in the infarction model, and the PZ served to delay propagation in subsequent beats. The presence of the PZ was found to be responsible for the increased vulnerability.

Bridging experiments, models and simulations: an integrative approach to validation in computational cardiac electrophysiology

Computational models in physiology often integrate functional and structural information from a large range of spatiotemporal scales from the ionic to the whole organ level. Their sophistication raises both expectations and skepticism concerning how computational methods can improve our understanding of living organisms and also how they can reduce, replace, and refine animal experiments. A fundamental requirement to fulfill these expectations and achieve the full potential of computational physiology is a clear understanding of what models represent and how they can be validated. The present study aims at informing strategies for validation by elucidating the complex interrelations among experiments, models, and simulations in cardiac electrophysiology. We describe the processes, data, and knowledge involved in the construction of whole ventricular multiscale models of cardiac electrophysiology. Our analysis reveals that models, simulations, and experiments are intertwined, in an assemblage that is a system itself, namely the model-simulation-experiment (MSE) system. We argue that validation is part of the whole MSE system and is contingent upon 1) understanding and coping with sources of biovariability; 2) testing and developing robust techniques and tools as a prerequisite to conducting physiological investigations; 3) defining and adopting standards to facilitate the interoperability of experiments, models, and simulations; 4) and understanding physiological validation as an iterative process that contributes to defining the specific aspects of cardiac electrophysiology the MSE system targets, rather than being only an external test, and that this is driven by advances in experimental and computational methods and the combination of both.

Optical imaging of voltage and calcium in cardiac cells & tissues

Cardiac optical mapping has proven to be a powerful technology for studying cardiovascular function and disease. The development and scientific impact of this methodology are well-documented. Because of its relevance in cardiac research, this imaging technology advances at a rapid pace. Here, we review technological and scientific developments during the past several years and look toward the future. First, we explore key components of a modern optical mapping set-up, focusing on: (1) new camera technologies; (2) powerful light-emitting-diodes (from ultraviolet to red) for illumination; (3) improved optical filter technology; (4) new synthetic and optogenetic fluorescent probes; (5) optical mapping with motion and contraction; (6) new multiparametric optical mapping techniques; and (7) photon scattering effects in thick tissue preparations. We then look at recent optical mapping studies in single cells, cardiomyocyte monolayers, atria, and whole hearts. Finally, we briefly look into the possible future roles of optical mapping in the development of regenerative cardiac research, cardiac cell therapies, and molecular genetic advances.

At the heart of computational modelling

The link between experimental data and biophysically based mathematical models is key to computational simulation meeting its potential to provide physiological insight. However, despite the importance of this link, scrutiny and analysis of the processes by which models are parameterised from data are currently lacking. While this situation is common to many areas of physiological modelling, to provide a concrete context, we use examples drawn from detailed models of cardiac electro-mechanics. Using this biophysically detailed cohort of models we highlight the specific issues of model parameterization and propose this process can be separated into three stages: observation, fitting and validation. Finally, future research challenges and directions in this area are discussed.

Investigating the role of the coronary vasculature in the mechanisms of defibrillation

BACKGROUND

The direct role of coronary vessels in defibrillation, although hypothesized to be important, remains to be elucidated. We investigated how vessel-induced virtual electrode polarizations assist reentry termination.

METHODS AND RESULTS

A highly anatomically detailed rabbit ventricular slice bidomain computer model was constructed from 25-μm magnetic resonance data, faithfully representing both structural and electric properties of blood vessels. For comparison, an equivalent simplified model with intramural cavities filled in was also built. Following fibrillation induction, 6 initial states were selected, and biphasic shocks (5-70 V) were applied using a realistic implanted cardioverter-defibrillator electrode configuration. A fundamental mechanism of biphasic defibrillation was uncovered in both models, involving successive break excitations (after each shock phase) emanating from opposing myocardial surfaces (in septum and left ventricle), which rapidly closed down excitable gaps. The presence of vessels accelerated this process, achieving more-rapid and successful defibrillation. Defibrillation failed in 5 cases (all because of initiation of new activity) compared with 8 with the simplified model (5/8 failures because of surviving activity). At stronger shocks, virtual electrodes formed around vessels, rapidly activating intramural tissue because of break excitations, assisting the main defibrillation mechanism, and eliminating all activity <15 ms of shock end in 60% of successful shocks (36% in simplified model). Subsequent analysis identified that only vessels >200 μm in diameter participated through this mechanism. Consequently, wavefronts could survive intramurally in the simplified model, leading to reentry and shock failure.

CONCLUSIONS

We provide new insight into defibrillation mechanisms by showing how intramural blood vessels facilitate more-effective elimination of existing wavefronts, rapid closing down of excitable gaps, and successful defibrillation and give guidance toward the required resolution of cardiac imaging and model generation endeavors for mechanistic defibrillation analysis.

The role of photon scattering in voltage-calcium fluorescent recordings of ventricular fibrillation

Recent optical mapping studies of cardiac tissue suggest that membrane voltage (V(m)) and intracellular calcium concentrations (Ca) become dissociated during ventricular fibrillation (VF), generating a proarrhythmic substrate. However, experimental methods used in these studies may accentuate measured dissociation due to differences in fluorescent emission wavelengths of optical voltage/calcium (V(opt)/Ca(opt)) signals. Here, we simulate dual voltage-calcium optical mapping experiments using a monodomain-Luo-Rudy ventricular-tissue model coupled to a photon-diffusion model. Dissociation of both electrical, V(m)/Ca, and optical, V(opt)/Ca(opt), signals is quantified by calculating mutual information (MI) for VF and rapid pacing protocols. We find that photon scattering decreases MI of V(opt)/Ca(opt) signals by 23% compared to unscattered V(m)/Ca signals during VF. Scattering effects are amplified by increasing wavelength separation between fluorescent voltage/calcium signals and respective measurement-location misalignment. In contrast, photon scattering does not affect MI during rapid pacing, but high calcium dye affinity can decrease MI by attenuating alternans in Ca(opt) but not in V(opt). We conclude that some dissociation exists between voltage and calcium at the cellular level during VF, but MI differences are amplified by current optical mapping methods.

Multiparametric optical mapping of the Langendorff-perfused rabbit heart

Optical imaging and fluorescent probes have significantly advanced research methodology in the field of cardiac electrophysiology in ways that could not have been accomplished by other approaches¹. With the use of the calcium- and voltage-sensitive dyes, optical mapping allows measurement of transmembrane action potentials and calcium transients with high spatial resolution without the physical contact with the tissue. This makes measurements of the cardiac electrical activity possible under many conditions where the use of electrodes is inconvenient or impossible¹. For example, optical recordings provide accurate morphological changes of membrane potential during and immediately after stimulation and defibrillation, while conventional electrode techniques suffer from stimulus-induced artifacts during and after stimuli due to electrode polarization¹.

The Langendorff-perfused rabbit heart is one of the most studied models of human heart physiology and pathophysiology. Many types of arrhythmias observed clinically could be recapitulated in the rabbit heart model. It was shown that wave patterns in the rabbit heart during ventricular arrhythmias, determined by effective size of the heart and the wavelength of reentry, are very similar to that in the human heart². It was also shown that critical aspects of excitation-contraction (EC) coupling in rabbit myocardium, such as the relative contribution of sarcoplasmic reticulum (SR), is very similar to human EC coupling³. Here we present the basic procedures of optical mapping experiments in Langendorff-perfused rabbit hearts, including the Langendorff perfusion system setup, the optical mapping systems setup, the isolation and cannulation of the heart, perfusion and dye-staining of the heart, excitation-contraction uncoupling, and collection of optical signals. These methods could be also applied to the heart from species other than rabbit with adjustments to flow rates, optics, solutions, etc.

Two optical mapping systems are described. The panoramic mapping system is used to map the entire epicardium of the rabbit heart^4-7. This system provides a global view of the evolution of reentrant circuits during arrhythmogenesis and defibrillation, and has been used to study the mechanisms of arrhythmias and antiarrhythmia therapy^8,9. The dual mapping system is used to map the action potential (AP) and calcium transient (CaT) simultaneously from the same field of view^10-13. This approach has enhanced our understanding of the important role of calcium in the electrical alternans and the induction of arrhythmia^14-16.

Single-sensor system for spatially resolved, continuous, and multiparametric optical mapping of cardiac tissue

Background

Simultaneous optical mapping of multiple electrophysiologically relevant parameters in living myocardium is desirable for integrative exploration of mechanisms underlying heart rhythm generation under normal and pathophysiologic conditions. Current multiparametric methods are technically challenging, usually involving multiple sensors and moving parts, which contributes to high logistic and economic thresholds that prevent easy application of the technique.

Objective

The purpose of this study was to develop a simple, affordable, and effective method for spatially resolved, continuous, simultaneous, and multiparametric optical mapping of the heart, using a single camera.

Methods

We present a new method to simultaneously monitor multiple parameters using inexpensive off-the-shelf electronic components and no moving parts. The system comprises a single camera, commercially available optical filters, and light-emitting diodes (LEDs), integrated via microcontroller-based electronics for frame-accurate illumination of the tissue. For proof of principle, we illustrate measurement of four parameters, suitable for ratiometric mapping of membrane potential (di-4-ANBDQPQ) and intracellular free calcium (fura-2), in an isolated Langendorff-perfused rat heart during sinus rhythm and ectopy, induced by local electrical or mechanical stimulation.

Results

The pilot application demonstrates suitability of this imaging approach for heart rhythm research in the isolated heart. In addition, locally induced excitation, whether stimulated electrically or mechanically, gives rise to similar ventricular propagation patterns.

Conclusion

Combining an affordable camera with suitable optical filters and microprocessor-controlled LEDs, single-sensor multiparametric optical mapping can be practically implemented in a simple yet powerful configuration and applied to heart rhythm research. The moderate system complexity and component cost is destined to lower the threshold to broader application of functional imaging and to ease implementation of more complex optical mapping approaches, such as multiparametric panoramic imaging. A proof-of-principle application confirmed that although electrically and mechanically induced excitation occur by different mechanisms, their electrophysiologic consequences downstream from the point of activation are not dissimilar.

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.

Whole-heart modeling: applications to cardiac electrophysiology and electromechanics

Recent developments in cardiac simulation have rendered the heart the most highly integrated example of a virtual organ. We are on the brink of a revolution in cardiac research, one in which computational modeling of proteins, cells, tissues, and the organ permit linking genomic and proteomic information to the integrated organ behavior, in the quest for a quantitative understanding of the functioning of the heart in health and disease. The goal of this review is to assess the existing state-of-the-art in whole-heart modeling and the plethora of its applications in cardiac research. General whole-heart modeling approaches are presented, and the applications of whole-heart models in cardiac electrophysiology and electromechanics research are reviewed. The article showcases the contributions that whole-heart modeling and simulation have made to our understanding of the functioning of the heart. A summary of the future developments envisioned for the field of cardiac simulation and modeling is also presented. Biophysically based computational modeling of the heart, applied to human heart physiology and the diagnosis and treatment of cardiac disease, has the potential to dramatically change 21st century cardiac research and the field of cardiology.

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.

Dual excitation wavelength epifluorescence imaging of transmural electrophysiological properties in intact hearts

BACKGROUND

Epifluorescence imaging using voltage-sensitive dyes has provided unique insights into cardiac electrical activity and arrhythmias. However, conventional dyes use blue-green excitation light, which has limited depth penetration.

OBJECTIVE

The aim of this study was to demonstrate that combining a short and a long excitation wavelength using near-infrared (NIR) dyes allows for epifluorescence imaging of transmural electrophysiological properties in intact hearts.

METHODS

Epifluorescence imaging was performed in rat hearts (N = 11) using DI-4-ANEPPS and the NIR dye DI-4-ANBDQBS. Activation and action potential duration (APD) patterns were investigated at 2 excitation wavelengths (530 and 660 nm) after epicardial stimulation at various cycle lengths (160 to 70 ms).

RESULTS

Optical action potential upstrokes acquired with 660-nm excitation of DI-4-ANBDQBS were significantly longer than upstrokes obtained with 530-nm excitation of DI-4-ANEPPS (P < .001). Comparison of activation maps showed counterclockwise rotation of isochrones consistent with a transmural rotation of myofibers. Pronounced APD modulation by the activation sequence was observed at both excitation wavelengths. Significantly prolonged APDs (P = .016) and steeper APD restitution curves were found with DI-4-ANBDQBS (660-nm excitation) when compared with DI-4-ANEPPS (530-nm excitation). Dual excitation wavelength experiments using solely DI-4-ANBDQBS yielded similar results. Monophasic action potential recordings showed prolonged APD and steeper APD restitution curves in the endocardium, indicating that 660-nm excitation provides a significant endocardial contribution to the signal. Three-dimensional computer simulations confirmed our findings.

CONCLUSION

Dual excitation wavelength epifluorescence allows detecting transmural heterogeneity in intact hearts. It therefore has the potential to become an important tool in experimental cardiac electrophysiology.