Search

Artificial intelligence (AI) and machine learning (ML) in medicine are currently areas of intense exploration, showing potential to automate human tasks and even perform tasks beyond human capabilities. Literacy and understanding of AI/ML methods are becoming increasingly important to researchers and clinicians. The first objective of this review is to provide the novice reader with literacy of AI/ML methods and provide a foundation for how one might conduct an ML study. We provide a technical overview of some of the most commonly used terms, techniques, and challenges in AI/ML studies, with reference to recent studies in cardiac electrophysiology to illustrate key points. The second objective of this review is to use examples from recent literature to discuss how AI and ML are changing clinical practice and research in cardiac electrophysiology, with emphasis on disease detection and diagnosis, prediction of patient outcomes, and novel characterization of disease. The final objective is to highlight important considerations and challenges for appropriate validation, adoption, and deployment of AI technologies into clinical practice.

Background: To facilitate ablation of ventricular tachycardia (VT), an automated localization system to identify the site of origin of left ventricular activation in real time using the 12-lead ECG was developed. The objective of this study was to prospectively assess its accuracy., Methods: The automated site of origin localization system consists of 3 steps: (1) localization of ventricular segment based on population templates, (2) population-based localization within a segment, and (3) patient-specific site localization. Localization error was assessed by the distance between the known reference site and the estimated site., Results: In 19 patients undergoing 21 catheter ablation procedures of scar-related VT, site of origin localization accuracy was estimated using 552 left ventricular endocardial pacing sites pooled together and 25 VT-exit sites identified by contact mapping. For the 25 VT-exit sites, localization error of the population-based localization steps was within 10 mm. Patient-specific site localization achieved accuracy of within 3.5 mm after including up to 11 pacing (training) sites. Using 3 remotes (67.8±17.0 mm from the reference VT-exit site), and then 5 close pacing sites, resulted in localization error of 7.2±4.1 mm for the 25 identified VT-exit sites. In 2 emulated clinical procedure with 2 induced VTs, the site of origin localization system achieved accuracy within 4 mm., Conclusions: In this prospective validation study, the automated localization system achieved estimated accuracy within 10 mm and could thus provide clinical utility.

Patients with myocardial infarction have an abundance of conduction channels (CC); however, only a small subset of these CCs sustain ventricular tachycardia (VT). Identifying these critical CCs (CCCs) in the clinic so that they can be targeted by ablation remains a significant challenge. The objective of this study is to use a personalized virtual-heart approach to conduct a three-dimensional (3D) assessment of CCCs sustaining VTs of different morphologies in these patients, to investigate their 3D structural features, and to determine the optimal ablation strategy for each VT. To achieve these goals, ventricular models were constructed from contrast enhanced magnetic resonance imagings of six postinfarction patients. Rapid pacing induced VTs in each model. CCCs that sustained different VT morphologies were identified. CCCs' 3D structure and type and the resulting rotational electrical activity were examined. Ablation was performed at the optimal part of each CCC, aiming to terminate each VT with a minimal lesion size. Predicted ablation locations were compared to clinical. Analyzing the simulation results, we found that the observed VTs in each patient model were sustained by a limited number (2.7 ± 1.2) of CCCs. Further, we identified three types of CCCs sustaining VTs: I-type and T-type channels, with all channel branches bounded by scar, and functional reentry channels, which were fully or partially bounded by conduction block surfaces. The different types of CCCs accounted for 43.8, 18.8, and 37.4% of all CCCs, respectively. The mean narrowest width of CCCs or a branch of CCC was 9.7 ± 3.6 mm. Ablation of the narrowest part of each CCC was sufficient to terminate VT. Our results demonstrate that a personalized virtual-heart approach can determine the possible VT morphologies in each patient and identify the CCCs that sustain reentry. The approach can aid clinicians in identifying accurately the optimal VT ablation targets in postinfarction patients., (Copyright © 2019 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Background: Mechanisms of ventricular tachycardia (VT) and ventricular fibrillation (VF) in patients with heart failure (HF) are undefined., Objective: The purpose of this study was to elucidate VT/VF mechanisms in HF by using a computational-clinical approach., Methods: In 53 patients with HF and 18 control patients, we established the relationship between low-amplitude action potential voltage alternans (APV-ALT) during ventricular pacing at near-resting heart rates and VT/VF on long-term follow-up. Mechanisms underlying the transition of APV-ALT to VT/VF, which cannot be ascertained in patients, were dissected with multiscale human ventricular models based on human electrophysiological and magnetic resonance imaging data (control and HF)., Results: For patients with APV-ALT k-score >1.7, complex action potential duration (APD) oscillations (≥2.3% of mean APD), rather than APD alternans, most accurately predicted VT/VF during long-term follow-up (+82%; -90% predictive values). In the failing human ventricular models, abnormal sarcoplasmic reticulum (SR) calcium handling caused APV-ALT (>1 mV) during pacing with a cycle length of 550 ms, which transitioned into large magnitude (>100 ms) discordant repolarization time alternans (RT-ALT) at faster rates. This initiated VT/VF (cycle length <400 ms) by steepening apicobasal repolarization (189 ms/mm) until unidirectional conduction block and reentry. Complex APD oscillations resulted from nonstationary discordant RT-ALT. Restoring SR calcium to control levels was antiarrhythmic by terminating electrical alternans., Conclusion: APV-ALT and complex APD oscillations at near-resting heart rates in patients with HF are linked to arrhythmogenic discordant RT-ALT. This may enable novel physiologically tailored, bioengineered indices to improve VT/VF risk stratification, where SR calcium handling and spatial apicobasal repolarization are potential therapeutic targets., (Copyright © 2016 Heart Rhythm Society. Published by Elsevier Inc. All rights reserved.)

