Comparison of Propagation Models and Forward Calculation Methods on Cellular, Tissue and Organ Scale Atrial Electrophysiology

The bidomain model and the finite element method are an established standard to mathematically describe cardiac electrophysiology, but are both suboptimal choices for fast and large-scale simulations due to high computational costs. We investigate to what extent simplified approaches for propagation models (monodomain, reaction-Eikonal and Eikonal) and forward calculation (boundary element and infinite volume conductor) deliver markedly accelerated, yet physiologically accurate simulation results in atrial electrophysiology. Methods: We compared action potential durations, local activation times (LATs), and electrocardiograms (ECGs) for sinus rhythm simulations on healthy and fibrotically infiltrated atrial models. Results: All simplified model solutions yielded LATs and P waves in accurate accordance with the bidomain results. Only for the Eikonal model with pre-computed action potential templates shifted in time to derive transmembrane voltages, repolarization behavior notably deviated from the bidomain results. ECGs calculated with the boundary element method were characterized by correlation coefficients >0.9 compared to the finite element method. The infinite volume conductor method led to lower correlation coefficients caused predominantly by systematic overestimations of P wave amplitudes in the precordial leads. Conclusion: Our results demonstrate that the Eikonal model yields accurate LATs and combined with the boundary element method precise ECGs compared to markedly more expensive full bidomain simulations. However, for an accurate representation of atrial repolarization dynamics, diffusion terms must be accounted for in simplified models. Significance: Simulations of atrial LATs and ECGs can be notably accelerated to clinically feasible time frames at high accuracy by resorting to the Eikonal and boundary element methods

A staggered-in-time and non-conforming-in-space numerical framework for realistic cardiac electrophysiology outputs

Computer-based simulations of non-invasive cardiac electrical outputs, such as electrocardiograms and body surface potential maps, usually entail severe computational costs due to the need of capturing fine-scale processes and to the complexity of the heart-torso morphology. In this work, we model cardiac electrical outputs by employing a coupled model consisting of a reaction-diffusion model - either the bidomain model or the most efficient pseudo-bidomain model - on the heart, and an elliptic model in the torso. We then solve the coupled problem with a segregated and staggered in-time numerical scheme, that allows for independent and infrequent solution in the torso region. To further reduce the computational load, main novelty of this work is in introduction of an interpolation method at the interface between the heart and torso domains, enabling the use of non-conforming meshes, and the numerical framework application to realistic cardiac and torso geometries. The reliability and efficiency of the proposed scheme is tested against the corresponding state-of-the-art bidomain-torso model. Furthermore, we explore the impact of torso spatial discretization and geometrical non-conformity on the model solution and the corresponding clinical outputs. The investigation of the interface interpolation method provides insights into the influence of torso spatial discretization and of the geometrical non-conformity on the simulation results and their clinical relevance.Comment: 26 pages,11 figures, 3 table

Multiscale Cohort Modeling of Atrial Electrophysiology : Risk Stratification for Atrial Fibrillation through Machine Learning on Electrocardiograms

