A discrete model to study reaction-diffusion-mechanics systems

Inverse mechano-electrical reconstruction of cardiac excitation wave patterns from mechanical deformation using deep learning

The inverse mechano-electrical problem in cardiac electrophysiology is the attempt to reconstruct electrical excitation or action potential wave patterns from the heart's mechanical deformation that occurs in response to electrical excitation. Because heart muscle cells contract upon electrical excitation due to the excitation-contraction coupling mechanism, the resulting deformation of the heart should reflect macroscopic action potential wave phenomena. However, whether the relationship between macroscopic electrical and mechanical phenomena is well-defined and unique enough to be utilized for an inverse imaging technique in which mechanical activation mapping is used as a surrogate for electrical mapping has yet to be determined. Here, we provide a numerical proof-of-principle that deep learning can be used to solve the inverse mechano-electrical problem in phenomenological two- and three-dimensional computer simulations of the contracting heart wall, or in elastic excitable media, with muscle fiber anisotropy. We trained a convolutional autoencoder neural network to learn the complex relationship between electrical excitation, active stress, and tissue deformation during both focal or reentrant chaotic wave activity and, consequently, used the network to successfully estimate or reconstruct electrical excitation wave patterns from mechanical deformation in sheets and bulk-shaped tissues, even in the presence of noise and at low spatial resolutions. We demonstrate that even complicated three-dimensional electrical excitation wave phenomena, such as scroll waves and their vortex filaments, can be computed with very high reconstruction accuracies of about 95% from mechanical deformation using autoencoder neural networks, and we provide a comparison with results that were obtained previously with a physics- or knowledge-based approach.

Discrete mechanical growth model for plant tissue

We present a discrete mechanical model to study plant development. The method is built up of mass points, springs and hinges mimicking the plant cell wall’s microstructure. To model plastic growth the resting lengths of springs are adjusted; when springs exceed a threshold length, new mass points, springs and hinges, are added. We formulate a stiffness tensor for the springs and hinges as a function of the fourth rank tensor of elasticity and the geometry of the mesh. This allows us to approximate the material law as a generalized orthotropic Hooke’s law, and control material properties during growth. The material properties of the model are illustrated in numerical simulations for finite strain and plastic growth. To solve the equations of motion of mass points we assume elastostatics and use Verlet integration. The method is demonstrated in simulations when anisotropic growth causes emergent residual strain fields in cell walls and a bending of tissue. The method can be used in multilevel models to study plant development, for example by coupling it to models for cytoskeletal, hormonal and gene regulatory processes.

Marker-Free Tracking for Motion Artifact Compensation and Deformation Measurements in Optical Mapping Videos of Contracting Hearts

Optical mapping is a high-resolution fluorescence imaging technique, which provides highly detailed visualizations of the electrophysiological wave phenomena, which trigger the beating of the heart. Recent advancements in optical mapping have demonstrated that the technique can now be performed with moving and contracting hearts and that motion and motion artifacts, once a major limitation, can now be overcome by numerically tracking and stabilizing the heart's motion. As a result, the optical measurement of electrical activity can be obtained from the moving heart surface in a co-moving frame of reference and motion artifacts can be reduced substantially. The aim of this study is to assess and validate the performance of a 2D marker-free motion tracking algorithm, which tracks motion and non-rigid deformations in video images. Because the tracking algorithm does not require markers to be attached to the tissue, it is necessary to verify that it accurately tracks the displacements of the cardiac tissue surface, which not only contracts and deforms, but also fluoresces and exhibits spatio-temporal physiology-related intensity changes. We used computer simulations to generate synthetic optical mapping videos, which show the contracting and fluorescing ventricular heart surface. The synthetic data reproduces experimental data as closely as possible and shows electrical waves propagating across the deforming tissue surface, as seen during voltage-sensitive imaging. We then tested the motion tracking and motion-stabilization algorithm on the synthetic as well as on experimental data. The motion tracking and motion-stabilization algorithm decreases motion artifacts approximately by 80% and achieves sub-pixel precision when tracking motion of 1–10 pixels (in a video image with 100 by 100 pixels), effectively inhibiting motion such that little residual motion remains after tracking and motion-stabilization. To demonstrate the performance of the algorithm, we present optical maps with a substantial reduction in motion artifacts showing action potential waves propagating across the moving and strongly deforming ventricular heart surface. The tracking algorithm reliably tracks motion if the tissue surface is illuminated homogeneously and shows sufficient contrast or texture which can be tracked or if the contrast is artificially or numerically enhanced. In this study, we also show how a reduction in dissociation-related motion artifacts can be quantified and linked to tracking precision. Our results can be used to advance optical mapping techniques, enabling them to image contracting hearts, with the ultimate goal of studying the mutual coupling of electrical and mechanical phenomena in healthy and diseased hearts.

Electromechanical vortex filaments during cardiac fibrillation

