Jason D. Bayer

A Computational Model to Predict the Effects of Class I Anti-Arrhythmic Drugs on Ventricular Rhythms

Two- and three-dimensional models of cardiac excitability based on sodium channel kinetics can predict the adverse effects of class I anti-arrhythmic drugs. Crowdsourcing the Heart for Drug Screening The old way: Consult a specialist to answer your question. The new way: Consult a crowd of generalists who in the aggregate can come up with a better answer. The old way—testing drugs on single cardiac cells in vitro—has not worked well for screening out potential anti-arrhythmia agents that can occasionally block conduction in the heart or exacerbate arrhythmia, serious problems that cause sudden death in treated patients. Instead, Moreno et al. have called on the crowd by building a model of heart tissue that includes many cardiac cells and their interactions. When anti-arrhythmia drugs are “applied” to the model’s beating heart tissue—but not when they are applied to the single cardiac cells that make up the model—the drugs that cause side effects, and the concentrations at which they do so, are revealed, results that the authors were able to validate experimentally. The model starts with the detailed kinetics of the heart’s sodium channels, first in the context of a single cell, then in two- and three-dimensional cardiac tissue. The authors compared the action of lidocaine, a class 1B anti-arrhythmic drug not known to cause conduction block, and flecainide, a prototypical class 1C drug that carries a warning from the Food and Drug Administration. In the modeled analyses of single cardiac cells, both drugs slowed excitability at concentrations that matched those used in patients, but the cells retained the ability to generate action potentials. But when the model incorporated coupled groups of cells, the behavior of the drugs diverged. Lidocaine lowered excitability without causing block, but at the higher concentrations (used clinically), flecainide caused serious conduction block when heart rates reached 160 beats per minute. Experiments in rabbit heart confirmed the results of the model. In scaled-up, 500 by 500 groups of cells, the authors’ model could also successfully predict the tendency of flecainide, but not lidocaine, to make the heart extra sensitive to heartbeats occurring too early or too late, an effect that causes even more severe arrhythmias in patients when they take anti-arrhythmia drugs. Again, experiments in rabbit hearts replicated the model’s predictions, as did simulations of anatomically accurate human hearts derived from magnetic resonance imaging images. The ability of this sophisticated model of living cardiac tissue to replicate the clinical adverse effects of lidocaine and flecainide is promising, but it will be necessary to validate its performance with other drugs to understand how to deploy it most effectively. Ideally, such models will be useful for screening out potential arrhythmic drugs that promote conduction block or exacerbate arrhythmias. Such a view of how drugs affect the collective activity of cardiac cells should help in these situations in which the cure proves more deadly than the disease. A long-sought, and thus far elusive, goal has been to develop drugs to manage diseases of excitability. One such disease that affects millions each year is cardiac arrhythmia, which occurs when electrical impulses in the heart become disordered, sometimes causing sudden death. Pharmacological management of cardiac arrhythmia has failed because it is not possible to predict how drugs that target cardiac ion channels, and have intrinsically complex dynamic interactions with ion channels, will alter the emergent electrical behavior generated in the heart. Here, we applied a computational model, which was informed and validated by experimental data, that defined key measurable parameters necessary to simulate the interaction kinetics of the anti-arrhythmic drugs flecainide and lidocaine with cardiac sodium channels. We then used the model to predict the effects of these drugs on normal human ventricular cellular and tissue electrical activity in the setting of a common arrhythmia trigger, spontaneous ventricular ectopy. The model forecasts the clinically relevant concentrations at which flecainide and lidocaine exacerbate, rather than ameliorate, arrhythmia. Experiments in rabbit hearts and simulations in human ventricles based on magnetic resonance images validated the model predictions. This computational framework initiates the first steps toward development of a virtual drug-screening system that models drug-channel interactions and predicts the effects of drugs on emergent electrical activity in the heart.

Action potential dynamics explain arrhythmic vulnerability in human heart failure: a clinical and modeling study implicating abnormal calcium handling.

