Dynamics of lattice spins as a model of arrhythmia

We consider evolution of initial disturbances in spatially extended systems with autonomous rhythmic activity, such as the heart. We consider the case when the activity is stable with respect to very smooth (changing little across the medium) disturbances and construct lattice models for description of not-so-smooth disturbances, in particular, topological defects; these models are modifications of the diffusive XY model. We find that when the activity on each lattice site is very rigid in maintaining its form, the topological defects - vortices or spirals - nucleate a transition to a disordered, turbulent state.Comment: 17 pages, revtex, 3 figure

Quantifying the effect of uncertainty in input parameters in a simplified bidomain model of partial thickness ischaemia

Reduced blood flow in the coronary arteries can lead to damaged heart tissue (myocardial ischaemia). Although one method for detecting myocardial ischaemia involves changes in the ST segment of the electrocardiogram, the relationship between these changes and subendocardial ischaemia is not fully understood. In this study, we modelled ST-segment epicardial potentials in a slab model of cardiac ventricular tissue, with a central ischaemic region, using the bidomain model, which considers conduction longitudinal, transverse and normal to the cardiac fibres. We systematically quantified the effect of uncertainty on the input parameters, fibre rotation angle, ischaemic depth, blood conductivity and six bidomain conductivities, on outputs that characterise the epicardial potential distribution. We found that three typical types of epicardial potential distributions (one minimum over the central ischaemic region, a tripole of minima, and two minima flanking a central maximum) could all occur for a wide range of ischaemic depths. In addition, the positions of the minima were affected by both the fibre rotation angle and the ischaemic depth, but not by changes in the conductivity values. We also showed that the magnitude of ST depression is affected only by changes in the longitudinal and normal conductivities, but not by the transverse conductivities

Nanoscale regulation of L-type calcium channels differentiates between ischemic and dilated cardiomyopathies.

Background Subcellular localization and function of L-type calcium channels (LTCCs) play an important role in regulating contraction of cardiomyocytes. Understanding how this is affected by the disruption of transverse tubules during heart failure could lead to new insights into the disease. Methods Cardiomyocytes were isolated from healthy donor hearts, as well as from patients with cardiomyopathies and with left ventricular assist devices. Scanning ion conductance and confocal microscopy was used to study membrane structures in the cells. Super-resolution scanning patch-clamp was used to examine LTCC function in different microdomains. Computational modeling predicted the impact of these changes to arrhythmogenesis at the whole-heart level. Findings We showed that loss of structural organization in failing myocytes leads to re-distribution of functional LTCCs from the T-tubules to the sarcolemma. In ischemic cardiomyopathy, the increased LTCC open probability in the T-tubules depends on the phosphorylation by protein kinase A, whereas in dilated cardiomyopathy, the increased LTCC opening probability in the sarcolemma results from enhanced phosphorylation by calcium-calmodulin kinase II. LVAD implantation corrected LTCCs pathophysiological activity, although it did not improve their distribution. Using computational modeling in a 3D anatomically-realistic human ventricular model, we showed how LTCC location and activity can trigger heart rhythm disorders of different severity. Interpretation Our findings demonstrate that LTCC redistribution and function differentiate between disease aetiologies. The subcellular changes observed in specific microdomains could be the consequence of the action of distinct protein kinases. Funding This work was supported by NIH grant (ROI-HL 126802 to NT-JG) and British Heart Foundation (grant RG/17/13/33173 to JG, project grant PG/16/17/32069 to RAC). Funders had no role in study design, data collection, data analysis, interpretation, writing of the repor

3D finite element electrical model of larval zebrafish ECG signals

Assessment of heart function in zebrafish larvae using electrocardiography (ECG) is a potentially useful tool in developing cardiac treatments and the assessment of drug therapies. In order to better understand how a measured ECG waveform is related to the structure of the heart, its position within the larva and the position of the electrodes, a 3D model of a 3 days post fertilisation (dpf) larval zebrafish was developed to simulate cardiac electrical activity and investigate the voltage distribution throughout the body. The geometry consisted of two main components; the zebrafish body was modelled as a homogeneous volume, while the heart was split into five distinct regions (sinoatrial region, atrial wall, atrioventricular band, ventricular wall and heart chambers). Similarly, the electrical model consisted of two parts with the body described by Laplace’s equation and the heart using a bidomain ionic model based upon the Fitzhugh-Nagumo equations. Each region of the heart was differentiated by action potential (AP) parameters and activation wave conduction velocities, which were fitted and scaled based on previously published experimental results. ECG measurements in vivo at different electrode recording positions were then compared to the model results. The model was able to simulate action potentials, wave propagation and all the major features (P wave, R wave, T wave) of the ECG, as well as polarity of the peaks observed at each position. This model was based upon our current understanding of the structure of the normal zebrafish larval heart. Further development would enable us to incorporate features associated with the diseased heart and hence assist in the interpretation of larval zebrafish ECGs in these conditions

