Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Whole-heart computational modelling provides further mechanistic insights into ST-elevation in Brugada syndrome

Background

Brugada syndrome (BrS) is characterized by dynamic ST-elevations in right precordial leads and increased risk of ventricular fibrillation and sudden cardiac death. As the mechanism underlying ST-elevation and malignant arrhythmias is controversial computational modeling can aid in exploring the disease mechanism. Thus we aim to test the main competing hypotheses ('delayed depolarization' vs. 'early repolarization') of BrS in a whole-heart computational model.

Methods

In a 3D whole-heart computational model, delayed epicardial RVOT activation with local conduction delay was simulated by reducing conductivity in the epicardial RVOT. Early repolarization was simulated by instead increasing the transient outward potassium current (Ito) in the same region. Additionally, a reduction in the fast sodium current (INa) was incorporated in both models.

Results

Delayed depolarization with local conduction delay in the computational model resulted in coved-type ST-elevation with negative T-waves in the precordial surface ECG leads. 'Saddleback'-shaped ST-elevation was obtained with reduced substrate extent or thickness. Increased Ito simulations showed early repolarization in the RVOT with a descending but not coved-type ST-elevation. Reduced INa did not show a significant effect on ECG morphology.

Conclusions

In this whole-heart BrS computational model of both major hypotheses, realistic coved-type ECG resulted only from delayed epicardial RVOT depolarization with local conduction delay but not early repolarizing ion channel modifications. These simulations provide further support for the depolarization hypothesis as electrophysiological mechanism underlying BrS.

Diffuse fibrosis and repolarization disorders explain ventricular arrhythmias in Brugada syndrome: a computational study

In this work, we reported a computational study to quantitatively determine the individual contributions of three candidate arrhythmic factors associated with Brugada Syndrome. In particular, we focused our analysis on the role of structural abnormalities, dispersion of repolarization, and size of the diseased region. We developed a human phenomenological model capable of replicating the action potential characteristics both in Brugada Syndrome and in healthy conditions. Inspired by physiological observations, we employed the phenomenological model in a 2D geometry resembling the pathological RVOT coupled with healthy epicardial tissue. We assessed the insurgence of sustained reentry as a function of electrophysiological and structural abnormalities. Our computational study indicates that both structural and repolarization abnormalities are essential to induce sustained reentry. Furthermore, our results suggest that neither dispersion of repolarization nor structural abnormalities are sufficient on their own to induce sustained reentry. It should be noted how our study seems to explain an arrhythmic mechanism that unifies the classic repolarization and depolarization hypotheses of the pathophysiology of the Brugada Syndrome. Finally, we believe that this work may offer a new perspective on the computational and clinical investigation of Brugada Syndrome and its arrhythmic behaviour.

Why Is Only Type 1 Electrocardiogram Diagnostic of Brugada Syndrome? Mechanistic Insights From Computer Modeling

BACKGROUND

Three types of characteristic ST-segment elevation are associated with Brugada syndrome but only type 1 is diagnostic. Why only type 1 ECG is diagnostic remains unanswered.

METHODS

Computer simulations were performed in single cells, 1-dimensional cables, and 2-dimensional tissues to investigate the effects of the peak and late components of the transient outward potassium current (Ito), sodium current, and L-type calcium current (ICa,L) as well as other potassium currents on the genesis of ECG morphologies and phase 2 reentry (P2R).

RESULTS

Although a sufficiently large peak Ito was required to result in the type 1 ECG pattern and P2R, increasing the late component of Ito converted type 1 ECG to type 2 ECG and suppressed P2R. Increasing the peak Ito promoted spiral wave breakup, potentiating the transition from tachycardia to fibrillation, but increasing the late Ito prevented spiral wave breakup by flattening the action potential duration restitution and preventing P2R. A sufficiently large ICa,L conductance was needed for P2R to occur, but once above the critical conductance, blocking ICa,L promoted P2R. However, selectively blocking the window and late components of ICa,L suppressed P2R, countering the effect of the late Ito. Blocking either the peak or late components of sodium current promoted P2R, with the late sodium current blockade having the larger effect. As expected, increasing other potassium currents potentiated P2R, with ATP-sensitive potassium current exhibiting a larger effect than rapid and slow component of the delayed rectifier potassium current.

