Effects of Electrical and Structural Remodeling on Atrial Fibrillation Maintenance: A Simulation Study

Atrial fibrillation, a common cardiac arrhythmia, often progresses unfavourably: in patients with long-term atrial fibrillation, fibrillatory episodes are typically of increased duration and frequency of occurrence relative to healthy controls. This is due to electrical, structural, and contractile remodeling processes. We investigated mechanisms of how electrical and structural remodeling contribute to perpetuation of simulated atrial fibrillation, using a mathematical model of the human atrial action potential incorporated into an anatomically realistic three-dimensional structural model of the human atria. Electrical and structural remodeling both shortened the atrial wavelength - electrical remodeling primarily through a decrease in action potential duration, while structural remodeling primarily slowed conduction. The decrease in wavelength correlates with an increase in the average duration of atrial fibrillation/flutter episodes. The dependence of reentry duration on wavelength was the same for electrical vs. structural remodeling. However, the dynamics during atrial reentry varied between electrical, structural, and combined electrical and structural remodeling in several ways, including: (i) with structural remodeling there were more occurrences of fragmented wavefronts and hence more filaments than during electrical remodeling; (ii) dominant waves anchored around different anatomical obstacles in electrical vs. structural remodeling; (iii) dominant waves were often not anchored in combined electrical and structural remodeling. We conclude that, in simulated atrial fibrillation, the wavelength dependence of reentry duration is similar for electrical and structural remodeling, despite major differences in overall dynamics, including maximal number of filaments, wave fragmentation, restitution properties, and whether dominant waves are anchored to anatomical obstacles or spiralling freely

Non-invasive localization of atrial ectopic beats by using simulated body surface P-wave integral maps

Non-invasive localization of continuous atrial ectopic beats remains a cornerstone for the treatment of atrial arrhythmias. The lack of accurate tools to guide electrophysiologists leads to an increase in the recurrence rate of ablation procedures. Existing approaches are based on the analysis of the P-waves main characteristics and the forward body surface potential maps (BSPMs) or on the inverse estimation of the electric activity of the heart from those BSPMs. These methods have not provided an efficient and systematic tool to localize ectopic triggers. In this work, we propose the use of machine learning techniques to spatially cluster and classify ectopic atrial foci into clearly differentiated atrial regions by using the body surface P-wave integral map (BSPiM) as a biomarker. Our simulated results show that ectopic foci with similar BSPiM naturally cluster into differentiated non-intersected atrial regions and that new patterns could be correctly classified with an accuracy of 97% when considering 2 clusters and 96% for 4 clusters. Our results also suggest that an increase in the number of clusters is feasible at the cost of decreasing accuracy.This work was partially supported by The "Programa Prometeu" from Conselleria d'Educacio Formacio I Ocupacio, Generalitat Valenciana (www.edu.gva.es/fio/index_es.asp) Award Number: PROMETEU/2016/088 to JS; The "Plan Estatal de Investigacion Cientifica y Tecnica y de Innovacion 2013-2016" from the Ministerio de Economia, Industria y Competitividad of Spain, Agencia Estatal de Investigacion (www.mineco.gob.es) and the European Commission (European Regional Development Funds - ERDF -FEDER) (ec.europa.eu/regional_policy/es/funding/erdf/) Award Number: DPI2016-75799-R to JS and The "Programa Estatal de Investigacion, Desarrollo e Innovacion Orientado a los Retos de la Sociedad" from the Ministerio de Economia y Competitividad of Spain, Agencia Estatal de Investigacion (www.mineco.gob.es) and the European Commission (European Regional Development Funds - ERDF -FEDER) (ec.europa.eu/regional_policy/es/funding/erdf/) Award Number: TIN2014-59932-JIN to AFA and RS. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Ferrer Albero, A.; Godoy, EJ.; Lozano, M.; Martínez Mateu, L.; Alonso Atienza, F.; Saiz Rodríguez, FJ.; Sebastián Aguilar, R. (2017). Non-invasive localization of atrial ectopic beats by using simulated body surface P-wave integral maps. PLoS ONE. 12(7):1-23. https://doi.org/10.1371/journal.pone.0181263S12312

Impact of Sarcoplasmic Reticulum Calcium Release on Calcium Dynamics and Action Potential Morphology in Human Atrial Myocytes: A Computational Study

