Search

Cardiac fibrosis stands as one of the most critical conditions leading to lethal cardiac arrhythmias. Identifying the precise location of cardiac fibrosis is crucial for planning clinical interventions in patients with various forms of ventricular and atrial arrhythmias. As fibrosis impedes and alters the path of electrical waves, detecting fibrosis in the heart can be achieved through analyzing electrical signals recorded from its surface. In current clinical practices, it has become feasible to record electrical activity from both the endocardial and epicardial surfaces of the heart. This paper presents a computational method for reconstructing 3D fibrosis using unipolar electrograms obtained from both surfaces of the ventricles. The proposed method calculates the percentage of fibrosis in various ventricular segments by analyzing the local activation times and peak-to-peak amplitudes of the electrograms. Initially, the method was tested using simulated data representing idealized fibrosis in a heart segment; subsequently, it was validated in the left ventricle with fibrosis obtained from a patient with nonischemic cardiomyopathy. The method successfully determined the location and extent of fibrosis in 204 segments of the left ventricle model with an average error of 0.0±4.3% (N = 204). Moreover, the method effectively detected fibrotic scars in the mid-myocardial region, a region known to present challenges in accurate detection using electrogram amplitude as the primary criterion., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Okenov et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Nonlinear waves of electrical excitation initiate cardiac contraction. Abnormal wave propagation in the heart, e.g., spiral waves, can lead to sudden cardiac arrest. This study analyzed the dynamics of spiral waves under the influence of an instability called negative filament tension, and examined how the spiral waves can be eliminated through high-frequency pacing. A generic anatomical model of the left ventricle of the human heart and the Aliev-Panfilov model for cardiac tissue were used. The study showed that the source of such arrhythmia is elongated filaments with lengths that can be 10-20 times greater than the characteristic thickness of the heart wall. In anisotropic tissue, the filament elongated before it was annihilated at the base of the heart. The spiral waves were eliminated through overdrive pacing with stimulation periods from 0.8 to 0.95 relative to the spiral wave period. The minimum time for the expulsion was about 10 s.

Cx43, a major cardiac connexin, forms precursor hemichannels that accrue at the intercalated disc to assemble as gap junctions. While gap junctions are crucial for electrical conduction in the heart, little is known about the potential roles of hemichannels. Recent evidence suggests that inhibiting Cx43 hemichannel opening with Gap19 has antiarrhythmic effects. Here, we used multiple electrophysiology, imaging, and super-resolution techniques to understand and define the conditions underlying Cx43 hemichannel activation in ventricular cardiomyocytes, their contribution to diastolic Ca2+ release from the sarcoplasmic reticulum, and their impact on electrical stability. We showed that Cx43 hemichannels were activated during diastolic Ca2+ release in single ventricular cardiomyocytes and cardiomyocyte cell pairs from mice and pigs. This activation involved Cx43 hemichannel Ca2+ entry and coupling to Ca2+ release microdomains at the intercalated disc, resulting in enhanced Ca2+ dynamics. Hemichannel opening furthermore contributed to delayed afterdepolarizations and triggered action potentials. In single cardiomyocytes, cardiomyocyte cell pairs, and arterially perfused tissue wedges from failing human hearts, increased hemichannel activity contributed to electrical instability compared with nonfailing rejected donor hearts. We conclude that microdomain coupling between Cx43 hemichannels and Ca2+ release is a potentially novel, targetable mechanism of cardiac arrhythmogenesis in heart failure.

High-voltage electrical defibrillation remains the only reliable method of quickly controlling life-threatening cardiac arrhythmias. This paper is devoted to studying an alternative approach, low-voltage cardioversion (LVC), which is based on ideas from non-linear dynamics and aims to remove sources of cardiac arrhythmias by applying high-frequency stimulation to cardiac tissue. We perform a detailed in-silico study of the elimination of arrhythmias caused by rotating spiral waves in a TP06 model of human cardiac tissue. We consider three parameter sets with slopes of the APD restitution curve of 0.7, 1.1 and 1.4, and we study LVC at the baseline and under the blocking of INa and ICaL and under the application of the drugs verapamil and amiodarone. We show that pacing can remove spiral waves; however, its efficiency can be substantially reduced by dynamic instabilities. We classify these instabilities and show that the blocking of INa and the application of amiodarone increase the efficiency of the method, while the blocking of ICaL and the application of verapamil decrease the efficiency. We discuss the mechanisms and the possible clinical applications resulting from our study.

