Calcium-Voltage Coupling in the Genesis of Early and Delayed Afterdepolarizations in Cardiac Myocytes

Early afterdepolarizations (EADs) and delayed afterdepolarizations (DADs) are voltage oscillations known to cause cardiac arrhythmias. EADs are mainly driven by voltage oscillations in the repolarizing phase of the action potential (AP), while DADs are driven by spontaneous calcium (Ca) release during diastole. Because voltage and Ca are bidirectionally coupled, they modulate each other’s behaviors, and new AP and Ca cycling dynamics can emerge from this coupling. In this study, we performed computer simulations using an AP model with detailed spatiotemporal Ca cycling incorporating stochastic openings of Ca channels and ryanodine receptors to investigate the effects of Ca-voltage coupling on EAD and DAD dynamics. Simulations were complemented by experiments in mouse ventricular myocytes. We show that: 1) alteration of the Ca transient due to increased ryanodine receptor leakiness and/or sarco/endoplasmic reticulum Ca ATPase activity can either promote or suppress EADs due to the complex effects of Ca on ionic current properties; 2) spontaneous Ca waves also exhibit complex effects on EADs, but cannot induce EADs of significant amplitude without the participation of I_Ca,L; 3) lengthening AP duration and the occurrence of EADs promote DADs by increasing intracellular Ca loading, and two mechanisms of DADs are identified, i.e., Ca-wave-dependent and Ca-wave-independent; and 4) Ca-voltage coupling promotes complex EAD patterns such as EAD alternans that are not observed for solely voltage-driven EADs. In conclusion, Ca-voltage coupling combined with the nonlinear dynamical behaviors of voltage and Ca cycling play a key role in generating complex EAD and DAD dynamics observed experimentally in cardiac myocytes, whose mechanisms are complex but analyzable.     

A Study of Early Afterdepolarizations in a Model for Human Ventricular Tissue

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.     

Mechanisms of Premature Ventricular Complexes Caused by QT Prolongation

QT prolongation, due to lengthening of the action potential duration in the ventricles, is a major risk factor of lethal ventricular arrhythmias. A widely known consequence of QT prolongation is the genesis of early afterdepolarizations (EADs), which are associated with arrhythmias through the generation of premature ventricular complexes (PVCs). However, the vast majority of the EADs observed experimentally in isolated ventricular myocytes are phase-2 EADs, and whether phase-2 EADs are mechanistically linked to PVCs in cardiac tissue remains an unanswered question. In this study, we investigate the genesis of PVCs using computer simulations with eight different ventricular action potential models of various species. Based on our results, we classify PVCs as arising from two distinct mechanisms: repolarization gradient (RG)-induced PVCs and phase-2 EAD-induced PVCs. The RG-induced PVCs are promoted by increasing RG and L-type calcium current and are insensitive to gap junction coupling. EADs are not required for this PVC mechanism. In a paced beat, a single or multiple PVCs can occur depending on the properties of the RG. In contrast, phase-2 EAD-induced PVCs occur only when the RG is small and are suppressed by increasing RG and more sensitive to gap junction coupling. Unlike with RG-induced PVCs, in each paced beat, only a single EAD-induced PVC can occur no matter how many EADs in an action potential. In the wide parameter ranges we explore, RG-induced PVCs can be observed in all models, but the EAD-induced PVCs can only be observed in five of the eight models. The links between these two distinct PVC mechanisms and arrhythmogenesis in animal experiments and clinical settings are discussed.     

Dynamics of early afterdepolarization-mediated triggered activity in cardiac monolayers

Early afterdepolarizations (EADs) are voltage oscillations that occur during the repolarizing phase of the cardiac action potential and cause cardiac arrhythmias in a variety of clinical settings. EADs occur in the setting of reduced repolarization reserve and increased inward-over-outward currents, which intuitively explains the repolarization delay but does not mechanistically explain the time-dependent voltage oscillations that are characteristic of EADs. In a recent theoretical study, we identified a dual Hopf-homoclinic bifurcation as a dynamical mechanism that causes voltage oscillations during EADs, depending on the amplitude and kinetics of the L-type Ca(2+) channel (LTCC) current relative to the repolarizing K(+) currents. Here we demonstrate this mechanism experimentally. We show that cardiac monolayers exposed to the LTCC agonists BayK8644 and isoproterenol produce EAD bursts that are suppressed by the LTCC blocker nitrendipine but not by the Na(+) current blocker tetrodoxin, depletion of intracellular Ca(2+) stores with thapsigargin and caffeine, or buffering of intracellular Ca(2+) with BAPTA-AM. These EAD bursts exhibited a key dynamical signature of the dual Hopf-homoclinic bifurcation mechanism, namely, a gradual slowing in the frequency of oscillations before burst termination. A detailed cardiac action potential model reproduced the experimental observations, and identified intracellular Na(+) accumulation as the likely mechanism for terminating EAD bursts. Our findings in cardiac monolayers provide direct support for the Hopf-homoclinic bifurcation mechanism of EAD-mediated triggered activity, and raise the possibility that this mechanism may also contribute to EAD formation in clinical settings such as long QT syndromes, heart failure, and increased sympathetic output.     