Local hyperactivation of L-type Ca2+ channels increases spontaneous Ca2+ release activity and cellular hypertrophy in right ventricular myocytes from heart failure rats

Right ventricle (RV) dysfunction is an independent predictor of patient survival in heart failure (HF). However, the mechanisms of RV progression towards failing are not well understood. We studied cellular mechanisms of RV remodelling in a rat model of left ventricle myocardial infarction (MI)-caused HF. RV myocytes from HF rats show significant cellular hypertrophy accompanied with a disruption of transverse-axial tubular network and surface flattening. Functionally these cells exhibit higher contractility with lower Ca2+ transients. The structural changes in HF RV myocytes correlate with more frequent spontaneous Ca2+ release activity than in control RV myocytes. This is accompanied by hyperactivated L-type Ca2+ channels (LTCCs) located specifically in the T-tubules of HF RV myocytes. The increased open probability of tubular LTCCs and Ca2+ sparks activation is linked to protein kinase A-mediated channel phosphorylation that occurs locally in T-tubules. Thus, our approach revealed that alterations in RV myocytes in heart failure are specifically localized in microdomains. Our findings may indicate the development of compensatory, though potentially arrhythmogenic, RV remodelling in the setting of LV failure. These data will foster better understanding of mechanisms of heart failure and it could promote an optimized treatment of patients

Comparison of diffusion tensor imaging by cardiovascular magnetic resonance and gadolinium enhanced 3D image intensity approaches to investigation of structural anisotropy in explanted rat hearts

Background: Cardiovascular magnetic resonance (CMR) can through the two methods 3D FLASH and diffusion tensor imaging (DTI) give complementary information on the local orientations of cardiomyocytes and their laminar arrays. Methods: Eight explanted rat hearts were perfused with Gd-DTPA contrast agent and fixative and imaged in a 9.4T magnet by two types of acquisition: 3D fast low angle shot (FLASH) imaging, voxels 50 × 50 × 50 μm, and 3D spin echo DTI with monopolar diffusion gradients of 3.6 ms duration at 11.5 ms separation, voxels 200 × 200 × 200 μm. The sensitivity of each approach to imaging parameters was explored. Results:The FLASH data showed laminar alignments of voxels with high signal, in keeping with the presumed predominance of contrast in the interstices between sheetlets. It was analysed, using structure-tensor (ST) analysis, to determine the most (v 1 ST ), intermediate (v 2 ST ) and least (v 3 ST ) extended orthogonal directions of signal continuity. The DTI data was analysed to determine the most (e 1 DTI ), intermediate (e 2 DTI ) and least (e 3 DTI ) orthogonal eigenvectors of extent of diffusion. The correspondence between the FLASH and DTI methods was measured and appraised. The most extended direction of FLASH signal (v 1 ST ) agreed well with that of diffusion (e 1 DTI ) throughout the left ventricle (representative discrepancy in the septum of 13.3 ± 6.7°: median ± absolute deviation) and both were in keeping with the expected local orientations of the long-axis of cardiomyocytes. However, the orientation of the least directions of FLASH signal continuity (v 3 ST ) and diffusion (e 3 ST ) showed greater discrepancies of up to 27.9 ± 17.4°. Both FLASH (v 3 ST ) and DTI (e 3 DTI ) where compared to directly measured laminar arrays in the FLASH images. For FLASH the discrepancy between the structure-tensor calculated v 3 ST and the directly measured FLASH laminar array normal was of 9 ± 7° for the lateral wall and 7 ± 9° for the septum (median ± inter quartile range), and for DTI the discrepancy between the calculated v 3 DTI and the directly measured FLASH laminar array normal was 22 ± 14° and 61 ± 53.4°. DTI was relatively insensitive to the number of diffusion directions and to time up to 72 hours post fixation, but was moderately affected by b-value (which was scaled by modifying diffusion gradient pulse strength with fixed gradient pulse separation). Optimal DTI parameters were b = 1000 mm/s2 and 12 diffusion directions. FLASH acquisitions were relatively insensitive to the image processing parameters explored. Conclusions: We show that ST analysis of FLASH is a useful and accurate tool in the measurement of cardiac microstructure. While both FLASH and the DTI approaches appear promising for mapping of the alignments of myocytes throughout myocardium, marked discrepancies between the cross myocyte anisotropies deduced from each method call for consideration of their respective limitations

A simulation study of the reaction of human heart to biphasic electrical shocks

