Dissolution of spiral wave's core using cardiac optogenetics.

Rotating spiral waves in the heart are associated with life-threatening cardiac arrhythmias such as ventricular tachycardia and fibrillation. These arrhythmias are treated by a process called defibrillation, which forces electrical resynchronization of the heart tissue by delivering a single global high-voltage shock directly to the heart. This method leads to immediate termination of spiral waves. However, this may not be the only mechanism underlying successful defibrillation, as certain scenarios have also been reported, where the arrhythmia terminated slowly, over a finite period of time. Here, we investigate the slow termination dynamics of an arrhythmia in optogenetically modified murine cardiac tissue both in silico and ex vivo during global illumination at low light intensities. Optical imaging of an intact mouse heart during a ventricular arrhythmia shows slow termination of the arrhythmia, which is due to action potential prolongation observed during the last rotation of the wave. Our numerical studies show that when the core of a spiral is illuminated, it begins to expand, pushing the spiral arm towards the inexcitable boundary of the domain, leading to termination of the spiral wave. We believe that these fundamental findings lead to a better understanding of arrhythmia dynamics during slow termination, which in turn has implications for the improvement and development of new cardiac defibrillation techniques.

In silico optical modulation of spiral wave trajectories in cardiac tissue.

Life-threatening cardiac arrhythmias such as ventricular tachycardia and fibrillation are common precursors to sudden cardiac death. They are associated with the occurrence of abnormal electrical spiral waves in the heart that rotate at a high frequency. In severe cases, arrhythmias are combated with a clinical method called defibrillation, which involves administering a single global high-voltage shock to the heart to reset all its activity and restore sinus rhythm. Despite its high efficiency in controlling arrhythmias, defibrillation is associated with several negative side effects that render the method suboptimal. The best approach to optimize this therapeutic technique is to deepen our understanding of the dynamics of spiral waves. Here, we use computational cardiac optogenetics to study and control the dynamics of a single spiral wave in a two-dimensional, electrophysiologically detailed, light-sensitive model of a mouse ventricle. First, we illuminate the domain globally by applying a sequence of periodic optical pulses with different frequencies in the sub-threshold regime where no excitation wave is induced. In doing so, we obtain epicycloidal, hypocycloidal, and resonant drift trajectories of the spiral wave core. Then, to effectively control the wave dynamics, we use a method called resonant feedback pacing. In this approach, each global optical pulse is applied when the measuring electrode positioned on the domain registers a predefined value of the membrane voltage. This enables us to steer the spiral wave in a desired direction determined by the position of the electrode. Our study thus provides valuable mechanistic insights into the success or failure of global optical stimulation in executing efficient arrhythmia control.

Optimising low-energy defibrillation in 2D cardiac tissue with a genetic algorithm.

Sequences of low-energy electrical pulses can effectively terminate ventricular fibrillation (VF) and avoid the side effects of conventional high-energy electrical defibrillation shocks, including tissue damage, traumatic pain, and worsening of prognosis. However, the systematic optimisation of sequences of low-energy pulses remains a major challenge. Using 2D simulations of homogeneous cardiac tissue and a genetic algorithm, we demonstrate the optimisation of sequences with non-uniform pulse energies and time intervals between consecutive pulses for efficient VF termination. We further identify model-dependent reductions of total pacing energy ranging from ∼4% to ∼80% compared to reference adaptive-deceleration pacing (ADP) protocols of equal success rate (100%).

Optical mapping of contracting hearts.

Optical mapping is a widely used tool to record and visualize the electrophysiological properties in a variety of myocardial preparations such as Langendorff-perfused isolated hearts, coronary-perfused wedge preparations, and cell culture monolayers. Motion artifact originating from the mechanical contraction of the myocardium creates a significant challenge to performing optical mapping of contracting hearts. Hence, to minimize the motion artifact, cardiac optical mapping studies are mostly performed on non-contracting hearts, where the mechanical contraction is removed using pharmacological excitation-contraction uncouplers. However, such experimental preparations eliminate the possibility of electromechanical interaction, and effects such as mechano-electric feedback cannot be studied. Recent developments in computer vision algorithms and ratiometric techniques have opened the possibility of performing optical mapping studies on isolated contracting hearts. In this review, we discuss the existing techniques and challenges of optical mapping of contracting hearts.