Mechanisms and determinants of ultralong action potential duration and slow rate-dependence in cardiac myocytes

In normal cardiac myocytes, the action potential duration (APD) is several hundred milliseconds. However, experimental studies showed that under certain conditions, APD could be excessively long (or ultralong), up to several seconds. Unlike the normal APD, the ultralong APD increases sensitively with pacing cycle length even when the pacing rate is very slow, exhibiting a sensitive slow rate-dependence. In addition, these long action potentials may or may not exhibit early afterdepolarizations (EADs). Although these phenomena are well known, the underlying mechanisms and ionic determinants remain incompletely understood. In this study, computer simulations were performed with a simplified action potential model. Modifications to the L-type calcium current (I(Ca,L)) kinetics and the activation time constant of the delayed rectifier K current were used to investigate their effects on APD. We show that: 1) the ultralong APD and its sensitive slow rate-dependence are determined by the steady-state window and pedestal I(Ca,L) currents and the activation speed and the recovery of the delayed rectifier K current; 2) whether an ultralong action potential exhibits EADs or not depends on the kinetics of I(Ca,L); 3) increasing inward currents elevates the plateau voltage, which in general prolongs APD, however, this can also shorten APD when the APD is already ultralong under certain conditions; and 4) APD alternans occurs at slow pacing rates due to the sensitive slow rate-dependence and the ionic determinants are different from the ones causing APD alternans at fast heart rates.     

Ca2+-activated chloride channel activity during Ca2+ alternans in ventricular myocytes

Cardiac alternans, defined beat-to-beat alternations in contraction, action potential (AP) morphology or cytosolic Ca transient (CaT) amplitude, is a high risk indicator for cardiac arrhythmias. We investigated mechanisms of cardiac alternans in single rabbit ventricular myocytes. CaTs were monitored simultaneously with membrane currents or APs recorded with the patch clamp technique. A strong correlation between beat-to-beat alternations of AP morphology and CaT alternans was observed. During CaT alternans application of voltage clamp protocols in form of pre-recorded APs revealed a prominent Ca²⁺-dependent membrane current consisting of a large outward component coinciding with AP phases 1 and 2, followed by an inward current during AP repolarization. Approximately 85% of the initial outward current was blocked by Cl^? channel blocker DIDS or lowering external Cl^? concentration identifying it as a Ca²⁺-activated Cl^? current (I_CaCC). The data suggest that I_CaCC plays a critical role in shaping beat-to-beat alternations in AP morphology during alternans.     

Maintenance of intercellular coupling by the antiarrhythmic peptide rotigaptide suppresses arrhythmogenic discordant alternans

Discordant action potential alternans creates large gradients of refractoriness, which are thought to be the mechanisms linking T-wave alternans to cardiac arrhythmogenesis. Since intercellular coupling acts to maintain synchronization of repolarization between cells, we hypothesized that intercellular uncoupling, such as during ischemia, would initiate discordant alternans and that restoration of intercellular coupling by the gap junction opener rotigaptide may provide a novel approach for suppressing arrhythmogenic discordant alternans. Optical mapping was used to record action potentials from ventricular epicardium of Langendorff-perfused guinea pig hearts. Threshold for spatially synchronized (i.e., concordant) alternans and discordant alternans was determined by increasing heart rate step-wise during 1) baseline, 2) treatment with rotigaptide or vehicle, and 3) global low-flow ischemia + rotigaptide or vehicle. Ischemia reduced the threshold for concordant alternans in both groups from 362 +/- 8 to 305 +/- 9 beats/min (P < 0.01) and for discordant alternans from 423 +/- 6 to 381 +/- 7 beats/min (P < 0.01). Interestingly, rotigaptide also increased the threshold for discordant alternans relative to vehicle both before (438 +/- 7 vs. 407 +/- 8 beats/min, P < 0.05) and during (394 +/- 7 vs. 364 +/- 9 beats/min, P < 0.05) ischemia. Rotigaptide increased conduction velocity and prevented conduction slowing and dispersion of repolarization during ischemia. Confocal immunofluorescence revealed that total connexin43 quantity and cellular distribution were unchanged before or after low-flow ischemia, with and without rotigaptide. However, connexin43 dephosphorylation in response to low-flow ischemia was significantly prevented by rotigaptide (15.9 +/- 7.0 vs. 0.3 +/- 6.4%, P < 0.001). These data suggest that intercellular uncoupling plays an important role in the transition from concordant to discordant alternans. By suppressing discordant alternans, repolarization gradients, and connexinx43 dephosphorylation, rotigaptide may protect against ischemia-induced arrhythmias. Drugs that selectively open gap junctions offer a novel strategy for antiarrhythmic therapy.     

