We numerically investigate the role of mechanical stress in modifying the conductivity properties of cardiac tissue, and also assess the impact of these effects in the solutions generated by ...computational models for cardiac electromechanics. We follow the recent theoretical framework from Cherubini et al. (2017), proposed in the context of general reaction-diffusion-mechanics systems emerging from multiphysics continuum mechanics and finite elasticity. In the present study, the adapted models are compared against preliminary experimental data of pig right ventricle fluorescence optical mapping. These data contribute to the characterization of the observed inhomogeneity and anisotropy properties that result from mechanical deformation. Our novel approach simultaneously incorporates two mechanisms for mechano-electric feedback (MEF): stretch-activated currents (SAC) and stress-assisted diffusion (SAD); and we also identify their influence into the nonlinear spatiotemporal dynamics. It is found that (i) only specific combinations of the two MEF effects allow proper conduction velocity measurement; (ii) expected heterogeneities and anisotropies are obtained via the novel stress-assisted diffusion mechanisms; (iii) spiral wave meandering and drifting is highly mediated by the applied mechanical loading. We provide an analysis of the intrinsic structure of the nonlinear coupling mechanisms using computational tests conducted with finite element methods. In particular, we compare static and dynamic deformation regimes in the onset of cardiac arrhythmias and address other potential biomedical applications.
A wide range of atrial arrythmias are caused by molecular defects in proteins that regulate calcium (Ca) cycling. In many cases, these defects promote the propagation of subcellular Ca waves in the ...cell, which can perturb the voltage time course and induce dangerous perturbations of the action potential (AP). However, subcellular Ca waves occur randomly in cells and, therefore, electrical coupling between cells substantially decreases their effect on the AP. In this study, we present evidence that Ca waves in atrial tissue can synchronize in-phase owing to an order-disorder phase transition. In particular, we show that, below a critical pacing rate, Ca waves are desynchronized and therefore do not induce substantial AP fluctuations in tissue. However, above this critical pacing rate, Ca waves gradually synchronize over millions of cells, which leads to a dramatic amplification of AP fluctuations. We exploit an underlying Ising symmetry of paced cardiac tissue to show that this transition exhibits universal properties common to a wide range of physical systems in nature. Finally, we show that in the heart, phase synchronization induces spatially out-of-phase AP duration alternans which drives wave break and reentry. These results suggest that cardiac tissue exhibits a phase transition that is required for subcellular Ca cycling defects to induce a life-threatening arrhythmia.
Most cardiac arrhythmias at the whole heart level result from alteration of cell membrane ionic channels and intracellular calcium concentration (Ca
) cycling with emerging spatiotemporal behavior ...through tissue-level coupling. For example, dynamically induced spatial dispersion of action potential duration, QT prolongation, and alternans are clinical markers for arrhythmia susceptibility in regular and heart-failure patients that originate due to changes of the transmembrane voltage (
) and Ca
. We present an optical-mapping methodology that permits simultaneous measurements of the
- Ca
signals using a single-camera without cross-talk, allowing quantitative characterization of favorable/adverse cell and tissue dynamical effects occurring from remodeling and/or drugs in heart failure. We demonstrate theoretically and experimentally in six different species the existence of a family of excitation wavelengths, we termed semasbestic, that give no change in signal for one dye, and thus can be used to record signals from another dye, guaranteeing zero cross-talk.
Customization of cardiac action potential models has become increasingly important with the recognition of patient-specific models and virtual patient cohorts as valuable predictive tools. ...Nevertheless, developing customized models by fitting parameters to data poses technical and methodological challenges: despite noise and variability associated with real-world datasets, traditional optimization methods produce a single “best-fit” set of parameter values. Bayesian estimation methods seek distributions of parameter values given the data by obtaining samples from the target distribution, but in practice widely known Bayesian algorithms like Markov chain Monte Carlo tend to be computationally inefficient and scale poorly with the dimensionality of parameter space. In this paper, we consider two computationally efficient Bayesian approaches: the Hamiltonian Monte Carlo (HMC) algorithm and the approximate Bayesian computation sequential Monte Carlo (ABC-SMC) algorithm. We find that both methods successfully identify distributions of model parameters for two cardiac action potential models using model-derived synthetic data and an experimental dataset from a zebrafish heart. Although both methods appear to converge to the same distribution family and are computationally efficient, HMC generally finds narrower marginal distributions, while ABC-SMC is less sensitive to the algorithmic settings including the prior distribution.
Graphical abstract
The monophasic action potential (MAP) is a near replica of the transmembrane potential recorded when an electrode is pushed firmly against cardiac tissue. Despite its many practical uses, the ...mechanism of MAP signal generation and the reason it is so different from unipolar recordings are not completely known and are a matter of controversy. In this work, we describe a method to simulate realistic MAP and intermediate forms, which are multiphasic electrograms different from an ideal MAP. The key ideas of our method are the formation of compressed zones and junctional spaces—regions of the extracellular and bath or blood pool directly in contact with electrodes that exhibit a pressure-induced reduction in electrical conductivity—and the presence of a complex network of passive components that acts as a high-pass filter to distort and attenuate the signal that reaches the recording amplifier. The network is formed by the interaction between the passive tissue properties and the double-layer capacitance of electrodes. The MAP and intermediate forms reside on a continuum of signals, which can be generated by the change of the model parameters. Our model helps to decipher the mechanisms of signal generation and can lead to a better design for electrodes, recording amplifiers, and experimental setups.
Finding appropriate values for parameters in mathematical models of cardiac cells is a challenging task. Here, we show that it is possible to obtain good parameterizations in as little as 30-40 s ...when as many as 27 parameters are fit simultaneously using a genetic algorithm and two flexible phenomenological models of cardiac action potentials. We demonstrate how our implementation works by considering cases of "model recovery" in which we attempt to find parameter values that match model-derived action potential data from several cycle lengths. We assess performance by evaluating the parameter values obtained, action potentials at fit and non-fit cycle lengths, and bifurcation plots for fidelity to the truth as well as consistency across different runs of the algorithm. We also fit the models to action potentials recorded experimentally using microelectrodes and analyze performance. We find that our implementation can efficiently obtain model parameterizations that are in good agreement with the dynamics exhibited by the underlying systems that are included in the fitting process. However, the parameter values obtained in good parameterizations can exhibit a significant amount of variability, raising issues of parameter identifiability and sensitivity. Along similar lines, we also find that the two models differ in terms of the ease of obtaining parameterizations that reproduce model dynamics accurately, most likely reflecting different levels of parameter identifiability for the two models.
Electrically based therapies for terminating atrial fibrillation (AF) currently fall into 2 categories: antitachycardia pacing and cardioversion. Antitachycardia pacing uses low-intensity pacing ...stimuli delivered via a single electrode and is effective for terminating slower tachycardias but is less effective for treating AF. In contrast, cardioversion uses a single high-voltage shock to terminate AF reliably, but the voltages required produce undesirable side effects, including tissue damage and pain. We propose a new method to terminate AF called far-field antifibrillation pacing, which delivers a short train of low-intensity electric pulses at the frequency of antitachycardia pacing but from field electrodes. Prior theoretical work has suggested that this approach can create a large number of activation sites ("virtual" electrodes) that emit propagating waves within the tissue without implanting physical electrodes and thereby may be more effective than point-source stimulation.
Using optical mapping in isolated perfused canine atrial preparations, we show that a series of pulses at low field strength (0.9 to 1.4 V/cm) is sufficient to entrain and subsequently extinguish AF with a success rate of 93% (69 of 74 trials in 8 preparations). We further demonstrate that the mechanism behind far-field antifibrillation pacing success is the generation of wave emission sites within the tissue by the applied electric field, which entrains the tissue as the field is pulsed.
AF in our model can be terminated by far-field antifibrillation pacing with only 13% of the energy required for cardioversion. Further studies are needed to determine whether this marked reduction in energy can increase the effectiveness and safety of terminating atrial tachyarrhythmias clinically.
Cardiac dynamics modeling has been useful for studying and treating arrhythmias. However, it is a multiscale problem requiring the solution of billions of differential equations describing the ...complex electrophysiology of interconnected cells. Therefore, large-scale cardiac modeling has been limited to groups with access to supercomputers and clusters. Many areas of computational science face similar problems where computational costs are too high for personal computers so that supercomputers or clusters currently are necessary. Here, we introduce a new approach that makes high-performance simulation of cardiac dynamics and other large-scale systems like fluid flow and crystal growth accessible to virtually anyone with a modest computer. For cardiac dynamics, this approach will allow not only scientists and students but also physicians to use physiologically accurate modeling and simulation tools that are interactive in real time, thereby making diagnostics, research, and education available to a broader audience and pushing the boundaries of cardiac science.