Background: Long QT syndromes (LQTS) arise from many genetic and nongenetic causes with certain characteristic ECG features preceding polymorphic ventricular tachyarrhythmias (PVTs). However, how the many molecular causes result in these characteristic ECG patterns and how these patterns are mechanistically linked to the spontaneous initiation of PVT remain poorly understood., Methods: Anatomic human ventricle and simplified tissue models were used to investigate the mechanisms of spontaneous initiation of PVT in LQTS., Results: Spontaneous initiation of PVT was elicited by gradually ramping up I Ca,L to simulate the initial phase of a sympathetic surge or by changing the heart rate, reproducing the different genotype-dependent clinical ECG features. In LQTS type 2 (LQT2) and LQTS type 3 (LQT3), T-wave alternans was observed followed by premature ventricular complexes (PVCs). Compensatory pauses occurred resulting in short-long-short sequences. As I Ca,L increased further, PVT episodes occurred, always preceded by a short-long-short sequence. However, in LQTS type 1 (LQT1), once a PVC occurred, it always immediately led to an episode of PVT. Arrhythmias in LQT2 and LQT3 were bradycardia dependent, whereas those in LQT1 were not. In all 3 genotypes, PVCs always originated spontaneously from the steep repolarization gradient region and manifested on ECG as R-on-T. We call this mechanism R-from-T, to distinguish it from the classic explanation of R-on-T arrhythmogenesis in which an exogenous PVC coincidentally encounters a repolarizing region. In R-from-T, the PVC and the T wave are causally related, where steep repolarization gradients combined with enhanced I Ca,L lead to PVCs emerging from the T wave. Since enhanced I Ca,L was required for R-from-T to occur, suppressing window I Ca,L effectively prevented arrhythmias in all 3 genotypes., Conclusions: Despite the complex molecular causes, these results suggest that R-from-T is likely a common mechanism for PVT initiation in LQTS. Targeting I Ca,L properties, such as suppressing window I Ca,L or preventing excessive I Ca,L increase, could be an effective unified therapy for arrhythmia prevention in LQTS.

Early Afterdepolarizations, EADs, are defined as the reversal of the action potential before completion of the repolarization phase, which can result in ectopic beats. However, the series of mechanisms of EADs leading to these ectopic beats and related cardiac arrhythmias are not well understood. Therefore, we aimed to investigate the influence of this single cell behavior on the whole heart level. For this study we used a modified version of the Ten Tusscher-Panfilov model of human ventricular cells (TP06) which we implemented in a 3D ventricle model including realistic fiber orientations. To increase the likelihood of EAD formation at the single cell level, we reduced the repolarization reserve (RR) by reducing the rapid delayed rectifier Potassium current and raising the L-type Calcium current. Varying these parameters defined a 2D parametric space where different excitation patterns could be classified. Depending on the initial conditions, by either exciting the ventricles with a spiral formation or burst pacing protocol, we found multiple different spatio-temporal excitation patterns. The spiral formation protocol resulted in the categorization of a stable spiral (S), a meandering spiral (MS), a spiral break-up regime (SB), spiral fibrillation type B (B), spiral fibrillation type A (A) and an oscillatory excitation type (O). The last three patterns are a 3D generalization of previously found patterns in 2D. First, the spiral fibrillation type B showed waves determined by a chaotic bi-excitable regime, i.e. mediated by both Sodium and Calcium waves at the same time and in same tissue settings. In the parameter region governed by the B pattern, single cells were able to repolarize completely and different (spiral) waves chaotically burst into each other without finishing a 360 degree rotation. Second, spiral fibrillation type A patterns consisted of multiple small rotating spirals. Single cells failed to repolarize to the resting membrane potential hence prohibiting the Sodium channel gates to recover. Accordingly, we found that Calcium waves mediated these patterns. Third, a further reduction of the RR resulted in a more exotic parameter regime whereby the individual cells behaved independently as oscillators. The patterns arose due to a phase-shift of different oscillators as disconnection of the cells resulted in continuation of the patterns. For all patterns, we computed realistic 9 lead ECGs by including a torso model. The B and A type pattern exposed the behavior of Ventricular Tachycardia (VT). We conclude that EADs at the single cell level can result in different types of cardiac fibrillation at the tissue and 3D ventricle level.

Fibrosis and altered gap junctional coupling are key features of ventricular remodelling and are associated with abnormal electrical impulse generation and propagation. Such abnormalities predispose to reentrant electrical activity in the heart. In the absence of tissue heterogeneity, high-frequency impulse generation can also induce dynamic electrical instabilities leading to reentrant arrhythmias. However, because of the complexity and stochastic nature of such arrhythmias, the combined effects of tissue heterogeneity and dynamical instabilities in these arrhythmias have not been explored in detail. Here, arrhythmogenesis was studied using in vitro and in silico monolayer models of neonatal rat ventricular tissue with 30% randomly distributed cardiac myofibroblasts and systematically lowered intercellular coupling achieved in vitro through graded knockdown of connexin43 expression. Arrhythmia incidence and complexity increased with decreasing intercellular coupling efficiency. This coincided with the onset of a specialized type of spatially discordant action potential duration alternans characterized by island-like areas of opposite alternans phase, which positively correlated with the degree of connexinx43 knockdown and arrhythmia complexity. At higher myofibroblast densities, more of these islands were formed and reentrant arrhythmias were more easily induced. This is the first study exploring the combinatorial effects of myocardial fibrosis and dynamic electrical instabilities on reentrant arrhythmia initiation and complexity.