Cardiac arrhythmia mechanisms in rats with heart failure induced by pulmonary hypertension

Pulmonary hypertension provokes right heart failure and arrhythmias. Better understanding of the mechanisms underlying these arrhythmias is needed to facilitate new therapeutic approaches for the hypertensive, failing right ventricle (RV). The aim of our study was to identify the mechanisms generating arrhythmias in a model of RV failure induced by pulmonary hypertension. Rats were injected with monocrotaline to induce either RV hypertrophy or failure or with saline (control). ECGs were measured in conscious, unrestrained animals by telemetry. In isolated hearts, electrical activity was measured by optical mapping and myofiber orientation by diffusion tensor-MRI. Sarcoplasmic reticular Ca(2+) handling was studied in single myocytes. Compared with control animals, the T-wave of the ECG was prolonged and in three of seven heart failure animals, prominent T-wave alternans occurred. Discordant action potential (AP) alternans occurred in isolated failing hearts and Ca(2+) transient alternans in failing myocytes. In failing hearts, AP duration and dispersion were increased; conduction velocity and AP restitution were steeper. The latter was intrinsic to failing single myocytes. Failing hearts had greater fiber angle disarray; this correlated with AP duration. Failing myocytes had reduced sarco(endo)plasmic reticular Ca(2+)-ATPase activity, increased sarcoplasmic reticular Ca(2+)-release fraction, and increased Ca(2+) spark leak. In hypertrophied hearts and myocytes, dysfunctional adaptation had begun, but alternans did not develop. We conclude that increased electrical and structural heterogeneity and dysfunctional sarcoplasmic reticular Ca(2+) handling increased the probability of alternans, a proarrhythmic predictor of sudden cardiac death. These mechanisms are potential therapeutic targets for the correction of arrhythmias in hypertensive, failing RVs.     

Nonlinear dynamics of two-dimensional cardiac action potential duration mapping model with memory

The aim of this work is the analysis of the nonlinear dynamics of two-dimensional mapping model of cardiac action potential duration (2D-map APD) with memory derived from one dimensional map (1D-map). Action potential duration (APD) restitution, which relates APD to the preceding diastolic interval (DI), is a useful tool for predicting cardiac arrhythmias. For a constant rate of stimulation the short action potential during alternans is followed by a longer DI and inversely. It has been suggested that these differences in DI are responsible for the occurrence and maintenance of APD alternans. We focus our attention on the observed bifurcations produced by a change in the stimulation period and a fixed value of a particular parameter in the model. This parameter provides new information about the dynamics of the APD with memory, such as the occurrence of bistabilities not previously described in the literature, as well as the fact that synchronization rhythms occur in different ways and in a new fashion as the stimulation frequency increases. Moreover, we show that this model is flexible enough as to accurately reflect the chaotic dynamics properties of the APD: we have highlighted the fractal structure of the strange attractor of the 2D-map APD, and we have characterized chaos by tools such as the calculation of the Lyapunov exponents, the fractal dimension and the Kolmogorov entropy, with the next objective of refining the study of the nonlinear dynamics of the duration of the action potential and to apply methods of controlling chaos.

     

Irregularly Appearing Early Afterdepolarizations in Cardiac Myocytes: Random Fluctuations or Dynamical Chaos?