OBJECTIVES The purpose of this study was to determine whether abnormalities of calcium cycling explain ventricular action potential (AP) oscillations and cause electrocardiogram T-wave alternans (TWA). BACKGROUND Mechanisms explaining why heart failure patients are at risk for malignant ventricular arrhythmias (ventricular tachycardia [VT]/ventricular fibrillation [VF]) are unclear. We studied whether oscillations in human ventricular AP explain TWA and predict VT/VF, and used computer modeling to suggest potential cellular mechanisms. METHODS We studied 53 patients with left ventricular ejection fraction 28 +/- 8% and 18 control subjects. Monophasic APs were recorded in the right ventricle (n = 62) and/or left ventricle (n = 9) at 109 beats/min. RESULTS Alternans of AP amplitude, computed spectrally, had higher magnitude in study patients than in controls (p = 0.03), particularly in AP phase II (p = 0.02) rather than phase III (p = 0.10). The AP duration and activation restitution (n = 11 patients) were flat at 109 beats/min and did not explain TWA. In computer simulations, only reduced sarcoplasmic reticulum calcium uptake explained our results, causing calcium oscillations, AP amplitude alternans, and TWA that were all abolished by calcium clamping. On prospective follow-up for 949 +/- 553 days, 17 patients had VT/VF. The AP amplitude alternans predicted VT/VF (p = 0.04), and was 78% concordant with simultaneous TWA (p = 0.003). CONCLUSIONS Patients with systolic dysfunction show ventricular AP amplitude alternans that prospectively predicted VT/VF. Alternans in AP amplitude, but not variations in AP duration or conduction, explained TWA at < or =109 beats/min. In computer models, these findings were best explained by reduced sarcoplasmic reticulum calcium uptake. Thus, in heart failure patients, in vivo AP alternans may reflect cellular calcium abnormalities and provide a mechanistic link with VT/VF.

Feasibility of image-based simulation to estimate ablation target in human ventricular arrhythmia

BACKGROUND Previous studies suggest that magnetic resonance imaging with late gadolinium enhancement (LGE) may identify slowly conducting tissues in scar-related ventricular tachycardia (VT). OBJECTIVE To test the feasibility of image-based simulation based on LGE to estimate ablation targets in VT. METHODS We conducted a retrospective study in 13 patients who had preablation magnetic resonance imaging for scar-related VT ablation. We used image-based simulation to induce VT and estimate target regions according to the simulated VT circuit. The estimated target regions were coregistered with the LGE scar map and the ablation sites from the electroanatomical map in the standard ablation approach. RESULTS In image-based simulation, VT was inducible in 12 (92.3%) patients. All VTs showed macroreentrant propagation patterns, and the narrowest width of estimated target region that an ablation line should span to prevent VT recurrence was 5.0 ± 3.4 mm. Of 11 patients who underwent ablation, the results of image-based simulation and the standard approach were consistent in 9 (82%) patients, where ablation within the estimated target region was associated with acute success (n = 8) and ablation outside the estimated target region was associated with failure (n = 1). In 1 (9%) case, the results of image-based simulation and the standard approach were inconsistent, where ablation outside the estimated target region was associated with acute success. CONCLUSIONS The image-based simulation can be used to estimate potential ablation targets of scar-related VT. The image-based simulation may be a powerful noninvasive tool for preprocedural planning of ablation procedures to potentially reduce the procedure time and complication rates.

Rate-dependent action potential alternans in human heart failure implicates abnormal intracellular calcium handling

BACKGROUND Alternans in action potential voltage (APV-ALT) at heart rates <110 bpm is a novel index to predict ventricular arrhythmias. However, the rate dependency of APV-ALT and its mechanisms in failing versus nonfailing human myocardium are poorly understood. It is hypothesized that APV-ALT in human heart failure (HF) reflects abnormal calcium handling. OBJECTIVE Using a modeling and clinical approach, our objectives were to (1) determine how APV-ALT varies with pacing rate and (2) ascertain whether abnormalities in calcium handling explain the rate dependence of APV-ALT in HF. METHODS APV-ALT was analyzed at several cycle lengths (CLs) using a dynamic pacing protocol applied to a human left ventricle wedge model with various alterations in calcium handling. Modeled APV-ALT was used to predict APV-ALT in left ventricle monophasic action potentials recorded from HF (n = 3) and control (n = 2) patients with the same pacing protocol. RESULTS Reducing the sarcoplasmic reticulum calcium uptake current < or =25%, the release current < or =11%, or the sarcolemmal L-type calcium channel current < or =43% of control predicted APV-ALT to arise at CL > or =600 ms and then increase in magnitude by >400% for CL <400 ms. In HF patients, APV-ALT arose at CL = 600 ms and then increased in magnitude by >500% at CL <350 ms. For all other model alterations and for control patients, APV-ALT occurred only at CL <500 ms. CONCLUSIONS APV-ALT shows differing rate dependence in HF versus control patients, arising at slower rates in HF and predicted by models with abnormal calcium handling. Future studies should investigate whether APV-ALT at slow rates identifies patients with deranged calcium handing, including HF patients before decompensation or at risk for arrhythmias.