Background: Left atrial flutter (LAFL) occurs in patients after atrial fibrillation ablation. Identification of optimal ablation targets to terminate LAFL remains challenging., Objective: The purpose of this study was to use patient-specific models to simulate LAFL and predict optimal ablation targets using a novel approach based on flow network theory., Methods: Late gadolinium-enhanced cardiac magnetic resonance scans from 10 patients with LAFL were used to construct atrial models incorporating fibrosis by investigators blinded to procedural findings. Rapid pacing was applied in silico to induce LAFL. In each LAFL, we represented reentrant wave propagation as an electric flow network and identified the "minimum cut" (MC), which was the smallest amount of tissue that separated the flow into 2 discontinuous components. In silico ablation was applied at MCs, and targets were compared to those that terminated LAFL during catheter ablation., Results: Patient-specific atrial models were successfully generated from patient scans. LAFL was induced in 7 of 10 models. Ablation of MCs terminated LAFL in 4 models and produced new, slower LAFL morphologies in the other 3. For the latter cases, flow analysis was repeated to identify MCs of emergent LAFLs. Ablation of these MCs terminated emergent LAFLs. The MC-based ablation lesions in simulations were similar in length and location to ablation targets that terminated LAFL during catheter ablation for these 7 patients., Conclusion: Personalized atrial simulations can predict ablation targets for LAFL. These simulations provide a powerful tool for planning ablation procedures and may reduce procedural times and complications., (Copyright © 2016 Heart Rhythm Society. Published by Elsevier Inc. All rights reserved.)

Background: The extent of left atrial (LA) late gadolinium enhancement (LGE), as a surrogate for fibrosis, has been associated with atrial fibrillation (AF) recurrence after catheter ablation. Furthermore, there is ex vivo evidence that islands of fibrosis may anchor fibrillatory rotors., Objective: The purpose of this study was to examine the anatomical association of AF rotors with LA and right atrial (RA) LGE on cardiac magnetic resonance., Methods: The cohort included 9 patients with persistent AF (mean age 61.1 ± 9.7 years) who underwent LGE cardiac magnetic resonance before AF ablation using the focal impulse and rotor modulation system. The extent of LA and RA LGE was quantified globally and in each of the 7 sectors: LA posterior/inferior wall, anterior wall, roof, left and right pulmonary vein antra, and RA lateral and septal regions. The multivariable association of rotor incidence with global and per sector LGE extent was examined using multivariable Bernoulli logistic regression estimated by generalized estimating equations., Results: The mean RA and LA volumes were 113.2 ± 37.31 and 143.03 ± 58.25 mL, respectively. The mean RA and LA LGE burden was 17.2% ± 11.0% and 17.4% ± 14.4%, respectively. A total of 18 LA rotors and 9 RA rotors were identified in all patients. No univariable or multivariable association was observed between global or per sector LGE extent and focal impulse and rotor modulation rotor incidence., Conclusion: In this cohort of patients, there was no association between AF rotor incidence and the global or regional extent of RA and LA LGE., (Copyright © 2016 Heart Rhythm Society. Published by Elsevier Inc. All rights reserved.)

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.

