Multiparametric analysis of geometric features of fibrotic textures leading to cardiac arrhythmias

Re-entry in models of cardiac ventricular tissue with scar represented as a Gaussian random field

Introduction: Fibrotic scar in the heart is known to act as a substrate for arrhythmias. Regions of fibrotic scar are associated with slowed or blocked conduction of the action potential, but the detailed mechanisms of arrhythmia formation are not well characterised and this can limit the effective diagnosis and treatment of scar in patients. The aim of this computational study was to evaluate different representations of fibrotic scar in models of 2D 10 × 10 cm ventricular tissue, where the region of scar was defined by sampling a Gaussian random field with an adjustable length scale of between 1.25 and 10.0 mm. Methods: Cellular electrophysiology was represented by the Ten Tusscher 2006 model for human ventricular cells. Fibrotic scar was represented as a spatially varying diffusion, with different models of the boundary between normal and fibrotic tissue. Dispersion of activation time and action potential duration (APD) dispersion was assessed in each sample by pacing at an S1 cycle length of 400 ms followed by a premature S2 beat with a coupling interval of 323 ms. Vulnerability to reentry was assessed with an aggressive pacing protocol. In all models, simulated fibrosis acted to delay activation, to increase the dispersion of APD, and to generate re-entry. Results: A higher incidence of re-entry was observed in models with simulated fibrotic scar at shorter length scale, but the type of model used to represent fibrotic scar had a much bigger influence on the incidence of reentry. Discussion: This study shows that in computational models of fibrotic scar the effects that lead to either block or propagation of the action potential are strongly influenced by the way that fibrotic scar is represented in the model, and so the results of computational studies involving fibrotic scar should be interpreted carefully.

Systematic in-silico evaluation of fibrosis effects on re-entrant wave dynamics in atrial tissue

Despite the key role of fibrosis in atrial fibrillation (AF), the effects of different spatial distributions and textures of fibrosis on wave propagation mechanisms in AF are not fully understood. To clarify these aspects, we performed a systematic computational study to assess fibrosis effects on the characteristics and stability of re-entrant waves in electrically-remodelled atrial tissues. A stochastic algorithm, which generated fibrotic distributions with controlled overall amount, average size, and orientation of fibrosis elements, was implemented on a monolayer spheric atrial model. 245 simulations were run at changing fibrosis parameters. The emerging propagation patterns were quantified in terms of rate, regularity, and coupling by frequency-domain analysis of correspondent synthetic bipolar electrograms. At the increase of fibrosis amount, the rate of reentrant waves significantly decreased and higher levels of regularity and coupling were observed (p < 0.0001). Higher spatial variability and pattern stochasticity over repetitions was observed for larger amount of fibrosis, especially in the presence of patchy and compact fibrosis. Overall, propagation slowing and organization led to higher stability of re-entrant waves. These results strengthen the evidence that the amount and spatial distribution of fibrosis concur in dictating re-entry dynamics in remodeled tissue and represent key factors in AF maintenance.

Computer based method for identification of fibrotic scars from electrograms and local activation times on the epi- and endocardial surfaces of the ventricles

Cardiac fibrosis stands as one of the most critical conditions leading to lethal cardiac arrhythmias. Identifying the precise location of cardiac fibrosis is crucial for planning clinical interventions in patients with various forms of ventricular and atrial arrhythmias. As fibrosis impedes and alters the path of electrical waves, detecting fibrosis in the heart can be achieved through analyzing electrical signals recorded from its surface. In current clinical practices, it has become feasible to record electrical activity from both the endocardial and epicardial surfaces of the heart. This paper presents a computational method for reconstructing 3D fibrosis using unipolar electrograms obtained from both surfaces of the ventricles. The proposed method calculates the percentage of fibrosis in various ventricular segments by analyzing the local activation times and peak-to-peak amplitudes of the electrograms. Initially, the method was tested using simulated data representing idealized fibrosis in a heart segment; subsequently, it was validated in the left ventricle with fibrosis obtained from a patient with nonischemic cardiomyopathy. The method successfully determined the location and extent of fibrosis in 204 segments of the left ventricle model with an average error of 0.0±4.3% (N = 204). Moreover, the method effectively detected fibrotic scars in the mid-myocardial region, a region known to present challenges in accurate detection using electrogram amplitude as the primary criterion.

