Computational prediction of proarrhythmogenic effect of the V241F KCNQ1 mutation in human atrium

Influence of genetic mutations to atria vulnerability to atrial fibrillation: An in-silico 3D human atria study

BACKGROUND AND OBJECTIVE

Personalized 3D computer models of atria have been extensively implemented in the last yearsas a tool to facilitate the understanding of the mechanisms underlying different forms of arrhythmia, such as atrial fibrillation (AF). Meanwhile, genetic mutations acting on potassium channel dynamics were demonstrated to induce fibrillatory episodes in asymptomatic patients. This research study aims at assessing the effects and the atrial susceptibility to AF of three gain-of-function mutations - namely, KCNH2 T895M, KCNH2 T436M, and KCNE3-V17M - associated with AF outbreaks, using highly detailed 3D atrial models with realistic wall thickness and heterogenous histological properties.

METHODS

The 3D atrial model was generated by reconstructing segmented anatomical structures from CT scans of an AF patient. Modified versions of the Courtemanche human atrial myocyte model were used to reproduce the electrophysiological activity of the WT and of the three mutant cells. Ectopic foci (EF) were simulated in sixteen locations across the atrial mesh using an S1-S2 protocol with two S2 basic cycle lengths (BCL) and eleven coupling intervals in order to induce arrhythmias.

RESULTS

The three genetic mutations at 3D level reduced the APD90. The KCNE3-V17M mutation provoked the highest shortening (55 % in RA and LA with respect to WT), followed by KCNH2 T895M (14 % in RA and 18 % LA with respect to WT)and KCNH2 T436M (7 % in RA and 9 % LA with respect to WT). The KCNE3-V17M mutation led to arrhythmia in 67 % of the cases simulated and in 94 % of ectopic foci considered, at S2 BCL equal to 100 ms. The KCNH2 T436M and KCNH2 T895M mutations increased the vulnerability to AF in a similar way, leading to arrhythmic episodes in 7 % of the simulated conditions, at S2 BCL set to 160 ms. Overall, 60 % of the arrhythmic events generated arise in the left atrium. Spiral waves, multiple rotors and disordered electrical pattern were elicited in the presence of the KCNE3-V17M mutation, exhibiting an instantaneous mean frequency of 7.6 Hz with a mean standard deviation of 1.12 Hz. The scroll waves induced in the presence of the KCNH2 T436M and KCNH2 T895M mutations showed steadiness and regularity with an instantaneous mean frequencies in the range of 4.9 - 5.1 Hz and a mean standard deviation within 0.19 - 0.53 Hz.

CONCLUSIONS

The pro-arrhythmogenicity of the KCNE3-V17M, KCNH2 T895M and KCNH2 T436M mutations was studied and proved on personalized 3D cardiac models. The three genetic mutations were demonstrated to increase the predisposition of atrial tissue to the formation of AF-susceptible substrate in different ways based on their effects on electrophysiological properties of the atria.

Arrhythmogenic Effects of Genetic Mutations Affecting Potassium Channels in Human Atrial Fibrillation: A Simulation Study

Genetic mutations in genes encoding for potassium channel protein structures have been recently associated with episodes of atrial fibrillation in asymptomatic patients. The aim of this study is to investigate the potential arrhythmogenicity of three gain-of-function mutations related to atrial fibrillation-namely, KCNH2 T895M, KCNH2 T436M, and KCNE3-V17M-using modeling and simulation of the electrophysiological activity of the heart. A genetic algorithm was used to tune the parameters' value of the original ionic currents to reproduce the alterations experimentally observed caused by the mutations. The effects on action potentials, ionic currents, and restitution properties were analyzed using versions of the Courtemanche human atrial myocyte model in different tissues: pulmonary vein, right, and left atrium. Atrial susceptibility of the tissues to spiral wave generation was also investigated studying the temporal vulnerability. The presence of the three mutations resulted in an overall more arrhythmogenic substrate. Higher current density, action potential duration shortening, and flattening of the restitution curves were the major effects of the three mutations at the single-cell level. The genetic mutations at the tissue level induced a higher temporal vulnerability to the rotor's initiation and progression, by sustaining spiral waves that perpetuate until the end of the simulation. The mutation with the highest pro-arrhythmic effects, exhibiting the widest sustained VW and the smallest meandering rotor's tip areas, was KCNE3-V17M. Moreover, the increased susceptibility to arrhythmias and rotor's stability was tissue-dependent. Pulmonary vein tissues were more prone to rotor's initiation, while in left atrium tissues rotors were more easily sustained. Re-entries were also progressively more stable in pulmonary vein tissue, followed by the left atrium, and finally the right atrium. The presence of the genetic mutations increased the susceptibility to arrhythmias by promoting the rotor's initiation and maintenance. The study provides useful insights into the mechanisms underlying fibrillatory events caused by KCNH2 T895M, KCNH2 T436M, and KCNE3-V17M and might aid the planning of patient-specific targeted therapies.