Electrophysiological studies of the human heart face the fundamental challenge that experimental data can be acquired only from patients with underlying heart disease. Regarding human atria, there exist sizable gaps in the understanding of the functional role of cellular Ca2+ dynamics, which differ crucially from that of ventricular cells, in the modulation of excitation-contraction coupling. Accordingly, the objective of this study was to develop a mathematical model of the human atrial myocyte that, in addition to the sarcolemmal (SL) ion currents, accounts for the heterogeneity of intracellular Ca2+ dynamics emerging from a structurally detailed sarcoplasmic reticulum (SR). Based on the simulation results, our model convincingly reproduces the principal characteristics of Ca2+ dynamics: 1) the biphasic increment during the upstroke of the Ca2+ transient resulting from the delay between the peripheral and central SR Ca2+ release, and 2) the relative contribution of SL Ca2+ current and SR Ca2+ release to the Ca2+ transient. In line with experimental findings, the model also replicates the strong impact of intracellular Ca2+ dynamics on the shape of the action potential. The simulation results suggest that the peripheral SR Ca2+ release sites define the interface between Ca2+ and AP, whereas the central release sites are important for the fire-diffuse-fire propagation of Ca2+ diffusion. Furthermore, our analysis predicts that the modulation of the action potential duration due to increasing heart rate is largely mediated by changes in the intracellular Na+ concentration. Finally, the results indicate that the SR Ca2+ release is a strong modulator of AP duration and, consequently, myocyte refractoriness/excitability. We conclude that the developed model is robust and reproduces many fundamental aspects of the tight coupling between SL ion currents and intracellular Ca2+ signaling. Thus, the model provides a useful framework for future studies of excitation-contraction coupling in human atrial myocytes

A three-dimensional human atrial model with fiber orientation. Electrograms and arrhythmic activation patterns relationship

The most common sustained cardiac arrhythmias in humans are atrial tachyarrhythmias, mainly atrial fibrillation. Areas of complex fractionated atrial electrograms and high dominant frequency have been proposed as critical regions for maintaining atrial fibrillation; however, there is a paucity of data on the relationship between the characteristics of electrograms and the propagation pattern underlying them. In this study, a realistic 3D computer model of the human atria has been developed to investigate this relationship. The model includes a realistic geometry with fiber orientation, anisotropic conductivity and electrophysiological heterogeneity. We simulated different tachyarrhythmic episodes applying both transient and continuous ectopic activity. Electrograms and their dominant frequency and organization index values were calculated over the entire atrial surface. Our simulations show electrograms with simple potentials, with little or no cycle length variations, narrow frequency peaks and high organization index values during stable and regular activity as the observed in atrial flutter, atrial tachycardia (except in areas of conduction block) and in areas closer to ectopic activity during focal atrial fibrillation. By contrast, cycle length variations and polymorphic electrograms with single, double and fragmented potentials were observed in areas of irregular and unstable activity during atrial fibrillation episodes. Our results also show: 1) electrograms with potentials without negative deflection related to spiral or curved wavefronts that pass over the recording point and move away, 2) potentials with a much greater proportion of positive deflection than negative in areas of wave collisions, 3) double potentials related with wave fragmentations or blocking lines and 4) fragmented electrograms associated with pivot points. Our model is the first human atrial model with realistic fiber orientation used to investigate the relationship between different atrial arrhythmic propagation patterns and the electrograms observed at more than 43000 points on the atrial surface.This work was partially supported by the Plan Nacional de Investigacion Cientifica, Desarrollo e Innovacion Tecnologica, Ministerio de Ciencia e Innovacion of Spain (TEC2008-02090), by the Plan Avanza (Accion Estrategica de Telecomunicaciones y Sociedad de la Informacion), Ministerio de Industria Turismo y Comercio of Spain (TSI-020100-2010-469), by the Programa Prometeo 2012 of the Generalitat Valenciana and by the Programa de Apoyo a la Investigacion y Desarrollo de la Universitat Politecnica de Valencia (PAID-06-11-2002). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Tobón Zuluaga, C.; Ruiz Villa, CA.; Heidenreich, E.; Romero Pérez, L.; Hornero, F.; Saiz Rodríguez, FJ. (2013). A three-dimensional human atrial model with fiber orientation. Electrograms and arrhythmic activation patterns relationship. PLoS ONE. 8(2):1-13. https://doi.org/10.1371/journal.pone.0050883S11382Ho SY, Sanchez-Quintana D, Anderson RH (1998) Can anatomy define electric pathways? In: International Workshop on Computer Simulation and Experimental Assessment of Electrical Cardiac Function, Lausanne, Switzerland. 77–86.Tobón C (2009) Evaluación de factores que provocan fibrilación auricular y de su tratamiento mediante técnicas quirúrgicas. Estudio de simulación. Master Thesis Universitat Politècnica de València.Ruiz C (2010) Estudio de la vulnerabilidad a reentradas a través de modelos matemáticos y simulación de la aurícula humana. Doctoral Thesis Universitat Politècnica de València.Tobón C (2010) Modelización y evaluación de factores que favorecen las arritmias auriculares y su tratamiento mediante técnicas quirúrgicas. Estudio de simulación. Doctoral Thesis Universitat Politècnica de València.Henriquez, C. S., & Papazoglou, A. A. (1996). Using computer models to understand the roles of tissue structure and membrane dynamics in arrhythmogenesis. Proceedings of the IEEE, 84(3), 334-354. doi:10.1109/5.486738Grimm, R. A., Chandra, S., Klein, A. L., Stewart, W. J., Black, I. W., Kidwell, G. A., & Thomas, J. D. (1996). Characterization of left atrial appendage Doppler flow in atrial fibrillation and flutter by Fourier analysis. American Heart Journal, 132(2), 286-296. doi:10.1016/s0002-8703(96)90424-xMaleckar, M. M., Greenstein, J. L., Giles, W. R., & Trayanova, N. A. (2009). K+ current changes account for the rate dependence of the action potential in the human atrial myocyte. American Journal of Physiology-Heart and Circulatory Physiology, 297(4), H1398-H1410. doi:10.1152/ajpheart.00411.200