Non-monotonous dose response function of the termination of spiral wave chaos.

The conventional termination technique of life threatening cardiac arrhythmia like ventricular fibrillation is the application of a high-energy electrical defibrillation shock, coming along with severe side-effects. In order to improve the current treatment reducing these side-effects, the application of pulse sequences of lower energy instead of a single high-energy pulse are promising candidates. In this study, we show that in numerical simulations the dose-response function of pulse sequences applied to two-dimensional spiral wave chaos is not necessarily monotonously increasing, but exhibits a non-trivial frequency dependence. This insight into crucial phenomena appearing during termination attempts provides a deeper understanding of the governing termination mechanisms in general, and therefore may open up the path towards an efficient termination of cardiac arrhythmia in the future.

A Mathematical Model for Electrical Activity in Pig Atrial Tissue.

State of the art mathematical models are currently used to bridge the gap between basic research conducted in the laboratory and preclinical research conducted on large animals, which ultimately paves the way for clinical translation. In this regard, there is a great need for models that can be used alongside experiments for in-depth investigation and validation. One such experimental model is the porcine atrium, which is commonly used to study the mechanisms of onset and control of atrial fibrillation in the context of its surgical management. However, a mathematical model of pig atria is lacking. In this paper, we present the first ionically detailed mathematical model of porcine atrial electrophysiology, at body temperature. The model includes 12 ionic currents, 4 of which were designed based on experimental patch-clamp data directly obtained from literature. The formulations for the other currents are adopted from the human atrial model, and modified for porcine specificity based on our measured restitution data for different action potential characteristics: resting membrane potential, action potential amplitude, maximum upstroke velocity and action potential duration and different levels of membrane voltage repolarization. The intracellular Ca 2+ dynamics follows the Luo-Rudy formulation for guinea pig ventricular cardiomyocytes. The resulting model represents "normal" cells which are formulated as a system of ordinary differential equations. We extend our model to two dimensions to obtain plane wave propagation in tissue with a velocity of 0.58 m/s and a wavelength of 8 cm. The wavelength reduces to 5 cm when the tissue is paced at 200 ms. Using S1-S2 cross-field protocol, we demonstrate in an 11.26 cm square simulation domain, the ability to initiate single spiral waves (rotation period ≃ 180 ms) that remain stable for more than 40 s. The spiral tip exhibits hypermeander. In agreement with previous experimental results using pig atria, our model shows that early repolarization is primarily driven by a calcium-mediated chloride current, I ClCa , which is completely inactivated at high pacing frequencies. This is a condition that occurs only in porcine atria. Furthermore, the model shows spatiotemporal chaos with reduced repolarization.

The role of pulse timing in cardiac defibrillation.

Life-threatening cardiac arrhythmias require immediate defibrillation. For state-of-the-art shock treatments, a high field strength is required to achieve a sufficient success rate for terminating the complex spiral wave (rotor) dynamics underlying cardiac fibrillation. However, such high energy shocks have many adverse side effects due to the large electric currents applied. In this study, we show, using 2D simulations based on the Fenton-Karma model, that also pulses of relatively low energy may terminate the chaotic activity if applied at the right moment in time. In our simplified model for defibrillation, complex spiral waves are terminated by local perturbations corresponding to conductance heterogeneities acting as virtual electrodes in the presence of an external electric field. We demonstrate that time series of the success rate for low energy shocks exhibit pronounced peaks which correspond to short intervals in time during which perturbations aiming at terminating the chaotic fibrillation state are (much) more successful. Thus, the low energy shock regime, although yielding very low temporal average success rates, exhibits moments in time for which success rates are significantly higher than the average value shown in dose-response curves. This feature might be exploited in future defibrillation protocols for achieving high termination success rates with low or medium pulse energies.

Taming cardiac arrhythmias: Terminating spiral wave chaos by adaptive deceleration pacing.

Sequences of weak electrical pulses are considered a promising alternative for terminating ventricular and atrial fibrillations while avoiding strong defibrillation shocks with adverse side effects. In this study, using numerical simulations of four different 2D excitable media, we show that pulse trains with increasing temporal intervals between successive pulses (deceleration pacing) provide high success rates at low energies. Furthermore, we propose a simple and robust approach to calculate inter-pulse spacing directly from the frequency spectrum of the dynamics (for instance, computed based on the electrocardiogram), which can be practically used in experiments and clinical applications.