Atrial fibrillation (AF) is the most common sustained arrhythmia in humans. The mechanisms that govern AF initiation and persistence are highly complex, of dynamic nature, and involve interactions across multiple temporal and spatial scales in the atria. This article aims to review the mathematical modeling and computer simulation approaches to understanding AF mechanisms and aiding in its management. Various atrial modeling approaches are presented, with descriptions of the methodological basis and advancements in both lower-dimensional and realistic geometry models. A review of the most significant mechanistic insights made by atrial simulations is provided. The article showcases the contributions that atrial modeling and simulation have made not only to our understanding of the pathophysiology of atrial arrhythmias, but also to the development of AF management approaches. A summary of the future developments envisioned for the field of atrial simulation and modeling is also presented. The review contends that computational models of the atria assembled with data from clinical imaging modalities that incorporate electrophysiological and structural remodeling could become a first line of screening for new AF therapies and approaches, new diagnostic developments, and new methods for arrhythmia prevention.

In addition to the left bundle branch block type of electrical activation, there are further remodeling aspects associated with dyssynchronous heart failure (HF) that affect the electromechanical behavior of the heart. Among the most important are altered ventricular structure (both geometry and fiber/sheet orientation), abnormal Ca(2+) handling, slowed conduction, and reduced wall stiffness. In dyssynchronous HF, the electromechanical delay (EMD), the time interval between local myocyte depolarization and myofiber shortening onset, is prolonged. However, the contributions of the four major HF remodeling aspects in extending EMD in the dyssynchronous failing heart remain unknown. The goal of this study was to determine the individual and combined contributions of HF-induced remodeling aspects to EMD prolongation. We used MRI-based models of dyssynchronous nonfailing and HF canine electromechanics and constructed additional models in which varying combinations of the four remodeling aspects were represented. A left bundle branch block electrical activation sequence was simulated in all models. The simulation results revealed that deranged Ca(2+) handling is the primary culprit in extending EMD in dyssynchronous HF, with the other aspects of remodeling contributing insignificantly. Mechanistically, we found that abnormal Ca(2+) handling in dyssynchronous HF slows myofiber shortening velocity at the early-activated septum and depresses both myofiber shortening and stretch rate at the late-activated lateral wall. These changes in myofiber dynamics delay the onset of myofiber shortening, thus giving rise to prolonged EMD in dyssynchronous HF.

Background: Experimental studies have demonstrated that unstable repolarization dynamics is a risk factor of arrhythmia. We have recently developed an algorithm to detect QT interval (QTI) instability from the clinical electrocardiogram (ECG)., Objective: To develop a clinical arrhythmia risk stratification index based on the detection of QTI instability., Methods: Intracardiac ECGs were recorded at rest in 114 patients with implanted implantable cardioverter-defibrillators (ICDs). Patients were followed up until appropriate implantable cardioverter-defibrillator therapy or death occurred, whichever came first. Each ECG recording was divided into 1-minute episodes (minECGs); the instability in QTI dynamics, if any, of each minECG was detected with our algorithm. An arrhythmia risk index termed QTI instability index (QTII) was defined as the number of minECGs with unstable QTI dynamics normalized by the number of minECGs with premature activations. The performance of QTII in arrhythmia risk stratification was examined with survival analysis and was compared with other risk indices, such as the mean RR interval (RRI), the standard deviation of the RRI and the QTI, and the frequency of premature activation. We hypothesized that the index QTII, which accounts for multiple risk factors and their interdependence, performs better than indices quantifying individual arrhythmia risk factors in the stratification of arrhythmia risk., Results: The results of survival analysis show that QTII outperformed all other studied indices in arrhythmia risk stratification and was the only independent indicator of arrhythmia propensity in a multivariate survival model., Conclusion: QTII is a promising arrhythmia risk stratification index., (Copyright © 2013 Heart Rhythm Society. Published by Elsevier Inc. All rights reserved.)