CONCLUSIONS

The peak Ito promotes type 1 ECG and P2R, whereas the late Ito converts type 1 ECG to type 2 ECG and suppresses P2R. Blocking the peak ICa,L and either the peak or the late sodium current promotes P2R, whereas blocking the window and late ICa,L suppresses P2R. These results provide important insights into the mechanisms of arrhythmogenesis and potential therapeutic targets for treatment of Brugada syndrome. Graphic Abstract: A graphic abstract is available for this article.

ESC working group on cardiac cellular electrophysiology position paper: relevance, opportunities, and limitations of experimental models for cardiac electrophysiology research

Cardiac arrhythmias are a major cause of death and disability. A large number of experimental cell and animal models have been developed to study arrhythmogenic diseases. These models have provided important insights into the underlying arrhythmia mechanisms and translational options for their therapeutic management. This position paper from the ESC Working Group on Cardiac Cellular Electrophysiology provides an overview of (i) currently available in vitro, ex vivo, and in vivo electrophysiological research methodologies, (ii) the most commonly used experimental (cellular and animal) models for cardiac arrhythmias including relevant species differences, (iii) the use of human cardiac tissue, induced pluripotent stem cell (hiPSC)-derived and in silico models to study cardiac arrhythmias, and (iv) the availability, relevance, limitations, and opportunities of these cellular and animal models to recapitulate specific acquired and inherited arrhythmogenic diseases, including atrial fibrillation, heart failure, cardiomyopathy, myocarditis, sinus node, and conduction disorders and channelopathies. By promoting a better understanding of these models and their limitations, this position paper aims to improve the quality of basic research in cardiac electrophysiology, with the ultimate goal to facilitate the clinical translation and application of basic electrophysiological research findings on arrhythmia mechanisms and therapies.

Basis for the Induction of Tissue-Level Phase-2 Reentry as a Repolarization Disorder in the Brugada Syndrome

Aims. Human action potentials in the Brugada syndrome have been characterized by delayed or even complete loss of dome formation, especially in the right ventricular epicardial layers. Such a repolarization pattern is believed to trigger phase-2 reentry (P2R); however, little is known about the conditions necessary for its initiation. This study aims to determine the specific mechanisms that facilitate P2R induction in Brugada-affected cardiac tissue in humans. Methods. Ionic models for Brugada syndrome in human epicardial cells were developed and used to study the induction of P2R in cables, sheets, and a three-dimensional model of the right ventricular free wall. Results. In one-dimensional cables, P2R can be induced by adjoining lost-dome and delayed-dome regions, as mediated by tissue excitability and transmembrane voltage profiles, and reduced coupling facilitates its induction. In two and three dimensions, sustained reentry can arise when three regions (delayed-dome, lost-dome, and normal epicardium) are present. Conclusions. Not only does P2R induction by Brugada syndrome require regions of action potential with delayed-dome and lost-dome, but in order to generate a sustained reentry from a triggered waveback multiple factors are necessary, including heterogeneity in action potential distribution, tissue coupling, direction of stimulation, the shape of the late plateau, the duration of lost-dome action potentials, and recovery of tissue excitability, which is predominantly modulated by tissue coupling.

Differential thermosensitivity in mixed syndrome cardiac sodium channel mutants