The self-organized dynamics of vortex-like rotating waves, which are also known as scroll waves, are the basis of the formation of complex spatiotemporal patterns in many excitable chemical and biological systems. In the heart, filament-like phase singularities that are associated with three-dimensional scroll waves are considered to be the organizing centres of life-threatening cardiac arrhythmias. The mechanisms that underlie the onset, maintenance and control of electromechanical turbulence in the heart are inherently three-dimensional phenomena. However, it has not previously been possible to visualize the three-dimensional spatiotemporal dynamics of scroll waves inside cardiac tissues. Here we show that three-dimensional mechanical scroll waves and filament-like phase singularities can be observed deep inside the contracting heart wall using high-resolution four-dimensional ultrasound-based strain imaging. We found that mechanical phase singularities co-exist with electrical phase singularities during cardiac fibrillation. We investigated the dynamics of electrical and mechanical phase singularities by simultaneously measuring the membrane potential, intracellular calcium concentration and mechanical contractions of the heart. We show that cardiac fibrillation can be characterized using the three-dimensional spatiotemporal dynamics of mechanical phase singularities, which arise inside the fibrillating contracting ventricular wall. We demonstrate that electrical and mechanical phase singularities show complex interactions and we characterize their dynamics in terms of trajectories, topological charge and lifetime. We anticipate that our findings will provide novel perspectives for non-invasive diagnostic imaging and therapeutic applications.

Mechanism for Mechanical Wave Break in the Heart Muscle

Using a reaction-diffusion-mechanics model we identify a mechanism for mechanical wave break in the heart muscle. For a wide range of strengths and durations an external mechanical load causes wave front dissipation leading to formation and breakup of spiral waves. We explain the mechanism, and discuss under which conditions it can cause or abolish cardiac arrhythmias.

Cardiac contraction induces discordant alternans and localized block

In this paper we use a simplified model of cardiac excitation-contraction coupling to study the effect of tissue deformation on the dynamics of alternans, i.e., alternations in the duration of the cardiac action potential, that occur at fast pacing rates and are known to be proarrhythmic. We show that small stretch-activated currents can produce large effects and cause a transition from in-phase to off-phase alternations (i.e., from concordant to discordant alternans) and to conduction blocks. We demonstrate numerically and analytically that this effect is the result of a generic change in the slope of the conduction velocity restitution curve due to electromechanical coupling. Thus, excitation-contraction coupling can potentially play a relevant role in the transition to reentry and fibrillation.

Higher dimensional Gaussian-type solitons of nonlinear Schrödinger equation with cubic and power-law nonlinearities in PT-symmetric potentials

Two families of Gaussian-type soliton solutions of the (n+1)-dimensional Schrödinger equation with cubic and power-law nonlinearities in -symmetric potentials are analytically derived. As an example, we discuss some dynamical behaviors of two dimensional soliton solutions. Their phase switches, powers and transverse power-flow densities are discussed. Results imply that the powers flow and exchange from the gain toward the loss regions in the cell. Moreover, the linear stability analysis and the direct numerical simulation are carried out, which indicates that spatial Gaussian-type soliton solutions are stable below some thresholds for the imaginary part of -symmetric potentials in the defocusing cubic and focusing power-law nonlinear medium, while they are always unstable for all parameters in other media.

Three-dimensional structures of the spatiotemporal nonlinear Schrödinger equation with power-law nonlinearity in PT-symmetric potentials

The spatiotemporal nonlinear Schrödinger equation with power-law nonlinearity in PT-symmetric potentials is investigated, and two families of analytical three-dimensional spatiotemporal structure solutions are obtained. The stability of these solutions is tested by the linear stability analysis and the direct numerical simulation. Results indicate that solutions are stable below some thresholds for the imaginary part of PT-symmetric potentials in the self-focusing medium, while they are always unstable for all parameters in the self-defocusing medium. Moreover, some dynamical properties of these solutions are discussed, such as the phase switch, power and transverse power-flow density. The span of phase switch gradually enlarges with the decrease of the competing parameter k in PT-symmetric potentials. The power and power-flow density are all positive, which implies that the power flow and exchange from the gain toward the loss domains in the PT cell.

Spiral-wave dynamics in ionically realistic mathematical models for human ventricular tissue: the effects of periodic deformation

We carry out an extensive numerical study of the dynamics of spiral waves of electrical activation, in the presence of periodic deformation (PD) in two-dimensional simulation domains, in the biophysically realistic mathematical models of human ventricular tissue due to (a) ten-Tusscher and Panfilov (the TP06 model) and (b) ten-Tusscher, Noble, Noble, and Panfilov (the TNNP04 model). We first consider simulations in cable-type domains, in which we calculate the conduction velocity θ and the wavelength λ of a plane wave; we show that PD leads to a periodic, spatial modulation of θ and a temporally periodic modulation of λ; both these modulations depend on the amplitude and frequency of the PD. We then examine three types of initial conditions for both TP06 and TNNP04 models and show that the imposition of PD leads to a rich variety of spatiotemporal patterns in the transmembrane potential including states with a single rotating spiral (RS) wave, a spiral-turbulence (ST) state with a single meandering spiral, an ST state with multiple broken spirals, and a state SA in which all spirals are absorbed at the boundaries of our simulation domain. We find, for both TP06 and TNNP04 models, that spiral-wave dynamics depends sensitively on the amplitude and frequency of PD and the initial condition. We examine how these different types of spiral-wave states can be eliminated in the presence of PD by the application of low-amplitude pulses by square- and rectangular-mesh suppression techniques. We suggest specific experiments that can test the results of our simulations.