Patienten mit Vorhofflimmern sind einem fünffach erhöhten Risiko für einen ischämischen Schlaganfall ausgesetzt. Eine frühzeitige Erkennung und Diagnose der Arrhythmie würde ein rechtzeitiges Eingreifen ermöglichen, um möglicherweise auftretende Begleiterkrankungen zu verhindern. Eine Vergrößerung des linken Vorhofs sowie fibrotisches Vorhofgewebe sind Risikomarker für Vorhofflimmern, da sie die notwendigen Voraussetzungen für die Aufrechterhaltung der chaotischen elektrischen Depolarisation im Vorhof erfüllen. Mithilfe von Techniken des maschinellen Lernens könnten Fibrose und eine Vergrößerung des linken Vorhofs basierend auf P Wellen des 12-Kanal Elektrokardiogramms im Sinusrhythmus automatisiert identifiziert werden. Dies könnte die Basis für eine nicht-invasive Risikostrat- ifizierung neu auftretender Vorhofflimmerepisoden bilden, um anfällige Patienten für ein präventives Screening auszuwählen. Zu diesem Zweck wurde untersucht, ob simulierte Vorhof-Elektrokardiogrammdaten, die dem klinischen Trainingssatz eines maschinellen Lernmodells hinzugefügt wurden, zu einer verbesserten Klassifizierung der oben genannten Krankheiten bei klinischen Daten beitra- gen könnten. Zwei virtuelle Kohorten, die durch anatomische und funktionelle Variabilität gekennzeichnet sind, wurden generiert und dienten als Grundlage für die Simulation großer P Wellen-Datensätze mit genau bestimmbaren Annotationen der zugrunde liegenden Patholo- gie. Auf diese Weise erfüllen die simulierten Daten die notwendigen Voraussetzungen für die Entwicklung eines Algorithmus für maschinelles Lernen, was sie von klinischen Daten unterscheidet, die normalerweise nicht in großer Zahl und in gleichmäßig verteilten Klassen vorliegen und deren Annotationen möglicherweise durch unzureichende Expertenannotierung beeinträchtigt sind. Für die Schätzung des Volumenanteils von linksatrialem fibrotischen Gewebe wurde ein merkmalsbasiertes neuronales Netz entwickelt. Im Vergleich zum Training des Modells mit nur klinischen Daten, führte das Training mit einem hybriden Datensatz zu einer Reduzierung des Fehlers von durchschnittlich 17,5 % fibrotischem Volumen auf 16,5 %, ausgewertet auf einem rein klinischen Testsatz. Ein Long Short-Term Memory Netzwerk, das für die Unterscheidung zwischen gesunden und P Wellen von vergrößerten linken Vorhöfen entwickelt wurde, lieferte eine Genauigkeit von 0,95 wenn es auf einem hybriden Datensatz trainiert wurde, von 0,91 wenn es nur auf klinischen Daten trainiert wurde, die alle mit 100 % Sicherheit annotiert wurden, und von 0,83 wenn es auf einem klinischen Datensatz trainiert wurde, der alle Signale unabhängig von der Sicherheit der Expertenannotation enthielt. In Anbetracht der Ergebnisse dieser Arbeit können Elektrokardiogrammdaten, die aus elektrophysiologischer Modellierung und Simulationen an virtuellen Patientenkohorten resul- tieren und relevante Variabilitätsaspekte abdecken, die mit realen Beobachtungen übereinstim- men, eine wertvolle Datenquelle zur Verbesserung der automatisierten Risikostratifizierung von Vorhofflimmern sein. Auf diese Weise kann den Nachteilen klinischer Datensätze für die Entwicklung von Modellen des maschinellen Lernens entgegengewirkt werden. Dies trägt letztendlich zu einer frühzeitigen Erkennung der Arrhythmie bei, was eine rechtzeitige Auswahl geeigneter Behandlungsstrategien ermöglicht und somit das Schlaganfallrisiko der betroffenen Patienten verringert

A Review of Healthy and Fibrotic Myocardium Microstructure Modeling and Corresponding Intracardiac Electrograms

Computational simulations of cardiac electrophysiology provide detailed information on the depolarization phenomena at different spatial and temporal scales. With the development of new hardware and software, in silico experiments have gained more importance in cardiac electrophysiology research. For plane waves in healthy tissue, in vivo and in silico electrograms at the surface of the tissue demonstrate symmetric morphology and high peak-to-peak amplitude. Simulations provided insight into the factors that alter the morphology and amplitude of the electrograms. The situation is more complex in remodeled tissue with fibrotic infiltrations. Clinically, different changes including fractionation of the signal, extended duration and reduced amplitude have been described. In silico, numerous approaches have been proposed to represent the pathological changes on different spatial and functional scales. Different modeling approaches can reproduce distinct subsets of the clinically observed electrogram phenomena. This review provides an overview of how different modeling approaches to incorporate fibrotic and structural remodeling affect the electrogram and highlights open challenges to be addressed in future research

Scalable and Accurate ECG Simulation for Reaction-Diffusion Models of the Human Heart