BACKGROUND: This article presents a study, which examines the effects of biphasic electrical shocks on human ventricular tissue. The effects of this type of shock are not yet fully understood. Animal experiments showed the superiority of biphasic shocks over monophasic ones in defibrillation. A mathematical computer simulation can increase the knowledge of human heart behavior. METHODS: The research presented in this article was done with different models representing a three-dimensional wedge of ventricular myocardium. The electrophysiology was described with Priebe-Beuckelmann model. The realistic fiber twist, which is specific to human myocardium was included. Planar electrodes were placed at the ends of the longest side of the virtual cardiac wedge, in a bath medium. They were sources of electrical shocks, which varied in magnitude from 0.1 to 5 V. In a second arrangement ring electrodes were placed directly on myocardium for getting a better view on secondary electrical sources. The electrical reaction of the tissue was generated with a bidomain model. RESULTS: The reaction of the tissue to the electrical shock was specific to the initial imposed characteristics. Depolarization appeared in the first 5 ms in different locations. A further study of the cardiac tissue behavior revealed, which features influence the response of the considered muscle. It was shown that the time needed by the tissue to be totally depolarized is much shorter when a biphasic shock is applied. Each simulation ended only after complete repolarization was achieved. This created the possibility of gathering information from all states corresponding to one cycle of the cardiac rhythm. CONCLUSIONS: The differences between the reaction of the homogeneous tissue and a tissue, which contains cleavage planes, reveals important aspects of superiority of biphasic pulses. ..

Heart Digital Twins Predict Features of Invasive Reentrant Circuits and Ablation Lesions in Scar-Dependent Ventricular Tachycardia

BACKGROUND: Catheter ablation of scar-dependent ventricular tachycardia (VT) is frequently hampered by hemodynamic instability, long procedure duration, and high recurrence rates. Magnetic resonance imaging-based personalized heart digital twins may overcome these challenges by noninvasively predicting VT circuits and optimum ablation lesion sites. In this combined clinical and digital twin study, we investigated the relationship between digital twin-predicted VTs and optimum ablation lesion sets with their invasively mapped counterparts during clinical VT ablation. METHODS: A total of 18 patients with scar-dependent VT underwent digital twin creation based on preprocedural, contrast-enhanced cardiac magnetic resonance imaging. Using rapid pacing protocols, VT was simulated and ablation targets were derived that would terminate all possible VTs in the models. Patients subsequently underwent invasive VT ablation, including targeting of diastolic activity and optimum entrainment sites. Digital twin-predicted VT circuits and ablation lesions were compared with their invasive clinical counterparts. RESULTS: Forty-three clinical VTs and 92 digital twin VTs were induced. Diastolic activity was seen in 16 of 43 (37.2%) clinical VTs. Sensitivity, specificity, positive predictive, and negative predictive values for the detection of critical VT sites by digital twins were 81.3%, 83.8%, 21.7%, and 98.8%, respectively. At an AHA-segment level, agreement between clinical VT critical sites and digital twin primary predicted sites was moderate, with a κ coefficient of 0.46 (±0.32; P≤0.001). Termination of VT with ablation was achieved at a digital twin-predicted site in 4 of 5 (80%) cases where attempted. A total of 426 of 709 (60.1%) lesions were within 5 mm of a predicted target site. In total, 54.0% (±28.9%) of the digital twin-predicted area was ablated per patient based on conventional mapping criteria. CONCLUSIONS: Heart digital twin VT circuits and ablation targets accurately predict many features of their respective clinical counterparts but have some limitations in spatial resolution. Our findings demonstrate the significant potential of digital twin technology in guiding catheter ablation for scar-dependent VT

Local hyperactivation of L-type Ca2+ channels increases spontaneous Ca2+ release activity and cellular hypertrophy in right ventricular myocytes from heart failure rats

Right ventricle (RV) dysfunction is an independent predictor of patient survival in heart failure (HF). However, the mechanisms of RV progression towards failing are not well understood. We studied cellular mechanisms of RV remodelling in a rat model of left ventricle myocardial infarction (MI)-caused HF. RV myocytes from HF rats show significant cellular hypertrophy accompanied with a disruption of transverse-axial tubular network and surface flattening. Functionally these cells exhibit higher contractility with lower Ca2+ transients. The structural changes in HF RV myocytes correlate with more frequent spontaneous Ca2+ release activity than in control RV myocytes. This is accompanied by hyperactivated L-type Ca2+ channels (LTCCs) located specifically in the T-tubules of HF RV myocytes. The increased open probability of tubular LTCCs and Ca2+ sparks activation is linked to protein kinase A-mediated channel phosphorylation that occurs locally in T-tubules. Thus, our approach revealed that alterations in RV myocytes in heart failure are specifically localized in microdomains. Our findings may indicate the development of compensatory, though potentially arrhythmogenic, RV remodelling in the setting of LV failure. These data will foster better understanding of mechanisms of heart failure and it could promote an optimized treatment of patients