Pulsed low-energy stimulation initiates electric turbulence in cardiac tissue.

Interruptions in nonlinear wave propagation, commonly referred to as wave breaks, are typical of many complex excitable systems. In the heart they lead to lethal rhythm disorders, the so-called arrhythmias, which are one of the main causes of sudden death in the industrialized world. Progress in the treatment and therapy of cardiac arrhythmias requires a detailed understanding of the triggers and dynamics of these wave breaks. In particular, two very important questions are: 1) What determines the potential of a wave break to initiate re-entry? and 2) How do these breaks evolve such that the system is able to maintain spatiotemporally chaotic electrical activity? Here we approach these questions numerically using optogenetics in an in silico model of human atrial tissue that has undergone chronic atrial fibrillation (cAF) remodelling. In the lesser studied sub-threshold illumination régime, we discover a new mechanism of wave break initiation in cardiac tissue that occurs for gentle slopes of the restitution characteristics. This mechanism involves the creation of conduction blocks through a combination of wavefront-waveback interaction, reshaping of the wave profile and heterogeneous recovery from the excitation of the spatially extended medium, leading to the creation of re-excitable windows for sustained re-entry. This finding is an important contribution to cardiac arrhythmia research as it identifies scenarios in which low-energy perturbations to cardiac rhythm can be potentially life-threatening.

Advanced Cardiac Rhythm Management by Applying Optogenetic Multi-Site Photostimulation in Murine Hearts.

Ventricular tachyarrhythmias are a major cause of mortality and morbidity worldwide. Electrical defibrillation using high-energy electric shocks is currently the only treatment for life-threatening ventricular fibrillation. However, defibrillation may have side-effects, including intolerable pain, tissue damage, and worsening of prognosis, indicating a significant medical need for the development of more gentle cardiac rhythm management strategies. Besides energy-reducing electrical approaches, cardiac optogenetics was introduced as a powerful tool to influence cardiac activity using light-sensitive membrane ion channels and light pulses. In the present study, a robust and valid method for successful photostimulation of Langendorff perfused intact murine hearts will be described based on multi-site pacing applying a 3 x 3 array of micro light-emitting diodes (micro-LED). Simultaneous optical mapping of epicardial membrane voltage waves allows the investigation of the effects of region-specific stimulation and evaluates the newly induced cardiac activity directly on-site. The obtained results show that the efficacy of defibrillation is strongly dependent on the parameters chosen for photostimulation during a cardiac arrhythmia. It will be demonstrated that the illuminated area of the heart plays a crucial role for termination success as well as how the targeted control of cardiac activity during illumination for modifying arrhythmia patterns can be achieved. In summary, this technique provides a possibility to optimize the on-site mechanism manipulation on the way to real-time feedback control of cardiac rhythm and, regarding the region specificity, new approaches in reducing the potential harm to the cardiac system compared to the usage of non-specific electrical shock applications.

Drift and termination of spiral waves in optogenetically modified cardiac tissue at sub-threshold illumination.

The development of new approaches to control cardiac arrhythmias requires a deep understanding of spiral wave dynamics. Optogenetics offers new possibilities for this. Preliminary experiments show that sub-threshold illumination affects electrical wave propagation in the mouse heart. However, a systematic exploration of these effects is technically challenging. Here, we use state-of-the-art computer models to study the dynamic control of spiral waves in a two-dimensional model of the adult mouse ventricle, using stationary and non-stationary patterns of sub-threshold illumination. Our results indicate a light-intensity-dependent increase in cellular resting membrane potentials, which together with diffusive cell-cell coupling leads to the development of spatial voltage gradients over differently illuminated areas. A spiral wave drifts along the positive gradient. These gradients can be strategically applied to ensure drift-induced termination of a spiral wave, both in optogenetics and in conventional methods of electrical defibrillation.

High-Resolution Optical Measurement of Cardiac Restitution, Contraction, and Fibrillation Dynamics in Beating vs. Blebbistatin-Uncoupled Isolated Rabbit Hearts.