Mechanisms of Human Atrial Fibrillation Initiation Clinical and Computational Studies of Repolarization Restitution and Activation Latency

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.

Computational cardiology: how computer simulations could be used to develop new therapies and advance existing ones.

This article reviews the latest developments in computational cardiology. It focuses on the contribution of cardiac modelling to the development of new therapies as well as the advancement of existing ones for cardiac arrhythmias and pump dysfunction. Reviewed are cardiac modelling efforts aimed at advancing and optimizing existent therapies for cardiac disease (defibrillation, ablation of ventricular tachycardia, and cardiac resynchronization therapy) and at suggesting novel treatments, including novel molecular targets, as well as efforts to use cardiac models in stratification of patients likely to benefit from a given therapy, and the use of models in diagnostic procedures.

Spatial Resolution Requirements for Accurate Identification of Drivers of Atrial Fibrillation

Background— Recent studies have demonstrated conflicting mechanisms underlying atrial fibrillation (AF), with the spatial resolution of data often cited as a potential reason for the disagreement. The purpose of this study was to investigate whether the variation in spatial resolution of mapping may lead to misinterpretation of the underlying mechanism in persistent AF. Methods and Results— Simulations of rotors and focal sources were performed to estimate the minimum number of recording points required to correctly identify the underlying AF mechanism. The effects of different data types (action potentials and unipolar or bipolar electrograms) and rotor stability on resolution requirements were investigated. We also determined the ability of clinically used endocardial catheters to identify AF mechanisms using clinically recorded and simulated data. The spatial resolution required for correct identification of rotors and focal sources is a linear function of spatial wavelength (the distance between wavefronts) of the arrhythmia. Rotor localization errors are larger for electrogram data than for action potential data. Stationary rotors are more reliably identified compared with meandering trajectories, for any given spatial resolution. All clinical high-resolution multipolar catheters are of sufficient resolution to accurately detect and track rotors when placed over the rotor core although the low-resolution basket catheter is prone to false detections and may incorrectly identify rotors that are not present. Conclusions— The spatial resolution of AF data can significantly affect the interpretation of the underlying AF mechanism. Therefore, the interpretation of human AF data must be taken in the context of the spatial resolution of the recordings.

Disrupted Calcium Release as a Mechanism for Atrial Alternans Associated with Human Atrial Fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhythmia, but our knowledge of the arrhythmogenic substrate is incomplete. Alternans, the beat-to-beat alternation in the shape of cardiac electrical signals, typically occurs at fast heart rates and leads to arrhythmia. However, atrial alternans have been observed at slower pacing rates in AF patients than in controls, suggesting that increased vulnerability to arrhythmia in AF patients may be due to the proarrythmic influence of alternans at these slower rates. As such, alternans may present a useful therapeutic target for the treatment and prevention of AF, but the mechanism underlying alternans occurrence in AF patients at heart rates near rest is unknown. The goal of this study was to determine how cellular changes that occur in human AF affect the appearance of alternans at heart rates near rest. To achieve this, we developed a computational model of human atrial tissue incorporating electrophysiological remodeling associated with chronic AF (cAF) and performed parameter sensitivity analysis of ionic model parameters to determine which cellular changes led to alternans. Of the 20 parameters tested, only decreasing the ryanodine receptor (RyR) inactivation rate constant (kiCa) produced action potential duration (APD) alternans seen clinically at slower pacing rates. Using single-cell clamps of voltage, fluxes, and state variables, we determined that alternans onset was Ca2+-driven rather than voltage-driven and occurred as a result of decreased RyR inactivation which led to increased steepness of the sarcoplasmic reticulum (SR) Ca2+ release slope. Iterated map analysis revealed that because SR Ca2+ uptake efficiency was much higher in control atrial cells than in cAF cells, drastic reductions in kiCa were required to produce alternans at comparable pacing rates in control atrial cells. These findings suggest that RyR kinetics may play a critical role in altered Ca2+ homeostasis which drives proarrhythmic APD alternans in patients with AF.