Background: Mechanisms of atrial fibrillation (AF) initiation are incompletely understood. We hypothesized that rate-dependent changes (restitution) in action potential duration (APD) and activation latency are central targets for clinical interventions that induce AF. We tested this hypothesis using clinical experiments and computer models., Methods and Results: In 50 patients (20 persistent, 23 paroxysmal AF, 7 controls), we used monophasic action potential catheters to define left atrial APD restitution, activation latency, and AF incidence from premature extrastimuli. Isoproterenol (n=14), adenosine (n=10), or rapid pacing (n=36) was then initiated to determine impact on these parameters. Compared with baseline in AF patients, isoproterenol and rapid pacing decreased activation latency (64±14 versus 31±13 versus 24±14 ms; P<0.05), steepened maximum APD restitution slope (0.8±0.7 versus 1.7±0.5 versus 1.1±0.5; P<0.05), and increased AF incidence (12% versus 64% versus 84%; P<0.05). Conversely, adenosine shortened APD (P<0.05), yet increased activation latency (86±27 ms; P=0.002) so that maximum APD restitution slope did not steepen (1.0±0.5; P=NS), and AF incidence was unchanged (10%; P=NS). In controls, no intervention steepened APD restitution or initiated AF. Computational modeling revealed that isoproterenol steepened APD restitution by increased L-type calcium current and decreased activation latency via enhanced rapid delayed potassium reactifier current inactivation, whereas rapid pacing steepened APD restitution via increased cardiac inward potassium rectifier current., Conclusions: Steep APD restitution is a common pathway for AF initiation by isoproterenol and tachycardia via reduced activation latency that enables engagement of steep APD restitution at rapid rates. Modeling suggests that AF initiation from each intervention uses distinct ionic mechanisms. This insight may help design interventions to prevent AF.

Personalized computational cardiac models are emerging as an important tool for studying cardiac arrhythmia mechanisms, and have the potential to become powerful instruments for guiding clinical anti-arrhythmia therapy. In this article, we present the methodology for constructing a patient-specific model of atrial fibrosis as a substrate for atrial fibrillation. The model is constructed from high-resolution late gadolinium-enhanced magnetic resonance imaging (LGE-MRI) images acquired in vivo from a patient suffering from persistent atrial fibrillation, accurately capturing both the patient's atrial geometry and the distribution of the fibrotic regions in the atria. Atrial fiber orientation is estimated using a novel image-based method, and fibrosis is represented in the patient-specific fibrotic regions as incorporating collagenous septa, gap junction remodeling, and myofibroblast proliferation. A proof-of-concept simulation result of reentrant circuits underlying atrial fibrillation in the model of the patient's fibrotic atrium is presented to demonstrate the completion of methodology development., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Background: Instability in ventricular repolarization in the presence of premature activations (PA) plays an important role in arrhythmogenesis. However, such instability cannot be detected clinically. This study developed a methodology for detecting QT interval (QTI) dynamics instability from the ECG and explored the contribution of PA and QTI instability to ventricular tachycardia (VT) onset., Methods and Results: To examine the contribution of PAs and QTI instability to VT onset, ECGs of 24 patients with acute myocardial infarction, 12 of whom had sustained VT (VT) and 12 nonsustained VT (NSVT), were used. From each patient ECG, 2 10-minute-long ECG recordings were extracted, 1 right before VT onset (onset epoch) and 1 at least 1 hour before it (control epoch). To ascertain how PA affects QTI dynamics stability, pseudo-ECGs were calculated from an MRI-based human ventricular model. Clinical and pseudo-ECGs were subdivided into 1-minute recordings (minECGs). QTI dynamics stability of each minECG was assessed with a novel approach. Frequency of PAs (f(PA)) and the number of minECGs with unstable QTI dynamics (N(us)) were determined for each patient. In the VT group, f(PA) and N(us) of the onset epoch were larger than in control. Positive regression relationships between f(PA) and N(us) were identified in both groups. The simulations showed that both f(PA) and the PA degree of prematurity contribute to QTI dynamics instability., Conclusions: Increased PA frequency and QTI dynamics instability precede VT onset in patients with acute myocardial infarction, as determined by novel methodology for detecting instability in QTI dynamics from clinical ECGs.

Electrical impulse propagation is an essential function in cardiac, skeletal muscle, and nervous tissue. Abnormalities in cardiac impulse propagation underlie lethal reentrant arrhythmias, including ventricular fibrillation. Temporary propagation block throughout the ventricular myocardium could possibly terminate these arrhythmias. Electrical stimulation has been applied to nervous tissue to cause reversible conduction block, but has not been explored sufficiently in cardiac tissue. We show that reversible propagation block can be achieved in cardiac tissue by holding myocardial cells in a refractory state for a designated period of time by applying a sustained sinusoidal high-frequency alternating current (HFAC); in doing so, reentrant arrhythmias are terminated. We demonstrate proof of concept using several models, including optically mapped monolayers of neonatal rat ventricular cardiomyocytes, Langendorff-perfused guinea pig and rabbit hearts, intact anesthetized adult rabbits, and computer simulations of whole-heart impulse propagation. HFAC may be an effective and potentially safer alternative to direct current application, currently used to treat ventricular fibrillation.