Cardiac arrhythmias are often associated with mutations in SCN5A the gene that encodes the cardiac paralogue of the voltage-gated sodium channel, NaV 1.5. The NaV 1.5 mutants R1193Q and E1784K give rise to both long QT and Brugada syndromes. Various environmental factors, including temperature, may unmask arrhythmia. We sought to determine whether temperature might be an arrhythmogenic trigger in these two mixed syndrome mutants. Whole-cell patch clamp was used to measure the biophysical properties of NaV 1.5 WT, E1784K and R1193Q mutants. Recordings were performed using Chinese hamster ovary (CHOk1) cells transiently transfected with the NaV 1.5 α subunit (WT, E1784K, or R1193Q), β1 subunit, and eGFP. The channels' voltage-dependent and kinetic properties were measured at three different temperatures: 10ºC, 22ºC, and 34ºC. The E1784K mutant is more thermosensitive than either WT or R1193Q channels. When temperature is elevated from 22°C to 34°C, there is a greater increase in late INa and use-dependent inactivation in E1784K than in WT or R1193Q. However, when temperature is lowered to 10°C, the two mutants show a decrease in channel availability. Action potential modelling using Q10 fit values, extrapolated to physiological and febrile temperatures, show a larger transmural voltage gradient in E1784K compared to R1193Q and WT with hyperthermia. The E1784K mutant is more thermosensitive than WT or R1193Q channels. This enhanced thermosensitivity may be a mechanism for arrhythmogenesis in patients with E1784K sodium channels.

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.

In silico investigation of the short QT syndrome, using human ventricle models incorporating electromechanical coupling

INTRODUCTION

Genetic forms of the Short QT Syndrome (SQTS) arise due to cardiac ion channel mutations leading to accelerated ventricular repolarization, arrhythmias and sudden cardiac death. Results from experimental and simulation studies suggest that changes to refractoriness and tissue vulnerability produce a substrate favorable to re-entry. Potential electromechanical consequences of the SQTS are less well-understood. The aim of this study was to utilize electromechanically coupled human ventricle models to explore electromechanical consequences of the SQTS.

METHODS AND RESULTS

The Rice et al. mechanical model was coupled to the ten Tusscher et al. ventricular cell model. Previously validated K(+) channel formulations for SQT variants 1 and 3 were incorporated. Functional effects of the SQTS mutations on [Ca(2+)] i transients, sarcomere length shortening and contractile force at the single cell level were evaluated with and without the consideration of stretch-activated channel current (I sac). Without I sac, at a stimulation frequency of 1Hz, the SQTS mutations produced dramatic reductions in the amplitude of [Ca(2+)] i transients, sarcomere length shortening and contractile force. When I sac was incorporated, there was a considerable attenuation of the effects of SQTS-associated action potential shortening on Ca(2+) transients, sarcomere shortening and contractile force. Single cell models were then incorporated into 3D human ventricular tissue models. The timing of maximum deformation was delayed in the SQTS setting compared to control.

CONCLUSION

The incorporation of I sac appears to be an important consideration in modeling functional effects of SQT 1 and 3 mutations on cardiac electro-mechanical coupling. Whilst there is little evidence of profoundly impaired cardiac contractile function in SQTS patients, our 3D simulations correlate qualitatively with reported evidence for dissociation between ventricular repolarization and the end of mechanical systole.

Altered sinoatrial node function and intra-atrial conduction in murine gain-of-function Scn5a+/ΔKPQ hearts suggest an overlap syndrome