Spiral-wave dynamics in a mathematical model of human ventricular tissue with myocytes and fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as , the fibroblast resting-membrane potential, the fibroblast conductance , and the MF gap-junctional coupling . Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as , and , and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity decreases as a function of , for zero-sided and one-sided couplings; however, for two-sided coupling, decreases initially and then increases as a function of , and, eventually, we observe that conduction failure occurs for low values of . In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling or . Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

Generalized and improved (G'/G)-expansion method for (3+1)-dimensional modified KdV-Zakharov-Kuznetsev equation

The generalized and improved -expansion method is a powerful and advantageous mathematical tool for establishing abundant new traveling wave solutions of nonlinear partial differential equations. In this article, we investigate the higher dimensional nonlinear evolution equation, namely, the (3+1)-dimensional modified KdV-Zakharov-Kuznetsev equation via this powerful method. The solutions are found in hyperbolic, trigonometric and rational function form involving more parameters and some of our constructed solutions are identical with results obtained by other authors if certain parameters take special values and some are new. The numerical results described in the figures were obtained with the aid of commercial software Maple.

A discrete electromechanical model for human cardiac tissue: effects of stretch-activated currents and stretch conditions on restitution properties and spiral wave dynamics

We introduce an electromechanical model for human cardiac tissue which couples a biophysical model of cardiac excitation (Tusscher, Noble, Noble, Panfilov, 2006) and tension development (adjusted Niederer, Hunter, Smith, 2006 model) with a discrete elastic mass-lattice model. The equations for the excitation processes are solved with a finite difference approach, and the equations of the mass-lattice model are solved using Verlet integration. This allows the coupled problem to be solved with high numerical resolution. Passive mechanical properties of the mass-lattice model are described by a generalized Hooke's law for finite deformations (Seth material). Active mechanical contraction is initiated by changes of the intracellular calcium concentration, which is a variable of the electrical model. Mechanical deformation feeds back on the electrophysiology via stretch-activated ion channels whose conductivity is controlled by the local stretch of the medium. We apply the model to study how stretch-activated currents affect the action potential shape, restitution properties, and dynamics of spiral waves, under constant stretch, and dynamic stretch caused by active mechanical contraction. We find that stretch conditions substantially affect these properties via stretch-activated currents. In constantly stretched medium, we observe a substantial decrease in conduction velocity, and an increase of action potential duration; whereas, with dynamic stretch, action potential duration is increased only slightly, and the conduction velocity restitution curve becomes biphasic. Moreover, in constantly stretched medium, we find an increase of the core size and period of a spiral wave, but no change in rotation dynamics; in contrast, in the dynamically stretching medium, we observe spiral drift. Our results may be important to understand how altered stretch conditions affect the heart's functioning.

Emergence of spiral wave activity in a mechanically heterogeneous reaction-diffusion-mechanics system

We perform a numerical study of emergent spiral wave activity in a two-dimensional reaction-diffusion-mechanics medium with a regional inhomogeneity in active and passive mechanical properties. We find that self-sustaining spiral wave activity emerges for a wide range of mechanical parameters of the inhomogeneity via five mechanisms. We classify these mechanisms, relate them to parameters of the inhomogeneity, and discuss how these results can be applied to understand the onset of cardiac arrhythmias due to regional mechanical heterogeneity.

New mechanism of spiral wave initiation in a reaction-diffusion-mechanics system

Spiral wave initiation in the heart muscle is a mechanism for the onset of dangerous cardiac arrhythmias. A standard protocol for spiral wave initiation is the application of a stimulus in the refractory tail of a propagating excitation wave, a region that we call the “classical vulnerable zone.” Previous studies of vulnerability to spiral wave initiation did not take the influence of deformation into account, which has been shown to have a substantial effect on the excitation process of cardiomyocytes via the mechano-electrical feedback phenomenon. In this work we study the effect of deformation on the vulnerability of excitable media in a discrete reaction-diffusion-mechanics (dRDM) model. The dRDM model combines FitzHugh-Nagumo type equations for cardiac excitation with a discrete mechanical description of a finite-elastic isotropic material (Seth material) to model cardiac excitation-contraction coupling and stretch activated depolarizing current. We show that deformation alters the “classical,” and forms a new vulnerable zone at longer coupling intervals. This mechanically caused vulnerable zone results in a new mechanism of spiral wave initiation, where unidirectional conduction block and rotation directions of the consequently initiated spiral waves are opposite compared to the mechanism of spiral wave initiation due to the “classical vulnerable zone.” We show that this new mechanism of spiral wave initiation can naturally occur in situations that involve wave fronts with curvature, and discuss its relation to supernormal excitability of cardiac tissue. The concept of mechanically induced vulnerability may lead to a better understanding about the onset of dangerous heart arrhythmias via mechano-electrical feedback.