Irregularly occurring early afterdepolarizations (EADs) in cardiac myocytes are traditionally hypothesized to be caused by random ion channel fluctuations. In this study, we combined 1), patch-clamp experiments in which action potentials were recorded at different pacing cycle lengths from isolated rabbit ventricular myocytes under several experimental conditions inducing EADs, including oxidative stress with hydrogen peroxide, calcium overload with BayK8644, and ionic stress with hypokalemia; 2), computer simulations using a physiologically detailed rabbit ventricular action potential model, in which repolarization reserve was reduced to generate EADs and random ion channel or path cycle length fluctuations were implemented; and 3), iterated maps with or without noise. By comparing experimental, modeling, and bifurcation analyses, we present evidence that noise-induced transitions between bistable states (i.e., between an action potential with and without an EAD) is not sufficient to account for the large variation in action potential duration fluctuations observed in experimental studies. We conclude that the irregular dynamics of EADs is intrinsically chaotic, with random fluctuations playing a nonessential, auxiliary role potentiating the complex dynamics.     

Synchronization of Early Afterdepolarizations and Arrhythmogenesis in Heterogeneous Cardiac Tissue Models

Early afterdepolarizations (EADs) are linked to both triggered arrhythmias and reentrant arrhythmias by causing premature ventricular complexes (PVCs), focal excitations, or heterogeneous tissue substrates for reentry formation. However, a critical number of cells that synchronously exhibit EADs are needed to result in arrhythmia triggers and substrates in tissue. In this study, we use mathematical modeling and computer simulations to investigate EAD synchronization and arrhythmia induction in tissue models with random cell-to-cell variations. Our major observations are as follows. Random cell-to-cell variations in action potential duration without EAD presence do not cause large dispersion of refractoriness in well-coupled tissue. In the presence of phase-2 EADs, the cells may synchronously exhibit the same number of EADs or no EADs with a very small dispersion of refractoriness, or synchronize regionally to result in large dispersion of refractoriness. In the presence of phase-3 EADs, regional synchronization leads to propagating EADs, forming PVCs in tissue. Interestingly, even though the uncoupled cells exhibit either no EAD or only a single EAD, when these cells are coupled to form a tissue, more than one PVC can occur. When the PVCs occur at different locations and time, multifocal arrhythmias are triggered, with the foci shifting in space and time in an irregular manner. The focal arrhythmias either spontaneously terminate or degenerate into reentrant arrhythmias due to heterogeneities and spatiotemporal chaotic dynamics of the foci.     

Effects of early afterdepolarizations on reentry in cardiac tissue: a simulation study

Early afterdepolarizations (EADs) are classically generated at slow heart rates when repolarization reserve is reduced by genetic diseases or drugs. However, EADs may also occur at rapid heart rates if repolarization reserve is sufficiently reduced. In this setting, spontaneous diastolic sarcoplasmic reticulum (SR) Ca release can facilitate cellular EAD formation by augmenting inward currents during the action potential plateau, allowing reactivation of the window L-type Ca current to reverse repolarization. Here, we investigated the effects of spontaneous SR Ca release-induced EADs on reentrant wave propagation in simulated one-, two-, and three-dimensional homogeneous cardiac tissue using a version of the Luo-Rudy dynamic ventricular action potential model modified to increase the likelihood of these EADs. We found: 1) during reentry, nonuniformity in spontaneous SR Ca release related to subtle differences in excitation history throughout the tissue created adjacent regions with and without EADs. This allowed EADs to initiate new wavefronts propagating into repolarized tissue; 2) EAD-generated wavefronts could propagate in either the original or opposite direction, as a single new wave or two new waves, depending on the refractoriness of tissue bordering the EAD region; 3) by suddenly prolonging local refractoriness, EADs caused rapid rotor displacement, shifting the electrical axis; and 4) rapid rotor displacement promoted self-termination by collision with tissue borders, but persistent EADs could regenerate single or multiple focal excitations that reinitiated reentry. These findings may explain many features of Torsades des pointes, such as perpetuation by focal excitations, rapidly changing electrical axis, frequent self-termination, and occasional degeneration to fibrillation.     

Evidence for direct arrhythmogenic action of endothelin.

We studied electrophysiological effects of endothelin on canine cardiac tissues. Endothelin prolonged action potential duration and decreased spontaneous firing rate of the right bundle branch cells. At a concentration of 2 x 10(-7)M the plateau phase of action potentials was flattened, followed by the abrupt occurrence of early afterdepolarizations (EADs). ET, at a concentration as low as 2 x 10(-9)M, was capable of inducing EADs although their incidence was low. The EADs were initiated from the membrane potential less negative than -30mV and were suppressed by nicardipine, suggesting the involvement of dihydropyridine-sensitive Ca2+ channels in the induction of EADs. Because EADs are considered to underlie certain types of arrhythmias endothelin per se may have arrhythmogenic action.     

