Nonlinear physics of electrical wave propagation in the heart: a review

The beating of the heart is a synchronized contraction of muscle cells (myocytes) that are triggered by a periodic sequence of electrical waves (action potentials) originating in the sino-atrial node and propagating over the atria and the ventricles. Cardiac arrhythmias like atrial and ventricular fibrillation (AF,VF) or ventricular tachycardia (VT) are caused by disruptions and instabilities of these electrical excitations, that lead to the emergence of rotating waves (VT) and turbulent wave patterns (AF,VF). Numerous simulation and experimental studies during the last 20 years have addressed these topics. In this review we focus on the nonlinear dynamics of wave propagation in the heart with an emphasis on the theory of pulses, spirals and scroll waves and their instabilities in excitable media and their application to cardiac modeling. After an introduction into electrophysiological models for action potential propagation, the modeling and analysis of spatiotemporal alternans, spiral and scroll meandering, spiral breakup and scroll wave instabilities like negative line tension and sproing are reviewed in depth and discussed with emphasis on their impact in cardiac arrhythmias.Peer ReviewedPreprin

Dynamics of cardiac re-entry in micro-CT and serial histological sections based models of mammalian hearts

Cardiac re-entry regime of self-organised abnormal synchronisation underlie dangerous arrhythmias and fatal fibrillation. Recent advances in the theory of dissipative vortices, experimental studies, and anatomically realistic computer simulations, elucidated the role of cardiac re-entry interaction with fine anatomical features in the heart, and anatomy induced drift. The fact that anatomy and structural anisotropy of the heart is consistent within a species suggested its possible functional effect on spontaneous drift of cardiac re-entry. A comparative study of the anatomy induced drift could be used in order to predict evolution of atrial arrhythmia, and improve low-voltage defibrillation protocols and ablation strategies. Here, in micro-CT based model of rat pulmonary vein wall, and in sheep atria models based on high resolution serial histological sections, we demonstrate effects of heart geometry and anisotropy on cardiac re-entry anatomy induced drift, its pinning to fluctuations of thickness in the layer. The data sets of sheep atria and rat pulmonary vein wall are incorporated into the BeatBox High Performance Computing simulation environment. Re-entry is initiated at prescribed locations in the spatially homogeneous mono-domain models of cardiac tissue. Excitation is described by FitzHugh-Nagumo kinetics. In the in-silico models, isotropic and anisotropic conduction show specific anatomy effects and the interplay between anatomy and anisotropy of the heart. The main objectives are to demonstrate the functional role of the species hearts geometry and anisotropy on cardiac re-entry anatomy induced drift. In case of the rat pulmonary vein wall with ~90 degree transmural fibre rotation, it is shown that the joint effect of the PV wall geometry and anisotropy turns a plane excitation wave into a re-entry pinned to a small fluctuation of thickness in the wall

Optogenetic Control of Cardiac Arrhythmias

The regular, coordinated contraction of the heart muscle is orchestrated by periodic waves generated by the heart’s natural pacemaker and transmitted through the heart’s electrical conduction system. Abnormalities occurring anywhere within the cardiac electrical conduction system can disrupt the propagation of these waves. Such dis- ruptions often lead to the development of high frequency spiral waves that override normal pacemaker activity and compromise cardiac function. The occurrence of high frequency spiral waves in the heart is associated with cardiac rhythm disorders such as tachycardia and fibrillation. While tachycardia may be terminated by rapid periodic stimulation known as anti-tachycardia pacing (ATP), life-threatening ventricular fibril- lation requires a single high-voltage electric shock that resets all the activity and restore the normal heart function. However, despite the high success rate of defibrillation, it is associated with significant side effects including tissue damage, intense pain and trauma. Thus, extensive research is conducted for developing low-energy alternatives to conventional defibrillation. An example of such an alternative is the low-energy anti-fibrillation pacing (LEAP). However, the clinical application of this technique, and other evolving techniques requires a detailed understanding of the dynamics of spiral waves that occur during arrhythmias. Optogenetics is a tool, that has recently gained popularity in the cardiac research, which serves as a probe to study biological processes. It involves genetically modifying cardiac muscle cells such that they become light sensitive, and then using light of specific wavelengths to control the electrical activity of these cells. Cardiac optogenetics opens up new ways of investigating the mechanisms underlying the onset, maintenance and control of cardiac arrhythmias. In this thesis, I employ optogenetics as a tool to control the dynamics of a spiral wave, in both computer simulations and in experiments.In the first study, I use optogenetics to investigate the mechanisms underlying de- fibrillation. Analogous to the conventional single electric-shock, I apply a single globally-illuminating light pulse to a two-dimensional cardiac tissue to study how wave termination occurs during defibrillation. My studies show a characteristic transient dynamics leading to the termination of the spiral wave at low light intensities, while at high intensities, the spiral waves terminate immediately. Next, I move on to explore the use of optogenetics to study spiral wave termina- tion via drift, theoretically well-known mechanism of arrhythmia termination in the context of electrical stimulation (e.g. ATP). I show that spiral wave drift can be induced by structured illumination patterns using lights of low intensity, that result in a spatial modulation of cardiac excitability. I observe that drift occurs in the positive direction of light intensity gradient, where the spiral also rotates with a longer period. I further show how modulation of the excitability in space can be used to control the dynamics of a spiral wave, resulting in the termination of the wave by collision with the domain boundary. Based on these observations, I propose a possible mechanism of optogenetic defibrillation. In the next chapter, I use optogenetics to demonstrate control over the dynamics of the spiral waves by periodic stimulation with light of different intensities and pacing frequencies resulting in a temporal modulation of cardiac excitability. I demonstrate how the temporal modulation of excitability leads to efficient termination of arrhythmia. In addition, I use computer simulations to identify mechanisms responsible for arrhyth- mia termination for sub- and supra-threshold light intensities. My numerical results are supported by experimental studies on intact hearts, extracted from transgenic mice. Finally, I demonstrate that cardiac optogenetics not only allows control of excita- tion waves, but also by generating new waves through the induction of wave breaks. We demonstrate the effects of high sub-threshold illumination on the morphology of the propagating wave, leading to the creation of new excitation windows in space that can serve as potential sites for re-entry initiation. In summary, this thesis investigates several approaches to control arrhythmia dy- namics using optogenetics. The experimental and numerical results demonstrate the potential of feedback-induced resonant pacing as a low-energy method to control arrhythmia.2022-01-1

