Search

Cardiac fibrosis occurs in many forms of heart disease and is considered to be one of the main arrhythmogenic factors. Regions with a high density of fibroblasts are likely to cause blocks of wave propagation that give rise to dangerous cardiac arrhythmias. Therefore, studies of the wave propagation through these regions are very important, yet the precise mechanisms leading to arrhythmia formation in fibrotic cardiac tissue remain poorly understood. Particularly, it is not clear how wave propagation is organized at the cellular level, as experiments show that the regions with a high percentage of fibroblasts (65-75%) are still conducting electrical signals, whereas geometric analysis of randomly distributed conducting and non-conducting cells predicts connectivity loss at 40% at the most (percolation threshold). To address this question, we used a joint in vitro-in silico approach, which combined experiments in neonatal rat cardiac monolayers with morphological and electrophysiological computer simulations. We have shown that the main reason for sustainable wave propagation in highly fibrotic samples is the formation of a branching network of cardiomyocytes. We have successfully reproduced the morphology of conductive pathways in computer modelling, assuming that cardiomyocytes align their cytoskeletons to fuse into cardiac syncytium. The electrophysiological properties of the monolayers, such as conduction velocity, conduction blocks and wave fractionation, were reproduced as well. In a virtual cardiac tissue, we have also examined the wave propagation at the subcellular level, detected wavebreaks formation and its relation to the structure of fibrosis and, thus, analysed the processes leading to the onset of arrhythmias. [ABSTRACT FROM AUTHOR]

Rotors are functional reentry sources identified in clinically relevant cardiac arrhythmias, such as ventricular and atrial fibrillation. Ablation targeting rotor sites has resulted in arrhythmia termination. Recent clinical, experimental and modelling studies demonstrate that rotors are often anchored around fibrotic scars or regions with increased fibrosis. However the mechanisms leading to abundance of rotors at these locations are not clear. The current study explores the hypothesis whether fibrotic scars just serve as anchoring sites for the rotors or whether there are other active processes which drive the rotors to these fibrotic regions. Rotors were induced at different distances from fibrotic scars of various sizes and degree of fibrosis. Simulations were performed in a 2D model of human ventricular tissue and in a patient-specific model of the left ventricle of a patient with remote myocardial infarction. In both the 2D and the patient-specific model we found that without fibrotic scars, the rotors were stable at the site of their initiation. However, in the presence of a scar, rotors were eventually dynamically anchored from large distances by the fibrotic scar via a process of dynamical reorganization of the excitation pattern. This process coalesces with a change from polymorphic to monomorphic ventricular tachycardia. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Atrial fibrillation (AF) is the most frequent form of arrhythmia occurring in the industrialized world. Because of its complex nature, each identified form of AF requires specialized treatment. Thus, an in-depth understanding of the bases of these arrhythmias is essential for therapeutic development. A variety of experimental studies aimed at understanding the mechanisms of AF are performed using primary cultures of neonatal rat atrial cardiomyocytes (NRAMs). Previously, we have shown that the distinct advantage of NRAM cultures is that they allow standardized, systematic, robust re-entry induction in the presence of a constitutively-active acetylcholine-mediated K+ current (IKACh-c). Experimental studies dedicated to mechanistic explorations of AF, using these cultures, often use computer models for detailed electrophysiological investigations. However, currently, no mathematical model for NRAMs is available. Therefore, in the present study we propose the first model for the action potential (AP) of a NRAM with constitutively-active acetylcholine-mediated K+ current (IKACh-c). The descriptions of the ionic currents were based on patch-clamp data obtained from neonatal rats. Our monolayer model closely mimics the action potential duration (APD) restitution and conduction velocity (CV) restitution curves presented in our previous in vitro studies. In addition, the model reproduces the experimentally observed dynamics of spiral wave rotation, in the absence and in the presence of drug interventions, and in the presence of localized myofibroblast heterogeneities. [ABSTRACT FROM AUTHOR]