One of the most important components of mechanoelectric coupling is stretch-activated channels, sarcolemmal channels that open upon mechanical stimuli. Uncovering the mechanisms by which stretch-activated channels contribute to ventricular arrhythmogenesis under a variety of pathologic conditions is hampered by the lack of experimental methodologies that can record the 3-dimensional electromechanical activity simultaneously at high spatiotemporal resolution. Computer modeling provides such an opportunity. The goal of this review is to illustrate the utility of sophisticated, physiologically realistic, whole heart computer simulations in determining the role of mechanoelectric coupling in ventricular arrhythmogenesis. We first present the various ways by which stretch-activated channels have been modeled and demonstrate how these channels affect cardiac electrophysiologic properties. Next, we use an electrophysiologic model of the rabbit ventricles to understand how so-called commotio cordis, the mechanical impact to the precordial region of the heart, can initiate ventricular tachycardia via the recruitment of stretch-activated channels. Using the same model, we also provide mechanistic insight into the termination of arrhythmias by precordial thump under normal and globally ischemic conditions. Lastly, we use a novel anatomically realistic dynamic 3-dimensional coupled electromechanical model of the rabbit ventricles to gain insight into the role of electromechanical dysfunction in arrhythmogenesis during acute regional ischemia., (Copyright © 2010 Elsevier Inc. All rights reserved.)

Background: After near-defibrillation threshold (DFT) shocks from an implantable cardioverter-defibrillator (ICD), the first postshock activation that leads to defibrillation failure arises focally after an isoelectric window (IW). The mechanisms underlying the IW remain incompletely understood., Objective: The goal of this study was to provide mechanistic insight into the origins of postshock activations and IW after ICD shocks, and to link shock outcome to the preshock state of the ventricles. We hypothesized that the nonuniform ICD field results in the formation of an intramural excitable area (tunnel) only in the left ventricular (LV) free wall, through which both pre-existing and new shock-induced wavefronts propagate during the IW., Methods: Simulations were conducted using a realistic three dimensional (3D) model of defibrillation in the rabbit ventricles. Biphasic ICD shocks of varying strengths were delivered to 27 different fibrillatory states., Results: After near-DFT shocks, regardless of preshock state, the main postshock excitable area was always located within LV free wall, creating an intramural tunnel. Either pre-existing fibrillatory or shock-induced wavefronts propagated during the IW (duration of up to 74 ms) in this tunnel and emerged as breakthroughs on LV epicardium. Preshock activity within the LV played a significant role in shock outcome: a large number of preshock filaments resulted in an IW associated with tunnel propagation of pre-existing rather than shock-induced wavefronts. Furthermore, shocks were more likely to succeed if the LV excitable area was smaller., Conclusion: The LV intramural excitable area is the primary reason for near-DFT failure. Any intervention that decreases the extent of this area will improve the likelihood of defibrillation success., (Copyright 2010 Heart Rhythm Society. Published by Elsevier Inc. All rights reserved.)

Background: During phase 1B of acute regional ischemia, the subepicardial and subendocardial layers coupled to the inexcitable midmyocardium remain viable., Objective: The purpose of this study was to examine how the degree of hyperkalemia in the surviving layers, the lateral width of border zone between the normal tissue and the central ischemic zone, and the degree of cellular uncoupling between the surviving layers and the midmyocardium contribute to initiation of reentry., Methods: Simulations were conducted on the state-of-the-art model of rabbit ventricles with realistic representation of the spatial distribution of the ischemic insult., Results: Hyperkalemia in the surviving layers led to induction of reentry by increasing refractoriness and slowing conduction in the layers. Such reentries were formed solely in the subepicardium. A minimal level of hyperkalemia was required for induction of reentry. Progress increase in hyperkalemia led to a biphasic change in vulnerability to reentry. For each level of hyperkalemia, increased cellular uncoupling between subepicardium and midmyocardium increased inducibility of reentry by restoring subepicardial tissue excitability via blocking midmyocardial electrotonic effect. In addition, increased lateral width of the border zone prevented inducibility of reentry as conduction block occurred in the central ischemic zone when the wave propagated across the border zone from the normal zone., Conclusion: The degree of hyperkalemia in the surviving subepicardium, the lateral width of border zone, and cellular uncoupling between the subepicardium and midmyocardium determine dispersion of refractoriness, conduction velocity, excitability, and, therefore, inducibility of reentry during phase 1B., (Copyright 2010 Heart Rhythm Society. Published by Elsevier Inc. All rights reserved.)

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.

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.