Optical mapping is a high-resolution fluorescence imaging technique, that uses voltage- or calcium-sensitive dyes to visualize electrical excitation waves on the heart surface. However, optical mapping is very susceptible to the motion of cardiac tissue, which results in so-called motion artifacts in the fluorescence signal. To avoid motion artifacts, contractions of the heart muscle are typically suppressed using pharmacological excitation-contraction uncoupling agents, such as Blebbistatin. The use of pharmacological agents, however, may influence cardiac electrophysiology. Recently, it has been shown that numerical motion tracking can significantly reduce motion-related artifacts in optical mapping, enabling the simultaneous optical measurement of cardiac electrophysiology and mechanics. Here, we combine ratiometric optical mapping with numerical motion tracking to further enhance the robustness and accuracy of these measurements. We evaluate the method's performance by imaging and comparing cardiac restitution and ventricular fibrillation (VF) dynamics in contracting, non-working vs. Blebbistatin-arrested Langendorff-perfused rabbit hearts (N = 10). We found action potential durations (APD) to be, on average, 25 ± 5% shorter in contracting hearts compared to hearts uncoupled with Blebbistatin. The relative shortening of the APD was found to be larger at higher frequencies. VF was found to be significantly accelerated in contracting hearts, i.e., 9 ± 2Hz with Blebbistatin and 15 ± 4Hz without Blebbistatin, and maintained a broader frequency spectrum. In contracting hearts, the average number of phase singularities was N PS = 11 ± 4 compared to N PS = 6 ± 3 with Blebbistatin during VF on the anterior ventricular surface. VF inducibility was reduced with Blebbistatin. We found the effect of Blebbistatin to be concentration-dependent and reversible by washout. Aside from the electrophysiological characterization, we also measured and analyzed cardiac motion. Our findings may have implications for the interpretation of optical mapping data, and highlight that physiological conditions, such as oxygenation and metabolic demand, must be carefully considered in ex vivo imaging experiments.

Excitable dynamics in neural and cardiac systems.

The application of mathematics, physics and engineering to medical research is continuously growing; interactions among these disciplines have become increasingly important and have contributed to an improved understanding of clinical and biological phenomena, with implications for disease prevention, diagnosis and treatment. This special issue presents examples of this synergy, with a particular focus on the investigation of cardiac and neural excitability. This issue includes 24 original research papers and covers a broad range of topics related to the physiological and pathophysiological function of the brain and the heart. Studies span scales from isolated neurons and small networks of neurons to whole-organ dynamics for the brain and from cardiac subcellular domains and cardiomyocytes to one-dimensional tissues for the heart. This preface is part of the Special Issue on "Excitable Dynamics in Neural and Cardiac Systems".

Extracting Robust Biomarkers From Multichannel EEG Time Series Using Nonlinear Dimensionality Reduction Applied to Ordinal Pattern Statistics and Spectral Quantities.

In this study, ordinal pattern analysis and classical frequency-based EEG analysis methods are used to differentiate between EEGs of different age groups as well as individuals. As characteristic features, functional connectivity as well as single-channel measures in both the time and frequency domain are considered. We compare the separation power of each feature set after nonlinear dimensionality reduction using t-distributed stochastic neighbor embedding and demonstrate that ordinal pattern-based measures yield results comparable to frequency-based measures applied to preprocessed data, and outperform them if applied to raw data. Our analysis yields no significant differences in performance between single-channel features and functional connectivity features regarding the question of age group separation.

Spiral wave unpinning facilitated by wave emitting sites in cardiac monolayers.

Rotating spiral waves of electrical activity in the heart can anchor to unexcitable tissue (an obstacle) and become stable pinned waves. A pinned rotating wave can be unpinned either by a local electrical stimulus applied close to the spiral core, or by an electric field pulse that excites the core of a pinned wave independently of its localization. The wave will be unpinned only when the pulse is delivered inside a narrow time interval called the unpinning window (UW) of the spiral. In experiments with cardiac monolayers, we found that other obstacles situated near the pinning centre of the spiral can facilitate unpinning. In numerical simulations, we found increasing or decreasing of the UW depending on the location, orientation and distance between the pinning centre and an obstacle. Our study indicates that multiple obstacles could contribute to unpinning in experiments with intact hearts.

Synchronization of viscoelastically coupled excitable oscillators.