Proarrhythmogenic Effect of the L532P and N588K KCNH2 Mutations in the Human Heart Using a 3D Electrophysiological Model

BACKGROUND

Atrial arrhythmia is a cardiac disorder caused by abnormal electrical signaling and transmission, which can result in atrial fibrillation and eventual death. Genetic defects in ion channels can cause myocardial repolarization disorders. Arrhythmia-associated gene mutations, including KCNH2 gene mutations, which are one of the most common genetic disorders, have been reported. This mutation causes abnormal QT intervals by a gain of function in the rapid delayed rectifier potassium channel (I_Kr). In this study, we demonstrated that mutations in the KCNH2 gene cause atrial arrhythmia.

METHODS

The N588K and L532P mutations were induced in the Courtemanche-Ramirez-Nattel (CRN) cell model, which was subjected to two-dimensional and three-dimensional simulations to compare the electrical conduction patterns of the wild-type and mutant-type genes.

RESULTS

In contrast to the early self-termination of the wild-type conduction waveforms, the conduction waveform of the mutant-type retained the reentrant wave (N588K) and caused a spiral break-up, resulting in irregular wave generation (L532P).

CONCLUSION

The present study confirmed that the KCNH2 gene mutation increases the vulnerability of the atrial tissue for arrhythmia.

Comparison of Electromechanical Delay during Ventricular Tachycardia and Fibrillation under Different Conductivity Conditions Using Computational Modeling

Electromechanical delay (EMD) is the time interval between local myocyte depolarization and the onset of myofiber shortening. Previously, researchers measured EMD during sinus rhythm and ectopic pacing in normal and heart failure conditions. However, to our knowledge, there are no reports regarding EMD during another type of rhythms or arrhythmia. The goal of this study was to quantify EMD during sinus rhythm, tachycardia, and ventricular fibrillation conditions. We hypothesized that EMD under sinus rhythm is longer due to isovolumetric contraction which is imprecise during arrhythmia. We used a realistic model of 3D electromechanical ventricles. During sinus rhythm, EMD was measured in the last cycle of cardiac systole under steady conditions. EMD under tachycardia and fibrillation conditions was measured during the entire simulation, resulting in multiple EMD values. We assessed EMD for the following 3 conduction velocities (CVs): 31 cm/s, 51 cm/s, and 69 cm/s. The average EMD during fibrillation condition was the shortest corresponding to 53.45 ms, 55.07 ms, and 50.77 ms, for the CVs of 31 cm/s, 51 cm/s, and 69 cm/s, respectively. The average EMD during tachycardia was 58.61 ms, 58.33 ms, and 52.50 ms for the three CVs. Under sinus rhythm with action potential duration restitution (APDR) slope 0.7, the average EMD was 66.35 ms, 66.41 ms, and 66.60 ms in line with the three CVs. This result supports our hypothesis that EMD under sinus rhythm is longer than that under tachyarrhythmia conditions. In conclusion, this study observed and quantified EMD under tachycardia and ventricular fibrillation conditions. This simulation study has widened our understanding of EMD in 3D ventricles under chaotic conditions.

Relationship Between Electrical Instability and Pumping Performance During Ventricular Tachyarrhythmia: Computational Study

