Search

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., (© 2022. International Federation for Medical and Biological Engineering.)

Purpose: Multi-polar diagnostic catheters are used to construct the 3D electro-anatomic mapping of the atrium during atrial fibrillation (AF) ablation procedures; however, it remains unclear how to use the electrograms recorded by these catheters to locate AF-driving sites known as focal and rotor source types. The purpose of this study is to present the first algorithm to iteratively navigate a circular multi-polar catheter to locate AF focal and rotor sources without the need to map the entire atria., Methods: Starting from an initial location, the algorithm, which was blinded to the location and type of the AF source, iteratively advanced a Lasso catheter based on its electrogram characteristics. The algorithm stopped the catheter when it located of an AF source and identified the type. The efficiency of the algorithm is validated using a set of simulated focal and rotor-driven arrhythmias in fibrotic human 2D and 3D atrial tissue., Results: Our study shows the feasibility of locating AF sources with a success rate of greater than 95.25% within average 7.56 ± 2.28 placements independently of the initial position of the catheter and the source type., Conclusions: The algorithm could play a critical role in clinical electrophysiology laboratories for mapping patient-specific ablation of AF sources located outside the pulmonary veins and improving the procedure success.

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.

Atrial fibrillation (AF) is an irregular heart rhythm that can lead to stroke and other heart-related complications. Catheter ablation has been commonly used to destroy triggering sources of AF in the atria and consequently terminate the arrhythmia. However, efficient and accurate localization of the AF sustaining sources known as rotors is a major challenge in catheter ablation. In this paper, we developed a novel probabilistic algorithm that can adaptively guide a Lasso diagnostic catheter to locate the center of a rotor. Our algorithm uses a Bayesian updating approach to search for and locate rotors based on the characteristics of electrogram signals collected at every catheter placement. The algorithm was evaluated using a 10 × 10 cm 2 D atrial tissue simulation of the Nygren human atrial cell model and was able to successfully guide the catheter to the rotor center in 3.37 ± 1.05 (mean±std) steps (including placement at the center) when starting from any location on the tissue. Our novel automated algorithm can potentially play a significant role in patient-specific ablation of AF sources and increase the success of AF elimination procedures.

Locating atrial fibrillation (AF) focal and rotor sources can help improve target ablation therapy for AF. However, it remains unclear how to use the information provided by multi-polar diagnostic catheters (MPDC) to locate AF sources. Our aim was to develop a catheter-guidance algorithm to locate AF focal and rotor sources using a conventional MPDC. We simulated a 10 cm × 10 cm atrial tissue with focal and rotor sources using the Nygren et al. ionic model. We modeled a Lasso MPDC with 20-unipole electrodes placed with a spacing of 4.5-1-4.5 mm (diameter, d=15 mm) along a circle to obtain 10-bipole electrograms. Staring from an initial location, the algorithm, which was blinded to the location and type of the AF source, iteratively advanced the MPDC by moving its center to the location of the first activated bipole (FAB). The algorithm located an AF source if a stopping condition for either source was satisfied using bipole electrogram characteristics extracted from the MPDC placement. We tested the algorithm for a single rotor and focal source for all possible initial positions on the simulated tissue and repeated it for a random placement with a maximum of 20 possible placements at every trial. The algorithm located the AF source for 100% of trials and on average required 5.99 ± 1.92 placements to an AF source. This algorithm may be used to iteratively direct an MPDC towards an AF source and allow the AF source to be localized for customized AF ablation.

We present the Spiral Classification Algorithm (SCA), a fast and accurate algorithm for classifying electrical spiral waves and their associated breakup in cardiac tissues. The classification performed by SCA is an essential component of the detection and analysis of various cardiac arrhythmic disorders, including ventricular tachycardia and fibrillation. Given a digitized frame of a propagating wave, SCA constructs a highly accurate representation of the front and the back of the wave, piecewise interpolates this representation with cubic splines, and subjects the result to an accurate curvature analysis. This analysis is more comprehensive than methods based on spiral-tip tracking, as it considers the entire wave front and back. To increase the smoothness of the resulting symbolic representation, the SCA uses weighted overlapping of adjacent segments which increases the smoothness at join points. SCA has been applied to a number of representative types of spiral waves, and, for each type, a distinct curvature evolution in time (signature) has been identified. Distinct signatures have also been identified for spiral breakup. These results represent a significant first step in automatically determining parameter ranges for which a computational cardiac-cell network accurately reproduces a particular kind of cardiac arrhythmia, such as ventricular fibrillation.