Mutations in SCN5A, the gene encoding the pore-forming subunit of cardiac Na⁺ channels, cause a spectrum of arrhythmic syndromes. Of these, sinoatrial node (SAN) dysfunction occurs in patients with both loss- and gain-of-function SCN5A mutations. We explored for corresponding alterations in SAN function and intracardiac conduction and clarified possible mechanisms underlying these in an established mouse long QT syndrome type 3 model carrying a mutation equivalent to human SCN5A-ΔKPQ. Electrophysiological characterizations of SAN function in living animals and in vitro sinoatrial preparations were compared with cellular SAN and two-dimensional tissue models exploring the consequences of Scn5a+/ΔKPQ mutations. Scn5a+/ΔKPQ mice showed prolonged electrocardiographic QT and corrected QT intervals confirming long QT phenotypes. They showed frequent episodes of sinus bradycardia, sinus pause/arrest, and significantly longer sinus node recovery times, suggesting compromised pacemaker activity compared with wild-type mice. Electrocardiographic waveforms suggested depressed intra-atrial, atrioventricular node, and intraventricular conduction in Scn5a+/ΔKPQ mice. Isolated Scn5a+/ΔKPQ sinoatrial preparations similarly showed lower mean intrinsic heart rates and overall slower conduction through the SAN to the surrounding atrium than did wild-type preparations. Computer simulations of both single SAN cells as well as two-dimensional SAN-atrial models could reproduce the experimental observations of impaired pacemaker and sinoatrial conduction in terms of changes produced by both augmented tail and reduced total Na⁺ currents, respectively. In conclusion, the gain-of-function long QT syndrome type 3 murine Scn5a+/ΔKPQ cardiac system, in overlap with corresponding features reported in loss-of-function Na⁺ channel mutations, shows compromised SAN pacemaker and conduction function explicable in modeling studies through a combination of augmented tail and reduced peak Na⁺ currents.

An efficient method for modeling kinetic behavior of channel proteins in cardiomyocytes

Characterization of the kinetic and conformational properties of channel proteins is a crucial element in the integrative study of congenital cardiac diseases. The proteins of the ion channels of cardiomyocytes represent an important family of biological components determining the physiology of the heart. Some computational studies aiming to understand the mechanisms of the ion channels of cardiomyocytes have concentrated on Markovian stochastic approaches. Mathematically, these approaches employ Chapman-Kolmogorov equations coupled with partial differential equations. As the scale and complexity of such subcellular and cellular models increases, the balance between efficiency and accuracy of algorithms becomes critical. We have developed a novel two-stage splitting algorithm to address efficiency and accuracy issues arising in such modeling and simulation scenarios. Numerical experiments were performed based on the incorporation of our newly developed conformational kinetic model for the rapid delayed rectifier potassium channel into the dynamic models of human ventricular myocytes. Our results show that the new algorithm significantly outperforms commonly adopted adaptive Runge-Kutta methods. Furthermore, our parallel simulations with coupled algorithms for multicellular cardiac tissue demonstrate a high linearity in the speedup of large-scale cardiac simulations.

Increased vulnerability of human ventricle to re-entrant excitation in hERG-linked variant 1 short QT syndrome

The short QT syndrome (SQTS) is a genetically heterogeneous condition characterized by abbreviated QT intervals and an increased susceptibility to arrhythmia and sudden death. This simulation study identifies arrhythmogenic mechanisms in the rapid-delayed rectifier K⁺ current (I_Kr)-linked SQT1 variant of the SQTS. Markov chain (MC) models were found to be superior to Hodgkin-Huxley (HH) models in reproducing experimental data regarding effects of the N588K mutation on KCNH2-encoded hERG. These ionic channel models were then incorporated into human ventricular action potential (AP) models and into 1D and 2D idealised and realistic transmural ventricular tissue simulations and into a 3D anatomical model. In single cell models, the N588K mutation abbreviated ventricular cell AP duration at 90% repolarization (APD₉₀) and decreased the maximal transmural voltage heterogeneity (δV) during APs. This resulted in decreased transmural heterogeneity of APD₉₀ and of the effective refractory period (ERP): effects that are anticipated to be anti-arrhythmic rather than pro-arrhythmic. However, with consideration of transmural heterogeneity of I_Kr density in the intact tissue model based on the ten Tusscher-Noble-Noble-Panfilov ventricular model, not only did the N588K mutation lead to QT-shortening and increases in T-wave amplitude, but δV was found to be augmented in some local regions of ventricle tissue, resulting in increased tissue vulnerability for uni-directional conduction block and predisposing to formation of re-entrant excitation waves. In 2D and 3D tissue models, the N588K mutation facilitated and maintained re-entrant excitation waves due to the reduced substrate size necessary for sustaining re-entry. Thus, in SQT1 the N588K-hERG mutation facilitates initiation and maintenance of ventricular re-entry, increasing the lifespan of re-entrant spiral waves and the stability of scroll waves in 3D tissue.