Dynamic origin of spatially discordant alternans in cardiac tissue

Alternans, a condition in which there is a beat-to-beat alternation in the electromechanical response of a periodically stimulated cardiac cell, has been linked to the genesis of life-threatening ventricular arrhythmias. Optical mapping of membrane voltage (V_m) and intracellular calcium (Ca_i) on the surface of animal hearts reveals complex spatial patterns of alternans. In particular, spatially discordant alternans has been observed in which regions with a large-small-large action potential duration (APD) alternate out-of-phase adjacent to regions of small-large-small APD. However, the underlying mechanisms that lead to the initiation of discordant alternans and govern its spatiotemporal properties are not well understood. Using mathematical modeling, we show that dynamic changes in the spatial distribution of discordant alternans can be used to pinpoint the underlying mechanisms. Optical mapping of V_m and Ca_i in paced rabbit hearts revealed that spatially discordant alternans induced by rapid pacing exhibits properties consistent with a purely dynamical mechanism as shown in theoretical studies. Our results support the viewpoint that spatially discordant alternans in the heart can be formed via a dynamical pattern formation process which does not require tissue heterogeneity.     

A Simulation Study of the Role of Mechanical Stretch in Arrhythmogenesis during Cardiac Alternans

The deformation of the heart tissue due to the contraction can modulate the excitation, a phenomenon referred to as mechanoelectrical feedback (MEF), via stretch-activated channels. The effects of MEF on the electrophysiology at high pacing rates are shown to be proarrhythmic in general. However, more studies need to be done to elucidate the underlying mechanism. In this work, we investigate the effects of MEF on cardiac alternans, which is an alternation in the width of the action potential that typically occurs when the heart is paced at high rates, using a biophysically detailed electromechanical model of cardiac tissue. We observe that the transition from spatially concordant alternans to spatially discordant alternans, which is more arrhythmogenic than concordant alternans, may occur in the presence of MEF and when its strength is sufficiently large. We show that this transition is due to the increase of the dispersion of conduction velocity. In addition, our results also show that the MEF effects, depending on the stretch-activated channels’ conductances and reversal potentials, can result in blocking action potential propagation.     

A rabbit ventricular action potential model replicating cardiac dynamics at rapid heart rates

Mathematical modeling of the cardiac action potential has proven to be a powerful tool for illuminating various aspects of cardiac function, including cardiac arrhythmias. However, no currently available detailed action potential model accurately reproduces the dynamics of the cardiac action potential and intracellular calcium (Ca_i) cycling at rapid heart rates relevant to ventricular tachycardia and fibrillation. The aim of this study was to develop such a model. Using an existing rabbit ventricular action potential model, we modified the L-type calcium (Ca) current (I_Ca,L) and Ca_i cycling formulations based on new experimental patch-clamp data obtained in isolated rabbit ventricular myocytes, using the perforated patch configuration at 35-37°C. Incorporating a minimal seven-state Markovian model of I_Ca,L that reproduced Ca- and voltage-dependent kinetics in combination with our previously published dynamic Ca_i cycling model, the new model replicates experimentally observed action potential duration and Ca_i transient alternans at rapid heart rates, and accurately reproduces experimental action potential duration restitution curves obtained by either dynamic or S1S2 pacing.     

Alternans Resonance and Propagation Block during Supernormal Conduction in Cardiac Tissue with Decreased [K]o

Cardiac restitution is an important factor in arrhythmogenesis. Steep positive action potential duration and conduction velocity (CV) restitution slopes promote alternans and reentrant arrhythmias. We examined the consequences of supernormal conduction (characterized by a negative CV restitution slope) on patterns of conduction and alternans in strands of Luo-Rudy model cells and in cultured cardiac cell strands. Interbeat intervals (IBIs) were analyzed as a function of distance during S1S2 protocols and during pacing at alternating cycle lengths. Supernormal conduction was induced by decreasing [K⁺]_o. In control [K⁺]_o simulations, S1S2 intervals converged toward basic cycle length with a length constant determined by both CV and the CV restitution slope. During alternant pacing, the amplitude of IBI alternans converged with a shorter length constant, determined also by the action potential duration restitution slope. In contrast, during supernormal conduction, S1S2 intervals and the amplitude of alternans diverged. This amplification (resonance) led to phase-locked or more complex alternans patterns, and then to distal conduction block. The convergence/divergence of IBIs was verified in the cultured strands, in which naturally occurring tissue heterogeneities resulted in prominent discontinuities of the spatial IBI profiles. We conclude that supernormal conduction potentiates alternans and spatial analysis of IBIs represents a powerful method to locate tissue heterogeneities.     