There are representative electrical parameters for understanding the mechanism of reentrant waves in studies on tachyarrhythmia, namely the action potential duration (APD), dominant frequency, phase singularity, and filament. However, there are no studies that have directly identified the correlation between these electrophysiological parameters and cardiac contractility. Therefore, we have identified individual and integrative correlations between these electrical phenomena and contractility during tachyarrhythmia by deriving regression equations and also investigated the electrophysiological parameters affecting cardiac contractility during tachyarrhythmia. We simulated ventricular tachyarrhythmia with 48 types of electrical patterns by applying four reentry generation methods and changing the electrical conductivity of the potassium channel, which has the greatest effect on ventricular tissue. The mechanical responses reflecting electrical complexity were obtained through deterministic simulations of excitation-contraction coupling. We used the stroke volume and amplitude of myocardial tension (ampTens) as the variables representing contractility. We derived stochastic models through single- and multivariable regression analyses to identify the electrical parameters affecting contractility during tachyarrhythmia. In single-variable regression analysis, the APD, dominant frequency, and filament, excluding phase singularity, have statistically significant correlations with the stroke volume and ampTens. Among them, the APD has the maximum influence on these two mechanical parameters (standard beta coefficient: 0.859 for stroke volume, 0.930 for ampTens). The stochastic model using all four electrical parameters fails to accurately predict contractility owing to the multicollinearity between the APD and dominant frequency. We have rederived the multi-variable stochastic model using three electrical parameters without the APD. The filament has the greatest effect on the stroke volume stochastically (standard beta coefficient: 0.853 and 0.752). The dominant frequency has the greatest effect on ampTens statistically (standard beta coefficient: -0.813). We conclude that among the electrical parameters, the APD has the highest individual influence on mechanical contraction, and the filament has the highest integrative influence in both statistical terms.

Computational prediction of the effect of D172N KCNJ2 mutation on ventricular pumping during sinus rhythm and reentry

The understanding of cardiac arrhythmia under genetic mutations has grown in interest among researchers. Previous studies focused on the effect of the D172N mutation on electrophysiological behavior. In this study, we analyzed not only the electrophysiological activity but also the mechanical responses during normal sinus rhythm and reentry conditions by using computational modeling. We simulated four different ventricular conditions including normal case of ten Tusscher model 2006 (TTM), wild-type (WT), heterozygous (WT/D172N), and homozygous D172N mutation. The 2D simulation result (in wire-shaped mesh) showed the WT/D172N and D172N mutation shortened the action potential duration by 14%, and by 23%, respectively. The 3D electrophysiological simulation results showed that the electrical wavelength between TTM and WT conditions were identical. Under sinus rhythm condition, the WT/D172N and D172N reduced the pumping efficacy with a lower left ventricle (LV) and aortic pressures, stroke volume, ejection fraction, and cardiac output. Under the reentry conditions, the WT condition has a small probability of reentry. However, in the event of reentry, WT has shown the most severe condition. Furthermore, we found that the position of the rotor or the scroll wave substantially influenced the ventricular pumping efficacy during arrhythmia. If the rotor stays in the LV, it will cause very poor pumping performance. Graphical Abstract A model of a ventricular electromechanical system. This whole model was established to observe the effect of D172N KCNJ2 mutation on ventricular pumping behavior during sinus rhythm and reentry conditions. The model consists of two components; electrical component and mechanical component. The electrophysiological model based on ten Tusscher et al. with the IK1 D172N KCNJ2 mutation, and the myofilament dynamic (cross-bridge) model based on Rice et al. study. The 3D electrical component is a ventricular geometry based on MRI which composed of nodes representing single-cell with electrophysiological activation. The 3D ventricular mechanic is a finite element mesh composed of single-cells myofilament dynamic model. Both components were coupled with Ca2+ concentration. We used Gaussian points for the calcium interpolation from the electrical mesh to the mechanical mesh.

The effect of myocardial action potential duration on cardiac pumping efficacy: a computational study

BACKGROUND AND AIMS

Although studies on the relation between arrhythmias and the action potential duration (APD) have been carried out, most of them are based only on electrophysiological factors of the heart and lack experiments that consider cardiac mechanical and electromechanical characteristics. Therefore, we conducted this study to clarify the relevance of the shortening of APD of a cell in relation to the mechanical contraction activity of the heart and the associated risk of arrhythmia.

METHODS

The human ventricular model used in this study has two dynamic characteristics: electrophysiological conduction and mechanical contraction. The model simulating electrophysiological characteristics was consisted of lumped parameter circuit that can mimic the phenomenon of ion exchange through the cell membrane of myocyte and consisted of 214,319 tetrahedral finite elements. In contrast, the model simulating mechanical contraction characteristics was constructed to mimic cardiac contraction by means of the crossbridge of a myofilament and consisted of 14,720 hermite-based finite elements to represent a natural 3D curve of the cardiac surface. First, we performed a single cell simulation and the electrophysiological simulation according to the change of the APD by changing the electrical conductivity of the I Ks channel. Thus, we confirmed the correlation between APD and intracellular Ca2+ concentration. Then, we compared mechanical response through mechanical simulation using Ca2+ data from electrical simulation.