Viscoelastically coupled excitable oscillators are used to model individually beating spatially separated cardiomyocytes surrounded by an extracellular matrix (ECM). We investigate how mechanical coupling via the ECM can synchronize two such oscillators with excitation contraction coupling and electromechanical feedback and how this synchronization depends on the rheological properties of the ECM. Extending our study to a linear chain of coupled oscillators we find a transition to synchronization as the ECM becomes stiffer. In the case of purely elastic coupling we observe antiphase chimera states.

Multiscale Modeling of Dyadic Structure-Function Relation in Ventricular Cardiac Myocytes.

Cardiovascular disease is often related to defects of subcellular components in cardiac myocytes, specifically in the dyadic cleft, which include changes in cleft geometry and channel placement. Modeling of these pathological changes requires both spatially resolved cleft as well as whole cell level descriptions. We use a multiscale model to create dyadic structure-function relationships to explore the impact of molecular changes on whole cell electrophysiology and calcium cycling. This multiscale model incorporates stochastic simulation of individual L-type calcium channels and ryanodine receptor channels, spatially detailed concentration dynamics in dyadic clefts, rabbit membrane potential dynamics, and a system of partial differential equations for myoplasmic and lumenal free Ca2+ and Ca2+-binding molecules in the bulk of the cell. We found action potential duration, systolic, and diastolic [Ca2+] to respond most sensitively to changes in L-type calcium channel current. The ryanodine receptor channel cluster structure inside dyadic clefts was found to affect all biomarkers investigated. The shape of clusters observed in experiments by Jayasinghe et al. and channel density within the cluster (characterized by mean occupancy) showed the strongest correlation to the effects on biomarkers.

Termination of Scroll Waves by Surface Impacts.

Three-dimensional scroll waves direct cell movement and gene expression, and induce chaos in the brain and heart. We found an approach to terminate multiple three-dimensional scrolls. A pulse of a properly configured electric field detaches scroll filaments from the surface. They shrink due to filament tension and disappear. Since wave emission from small heterogeneities is not used, this approach requires a much lower electric field. It is not sensitive to the details of the excitable medium. It may affect future studies of low-energy chaos termination in the heart.

Spontaneous termination of chaotic spiral wave dynamics in human cardiac ion channel models.

Chaotic spiral or scroll wave dynamics can be found in diverse systems. In cardiac dynamics, spiral or scroll waves of electrical excitation determine the dynamics during life-threatening arrhythmias like ventricular fibrillation. In numerical studies it was found that chaotic episodes of spiral and scroll waves can be transient, thus they terminate spontaneously. We show in this study that this behavior can also be observed using models which describe the ion channel dynamics of human cardiomyocytes (Bueno-Orovio-Cherry-Fenton model and the Ten Tusscher-Noble-Noble-Panfilov model). For both models we find that the average lifetime of the chaotic transients grows exponentially with the system size. With this behavior, we classify the systems into the group of type-II supertransients. We observe a significant difference of the breakup behavior between the models, which results in a distinct dynamics during the final phase just before the termination. The observation of a (temporally) stable single-spiral state affects the prevailing description of the dynamics of type-II supertransients as being "quasi-stationary" and also the feasibility of predicting the spontaneous termination of the spiral wave dynamics. In the long term, the relation between the breakup behavior of spiral waves and properties of chaotic transients like predictability or average transient lifetime may contribute to an improved understanding and classification of cardiac arrhythmias.

Simultaneous unpinning of multiple vortices in two-dimensional excitable media.

There are many examples of excitable media, such as the heart, that can show complex dynamics and where control is a challenging task. Heavy means like a strong electric shock are nowadays still necessary to control and terminate ventricular fibrillation (VF). It is known that heterogeneities in an excitable medium can stabilize the activity, e.g., spiral waves can pin to such obstacles. This might also be a reason for the persistence of VF and the difficulty to control it. Previous studies investigated systems with a single pinned spiral wave and demonstrated how the spiral can be unpinned. In this article, we extend this case and investigate a generic excitable system with multiple pinned spiral waves. We describe a control technique that allows the simultaneous unpinning of pinned spiral waves. Apart from theoretical considerations, we provide numerical evidence that the proposed technique is superior to the underdrive pacing method that has reportedly high success rates when applied to a single pinned spiral.