Early Afterdepolarisations Induced by an Enhancement in the Calcium Current

Bifurcations and Proarrhythmic Behaviors in Cardiac Electrical Excitations

The heart is a hierarchical dynamic system consisting of molecules, cells, and tissues, and acts as a pump for blood circulation. The pumping function depends critically on the preceding electrical activity, and disturbances in the pattern of excitation propagation lead to cardiac arrhythmia and pump failure. Excitation phenomena in cardiomyocytes have been modeled as a nonlinear dynamical system. Because of the nonlinearity of excitation phenomena, the system dynamics could be complex, and various analyses have been performed to understand the complex dynamics. Understanding the mechanisms underlying proarrhythmic responses in the heart is crucial for developing new ways to prevent and control cardiac arrhythmias and resulting contractile dysfunction. When the heart changes to a pathological state over time, the action potential (AP) in cardiomyocytes may also change to a different state in shape and duration, often undergoing a qualitative change in behavior. Such a dynamic change is called bifurcation. In this review, we first summarize the contribution of ion channels and transporters to AP formation and our knowledge of ion-transport molecules, then briefly describe bifurcation theory for nonlinear dynamical systems, and finally detail its recent progress, focusing on the research that attempts to understand the developing mechanisms of abnormal excitations in cardiomyocytes from the perspective of bifurcation phenomena.

Dynamical analysis of early afterdepolarization patterns in a biophysically detailed cardiac model

Arrhythmogenic early afterdepolarizations (EADs) are investigated in a biophysically detailed mathematical model of a rabbit ventricular myocyte, providing their location in the parameter phase space and describing their dynamical mechanisms. Simulations using the Sato model, defined by 27 state variables and 177 parameters, are conducted to generate electrical action potentials (APs) for different values of the pacing cycle length and other parameters related to sodium and calcium concentrations. A detailed study of the different AP patterns with or without EADs is carried out, showing the presence of a high variety of temporal AP configurations with chaotic and quasiperiodic behaviors. Regions of bistability are identified and, importantly, linked to transitions between different behaviors. Using sweeping techniques, one-, two-, and three-parameter phase spaces are provided, allowing ascertainment of the role of the selected parameters as well as location of the transition regions. A Devil's staircase, with symbolic sequence analysis, is proposed to describe transitions in the ratio between the number of voltage (EAD and AP) peaks and the number of APs. To conclude, the obtained results are linked to recent studies for low-dimensional models and a conjecture is made for the internal dynamical structure of the transition region from non-EAD to EAD behavior using fold and cusp bifurcations and maximal canards.

Dynamics of a neuron-glia system: the occurrence of seizures and the influence of electroconvulsive stimuli : A mathematical and numerical study

In this paper, we investigate the dynamics of a neuron-glia cell system and the underlying mechanism for the occurrence of seizures. For our mathematical and numerical investigation of the cell model we will use bifurcation analysis and some computational methods. It turns out that an increase of the potassium concentration in the reservoir is one trigger for seizures and is related to a torus bifurcation. In addition, we will study potassium dynamics of the model by considering a reduced version and we will show how both mechanisms are linked to each other. Moreover, the reduction of the potassium leak current will also induce seizures. Our study will show that an enhancement of the extracellular potassium concentration, which influences the Nernst potential of the potassium current, may lead to seizures. Furthermore, we will show that an external forcing term (e.g. electroshocks as unidirectional rectangular pulses also known as electroconvulsive therapy) will establish seizures similar to the unforced system with the increased extracellular potassium concentration. To this end, we describe the unidirectional rectangular pulses as an autonomous system of ordinary differential equations. These approaches will explain the appearance of seizures in the cellular model. Moreover, seizures, as they are measured by electroencephalography (EEG), spread on the macro-scale (cm). Therefore, we extend the cell model with a suitable homogenised monodomain model, propose a set of (numerical) experiment to complement the bifurcation analysis performed on the single-cell model. Based on these experiments, we introduce a bidomain model for a more realistic modelling of white and grey matter of the brain. Performing similar (numerical) experiment as for the monodomain model leads to a suitable comparison of both models. The individual cell model, with its seizures explained in terms of a torus bifurcation, extends directly to corresponding results in both the monodomain and bidomain models where the neural firing spreads almost synchronous through the domain as fast traveling waves, for physiologically relevant paramenters.

Special Issue: Methods in Computational Biology

Biological systems are multiscale with respect to time and space, exist at the interface of biological and physical constraints, and their interactions with the environment are often nonlinear [...]