International audienceRealistic electrocardiogram (ECG) simulation with numerical models is important for research linking cellular and molecular physiology to clinically observable signals, and crucial for patient tailoring of numerical heart models. However, ECG simulation with a realistic torso model is computationally much harder than simulation of cardiac activity itself, so that many studies with sophisticated heart models have resorted to crude approximations of the ECG. This paper shows how the classical concept of electrocardiographic lead fields can be used for an ECG simulation method that matches the realism of modern heart models. The accuracy and resource requirements were compared to those of a full-torso solution for the potential and scaling was tested up to 14,336 cores with a heart model consisting of 11 million nodes. Reference ECGs were computed on a 3.3 billion-node heart-torso mesh at 0.2 mm resolution. The results show that the lead-field method is more efficient than a full-torso solution when the number of simulated samples is larger than the number of computed ECG leads. While the initial computation of the lead fields remains a hard and poorly scalable problem, the ECG computation itself scales almost perfectly and, even for several hundreds of ECG leads, takes much less time than the underlying simulation of cardiac activity

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

Three-dimensional Multiscale Modelling and Simulation of Atria and Torso Electrophysiology

A better understanding of the electrical activity of the heart under physiological and pathological conditions has always been key for clinicians and researchers. Over the last years, the information in the P-wave signals has been extensively analysed to un-cover the mechanisms underlying atrial arrhythmias by 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 P-wave signals or body surface poten-tial maps are still far from being completely understood. Multiscale anatomical and functional models of the heart are a new technological framework that can enable the investigation of the heart as a complex system. This thesis is centred in the construction of a multiscale framework that allows the realistic simulation of atrial and torso electrophysiology and integrates all the anatom-ical and functional descriptions described in the literature. The construction of such model involves the development of heterogeneous cellular and tissue electrophysiolo-gy models fitted to empirical data. It also requires an accurate 3D representation of the atrial anatomy, including tissue fibre arrangement, and preferential conduction axes. This multiscale model aims to reproduce faithfully the activation of the atria under physiological and pathological conditions. We use the model for two main applica-tions. First, to study the relationship between atrial activation and surface signals in sinus rhythm. This study should reveal the best places for recording P-waves signals in the torso, and which are the regions of the atria that make the most significant contri-bution to the body surface potential maps and determine the main P-wave characteris-tics. Second, to spatially cluster and classify ectopic atrial foci into clearly differenti-ated atrial regions by using the body surface P-wave integral map (BSPiM) as a bi-omarker. We develop a machine-learning pipeline trained from simulations obtained from the atria-torso model aiming to validate whether ectopic foci with similar BSPiM naturally cluster into differentiated non-intersected atrial regions, and whether new BSPiM could be correctly classified with high accuracy.En la actualidad, una mejor compresión de la actividad eléctrica del corazón en condi-ciones fisiológicas y patológicas es clave para médicos e investigadores. A lo largo de los últimos años, la información derivada de la onda P se ha utilizado para intentar descubrir los mecanismos subyacentes a las arritmias auriculares mediante la localiza-ción de focos ectópicos y rotores de alta frecuencia. Sin embargo, la relación entre la activación de distintas regiones auriculares y las características tanto de las ondas P como de la distribución de potencial en la superficie del torso está lejos de entenderse completamente. Los modelos cardíacos funcionales y anatómicos son una nueva he-rramienta que puede facilitar la investigación relativa al corazón entendido como sis-tema complejo. La presente tesis se centra en la construcción de un modelo multiescala para la simula-ción realista de la electrofisiología cardíaca tanto a nivel auricular como de torso, integrando toda la información anatómica y funcional disponible en la literatura. La construcción de este modelo implica el desarrollo, en base a datos experimentales, de modelos electrofisiológicos heterogéneos tanto celulares como tisulares. Así mismo, es imprescindible una representación tridimensional precisa de la anatomía auricular, incluyendo la dirección de fibras y los haces de conducción preferentes. Este modelo multiescala busca reproducir fielmente la activación auricular en condiciones fisiológi-cas y patológicas. Su uso se ha centrado fundamentalmente en dos aplicaciones. En primer lugar, estudiar la relación entre la activación auricular en ritmo sinusal y las señales en la superficie del torso. Este estudio busca definir la mejor ubicación para el registro de las ondas P en el torso así como determinar aquellas regiones auriculares que contribuyen fundamentalmente a la formación y distribución de potenciales super-ficiales así como a las características de las ondas P. En segundo lugar, agrupar y cla-sificar espacialmente los focos ectópicos en regiones auriculares claramente diferen-ciables empleando como biomarcador los mapas superficiales de integral de la onda P (BSPiM). Se ha desarrollado para ello una metodología de aprendizaje automático en la que las simulaciones obtenidas con el modelo multiescala aurícula-torso sirven de entrenamiento, permitiendo validar si los focos ectópicos cuyos BSPiMs son similares se agrupan de forma natural en regiones auriculares no intersectadas y si BSPiMs nue-vos podrían ser clasificados prospectivamente con gran precisión.Avui en dia, una millor comprenssió de l'activitat elèctrica del cor en condicions fisio-lògiques i patològiques és clau per a metges i investigadors. Al llarg dels últims anys, la informació derivada de l'ona P s'ha utilitzat per intentar descobrir els mecanismes subjacents a les arítmies auriculars mitjançant la localització de focus ectòpics i rotors d'alta freqüència. No obstant això, la relació entre l'activació de diferents regions auri-culars i les característiques tant de les ones P com de la distribució de potencial en la superfície del tors està lluny d'entendre's completament. Els models cardíacs funcionals i anatòmics són una nova eina que pot facilitar la recerca relativa al cor entès com a sistema complex. La present tesi es centra en la construcció d'un model multiescala per a la simulació realista de la electrofisiologia cardíaca tant a nivell auricular com de tors, integrant tota la informació anatòmica i funcional disponible en la literatura. La construcció d'aquest model implica el desenvolupament, sobre la base de dades experimentals, de models electrofisiològics heterogenis, tant cel·lulars com tissulars. Així mateix, és imprescindible una representació tridimensional precisa de l'anatomia auricular, in-cloent la direcció de fibres i els feixos de conducció preferents. Aquest model multies-cala busca reproduir fidelment l'activació auricular en condicions fisiològiques i pa-tològiques. El seu ús s'ha centrat fonamentalment en dues aplicacions. En primer lloc, estudiar la relació entre l'activació auricular en ritme sinusal i els senyals en la superfí-cie del tors. A més a més, amb aquest estudi també es busca definir la millor ubicació per al registre de les ones P en el tors, així com, determinar aquelles regions auriculars que contribueixen fonamentalment a la formació i distribució de potencials superfi-cials a l'hora que es caracteritzen les ones P. En segon lloc, agrupar i classificar espa-cialment els focus ectòpics en regions auriculars clarament diferenciables emprant com a biomarcador els mapes superficials d'integral de l'ona P (BSPiM). És per això que s'ha desenvolupat una metodologia d'aprenentatge automàtic en la qual les simulacions obtingudes amb el model multiescala aurícula-tors serveixen d'entrenament, la qual cosa permet validar si els focus ectòpics, llurs BSPiMs són similars, s'agrupen de for-ma natural en regions auriculars no intersectades i si BSPiMs nous podrien ser classifi-cats de manera prospectiva amb precisió.Ferrer Albero, A. (2017). Three-dimensional Multiscale Modelling and Simulation of Atria and Torso Electrophysiology [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/88402TESI

Arrhythmic risk biomarkers for the assessment of drug cardiotoxicity: from experiments to computer simulations

In this paper, we illustrate how advanced computational modelling and simulation can be used to investigate drug-induced effects on cardiac electrophysiology and on specific biomarkers of pro-arrhythmic risk. To do so, we first perform a thorough literature review of proposed arrhythmic risk biomarkers from the ionic to the electrocardiogram levels. The review highlights the variety of proposed biomarkers, the complexity of the mechanisms of drug-induced pro-arrhythmia and the existence of significant animal species differences in drug-induced effects on cardiac electrophysiology. Predicting drug-induced pro-arrhythmic risk solely using experiments is challenging both preclinically and clinically, as attested by the rise in the cost of releasing new compounds to the market. Computational modelling and simulation has significantly contributed to the understanding of cardiac electrophysiology and arrhythmias over the last 40 years. In the second part of this paper, we illustrate how state-of-the-art open source computational modelling and simulation tools can be used to simulate multi-scale effects of drug-induced ion channel block in ventricular electrophysiology at the cellular, tissue and whole ventricular levels for different animal species. We believe that the use of computational modelling and simulation in combination with experimental techniques could be a powerful tool for the assessment of drug safety pharmacology