Experimental and clinical studies have shown that ventricular dilatation is associated with increased arrhythmogenesis and elevated defibrillation threshold; however, the underlying mechanisms remain poorly understood. The goal of the present study was to test the hypothesis that (1) stretch-activated channel (SAC) recruitment and (2) geometrical deformations in organ shape and fiber architecture lead to increased arrhythmogenesis by electric shocks following acute ventricular dilatation. To elucidate the contribution of these two factors, the study employed, for the first time, a combined electro-mechanical simulation approach. Acute dilatation was simulated in a model of rabbit ventricular mechanics by raising the LV end-diastolic pressure from 0.6 (control) to 4.2 kPa (dilated). The output of the mechanics model was used in the electrophysiological model. Vulnerability to shocks was examined in the control, the dilated ventricles, and in the dilated ventricles that also incorporated currents through SAC as a function of local strain, by constructing vulnerability grids. Results showed that dilatation-induced deformation alone decreased upper limit of vulnerability (ULV) slightly and did not result in increased vulnerability. With SAC recruitment in the dilated ventricles, the number of shock-induced arrhythmia episodes increased by 37% (from 41 to 56) and the lower limit of vulnerability (LLV) decreased from 9 to 7 V/cm, while ULV did not change. The heterogeneous activation of SAC caused by the heterogeneous fiber strain in the ventricular walls was the main reason for increased vulnerability to electric shocks since it caused dispersion of electrophysiological properties in the tissue, resulting in postshock unidirectional block and establishment of reentry.

Comprehensive understanding of the ventricular response to shocks is the approach most likely to succeed in reducing defibrillation threshold. We propose a new theory of shock-induced arrhythmogenesis that unifies all known aspects of the response of the heart to monophasic (MS) and biphasic (BS) shocks. The central hypothesis is that submerged "tunnel" propagation of postshock activations through shock-induced intramural excitable areas underlies fibrillation induction and the existence of isoelectric window. We conducted simulations of fibrillation induction using a realistic bidomain model of rabbit ventricles. Following pacing, MS and BS of various strengths/timings were delivered. The results demonstrated that, during the isoelectric window, an activation originated deep within the ventricular wall, arising from virtual electrodes; it then propagated fully intramurally through an excitable tunnel induced by the shock, until it emerged onto the epicardium, becoming the earliest-propagated postshock activation. Differences in shock outcomes for MS and BS were found to stem from the narrower BS intramural postshock excitable area, often resulting in conduction block, and the difference in the mechanisms of origin of the postshock activations, namely intramural virtual electrode-induced phase singularity for MS and virtual electrode-induced propagated graded response for BS. This study provides a novel analysis of the 3D mechanisms underlying the origin of postshock activations in the process of fibrillation induction by MS and BS and the existence of isoelectric window. The tunnel propagation hypothesis could open a new avenue for interventions exploration to achieve significantly lower defibrillation threshold.

Background: Optical mapping is a widely used experimental tool providing high-resolution recordings of cardiac electrical activity. However, the technique is limited by signal distortion due to photon scattering in the tissue. Computational models of the fluorescence recording are capable of assessing these distortion effects, providing important insight to assist experimental data interpretation., Methods: We present results from a new panoramic optical mapping model, which is used to assess distortion in ventricular optical mapping signals during pacing and arrhythmogenesis arising from 3-dimensional photon scattering., Results/conclusions: We demonstrate that accurate consideration of wavefront propagation within the complex ventricular structure, along with accurate representation of photon scattering in 3 dimensions, is essential to faithfully assess distortion effects arising during optical mapping. In this article, examined effects include (1) the specific morphology of the optical action potential upstroke during pacing and (2) the shift in the location of epicardial phase singularities obtained from fluorescent maps.

The goal of this study is to investigate the mechanisms responsible for the increase in the upper limit of vulnerability (ULV; highest shock strength that induces arrhythmia) following the increase in pacing rate. To accomplish this goal, the study employs a three-dimensional bidomain finite element model of a slice through the canine ventricles. The preparation was paced eight times at a basic cycle length (BCL) of either 80 or 150ms followed by delivery of shocks of various strengths and timings. Our results demonstrate that the shock strength, which induced an arrhythmia 50% of the time, increased 20% for the faster pacing compared to the slower pacing. Analysis of the mechanisms underlying the increased vulnerability revealed that delayed post-shock activations originating in the tissue depths appear as breakthrough activations on the surfaces of the preparation following an isoelectric window (IW). However, the IW duration was consistently shorter in the faster-paced preparation. Consequently, breakthrough activations appeared on the surfaces of this preparation earlier, when the tissue was less recovered, resulting in higher probability of unidirectional block and reentry. This explains why shocks of the same strength were more likely to result in arrhythmia induction when delivered to a preparation that was rapidly paced.