Effects of Mechano-Electric Feedback on Scroll Wave Stability in Human Ventricular Fibrillation

Recruitment of stretch-activated channels, one of the mechanisms of mechano-electric feedback, has been shown to influence the stability of scroll waves, the waves that underlie reentrant arrhythmias. However, a comprehensive study to examine the effects of recruitment of stretch-activated channels with different reversal potentials and conductances on scroll wave stability has not been undertaken; the mechanisms by which stretch-activated channel opening alters scroll wave stability are also not well understood. The goals of this study were to test the hypothesis that recruitment of stretch-activated channels affects scroll wave stability differently depending on stretch-activated channel reversal potential and channel conductance, and to uncover the relevant mechanisms underlying the observed behaviors. We developed a strongly-coupled model of human ventricular electromechanics that incorporated human ventricular geometry and fiber and sheet orientation reconstructed from MR and diffusion tensor MR images. Since a wide variety of reversal potentials and channel conductances have been reported for stretch-activated channels, two reversal potentials, −60 mV and −10 mV, and a range of channel conductances (0 to 0.07 mS/µF) were implemented. Opening of stretch-activated channels with a reversal potential of −60 mV diminished scroll wave breakup for all values of conductances by flattening heterogeneously the action potential duration restitution curve. Opening of stretch-activated channels with a reversal potential of −10 mV inhibited partially scroll wave breakup at low conductance values (from 0.02 to 0.04 mS/µF) by flattening heterogeneously the conduction velocity restitution relation. For large conductance values (>0.05 mS/µF), recruitment of stretch-activated channels with a reversal potential of −10 mV did not reduce the likelihood of scroll wave breakup because Na channel inactivation in regions of large stretch led to conduction block, which counteracted the increased scroll wave stability due to an overall flatter conduction velocity restitution.

Novel Radiofrequency Ablation Strategies for Terminating Atrial Fibrillation in the Left Atrium: A Simulation Study

Pulmonary vein isolation (PVI) with radiofrequency ablation (RFA) is the cornerstone of atrial fibrillation (AF) therapy, but few strategies exist for when it fails. To guide RFA, phase singularity (PS) mapping locates reentrant electrical waves (rotors) that perpetuate AF. The goal of this study was to test existing and develop new RFA strategies for terminating rotors identified with PS mapping. It is unsafe to test experimental RFA strategies in patients, so they were evaluated in silico using a bilayer computer model of the human atria with persistent AF (pAF) electrical (ionic) and structural (fibrosis) remodeling. pAF was initiated by rapidly pacing the right (RSPV) and left (LSPV) superior pulmonary veins during sinus rhythm, and rotor dynamics quantified by PS analysis. Three RFA strategies were studied: (i) PVI, roof, and mitral lines; (ii) circles, perforated circles, lines, and crosses 0.5–1.5 cm in diameter/length administered near rotor locations/pathways identified by PS mapping; and (iii) 4–8 lines streamlining the sequence of electrical activation during sinus rhythm. As in pAF patients, 2 ± 1 rotors with cycle length 185 ± 4 ms and short PS duration 452 ± 401 ms perpetuated simulated pAF. Spatially, PS density had weak to moderate positive correlations with fibrosis density (RSPV: r = 0.38, p = 0.35, LSPV: r = 0.77, p = 0.02). RFA PVI, mitral, and roof lines failed to terminate pAF, but RFA perforated circles and lines 1.5 cm in diameter/length terminated meandering rotors from RSPV pacing when placed at locations with high PS density. Similarly, RFA circles, perforated circles, and crosses 1.5 cm in diameter/length terminated stationary rotors from LSPV pacing. The most effective strategy for terminating pAF was to streamline the sequence of activation during sinus rhythm with >4 RFA lines. These results demonstrate that co-localizing 1.5 cm RFA lesions with locations of high PS density is a promising strategy for terminating pAF rotors. For patients immune to PVI, roof, mitral, and PS guided RFA strategies, streamlining patient-specific activation sequences during sinus rhythm is a robust but challenging alternative.