Slow Conduction in the Border Zones of Patchy Fibrosis Stabilizes the Drivers for Atrial Fibrillation: Insights from Multi-Scale Human Atrial Modeling

Introduction: The genesis of atrial fibrillation (AF) and success of AF ablation therapy have been strongly linked with atrial fibrosis. Increasing evidence suggests that patient-specific distributions of fibrosis may determine the locations of electrical drivers (rotors) sustaining AF, but the underlying mechanisms are incompletely understood. This study aims to elucidate a missing mechanistic link between patient-specific fibrosis distributions and AF drivers. Methods: 3D atrial models integrated human atrial geometry, rule-based fiber orientation, region-specific electrophysiology, and AF-induced ionic remodeling. A novel detailed model for an atrial fibroblast was developed, and effects of myocyte-fibroblast (M-F) coupling were explored at single-cell, 1D tissue and 3D atria levels. Left atrial LGE MRI datasets from 3 chronic AF patients were segmented to provide the patient-specific distributions of fibrosis. The data was non-linearly registered and mapped to the 3D atria model. Six distinctive fibrosis levels (0–healthy tissue, 5–dense fibrosis) were identified based on LGE MRI intensity and modeled as progressively increasing M-F coupling and decreasing atrial tissue coupling. Uniform 3D atrial model with diffuse (level 2) fibrosis was considered for comparison. Results: In single cells and tissue, the largest effect of atrial M-F coupling was on the myocyte resting membrane potential, leading to partial inactivation of sodium current and reduction of conduction velocity (CV). In the 3D atria, further to the M-F coupling, effects of fibrosis on tissue coupling greatly reduce atrial CV. AF was initiated by fast pacing in each 3D model with either uniform or patient-specific fibrosis. High variation in fibrosis distributions between the models resulted in varying complexity of AF, with several drivers emerging. In the diffuse fibrosis models, waves randomly meandered through the atria, whereas in each the patient-specific models, rotors stabilized in fibrotic regions. The rotors propagated slowly around the border zones of patchy fibrosis (levels 3–4), failing to spread into inner areas of dense fibrosis. Conclusion: Rotors stabilize in the border zones of patchy fibrosis in 3D atria, where slow conduction enable the development of circuits within relatively small regions. Our results can provide a mechanistic explanation for the clinical efficacy of ablation around fibrotic regions

Investigating Calcium-Mediated Arrhythmias via a Computational Model of a Rabbit Atrial Myocyte

The system of transverse and longitudinal sarcolemmal tubules (T-system) is observed to remodel in atrial fibrillation (AF) and heart failure (HF). The resulting calcium dysregulation has been suggested to underlie disruptions in excitation-contraction coupling, and increase the frequency of arrhythmic events at the cellular scale; however, these mechanisms and their importance are yet to be fully described. A stochastic, 3D, spatiotemporal model of the rabbit atrial myocyte was developed in order to study calcium-mediated arrhythmic phenomena at the cellular scale. Preliminary findings suggest a relationship between the severity of detubulation, and the promotion of spontaneous activity, and provides insight into the conditions required for the emergence of spontaneous activity within atrial myocytes in HF

High resolution 3-Dimensional imaging of the human cardiac conduction system from microanatomy to mathematical modeling

Cardiac arrhythmias and conduction disturbances are accompanied by structural remodelling of the specialised cardiomyocytes known collectively as the cardiac conduction system. Here, using contrast enhanced micro-computed tomography, we present, in attitudinally appropriate fashion, the first 3-dimensional representations of the cardiac conduction system within the intact human heart. We show that cardiomyocyte orientation can be extracted from these datasets at spatial resolutions approaching the single cell. These data show that commonly accepted anatomical representations are oversimplified. We have incorporated the high-resolution anatomical data into mathematical simulations of cardiac electrical depolarisation. The data presented should have multidisciplinary impact. Since the rate of depolarisation is dictated by cardiac microstructure, and the precise orientation of the cardiomyocytes, our data should improve the fidelity of mathematical models. By showing the precise 3-dimensional relationships between the cardiac conduction system and surrounding structures, we provide new insights relevant to valvar replacement surgery and ablation therapies. We also offer a practical method for investigation of remodelling in disease, and thus, virtual pathology and archiving. Such data presented as 3D images or 3D printed models, will inform discussions between medical teams and their patients, and aid the education of medical and surgical trainees

