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

Mathematical Modeling and Simulation of Ventricular Activation Sequences: Implications for Cardiac Resynchronization Therapy

Next to clinical and experimental research, mathematical modeling plays a crucial role in medicine. Biomedical research takes place on many different levels, from molecules to the whole organism. Due to the complexity of biological systems, the interactions between components are often difficult or impossible to understand without the help of mathematical models. Mathematical models of cardiac electrophysiology have made a tremendous progress since the first numerical ECG simulations in the 1960s. This paper briefly reviews the development of this field and discusses some example cases where models have helped us forward, emphasizing applications that are relevant for the study of heart failure and cardiac resynchronization therapy

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

2017 HRS/EHRA/ECAS/APHRS/SOLAECE expert consensus statement on catheter and surgical ablation of atrial fibrillation: executive summary.

S

In-silico study of the ionic current gradients determining left-to-right atrial frequencies during paroxysmal atrial fibrillation

Atrial fibrillation is the most prevalent cardiac arrhythmia. Paroxysmal atrial fibrillation (pAF) may occur in episodes lasting from minutes to days. Recent studies suggest that some pAF episodes present a left-to-right dominant frequency gradient caused by ionic current gradients. However, how each ionic current gradient affects the left-to-right dominant frequency gradient during pAF has not been studied. In this work, we use a 3D model of human atria to study how the ionic current gradients affect the dominant frequency gradient during pAF induced by continuous ectopic activity. The role of the specific gradients of acetylcholine-activated potassium current (I KACh ) and inward-rectifier potassium current (I K1 ) on determining the left-to-right dominant frequency gradient was assessed. The main outcome of this study is that either or both of the I KACh or I K1 gradients are necessary to induce a left-to-right dominant frequency gradient during pAF. However, both gradients are necessary to the left atrium maintaining, by itself, the pAF episode. These findings have potentially important implications for the development of atrial-selective therapeutic approaches. © The Author(s) 2019