Early afterdepolarizations (EADs), which are abnormal oscillations of the membrane potential at the plateau phase of an action potential, are implicated in the development of cardiac arrhythmias like Torsade de Pointes. We carry out extensive numerical simulations of the TP06 and ORd mathematical models for human ventricular cells with EADs. We investigate the different regimes in both these models, namely, the parameter regimes where they exhibit (1) a normal action potential (AP) with no EADs, (2) an AP with EADs, and (3) an AP with EADs that does not go back to the resting potential. We also study the dependence of EADs on the rate of at which we pace a cell, with the specific goal of elucidating EADs that are induced by slow or fast rate pacing. In our simulations in two- and three-dimensional domains, in the presence of EADs, we find the following wave types: (A) waves driven by the fast sodium current and the L-type calcium current (Na-Ca-mediated waves); (B) waves driven only by the L-type calcium current (Ca-mediated waves); (C) phase waves, which are pseudo-travelling waves. Furthermore, we compare the wave patterns of the various wave-types (Na-Ca-mediated, Ca-mediated, and phase waves) in both these models. We find that the two models produce qualitatively similar results in terms of exhibiting Na-Ca-mediated wave patterns that are more chaotic than those for the Ca-mediated and phase waves. However, there are quantitative differences in the wave patterns of each wave type. The Na-Ca-mediated waves in the ORd model show short-lived spirals but the TP06 model does not. The TP06 model supports more Ca-mediated spirals than those in the ORd model, and the TP06 model exhibits more phase-wave patterns than does the ORd model.

Scroll waves are three-dimensional vortices which occur in excitable media. Their formation in the heart results in the onset of cardiac arrhythmias, and the dynamics of their filaments determine the arrhythmia type. Most studies of filament dynamics were performed in domains with simple geometries and generic description of the anisotropy of cardiac tissue. Recently, we developed an analytical model of fibre structure and anatomy of the left ventricle (LV) of the human heart. Here, we perform a systematic study of the dynamics of scroll wave filaments for the cases of positive and negative tension in this anatomical model. We study the various possible shapes of LV and different degree of anisotropy of cardiac tissue. We show that, for positive filament tension, the final position of scroll wave filament is mainly determined by the thickness of the myocardial wall but, however, anisotropy attracts the filament to the LV apex. For negative filament tension, the filament buckles, and for most cases, tends to the apex of the heart with no or slight dependency on the thickness of the LV. We discuss the mechanisms of the observed phenomena and their implications for cardiac arrhythmias.

Rotors occurring in the heart underlie the mechanisms of cardiac arrhythmias. Answering the question whether or not the location of rotors is related to local properties of cardiac tissue has important practical applications. This is because ablation of rotors has been shown to be an effective way to fight cardiac arrhythmias. In this study, we investigate, in silico, the dynamics of rotors in two-dimensional and in an anatomical model of human ventricles using a Ten Tusscher-Noble-Noble-Panfilov (TNNP) model for ventricular cells. We study the effect of small size ionic heterogeneities, similar to those measured experimentally. It is shown that such heterogeneities cannot only anchor, but can also attract, rotors rotating at a substantial distance from the heterogeneity. This attraction distance depends on the extent of the heterogeneities and can be as large as 5-6 cm in realistic conditions. We conclude that small size ionic heterogeneities can be preferred localization points for rotors and discuss their possible mechanism and value for applications., (Copyright © 2014 the American Physiological Society.)

Sudden cardiac death is often caused by cardiac arrhythmias. Recently, special attention has been given to a certain arrhythmogenic condition, the long-QT syndrome, which occurs as a result of genetic mutations or drug toxicity. The underlying mechanisms of arrhythmias, caused by the long-QT syndrome, are not fully understood. However, arrhythmias are often connected to special excitations of cardiac cells, called early afterdepolarizations (EADs), which are depolarizations during the repolarizing phase of the action potential. So far, EADs have been studied mainly in isolated cardiac cells. However, the question on how EADs at the single-cell level can result in fibrillation at the tissue level, especially in human cell models, has not been widely studied yet. In this paper, we study wave patterns that result from single-cell EAD dynamics in a mathematical model for human ventricular cardiac tissue. We induce EADs by modeling experimental conditions which have been shown to evoke EADs at a single-cell level: by an increase of L-type Ca currents and a decrease of the delayed rectifier potassium currents. We show that, at the tissue level and depending on these parameters, three types of abnormal wave patterns emerge. We classify them into two types of spiral fibrillation and one type of oscillatory dynamics. Moreover, we find that the emergent wave patterns can be driven by calcium or sodium currents and we find phase waves in the oscillatory excitation regime. From our simulations we predict that arrhythmias caused by EADs can occur during normal wave propagation and do not require tissue heterogeneities. Experimental verification of our results is possible for experiments at the cell-culture level, where EADs can be induced by an increase of the L-type calcium conductance and by the application of I[Formula: see text] blockers, and the properties of the emergent patterns can be studied by optical mapping of the voltage and calcium.