Sudden cardiac death may arise in individuals with diseased heart tissue, or in apparently healthy subjects who suffer from genetic defects in ‘ion channel’ proteins, which increase cardiac arrhythmia risk and are associated with significant morbidity and mortality. One rare, though serious, genetic condition is the ‘short QT syndrome’ (SQTS). Although it is now known that the KCNH2-encoded N588K-hERG mutation is associated with the main (SQT1) variant of the SQTS, the mechanisms by which ventricular arrhythmia is initiated and sustained are still unclear due to lack of genotypically accurate experimental models. In this study, we used sophisticated multi-scale computer models of human ventricles in order to investigate the pro-arrhythmic effects of the N588K hERG mutation. It was found that the mutation accelerated the ventricular repolarization process, produced augmented electrical heterogeneity in some local regions of the tissue, leading to increased risk of arrhythmia genesis. It was also found that accelerated ventricular repolarization reduced the substrate size of the tissue required to sustain re-entrant circuits in both two and three dimensions. This study provides new mechanistic insight into understanding of how changes to hERG channel function in SQT1 lead to exacerbated ventricular arrhythmia risk in this inherited arrhythmia syndrome.

Mathematical models of canine right and left atria cardiomyocytes

The aim of this study is to build two mathematical models of canine ionic currents specific to right atria and left atria. The canine left atria mathematical model was firstly modified from the Ramirez-Nattel-Courtemanche (RNC) model using the recently available experimental data of ionic currents and was further developed based on our own experimental data. A model of right atria was then built by considering the differences between right atria and left atria. The two developed models well reproduced the experimental data on action potential morphology, the rate dependence, and action potential duration restitution. They are useful for investigating the mechanisms underlying the heterogeneity of canine regional action potentials and would help the simulation of whole heart excitation propagation and cardiac arrhythmia in the near future.

Peroxynitrite formation mediates LPC-induced augmentation of cardiac late sodium currents

Lysophosphatidylcholine (LPC) accumulates in the ischaemic myocardium and is arrhythmogenic. We have examined the mechanisms underlying the effects of LPC on the late cardiac Na(+) current (I(L)Na). Na(+) currents were recorded in HEK293 cells expressing Na(V)1.5 and isolated rat ventricular myocytes. LPC enhanced recombinant I(L)Na, while it reduced peak Na(+) current. Computer modeling of human ventricular myocyte action potentials predicted a marked duration prolonging effect and arrhythmogenic potential due to these effects of LPC on peak and late currents. Enhancement of recombinant I(L)Na was suppressed by the antioxidant ascorbic acid and by the NADPH oxidase inhibitor DPI. Inhibitors of the mitochondrial electron transport chain (rotenone, TTFA and myxothiazol) were without effect on LPC responses. The superoxide donor pyrogallol was without effect on I(L)Na. Enhancement of I(L)Na was abrogated by the NOS inhibitors l-NAME and 7-nitroindazole, while LPC induced an l-NAME-sensitive production of NO, measured as enhanced DAF-FM fluorescence, in both HEK293 cells and ventricular myocytes. Despite this, the NO donors SNAP and SNP caused no change in I(L)Na. However, SNAP enhanced TTX-sensitive recombinant and native I(L)Na in the presence of pyrogallol, suggesting peroxynitrite formation as a mediator of the response to LPC. In support of this, the peroxynitrite scavenger FeTPPS prevented the response of I(L)Na to LPC. Peroxynitrite formation provides a novel mechanism by which LPC regulates the late cardiac Na(+) current.