Atrial fibrillation (AF) is believed to be perpetuated by recirculating spiral waves. Atrial structures are often characterized with action potentials of varying morphologies; however, the role of the structure-dependent atrial electrophysiological heterogeneity in spiral wave behaviour is not well understood. The purpose of this study is to determine the effect of action potential morphology heterogeneity associated with the major atrial structures in spiral wave maintenance. The present study also focuses on how this effect is further modulated by the presence of the inherent periodicity in atrial structure. The goals of the study are achieved through the simulation of electrical behaviour in a two-dimensional atrial tissue model that incorporates the representation of action potentials in various structurally distinct regions in the right atrium. Periodic boundary conditions are then imposed to form a cylinder (quasi three-dimensional), thus allowing exploration of the additional effect of structure periodicity on spiral wave behaviour. Transmembrane potential maps and phase singularity traces are analysed to determine effects on spiral wave behaviour. Results demonstrate that the prolonged refractoriness of the crista terminalis (CT) affects the pattern of spiral wave reentry, while the variation in action potential morphology of the other structures does not. The CT anchors the spiral waves, preventing them from drifting away. Spiral wave dynamics is altered when the ends of the sheet are spliced together to form a cylinder. The main effect of the continuous surface is the generation of secondary spiral waves which influences the primary rotors. The interaction of the primary and secondary spiral waves decreased as cylinder diameter increased.

Background: Studies have demonstrated that failed defibrillation shocks often are followed by an electrically quiescent period (isoelectric window); however, the underlying mechanisms remain incompletely understood. We recently suggested a new mechanism termed "virtual electrode polarization-induced propagated graded responses" (VEPiPGRs) that might play a role in the origin of the global postshock activation following the isoelectric window., Objectives: The purpose of this study to elucidate the circumstances under which VEPiPGR activations originate for shocks given to paced right ventricular preparations. Specifically, we examined the dependence of VEPiPGRs on coupling interval (CI) and shock polarity and whether VEPiPGRs emerge preferentially on the epicardium or the endocardium., Methods: Simultaneous endocardial and epicardial activity in isolated right ventricular preparations (n = 4) was imaged optically following shocks of strength +/-5A. All VEPiPGRs were analyzed, and the time T from shock end to activation onset was recorded (isoelectric window is the smallest T among activations that propagated globally)., Results: VEPiPGR activations occurred for CIs in the range from 80 to 150 ms. Average duration of T was 64.5 +/- 18.15 ms, with T decreasing as CI increased (Tmax = 82 ms, Tmin = 46 ms, linear-fit slope = -0.675). The average earliest CI at which cathodal (+5A) shocks resulted in VEPiPGRs was 87 ms compared with 116 ms for anodal (-5A) shocks. All VEPiPGR activations emerged first on the epicardium in a focal pattern, and all induced ventricular fibrillation., Conclusion: The global activation that terminates the isoelectric window could result from VEPiPGRs that find an exit pathway. VEPiPGRs originate at the sites of maximum action potential abbreviation by the shock, always on the epicardium for the preparation used here.