Detailed Anatomical and Electrophysiological Models of Human Atria and Torso for the Simulation of Atrial Activation

Atrial arrhythmias, and specifically atrial fibrillation (AF), induce rapid and irregular activation patterns that appear on the torso surface as abnormal P-waves in electrocardiograms and body surface potential maps (BSPM). In recent years both P-waves and the BSPM have been used to identify the mechanisms underlying AF, such as localizing ectopic foci or high-frequency rotors. However, the relationship between the activation of the different areas of the atria and the characteristics of the BSPM and P-wave signals are still far from being completely understood. In this work we developed a multi-scale framework, which combines a highly-detailed 3D atrial model and a torso model to study the relationship between atrial activation and surface signals in sinus rhythm. Using this multi scale model, it was revealed that the best places for recording P-waves are the frontal upper right and the frontal and rear left quadrants of the torso. Our results also suggest that only nine regions (of the twenty-one structures in which the atrial surface was divided) make a significant contribution to the BSPM and determine the main P-wave characteristics.This work was partially supported by the "VI Plan Nacional de Investigacion Cientifica, Desarrollo e Innovacion Tecnologica" from the Ministerio de Economia y Competitividad of Spain and the European Commission (European Regional Development Funds - ERDF - FEDER), Award Number: TIN2012-37546-C03-01 (Recipient: Ana Ferrer); the "Programa Estatal de Investigacion, Desarrollo e Innovacion Orientado a los Retos de la Sociedad" from the Ministerio de Economia y Competitividad and the European Commission (European Regional Development Funds - ERDF - FEDER), Award Number: TIN2014-59932-JIN (Recipient: Rafael Sebastion); and the "Programa Prometeo" from the Generalitat Valenciana, Award Number: 2012/030 (Recipient: Laura Martinez).Ferrer Albero, A.; Sebastián Aguilar, R.; Sánchez Quintana, D.; Rodriguez, JF.; Godoy, EJ.; Martinez, L.; Saiz Rodríguez, FJ. (2015). Detailed Anatomical and Electrophysiological Models of Human Atria and Torso for the Simulation of Atrial Activation. PLoS ONE. 10(11):1-29. https://doi.org/10.1371/journal.pone.0141573S129101

Human Atrial Electrophysiological Models Under Fractional Derivative: Depolarization and Repolarization Dynamics During Normal and Fibrillation Conditions

Atrial fibrillation (AF) is the most common arrhythmia within the clinical context. Advanced stages of the AF involve several difficulties in its management and treatment. This occurs mostly because the initiation and perpetuation mechanisms of the AF are still not fully understood. Cardiac scientific computation has become an important tool in researching the underlying mechanisms of the AF. In this work, an equation of action potential propagation that implements fractional order derivatives is used to model the atrial dynamics. The fractional derivative order represents the structural heterogeneities of the atrial myocardium. Using such mathematical operator, the Courtemanche and Maleckar human atrial electrophysiological models, during healthy and AF conditions, are assessed. The results indicate that, through the fractional order variations, there are electrophysiological properties whose behavior do not depend on the cellular model or physiological conditions. On the other hand, there are properties whose behavior under distinct values of the fractional order, are specific to the cellular model and to the physiological condition and they can be characterized quantitatively and qualitatively. Therefore, the fractional atrial propagation model can be a useful tool for modeling a wide range of electrophysiological scenarios in both healthy and AF conditions. © 2019, Springer Nature Switzerland AG

Genesis of Atrial Fibrillation Under Different Diffuse Fibrosis Density Related with Atmospheric Pollution. In-Silico Study

Atrial remodeling is a widely acknowledged process that accelerates the susceptibility to and progression of atrial fibrillation. An increasingly recognized structural component is atrial fibrosis. Recent studies have shown that air pollution increases the risk of heart arrhythmias, where the exposure to particulate matter (PM) contributes to the generation of myocardial fibrosis, increasing the cardiovascular risk. The density and patterns of fibrosis (interstitial, compact and diffuse) are relevant in abnormal conduction and vulnerability to cardiac arrhythmias. Taking into account that fibrosis has been widely reported as one of the consequences of PM exposure, in this work, we evaluated the effects of low and high diffuse fibrosis density on conduction velocity and arrhythmic propagation patterns. For this purpose, cellular models of atrial myocyte and fibroblast were implemented in a 3D model of the human atria. Low (6.25%) and high (25%) fibrosis densities were simulated in the left atrium and its effect on conduction velocity and fibrillatory dynamics was evaluated. Results showed a conduction velocity reduction of 71% associated with a high fibrosis density. At low fibrosis density, few reentries were observed. On the other hand, at high fibrosis density, irregular propagation patterns, characterized by multiple wavelets and rotors, were observed. Our results suggest that high diffuse fibrosis density is associated with a significant conduction velocity reduction and with chaotic propagation patterns during atrial fibrillation. © 2020, Springer Nature Switzerland AG