Background: One of the main factors affecting propagation of electrical waves and contraction in ventricles of the heart is anisotropy of cardiac tissue. Anisotropy is determined by orientation of myocardial fibres. Determining fibre orientation field and shape of the heart is important for anatomically accurate modelling of electrical and mechanical function of the heart. The aim of this paper is to introduce a theoretical rule-based model for anatomy and fibre orientation of the left ventricle (LV) of the heart and to compare it with experimental data. We suggest explicit analytical formulae that allow us to obtain the left ventricle form and its fibre direction field. The ventricle band concept of cardiac architecture given by Torrent-Guasp is chosen as the model postulate., Methods: In our approach, anisotropy of the heart is derived from some general principles. The LV is considered as a set of identical spiral surfaces, each of which can be produced from the other by rotation around one vertical axis. Each spiral surface is filled with non-intersecting curves which represent myocardial fibres.For model verification, we use experimental data on fibre orientation in human and canine hearts., Results: LV shape and anisotropy are represented by explicit analytical expressions in a curvilinear 3-D coordinate system. The derived fibre orientation field shows good qualitative agreement with experimental data. The model reveals the most thorough quantitative simulation of fibre angles at the LV middle zone., Conclusions: Our analysis shows that the band concept can generate realistic anisotropy of the LV. Our model shows good qualitative agreement between the simulated fibre orientation field and the experimental data on LV anisotropy, and the model can be used for various numerical simulations to study the effects of anisotropy on cardiac excitation and mechanical function.

Aims: During ageing, after infarction, in cardiomyopathies and other cardiac diseases, the percentage of fibrotic (connective) tissue may increase from 6% up to 10-35%. The presence of increased amounts of connective tissue is strongly correlated with the occurrence of arrhythmias and sudden cardiac death., Methods and Results: In this article, we investigate the role of diffuse fibrosis on wave propagation, arrhythmogenesis, and arrhythmia mechanism in human ventricular tissue using computer modelling. We show that diffuse fibrosis slows down wave propagation and increases tissue vulnerability to wave break and spiral wave formation. We also demonstrate that diffuse fibrosis increases the period of re-entrant arrhythmias and can suppress the restitution-induced transition from tachycardia to fibrillation., Conclusion: The latter suggests that mechanisms different from restitution-induced spiral break-up might be more likely to account for the onset of fibrillation in the presence of large amounts of diffuse fibrotic tissue.

Sudden cardiac death is a major cause of death in the industrialized world, claiming approximately 300,000 victims annually in the United States alone. In most cases, sudden cardiac death is caused by ventricular fibrillation (VF). Experimental studies in large animal hearts have shown that the uncoordinated contractions during VF are caused by large numbers of chaotically wandering reentrant waves of electrical activity. However, recent clinical data on VF in the human heart seem to suggest that human VF may have a markedly different organization. Here, we use a detailed model of the human ventricles, including a detailed description of cell electrophysiology, ventricular anatomy, and fiber direction anisotropy, to study the organization of human VF. We show that characteristics of our simulated VF are qualitatively similar to the clinical data. Furthermore, we find that human VF is driven by only approximately 10 reentrant sources and thus is much more organized than VF in animal hearts of comparable size, where VF is driven by approximately 50 sources. We investigate the influence of anisotropy ratio, tissue excitability, and restitution properties on the number of reentrant sources driving VF. We find that the number of rotors depends strongest on minimum action potential duration, a property that differs significantly between human and large animal hearts. Based on these findings, we suggest that the simpler spatial organization of human VF relative to VF in large animal hearts may be caused by differences in minimum action potential duration. Both the simpler spatial organization of human VF and its suggested cause may have important implications for treating and preventing this dangerous arrhythmia in humans.

Cardiovascular disease remains the leading cause of death in otherwise healthy humans. In particular, most cases of sudden cardiac death occur as a result of failure of the mechanical function of the heart which is triggered by a turbulent pattern of electrical excitation of the heart e.g., ventricular fibrillation (VF). Although the exact mechanisms of VF remain unknown, increasing evidence indicates that it is organized by multiple reentrant sources (wavelets).