Relevance of Positron-Emission Therapy for Optimization of Treatment of Advanced Hodgkin’s Lymphoma Using Intensive ЕАСОРР-14 Program

Aim. To evaluate the relevance of the positron-emission therapy (PET) for optimization of the therapy of advanced Hodgkin’s lymphoma (HL) using the intensive EACOPP-14 program. Materials & Methods. 91 patients with advanced HL (IIX–IIE, III–IV) received the treatment according to the “ЛХМосква1-3” protocol over the period from November 2009 to February 2015, and then the treatment was analyzed. The median age was 29 years (range: 17–50); there were 42 men (46.3 %) and 49 (53.7 %) women. The treatment included 6 cycles of polychemotherapy according to the regimen ЕА(50)СОРР-14 ± radiation therapy. The radiation therapy was performed in 66 patients (72.5 %) after the completion of the chemotherapy. The cumulative focal dose was 30 Gy onto the areas of residual lesions and/or initially large tumor masses. Results. PET performed during the initial HL diagnosing permited to identify new areas of neoplastic lesions without changes in staging and treatment scheme, as well as specify areas and field size of planned radiation consolidation. The paper confirms the prognostic value of the intermediate PET in patients with advanced HL during the intensive first-line chemotherapy. The intensive therapy at the beginning of the treatment program is associated with higher chances for survival for patients with extremely unfavorable prognosis. After completion of the drug therapy, negative PET findings had a higher prognostic value, than the positive ones. The analysis of the relevance of residual tumor dimensions in the PET negative group demonstrated that the relapses were more common, if the residual tumor was more than 4.5 cm (according to CT findings). Conclusion. This study confirmed that it reasonable to discuss the discontinuation of the radiation therapy in patients with advanced HL, negative PET findings and small (< 2.5 cm) residual tumor after the intensive ЕАСОРР-14 program. This tactics permits avoiding a number of delayed complications

Atrial arrhythmogenicity of KCNJ2 mutations in short QT syndrome: Insights from virtual human atria

Gain-of-function mutations in KCNJ2-encoded Kir2.1 channels underlie variant 3 (SQT3) of the short QT syndrome, which is associated with atrial fibrillation (AF). Using biophysically-detailed human atria computer models, this study investigated the mechanistic link between SQT3 mutations and atrial arrhythmogenesis, and potential ion channel targets for treatment of SQT3. A contemporary model of the human atrial action potential (AP) was modified to recapitulate functional changes in IK1 due to heterozygous and homozygous forms of the D172N and E299V Kir2.1 mutations. Wild-type (WT) and mutant formulations were incorporated into multi-scale homogeneous and heterogeneous tissue models. Effects of mutations on AP duration (APD), conduction velocity (CV), effective refractory period (ERP), tissue excitation threshold and their rate-dependence, as well as the wavelength of re-entry (WL) were quantified. The D172N and E299V Kir2.1 mutations produced distinct effects on IK1 and APD shortening. Both mutations decreased WL for re-entry through a reduction in ERP and CV. Stability of re-entrant excitation waves in 2D and 3D tissue models was mediated by changes to tissue excitability and dispersion of APD in mutation conditions. Combined block of IK1 and IKr was effective in terminating re-entry associated with heterozygous D172N conditions, whereas IKr block alone may be a safer alternative for the E299V mutation. Combined inhibition of IKr and IKur produced a synergistic anti-arrhythmic effect in both forms of SQT3. In conclusion, this study provides mechanistic insights into atrial proarrhythmia with SQT3 Kir2.1 mutations and highlights possible pharmacological strategies for management of SQT3-linked AF

Novel non-invasive algorithm to identify the origins of re-entry and ectopic foci in the atria from 64-lead ECGs: A computational study.