Aim: Existing models of myocardial membrane kinetics have not been able to reproduce the experimentally-observed negative bias in the asymmetry of transmembrane potential changes (DeltaV(m)) induced by strong electric shocks. The goals of this study are (1) to demonstrate that this negative bias could be reproduced by the addition, to the membrane model, of electroporation and an outward current, I(a), part of the K(+) flow through the L-type Ca(2+)-channel, and (2) to determine how such modifications in the membrane model affect shock-induced break excitation in a 2D preparation., Methods and Results: We conducted simulations of shocks in bidomain fibres and sheets with membrane dynamics represented by the Luo-Rudy dynamic model (LRd'2000), to which electroporation (LRd + EP model) and the outward current, I(a), activated upon strong shock-induced depolarization (aLRd model) was added. Assuming I(a) is a part of K(+) flow through the L-type Ca(2+)-channel enabled us to reproduce both the experimentally observed rectangularly-shaped positive DeltaV(m) and the value of near 2 of the negative-to-positive DeltaV(m) ratio. In the sheet, I(a) not only contributed to the negative bias in DeltaV(m) asymmetry at sites polarized by physical and virtual electrodes, but also restricted positive DeltaV(m). Electroporation, in its turn, was responsible for the decrease in cathode-break excitation threshold in the aLRd sheet, compared with the other two cases, as well as for the occurrence of the excitation after the shock-end rather than during the shock., Conclusions: The incorporation of electroporation and I(a) in a membrane model ensures match between simulation results and experimental data. The use of the aLRd model results in a lower threshold for shock-induced break excitation.

Models of myocardial membrane dynamics have not been able to reproduce the experimentally observed negative bias in the asymmetry of transmembrane potential changes (DeltaVm) induced by strong electric shocks delivered during the action potential plateau. The goal of this study is to determine what membrane model modifications can bridge this gap between simulation and experiment. We conducted simulations of shocks in bidomain fibers and sheets with membrane dynamics represented by the LRd'2000 model. We found that in the fiber, the negative bias in DeltaVm asymmetry could not be reproduced by addition of electroporation only, but by further addition of hypothetical outward current, Ia, activated upon strong shock-induced depolarization. Furthermore, the experimentally observed rectangularly shaped positive DeltaVm, negative-to-positive DeltaVm ratio (asymmetry ratio) = approximately 2, electroporation occurring at the anode only, and the increase in positive DeltaVm caused by L-type Ca2+-channel blockade were reproduced in the strand only if Ia was assumed to be a part of K+ flow through the L-type Ca2+-channel. In the sheet, Ia not only contributed to the negative bias in DeltaVm asymmetry at sites polarized by physical and virtual electrodes, but also restricted positive DeltaVm. Inclusion of Ia and electroporation is thus the bridge between experiment and simulation., (Copyright 2004 Biophysical Society)

Introduction: Our previous article examined cardiac vulnerability to reentry in the near field within the framework of the virtual electrode polarization (VEP) concept. The present study extends this examination to the far field and compares its predictions to the critical point hypothesis., Methods and Results: We simulate the electrical behavior of a sheet of myocardium using a two-dimensional bidomain model. The fiber field is extrapolated from a set of rabbit heart fiber directions obtained experimentally. An S1 stimulus is applied along the top or left border. An extracellular line electrode on the top delivers a cathodal or anodal S2 stimulus. A VEP pattern matching that seen experimentally is observed and covers the entire sheet, thus constituting a far-field effect. Reentry arises from break excitation, make excitation, or a combination of both, and subsequent propagation through deexcited and recovered areas. Reentry occurs in cross-field, parallel-field, and uniform refractoriness protocols. For long coupling intervals (CIs) above CImake(min) (defined as the shortest CI at which make excitation can take place), rotors move away from the cathodal electrode and the S1 site for increases in S2 strength and CI, respectively. For cathodal S2 stimuli, findings are consistent with the critical point hypothesis. For CIs below CImake(min), reentry is initiated by break excitation only, and the resulting reentrant patterns are no longer consistent with those predicted by the critical point hypothesis., Conclusion: Shock-induced VEP can explain vulnerability in the far field. The VEP theory of vulnerability encompasses the critical point hypothesis for cathodal S2 shocks at long CIs.

Introduction: Recent experimental evidence demonstrates that a point stimulus generates a nonuniform distribution of transmembrane potential (virtual electrode pattern) consisting of large adjacent areas of depolarization and hyperpolarization. This simulation study focuses on the role of virtual electrodes in reentry induction., Methods and Results: We simulated the electrical behavior of a sheet of myocardium using a two-dimensional bidomain model with straight fibers. Membrane kinetics were represented by the Beeler-Reuter Drouhard-Roberge model. Simulations were conducted for equal and unequal anisotropy ratios. S1 wavefront was planar and propagated parallel or perpendicular to the fibers. S2 unipolar stimulus was cathodal or anodal. With regard to unequal anisotropy, for both cathodal and anodal stimuli, the S2 stimulus negatively polarizes some portion of membrane, deexciting it and opening an excitable pathway in a region of otherwise unexcitable tissue. Reentry is generated by break excitation of this tissue and subsequent propagation through deexcited and recovered areas of myocardium. Figure-of-eight and quatrefoil reentry are observed, with figure-of-eight most common. Figure-of-eight rotation is seen in the direction predicted by the critical point hypothesis. With regard to equal anisotropy, reentry was observed for cathodal stimuli only at strengths > -95 A/m., Conclusion: The key to reentry induction is the close proximity of S2-induced excited and deexcited areas, with adjacent nonexcited areas available for propagation.