Waveblock formation is the main cause of reentry. We have performed a comprehensive numerical modeling study of block formation due to anisotropy in Ten Tusscher and Panfilov (2006) ionic model for human ventricular tissue. We have examined the border between different areas of myocardial fiber alignment and have shown that blockage can occur for a wave traveling from a transverse fiber area to a longitudinal one. Such blockage occurs for reasonable values of the anisotropy ratio (AR): from 2.4 to 6.2 with respect to propagation velocities. This critical AR decreases by the suppression of INa and ICa, slightly decreases by the suppression of IKr and IKs, and substantially increases by the suppression of IK 1. Hyperkalemia affects the block formation in a complex, biphasic way. We provide examples of reentry formation due to the studied effects and have concluded that the suppression of IK 1 should be the most effective way to prevent waveblock at the areas of abrupt change in anisotropy. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Acute regional ischemia in the heart can lead to cardiac arrhythmias such as ventricular fibrillation (VF), which in turn compromise cardiac output and result in secondary global cardiac ischemia. The secondary ischemia may influence the underlying arrhythmia mechanism. A recent clinical study documents the effect of global cardiac ischaemia on the mechanisms of VF. During 150 seconds of global ischemia the dominant frequency of activation decreased, while after reperfusion it increased rapidly. At the same time the complexity of epicardial excitation, measured as the number of epicardical phase singularity points, remained approximately constant during ischemia. Here we perform numerical studies based on these clinical data and propose explanations for the observed dynamics of the period and complexity of activation patterns. In particular, we study the effects on ischemia in pseudo-1D and 2D cardiac tissue models as well as in an anatomically accurate model of human heart ventricles. We demonstrate that the fall of dominant frequency in VF during secondary ischemia can be explained by an increase in extracellular potassium, while the increase during reperfusion is consistent with washout of potassium and continued activation of the ATP-dependent potassium channels. We also suggest that memory effects are responsible for the observed complexity dynamics. In addition, we present unpublished clinical results of individual patient recordings and propose a way of estimating extracellular potassium and activation of ATP-dependent potassium channels from these measurements. [ABSTRACT FROM AUTHOR]

We develop a numerical approach based on our recent analytical model of fiber structure in the left ventricle of the human heart. A special curvilinear coordinate system is proposed to analytically include realistic ventricular shape and myofiber directions. With this anatomical model, electrophysiological simulations can be performed on a rectangular coordinate grid. We apply our method to study the effect of fiber rotation and electrical anisotropy of cardiac tissue (i.e., the ratio of the conductivity coefficients along and across the myocardial fibers) on wave propagation using the ten Tusscher–Panfilov (2006) ionic model for human ventricular cells. We show that fiber rotation increases the speed of cardiac activation and attenuates the effects of anisotropy. Our results show that the fiber rotation in the heart is an important factor underlying cardiac excitation. We also study scroll wave dynamics in our model and show the drift of a scroll wave filament whose velocity depends non-monotonically on the fiber rotation angle; the period of scroll wave rotation decreases with an increase of the fiber rotation angle; an increase in anisotropy may cause the breakup of a scroll wave, similar to the mother rotor mechanism of ventricular fibrillation. [ABSTRACT FROM AUTHOR]

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 blockers, and the properties of the emergent patterns can be studied by optical mapping of the voltage and calcium. [ABSTRACT FROM AUTHOR]

Action potential duration (APD) heterogeneity of cardiac tissue is one of the most important factors underlying initiation of deadly cardiac arrhythmias. In many cases such heterogeneity can be measured at tissue level only, while it originates from differences between the individual cardiac cells. The extent of heterogeneity at tissue and single cell level can differ substantially and in many cases it is important to know the relation between them. Here we study effects from cell coupling on APD heterogeneity in cardiac tissue in numerical simulations using the ionic TP06 model for human cardiac tissue. We show that the effect of cell coupling on APD heterogeneity can be described mathematically using a Gaussian Green's function approach. This relates the problem of electrotonic interactions to a wide range of classical problems in physics, chemistry and biology, for which robust methods exist. We show that, both for determining effects of tissue heterogeneity from cell heterogeneity (forward problem) as well as for determining cell properties from tissue level measurements (inverse problem), this approach is promising. We illustrate the solution of the forward and inverse problem on several examples of 1D and 2D systems. [ABSTRACT FROM AUTHOR]

We introduce an electromechanical model for human cardiac tissue which couples a biophysical model of cardiac excitation (Tusscher, Noble, Noble, Panfilov, 2006) and tension development (adjusted Niederer, Hunter, Smith, 2006 model) with a discrete elastic mass-lattice model. The equations for the excitation processes are solved with a finite difference approach, and the equations of the mass-lattice model are solved using Verlet integration. This allows the coupled problem to be solved with high numerical resolution. Passive mechanical properties of the mass-lattice model are described by a generalized Hooke's law for finite deformations (Seth material). Active mechanical contraction is initiated by changes of the intracellular calcium concentration, which is a variable of the electrical model. Mechanical deformation feeds back on the electrophysiology via stretch-activated ion channels whose conductivity is controlled by the local stretch of the medium. We apply the model to study how stretch-activated currents affect the action potential shape, restitution properties, and dynamics of spiral waves, under constant stretch, and dynamic stretch caused by active mechanical contraction. We find that stretch conditions substantially affect these properties via stretch-activated currents. In constantly stretched medium, we observe a substantial decrease in conduction velocity, and an increase of action potential duration; whereas, with dynamic stretch, action potential duration is increased only slightly, and the conduction velocity restitution curve becomes biphasic. Moreover, in constantly stretched medium, we find an increase of the core size and period of a spiral wave, but no change in rotation dynamics; in contrast, in the dynamically stretching medium, we observe spiral drift. Our results may be important to understand how altered stretch conditions affect the heart's functioning. [ABSTRACT FROM AUTHOR]