The reconstruction, from a set of points, and analysis of the interior surface of the heart chamber

Adequate description of heart muscle electrical activity is essential for the proper treatment of cardiac arrhythmias. Contemporary mapping and ablating systems allow a physician to introduce an electrode (catheter) into the human heart, to measure the position of the electrode in space and, simultaneously, the electrical activity timing and the bipolar and unipolar signal amplitudes--which correspond to the electrical viability of the heart muscle. If enough data points are collected, an approximate reconstruction of the heart chamber geometry (anatomy) is possible using also surface data such as the viability and local activity isochrones. Myocardial viability in patients after myocardial infarction is crucial for understanding and treating life threatening arrhythmias. Although there are commercial tools for heart chamber reconstruction, they lack the ability to quantitatively analyse the reconstructed data. Here, we show a method of reconstruction of the left ventricle of the heart from a measured set of data points and perform an interpolation of the measured voltages over the reconstructed surface. Next, we detect regions with voltage in a specified range and compute their areas and circumferences. Our methods allowed us to quantitatively describe the 'normal' muscle, the damaged or scar areas and the border zones between healthy muscle and the scars. In particular, we are able to find geometries of the damaged muscle areas that may be dangerous, e.g. when two such areas lie close to each other creating an isthmus--a macroreentry arrhythmia substrate. This work was inspired by a clinical hypothesis that the size of the border zone corresponds to the rate of occurrence of ventricular arrhythmia in patients after myocardial infarction.

Nonlinear oscillator model reproducing various phenomena in the dynamics of the conduction system of the heart

A dedicated nonlinear oscillator model able to reproduce the pulse shape, refractory time, and phase sensitivity of the action potential of a natural pacemaker of the heart is developed. The phase space of the oscillator contains a stable node, a hyperbolic saddle, and an unstable focus. The model reproduces several phenomena well known in cardiology, such as certain properties of the sinus rhythm and heart block. In particular, the model reproduces the decrease of heart rate variability with an increase in sympathetic activity. A sinus pause occurs in the model due to a single, well-timed, external pulse just as it occurs in the heart, for example due to a single supraventricular ectopy. Several ways by which the oscillations cease in the system are obtained (models of the asystole). The model simulates properly the way vagal activity modulates the heart rate and reproduces the vagal paradox. Two such oscillators, coupled unidirectionally and asymmetrically, allow us to reproduce the properties of heart rate variability obtained from patients with different kinds of heart block including sino-atrial blocks of different degree and a complete AV block (third degree). Finally, we demonstrate the possibility of introducing into the model a spatial dimension that creates exciting possibilities of simulating in the future the SA the AV nodes and the atrium including their true anatomical structure.

Influence of periprocedural anticoagulation strategies on complication rate and hospital stay in patients undergoing catheter ablation for persistent atrial fibrillation

BACKGROUND

The use of non-vitamin K antagonists (NOACs), uninterrupted (uVKA) and interrupted vitamin K antagonists (iVKA) are common periprocedural oral anticoagulation (OAC) strategies for atrial fibrillation (AF) ablation. Comparative data on complication rates resulting from OAC strategies for solely persistent AF (persAF) undergoing ablation are sparse. Thus, we sought to determine the impact of these OAC strategies on complication rates among patients with persAF undergoing catheter ablation.

METHODS

Consecutive patients undergoing persAF ablation were included. Depending on preprocedural OAC, three groups were defined: (1) NOACs (paused 48 h preablation), (2) uVKA, and (3) iVKA with heparin bridging. A combined complication endpoint (CCE) composed of bleeding and thromboembolic events was analyzed.

RESULTS