Characterizing the contribution of voltage- and calcium-dependent coupling to action potential stability: implications for repolarization alternans

Experiments have provided suggestive but inconclusive insights into the relative contributions of transmembrane voltage and intracellular calcium handling to the development of cardiac electrical instabilities such as repolarization alternans. In this study, we applied a novel combination of techniques (action potential voltage clamping, calcium-transient clamping, and stability analysis) to cardiac cell models to more clearly determine the roles that voltage- and calcium-dependent coupling play in regulating action potential stability and the development of alternans subsequent to the loss of stability. Using these techniques, we are able to demonstrate that voltage- and calcium-dependent coupling exhibit varying degrees of influence on action potential stability across models. Our results indicate that cellular dynamic instabilities such as alternans may be initiated by either voltage- or calcium-dependent mechanisms or by some combination of the two. Based on these modeling results, we propose novel single-cell experiments that incorporate action-potential voltage clamping, calcium imaging, and real-time measurement of action potential stability. These experiments will make it possible to experimentally determine the relative contribution of voltage coupling to the regulation of action potential stability in real cardiac myocytes, thereby providing further insights into the mechanism of alternans.     

Circadian Rhythms of Early Afterdepolarizations and Ventricular Arrhythmias in a Cardiomyocyte Model

Sudden cardiac arrest is a malfunction of the heart’s electrical system, typically caused by ventricular arrhythmias, that can lead to sudden cardiac death (SCD) within minutes. Epidemiological studies have shown that SCD and ventricular arrhythmias are more likely to occur in the morning than in the evening, and laboratory studies indicate that these daily rhythms in adverse cardiovascular events are at least partially under the control of the endogenous circadian timekeeping system. However, the biophysical mechanisms linking molecular circadian clocks to cardiac arrhythmogenesis are not fully understood. Recent experiments have shown that L-type calcium channels exhibit circadian rhythms in both expression and function in guinea pig ventricular cardiomyocytes. We developed an electrophysiological model of these cells to simulate the effect of circadian variation in L-type calcium conductance. In our simulations, we found that there is a circadian pattern in the occurrence of early afterdepolarizations (EADs), which are abnormal depolarizations during the repolarization phase of a cardiac action potential that can trigger fatal ventricular arrhythmias. Specifically, the model produces EADs in the morning, but not at other times of day. We show that the model exhibits a codimension-2 Takens-Bogdanov bifurcation that serves as an organizing center for different types of EAD dynamics. We also simulated a two-dimensional spatial version of this model across a circadian cycle. We found that there is a circadian pattern in the breakup of spiral waves, which represents ventricular fibrillation in cardiac tissue. Specifically, the model produces spiral wave breakup in the morning, but not in the evening. Our computational study is the first, to our knowledge, to propose a link between circadian rhythms and EAD formation and suggests that the efficacy of drugs targeting EAD-mediated arrhythmias may depend on the time of day that they are administered.     

Spatially Discordant Repolarization Alternans in the Absence of Conduction Velocity Restitution

Spatially discordant alternans (SDA) of action potential duration (APD) has been widely observed in cardiac tissue and is linked to cardiac arrhythmogenesis. Theoretical studies have shown that conduction velocity restitution (CVR) is required for the formation of SDA. However, this theory is not completely supported by experiments, indicating that other mechanisms may exist. In this study, we carried out computer simulations using mathematical models of action potentials to investigate the mechanisms of SDA in cardiac tissue. We show that when CVR is present and engaged, such as fast pacing from one side of the tissue, the spatial pattern of APD in the tissue undergoes either spatially concordant alternans or SDA, independent of initial conditions or tissue heterogeneities. When CVR is not engaged, such as simultaneous pacing of the whole tissue or under normal/slow heart rates, the spatial pattern of APD in the tissue can have multiple solutions, including spatially concordant alternans and different SDA patterns, depending on heterogeneous initial conditions or pre-existing repolarization heterogeneities. In homogeneous tissue, curved nodal lines are not stable, which either evolve into straight lines or disappear. However, in heterogeneous itssue, curved nodal lines can be stable, depending on their initial locations and shapes relative to the structure of the heterogeneity. Therefore, CVR-induced SDA and non-CVR-induced SDA exhibit different dynamical properties, which may be responsible for the different SDA properties observed in experimental studies and arrhythmogenesis in different clinical settings.