Cardiac re-entry dynamics and self-termination in DT-MRI based model of Human Foetal Heart

The effect of human fetal heart geometry and anisotropy on anatomy induced drift and self-termination of cardiac re-entry is studied here in MRI based 2D slice and 3D whole heart computer simulations. Isotropic and anisotropic models of 20 weeks of gestational age human fetal heart obtained from 100 μm voxel diffusion tensor MRI data sets were used in the computer simulations. The fiber orientation angles of the heart were obtained from the orientation of the DT-MRI primary eigenvectors. In a spatially homogeneous electrophysiological monodomain model with the DT-MRI based heart geometries, cardiac re-entry was initiated at a prescribed location in a 2D slice, and in the 3D whole heart anatomy models. Excitation was described by simplified FitzHugh-Nagumo kinetics. In a slice of the heart, with propagation velocity twice as fast along the fibers than across the fibers, DT-MRI based fiber anisotropy changes the re-entry dynamics from pinned to an anatomical re-entry. In the 3D whole heart models, the fiber anisotropy changes cardiac re-entry dynamics from a persistent re-entry to the re-entry self-termination. The self-termination time depends on the re-entry's initial position. In all the simulations with the DT-MRI based cardiac geometry, the anisotropy of the myocardial tissue shortens the time to re-entry self-termination several folds. The numerical simulations depend on the validity of the DT-MRI data set used. The ventricular wall showed the characteristic transmural rotation of the helix angle of the developed mammalian heart, while the fiber orientation in the atria was irregula

BeatBox - HPC simulation environment for biophysically and anatomically realistic cardiac electrophysiology

The BeatBox simulation environment combines flexible script language user interface with the robust computational tools, in order to setup cardiac electrophysiology in-silico experiments without re-coding at low-level, so that cell excitation, tissue/anatomy models, stimulation protocols may be included into a BeatBox script, and simulation run either sequentially or in parallel (MPI) without re-compilation. BeatBox is a free software written in C language to be run on a Unix-based platform. It provides the whole spectrum of multi scale tissue modelling from 0-dimensional individual cell simulation, 1-dimensional fibre, 2-dimensional sheet and 3-dimensional slab of tissue, up to anatomically realistic whole heart simulations, with run time measurements including cardiac re-entry tip/filament tracing, ECG, local/global samples of any variables, etc. BeatBox solvers, cell, and tissue/anatomy models repositories are extended via robust and flexible interfaces, thus providing an open framework for new developments in the field. In this paper we give an overview of the BeatBox current state, together with a description of the main computational methods and MPI parallelisation approaches.Comment: 37 pages, 10 figures, last version submitted to PLOS ON

A Biomagnetic Field Mapping System for Detection of Heart Disease in a Clinical Environment

This PhD was inspired by the clinical demand for a system to triage chest pain and the untapped diagnostic potential of magnetocardiography (MCG) to reliably detect silent ischaemic heart disease, which is responsible for the highest mortality rate of any single disease category. The aim was to develop a low cost and portable biomagnetic field mapping system capable of differentiating between healthy and diseased hearts within an unshielded hospital environment. This entailed the development of a system based on an array of magnetometers with sufficient sensitivity (104fT/ Hz at 10Hz) and noise rejection performance (68.4 ± 3.9 dB) to measure the small magnetic field associated with the heart beat, the magneocardiogram, within a much larger background noise. The array of induction coil magnetometers (ICM) we developed had sufficient sensitivity and were robust to high amplitude noise. These sensors were also cheap to manufacture and capable of operating on battery power, allowing a low cost, portable device to be developed. The key element that allowed us to achieve unshielded operation was the development of an algorithmic spatial filter, used as a substitute to operation within a magnetically shielded room. This coherent noise rejection (CNR) algorithm exploits the difference in spatial coherence between the local cardiac signals and the distant background noise sources. The observed coherence width during a clinical trial of the system within a hospital ward was 2.8 ± 0.9 × 10 6 mm 2 . This allowed us to capture MCG signals with a signal to noise ratio of SNR QRS = 0.93 ± 4.43dB. The performance of CNR was found to improve by 9dB per order of magnitude increase in environmental spatial coherence width. The coherence width can be increased by changes to hospital architecture, electromagnetic field regulation and device design optimisation. The thesis also explores a variety of approaches to obtain binary diagnostic information from MCG, from traditional statistical learning on manually engineered features to machine learning. I found that machine learning techniques, in particular convolutional neural networks (CNN), were able to capture more diagnostic information than traditional techniques and achieved world class prediction accuracy of 88% on the clinical trial dataset