RESULTS

The APD and the sum of the intracellular Ca2+ concentrations showed a positive correlation. The shortened APD reduced the conduction wavelength of ventricular cells by shortening the plateau and early repolarization in myocardial cells. The decrease in APD reduced ventricular pumping efficiency by more than 60% as compared with the normal group (normal conditions). This change is caused by the decline of ventricular output owing to reduced ATP consumption during the crossbridge of myofilaments and decreased tension.

CONCLUSION

The shortening of APD owing to increased electrical conductivity of a protein channel on myocardial cells likely decreases the wavelength and the pumping efficiency of the ventricles. Additionally, it may increase tissue sensitivity to ventricular fibrillation, including reentry, and cause symptoms such as dyspnea and dizziness.

Effect of KCNQ1 G229D mutation on cardiac pumping efficacy and reentrant dynamics in ventricles: Computational study

There is growing interest in genetic arrhythmia since mutations in gene which encodes the ion channel underlie numerous arrhythmias. Hasegawa et al reported that G229D mutation in KCNQ1 underlies atrial fibrillation due to significant shortening of action potential duration (APD) in atrial cells. Here, we predicted whether KCNQ1 G229D mutation affects ventricular fibrillation generation, although it shortens APD slightly compared with the atrial cell. We analyzed the effects of G229D mutation on electrical and mechanical ventricle behaviors (not considered in previous studies). We compared action potential shapes under wild-type and mutant conditions. Electrical wave propagations through ventricles were analyzed during sinus rhythm and reentrant conditions. I_Ks enhancement due to G229D mutation shortened the APD in the ventricular cells (6%, 0.3%, and 8% for endo, M, and epi-cells, respectively). The shortened APD contributed to 7% shortening of QT intervals, 29% shortening of wavelengths, 20% decrease in intraventricular pressure, and increase in end-systolic volume 17%, end-diastolic volume 7%, and end-diastolic pressure 11%, which further resulted in reduction in stroke volume as well as cardiac output (28%), ejection fraction 33% stroke work 44%, and ATP consumption 28%. In short, using computational model of the ventricle, we predicted that G229D mutation decreased cardiac pumping efficacy and increased the vulnerability of ventricular fibrillation.

Signal Reconstruction of Pulmonary Vein Recordings Using a Phenomenological Mathematical Model: Application to Pulmonary Vein Isolation Therapy

Atrial fibrillation (AF), the most prevalent cardiac arrhythmia, is commonly initiated by ectopic beats originating from a small myocardial sleeve extending over the pulmonary veins. Pulmonary vein isolation therapy attempts to isolate the pulmonary veins from the left atrium by ablating tissue, commonly by using radiofrequency ablation. During this procedure, the cardiologist records electrical activity using a lasso catheter, and the activation pattern recorded is used as a guide toward which regions to ablate. However, poor contact between electrode and tissue can lead to important regions of electrical activity not being recorded in clinic. We reproduce these signals through the use of a phenomenological model of the cardiac action potential on a cylinder, which we fit to post-AF atrial cells, and model the bipolar electrodes of the lasso catheter by an approximation of the surface potential. The resulting activation pattern is validated by direct comparison with those of clinical recordings. A potential application of the model is to reconstruct the missing electrical activity, minimizing the impact of the information loss on the clinical procedure, and we present results to demonstrate this.

Computational simulations of the effects of the G229D KCNQ1 mutation on human atrial fibrillation

Atrial fibrillation (AF) is related to mutations at the genetic level. This includes mutations in genes that encode KCNQ1, a subunit of the I Ks channel. Here, we investigate the mechanism of gain-of-function in I Ks towards the occurrence of AF. We used the Courtemanche-Ramirez-Nattel (CRN) human atrial cell model (Am J Physiol Heart Circ Physiol 275:H301-H321, 1998) and applied the modification proposed by Hasegawa et al. (Heart Rhythm 11:67-75, 2014) to fit the behavior of I Ks due to the G229D mutation in KCNQ1 under a heterozygous mutant form. This was incorporated into two-(2D) and three-dimensional (3D) tissue models, where the mutation sustained a reentrant wave. However, under the wild-type condition, the reentrant wave terminated before the end of our simulations (in 2D, the spiral wave terminated before 10 s, while in 3D, the spiral wave terminated before 13 s). Sustained reentry under the mutation conditions also resulted in a spiral wave breakup in the 3D model, which was sustained until the end of the simulation (20 s), indicating AF.