Atrial tachy-arrhytmias, such as atrial fibrillation (AF), are characterised by irregular electrical activity in the atria, generally associated with erratic excitation underlain by re-entrant scroll waves, fibrillatory conduction of multiple wavelets or rapid focal activity. Epidemiological studies have shown an increase in AF prevalence in the developed world associated with an ageing society, highlighting the need for effective treatment options. Catheter ablation therapy, commonly used in the treatment of AF, requires spatial information on atrial electrical excitation. The standard 12-lead electrocardiogram (ECG) provides a method for non-invasive identification of the presence of arrhythmia, due to irregularity in the ECG signal associated with atrial activation compared to sinus rhythm, but has limitations in providing specific spatial information. There is therefore a pressing need to develop novel methods to identify and locate the origin of arrhythmic excitation. Invasive methods provide direct information on atrial activity, but may induce clinical complications. Non-invasive methods avoid such complications, but their development presents a greater challenge due to the non-direct nature of monitoring. Algorithms based on the ECG signals in multiple leads (e.g. a 64-lead vest) may provide a viable approach. In this study, we used a biophysically detailed model of the human atria and torso to investigate the correlation between the morphology of the ECG signals from a 64-lead vest and the location of the origin of rapid atrial excitation arising from rapid focal activity and/or re-entrant scroll waves. A focus-location algorithm was then constructed from this correlation. The algorithm had success rates of 93% and 76% for correctly identifying the origin of focal and re-entrant excitation with a spatial resolution of 40 mm, respectively. The general approach allows its application to any multi-lead ECG system. This represents a significant extension to our previously developed algorithms to predict the AF origins in association with focal activities

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

In silico assessment of the effects of quinidine, disopyramide and E-4031 on short QT syndrome variant 1 in the human ventricles.

Aims Short QT syndrome (SQTS) is an inherited disorder associated with abnormally abbreviated QT intervals and an increased incidence of atrial and ventricular arrhythmias. SQT1 variant (linked to the rapid delayed rectifier potassium channel current, IKr) of SQTS, results from an inactivation-attenuated, gain-of-function mutation (N588K) in the KCNH2-encoded potassium channels. Pro-arrhythmogenic effects of SQT1 have been well characterized, but less is known about the possible pharmacological antiarrhythmic treatment of SQT1. Therefore, this study aimed to assess the potential effects of E-4031, disopyramide and quinidine on SQT1 using a mathematical model of human ventricular electrophysiology. Methods The ten Tusscher et al. biophysically detailed model of the human ventricular action potential (AP) was modified to incorporate IKr Markov chain (MC) formulations based on experimental data of the kinetics of the N588K mutation of the KCNH2-encoded subunit of the IKr channels. The modified ventricular cell model was then integrated into one-dimensional (1D) strand, 2D regular and realistic tissues with transmural heterogeneities. The channel-blocking effect of the drugs on ion currents in healthy and SQT1 cells was modeled using half-maximal inhibitory concentration (IC50) and Hill coefficient (nH) values from literatures. Effects of drugs on cell AP duration (APD), effective refractory period (ERP) and pseudo-ECG traces were calculated. Effects of drugs on the ventricular temporal and spatial vulnerability to re-entrant excitation waves were measured. Re-entry was simulated in both 2D regular and realistic ventricular tissue. Results At the single cell level, the drugs E-4031 and disopyramide had hardly noticeable effects on the ventricular cell APD at 90% repolarization (APD90), whereas quinidine caused a significant prolongation of APD90. Quinidine prolonged and decreased the maximal transmural AP heterogeneity (δV); this led to the decreased transmural heterogeneity of APD across the 1D strand. Quinidine caused QT prolongation and a decrease in the T-wave amplitude, and increased ERP and decreased temporal susceptibility of the tissue to the initiation of re-entry and increased the minimum substrate size necessary to prevent re-entry in the 2D regular model, and further terminated re-entrant waves in the 2D realistic model. Quinidine exhibited significantly better therapeutic effects on SQT1 than E-4031 and disopyramide. Conclusions The simulated pharmacological actions of quinidine exhibited antiarrhythmic effects on SQT1. This study substantiates a causal link between quinidine and QT interval prolongation in SQT1 and suggests that quinidine may be a potential pharmacological agent for treating SQT1 patients

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