Between 2011 and 2014, 1440 persAF ablation procedures were performed in 1092 patients. NOACs were given in 441 procedures (31 %; rivaroxaban 57 %, dabigatran 33 %, and apixaban 10 %), uVKA in 488 (34 %), and iVKA in 511 (35 %). Adjusted CCE rates were 5.5 % [95 % confidence interval (CI) (3.1-7.8)] in group 1 (NOACs), 7.5 % [95 % CI (5.0-10.1)] in group 2 (uVKA), and 9.9 % [95 % CI (6.6-13.2)] in group 3. Compared to group 1, the combined complication risk was almost twice as high in group 3 [odd's ratio (OR) 1.9, 95 % CI (1.0-3.7), p = 0.049)]. The major complication rate was low (0.9 %). Bleeding complications, driven by minor groin complications, are more frequent than thromboembolic events (n = 112 vs. 1, p < 0.0001).

CONCLUSIONS

Patients undergoing persAF ablation with iVKA anticoagulation have an increased risk of complications compared to NOACs. Major complications, such as thromboembolic events, are generally rare and are exceeded by minor bleedings.

Spiral wave breakup in excitable media with an inhomogeneity of conduction anisotropy

Many conditions remodel the heart muscle such that it results in a perturbation of cells coupling. The effect of this perturbation on the stability of the spiral waves of electrochemical activity is not clear. We used the FitzHugh-Nagumo model of an excitable medium to model the conduction of the activation waves in a two-dimensional system with inhomogeneous anisotropy level. Inhomogeneity of the anisotropy level was modeled by adding Gaussian noise to diffusion coefficients corresponding with lateral coupling of the cells. Low noise levels resulted in a stable propagation of the spiral wave. For large noise level conduction was not possible due to insufficient coupling in direction perpendicular to fibers. For intermediate noise intensities, the initial wave broke up into several independent spiral waves or waves circulating around conduction obstacles. At an optimal noise intensity, the number of wavelets was maximized-a form of anti-coherent resonance was obtained. Our results suggest that the inhomogeneity of conduction anisotropy may promote wave breakup and hence play an important role in the initiation and perpetuation of the cardiac arrhythmias.

Reentry wave formation in excitable media with stochastically generated inhomogeneities

Clinical research shows that the frequency of arrhythmia events depends on the number and area of the border zones of infarct scars. We investigate the possibility that arrhythmia is initiated by reentry waves generated by the inhomogeneity of conduction velocity at the border zone. The interaction of a plane wave with a spatially extended inhomogeneity is simulated in the FitzHugh- Nagumo model. The inhomogeneity is introduced into the model by modifying the spatial dependence of the diffusion coefficient in a stochastic manner. This results in a rich variety of spatial distributions of conductivity. A plane wave propagating through such a system may break up on the regions with low conductivity and produce numerous spiral waves. The frequency of reentry wave formation is studied as a function of the parameters of the inhomogeneity generation algorithm. Three main scenarios of reentry wave formation were found: unidirectional block, main wave-wavelet collision, and wave break up during collision, on a region in which a conduction velocity gradient occurs. These scenarios are likely candidates for the mechanisms of arrhythmia initiation in a damaged tissue, e.g., the border zone of an infarct scar.

Quantitative description of the 3D regional mechanics of the left atrium using cardiac magnetic resonance imaging

The left atrium (LA) plays an important role in the maintenance of hemodynamic and electrical stability of the heart. One of the conditions altering the atrial mechanical function is atrial fibrillation (AF), leading to an increased thromboembolic risk due to impaired mechanical function. Preserving the regions of the LA that contribute the greatest to atrial mechanical function during curative strategies for AF is important. The purpose of this study is to introduce a novel method of regional assessment of mechanical function of the LA. We used cardiac MRI to reconstruct the 3D geometry of the LA in nine control and nine patients with paroxysmal atrial fibrillation (PAF). Regional mechanical function of the LA in pre-defined segments of the atrium was calculated using regional ejection fraction and wall velocity. We found significantly greater mechanical function in anterior, septal and lateral segments as opposed to roof and posterior segments, as well as a significant decrease of mechanical function in the PAF group. We suggest that in order to minimize the impact of the AF treatment on global atrial mechanical function, damage related to therapeutic intervention, such as catheter ablation, in those areas should be minimized.

Complex activity patterns in arterial wall: results from a model of calcium dynamics

Using a dynamical model of smooth muscle cells in an arterial wall, defined as a system of coupled five-dimensional nonlinear oscillators, on a grid with cylindrical symmetry, we compare the admissible activity patterns with those known from the heart tissue. We postulate on numerical basis the possibility to induce a stable spiral wave in the arterial wall. Such a spiral wave can inhibit the propagation of the axial calcium wave and effectively stop the vasomotion. We also discuss the dynamics of the circumferential calcium wave in comparison to rotors in venous ostia that are a common source of supraventricular ectopy. We show that the velocity and in consequence the frequency range of the circumferential calcium wave is by orders of magnitude too small compared to that of the rotors. The mechanism of the rotor is not likely to involve the calcium-related dynamics of the smooth muscle cells. The calcium-related dynamics which is voltage-independent and hard to be reset seems to actually protect the blood vessels against the electric activity of the atria. We also discuss the microreentry phenomenon, which was found in numerical experiments in the studied model.