title: Period doubling cascades of limit cycles in cardiac action potential models as precursors to chaotic early Afterdepolarizations.
abstract: Early afterdepolarizations (EADs) are pathological voltage oscillations during the repolarization phase of cardiac action potentials (APs). EADs are caused by drugs, oxidative stress or ion channel disease, and they are considered as potential precursors to cardiac arrhythmias in recent attempts to redefine the cardiac drug safety paradigm. The irregular behaviour of EADs observed in experiments has been previously attributed to chaotic EAD dynamics under periodic pacing, made possible by a homoclinic bifurcation in the fast subsystem of the deterministic AP system of differential equations.
  
  
  
  In this article we demonstrate that a homoclinic bifurcation in the fast subsystem of the action potential model is neither a necessary nor a sufficient condition for the genesis of chaotic EADs. We rather argue that a cascade of period doubling (PD) bifurcations of limit cycles in the full AP system paves the way to chaotic EAD dynamics across a variety of models including a) periodically paced and spontaneously active cardiomyocytes, b) periodically paced and non-active cardiomyocytes as well as c) unpaced and spontaneously active cardiomyocytes. Furthermore, our bifurcation analysis reveals that chaotic EAD dynamics may coexist in a stable manner with fully regular AP dynamics, where only the initial conditions decide which type of dynamics is displayed.
  
  
  
  EADs are a potential source of cardiac arrhythmias and hence are of relevance both from the viewpoint of drug cardiotoxicity testing and the treatment of cardiomyopathies. The model-independent association of chaotic EADs with period doubling cascades of limit cycles introduced in this article opens novel opportunities to study chaotic EADs by means of bifurcation control theory and inverse bifurcation analysis. Furthermore, our results may shed new light on the synchronization and propagation of chaotic EADs in homogeneous and heterogeneous multicellular and cardiac tissue preparations.
authors.affiliation: Institute of Applied Mathematics and Statistics, University of Hohenheim, Schloss 1, Stuttgart, 70599, Germany. philipp.kuegler@uni-hohenheim.de.|Radon Institute for Computational and Applied Mathematics, Austrian Academy of Sciences, Altenbergerstrasse 69, Linz, 4040, Austria. philipp.kuegler@uni-hohenheim.de.
authors.forename: Philipp
authors.initials: P
authors.lastname: Kügler
authors.affiliation: Radon Institute for Computational and Applied Mathematics, Austrian Academy of Sciences, Altenbergerstrasse 69, Linz, 4040, Austria.
authors.forename: M A K
authors.initials: MA
authors.lastname: Bulelzai
authors.affiliation: Institute of Applied Mathematics and Statistics, University of Hohenheim, Schloss 1, Stuttgart, 70599, Germany.
authors.forename: André H
authors.initials: AH
authors.lastname: Erhardt
mesh_terms: D000200:Action Potentials; D000328:Adult; D000818:Animals; D001683:Biological Clocks; D006801:Humans; D008955:Models, Cardiovascular; D032383:Myocytes, Cardiac / Q000166:cytology*
keywords: Bifurcation theory; Cardiac action potential; Chaos; Early afterdepolarizations; Nonlinear dynamics