In silico analysis of the contribution of cardiomyocyte-fibroblast electromechanical interaction to the arrhythmia

Although fibroblasts are about 5–10 times smaller than cardiomyocytes, their number in the ventricle is about twice that of cardiomyocytes. The high density of fibroblasts in myocardial tissue leads to a noticeable effect of their electromechanical interaction with cardiomyocytes on the electrical and mechanical functions of the latter. Our work focuses on the analysis of the mechanisms of spontaneous electrical and mechanical activity of the fibroblast-coupled cardiomyocyte during its calcium overload, which occurs in a variety of pathologies, including acute ischemia. For this study, we developed a mathematical model of the electromechanical interaction between cardiomyocyte and fibroblasts and used it to simulate the impact of overloading cardiomyocytes. In contrast to modeling only the electrical interaction between cardiomyocyte and fibroblasts, the following new features emerge in simulations with the model that accounts for both electrical and mechanical coupling and mechano-electrical feedback loops in the interacting cells. First, the activity of mechanosensitive ion channels in the coupled fibroblasts depolarizes their resting potential. Second, this additional depolarization increases the resting potential of the coupled myocyte, thus augmenting its susceptibility to triggered activity. The triggered activity associated with the cardiomyocyte calcium overload manifests itself in the model either as early afterdepolarizations or as extrasystoles, i.e., extra action potentials and extra contractions. Analysis of the model simulations showed that mechanics contribute significantly to the proarrhythmic effects in the cardiomyocyte overloaded with calcium and coupled with fibroblasts, and that mechano-electrical feedback loops in both the cardiomyocyte and fibroblasts play a key role in this phenomenon.

Role of Scar and Border Zone Geometry on the Genesis and Maintenance of Re-Entrant Ventricular Tachycardia in Patients With Previous Myocardial Infarction

In patients with healed myocardial infarction, the left ventricular ejection fraction is characterized by low sensitivity and specificity in the prediction of future malignant arrhythmias. Thus, there is the need for new parameters in daily practice to perform arrhythmic risk stratification. The aim of this study is to identify some features of proarrhythmic geometric configurations of scars and border zones (BZ), by means of numerical simulations based on left ventricular models derived from post myocardial infarction patients. Two patients with similar clinical characteristics were included in this study. Both patients exhibited left ventricular scars characterized by subendo- and subepicardial BZ and a transmural BZ isthmus. The scar of patient #1 was significantly larger than that of patient #2, whereas the transmural BZ isthmus and the subdendo- and subepicardial BZs of patient #2 were thicker than those of patient #1. Patient #1 was positive at electrophysiologic testing, whereas patient #2 was negative. Based on the cardiac magnetic resonance (CMR) data, we developed a geometric model of the left ventricles of the two patients, taking into account the position, extent, and topological features of scars and BZ. The numerical simulations were based on the anisotropic monodomain model of electrocardiology. In the model of patient #1, sustained ventricular tachycardia (VT) was inducible by an S2 stimulus delivered at any of the six stimulation sites considered, while in the model of patient #2 we were not able to induce sustained VT. In the model of patient #1, making the subendo- and subepicardial BZs as thick as those of patient #2 did not affect the inducibility and maintenance of VT. On the other hand, in the model of patient #2, making the subendo- and subepicardial BZs as thin as those of patient #1 yielded sustained VT. In conclusion, the results show that the numerical simulations have an effective predictive capability in discriminating patients at high arrhythmic risk. The extent of the infarct scar and the presence of transmural BZ isthmuses and thin subendo- and subepicardial BZs promote sustained VT.