ECG modeling for simulation of arrhythmias in time-varying conditions

The present paper proposes an ECG simulator that advances modeling of arrhythmias and noise by introducing time-varying signal characteristics. The simulator is built around a discrete-time Markov chain model for simulating atrial and ventricular arrhythmias of particular relevance when analyzing atrial fibrillation (AF). Each state is associated with statistical information on episode duration and heartbeat characteristics. Statistical, time-varying modeling of muscle noise, motion artifacts, and the influence of respiration is introduced to increase the complexity of simulated ECGs, making the simulator well suited for data augmentation in machine learning. Modeling of how the PQ and QT intervals depend on heart rate is also introduced. The realism of simulated ECGs is assessed by three experienced doctors, showing that simulated ECGs are difficult to distinguish from real ECGs. Simulator usefulness is illustrated in terms of AF detection performance when either simulated or real ECGs are used to train a neural network for signal quality control. The results show that both types of training lead to similar performance

Atrial fibrillation dynamics and ionic block effects in six heterogeneous human 3D virtual atria with distinct repolarization dynamics

Atrial fibrillation (AF) usually manifests as reentrant circuits propagating through the whole atria creating chaotic activation patterns. Little is yet known about how differences in electrophysiological and ionic properties between patients modulate reentrant patterns in AF. The goal of this study is to quantify how variability in action potential duration (APD) at different stages of repolarization determines AF dynamics and their modulation by ionic block using a set of virtual whole-atria human models. Six human whole-atria models are constructed based on the same anatomical structure and fiber orientation, but with different electrophysiological phenotypes. Membrane kinetics for each whole-atria model are selected with distinct APD characteristics at 20, 50, and 90% repolarization, from an experimentally calibrated population of human atrial action potential models, including AF remodeling and acetylcholine parasympathetic effects. Our simulations show that in all whole-atria models, reentrant circuits tend to organize around the pulmonary veins and the right atrial appendage, thus leading to higher dominant frequency (DF) and more organized activation in the left atrium than in the right atrium. Differences in APD in all phases of repolarization (not only APD90) yielded quantitative differences in fibrillation patterns with long APD associated with slower and more regular dynamics. Long APD50 and APD20 were associated with increased interatrial conduction block and interatrial differences in DF and organization index, creating reentry instability and self-termination in some cases. Specific inhibitions of IK1, INaK, or INa reduce DF and organization of the arrhythmia by enlarging wave meandering, reducing the number of secondary wavelets, and promoting interatrial block in all six virtual patients, especially for the phenotypes with short APD at 20, 50, and/or 90% repolarization. This suggests that therapies aiming at prolonging the early phase of repolarization might constitute effective antiarrhythmic strategies for the pharmacological management of AF. In summary, simulations report significant differences in atrial fibrillatory dynamics resulting from differences in APD at all phases of repolarization

Influence of heart rate in non-linear HRV indices as a sampling rate effect evaluated on supine and standing

The purpose of this study is to characterize and attenuate the influence of mean heart rate (HR) on nonlinear heart rate variability (HRV) indices (correlation dimension, sample, and approximate entropy) as a consequence of being the HR the intrinsic sampling rate of HRV signal. This influence can notably alter nonlinear HRV indices and lead to biased information regarding autonomic nervous system (ANS) modulation. First, a simulation study was carried out to characterize the dependence of nonlinear HRV indices on HR assuming similar ANS modulation. Second, two HR-correction approaches were proposed: one based on regression formulas and another one based on interpolating RR time series. Finally, standard and HR-corrected HRV indices were studied in a body position change database. The simulation study showed the HR-dependence of non-linear indices as a sampling rate effect, as well as the ability of the proposed HR-corrections to attenuate mean HR influence. Analysis in a body position changes database shows that correlation dimension was reduced around 21% in median values in standing with respect to supine position (p < 0.05), concomitant with a 28% increase in mean HR (p < 0.05). After HR-correction, correlation dimension decreased around 18% in standing with respect to supine position, being the decrease still significant. Sample and approximate entropy showed similar trends. HR-corrected nonlinear HRV indices could represent an improvement in their applicability as markers of ANS modulation when mean HR changes

The Application of Computer Techniques to ECG Interpretation

This book presents some of the latest available information on automated ECG analysis written by many of the leading researchers in the field. It contains a historical introduction, an outline of the latest international standards for signal processing and communications and then an exciting variety of studies on electrophysiological modelling, ECG Imaging, artificial intelligence applied to resting and ambulatory ECGs, body surface mapping, big data in ECG based prediction, enhanced reliability of patient monitoring, and atrial abnormalities on the ECG. It provides an extremely valuable contribution to the field

Liikeartefaktat elektrokardiografiassa

Movement of the patient during electrocardiograph (ECG) recording is a severe source of artifacts. Recent technical developments have enabled ECG recording without continuous supervision by experts. However, ECG recording outside of hospitals is prone to poor quality and movement artifacts. Therefore, it is important to study how and how much ECG recordings are affected by movement. Movement artifacts can hide signal components or mimic them, which causes false negative or false positive detections. Methods to manage movement artifacts include both computational and non-computational approaches. Computational approaches include, for example, adaptive filtering and machine learning methods. Additional variables that correlate with the artifact sources can be utilized in artifact recognition. For example, acceleration, impedance, and pressure signals have been studied as possible movement references. These additional signals are recorded by sensors that are placed on the ECG electrodes or on the patient’s body. In this thesis, the effect of movement artifacts is quantified using a simulation. The simulation makes use of open ECG databases. This study investigates how automated ECG analysis is affected by incremental increase in the movement artifact level. According to the results QRS detection statistics worsen with increased artifact levels. Capturing a movement reference for ECG is studied by experimental research. ECG and inertial measurement unit signals were recorded during different movements in order to analyze the creation of movement artifacts and movement reference signals. According to the results, placement of the movement reference signal sensor has a significant effect on the results. Different movements are captured better by different sensors and affect different ECG leads with different strengths.Potilaan liike sydänsähkökäyrämittauksen (EKG) aikana on merkittävä artefaktien lähde. Viimeaikainen teknologinen kehitys on mahdollistanut EKG-mittauksen ilman asiantuntijoiden jatkuvaa valvontaa. EKG-mittaukset sairaalaolosuhteiden ulkopuolella ovat kuitenkin erityisen alttiita huonolle signaalilaadulle ja liikeartefaktoille. Tämän vuoksi on tärkeää tutkia, miten ja kuinka paljon liike vaikuttaa EKG-mittauksiin. Liikeartefaktat voivat joko peittää tai jäljitellä EKG-signaalin eri osia, aiheuttaen vääriä negatiivisia tai vääriä positiivisia havaintoja. Liikeartefaktojen vaikutusta voidaan vähentää sekä laskennallisten että muiden menetelmien avulla. Laskennallisia menetelmiä ovat esimerkiksi adaptiivinen suodatus ja koneoppimismenetelmät. Artefaktojen lähteen kanssa korreloivia muuttujia mittaamalla voidaan edistää artefaktojen tunnistusta EKG-signaalista. Esimerkiksi kiihtyvyys-, impedanssi- ja painesignaalien käyttöä liikereferensseinä on tutkittu. Kyseisiä referenssisignaaleja voidaan mitata EKG-elektrodeihin tai potilaan kehoon kiinnitettävillä sensoreilla. Liikeartefaktojen vaikutuksen suuruutta tutkitaan tässä työssä simulaation avulla. Simulaatiossa hyödynnetään avoimia EKG-tietokantoja. Tutkimuksessa tarkastellaan sitä, miten vähittäinen liikeartefaktatason kasvu vaikuttaa automaattiseen EKG-analyysiin. Tulosten mukaan QRS-detektioon liittyvät tilastot huononevat artefaktatason kasvaessa. Liikereferenssin luomista tarkastellaan kokeellisen tutkimuksen avulla. EKG- ja inertiamittausyksikkö-signaaleja mitattiin erilaisten liikkeiden aikana, jotta voitaisiin havainnoida liikeartefaktojen ja liikesignaalin syntymistä. Tulosten mukaan liikereferenssiä mittaavan sensorin sijoituspaikalla on merkittävä vaikutus tuloksiin. Tietyt liikkeet saadaan paremmin mitattua eri tavoin sijoitettujen sensorien avulla. Lisäksi liikkeet vaikuttavat eri vahvuuksilla eri EKG-kytkentöihin

Noninvasive autonomic nervous system assessment in respiratory disorders and sport sciences applications

La presente tesis está centrada en el análisis no invasivo de señales cardíacas y respiratorias, con el objetivo de evaluar la actividad del sistema nervioso autónomo (ANS) en diferentes escenarios, tanto clínicos como no clínicos. El documento está estructurado en tres partes principales. La primera parte consiste en una introducción a los aspectos fisiológicos y metodológicos que serán cubiertos en el resto de la tesis. En la segunda parte, se analiza la variabilidad del ritmo cardiaco (HRV) en el contexto de enfermedades respiratorias, concretamente asma (tanto en niños como en adultos) y apnea del sueño. En la tercera parte, se estudian algunas aplicaciones novedosas del análisis de señales cardiorespiratorias en el campo de las ciencias del deporte. La primera parte está compuesta por los capítulos 1 y 2. El capítulo 1 consiste en una extensa introducción al funcionamiento del sistema nervioso autónomo y las características de las bioseñales analizadas a lo largo de la tesis. Por otro lado, se aborda la patofisiología del asma y la apnea del sueño, su relación con el funcionamiento del ANS y las estrategias de diagnóstico y tratamiento de lasmismas. El capítulo concluye con una introducción a la fisiología del ejercicio, así como al interés en la estimación del volumen tidal y del umbral anaeróbico en el campo de las ciencias del deporte.En cuanto al capítulo 2, se presenta un marco de trabajo para el análisis contextualizado de la HRV. Después de una descripción de las técnicas de evaluación y acondicionamiento de la señal de HRV, el capítulo se centra en el efecto de los latidos ectópicos, la arritmia sinusal respiratoria y la frecuencia respiratoria en el análisis de la HRV.Además, se discute el uso de un índice para la evaluación de la distribución de la potencia en los espectros de HRV, así como diferentes medidas de acoplo cardiorespiratorio.La segunda parte está compuesta por los capítulos 3, 4 y 5, todos ellos relacionados con el análisis de la HRV en enfermedades respiratorias. Mientras que los capítulos 3 y 4 están centrados en asma infantil y en adultos respectivamente, el capítulo 5 aborda la apnea del sueño. El asma es una enfermedad respiratoria crónica que aparece habitualmente acompañada por una inflamación de las vías respiratorias. Aunque afecta a personas detodas las edades, normalmente se inicia en edades tempranas, y ha llegado a constituir una de las enfermedades crónicasmás comunes durante la infancia. Sin embargo, todavía no existe un método adecuado para el diagnóstico de asma en niños pequeños. Por otro lado, el rol fundamental que desempeña el sistema nervioso parasimpático en el control del tono bronco-motor y la bronco-dilatación sugiere que la rama parasimpática del ANS podría estar implicada en la patogénesis del asma. De estemodo, en el capítulo 3 se evalúa el ANS mediante el análisis de la HRV en dos bases de datos diferentes, compuestas por niños en edad pre-escolar clasificados en función de su riesgo de desarrollar asma, o de su condición asmática actual. Los resultados del análisis revelaron un balance simpáticovagal reducido y una componente espectral de alta frecuencia más picuda en aquellos niños con un mayor riesgo de desarrollar asma. Además, la actividad parasimpática y el acoplo cardiorespiratorio se redujeron en un grupo de niños con bajo riesgo de asma al finalizar un tratamiento para bronquitis obstructiva, mientras que estos permanecieron inalterados en aquellos niños con una peor prógnosis.A diferencia de los niños pequeños, en el caso de adultos el diagnóstico de asma se realiza a través de una rutina clínica bien definida. Sin embargo, la estratificación de los pacientes en función de su grado de control de los síntomas se basa generalmente en el uso de cuestionarios auto-aplicados, que pueden tener un carácter subjetivo. Por otro lado, la evaluación de la severidad del asma requiere de una visita hospitalaria y de incómodas pruebas, que no pueden aplicarse de una forma continua en el tiempo. De este modo, en el capítulo 4 se estudia el valor de la evaluación del ANS para la estratificación de adultos asmáticos. Para ello, se emplearon diferentes características extraídas de la HRV y la respiración, junto con varios parámetros clínicos, para entrenar un conjunto de algoritmos de clasificación. La inclusión de características relacionadas con el ANS para clasificar los sujetos atendiendo a la severidad del asma derivó en resultados similares al caso de utilizar únicamente parámetros clínicos, superando el desempeño de estos últimos en algunos casos. Por lo tanto, la evaluación del ANS podría representar un potencial complemento para la mejora de la monitorización de sujetos asmáticos.En el capítulo 5, se analiza la HRV en sujetos que padecen el síndrome de apnea del sueño (SAS) y comorbididades cardíacas asociadas. El SAS se ha relacionado con un incremento de 5 veces en el riesgo de desarrollar enfermedades cardiovasculares (CVD), que podría aumentar hasta 11 veces si no se trata convenientemente. Por otro lado, una HRV alterada se ha relacionado independientemente con el SAS y con numerosos factores de riesgo para el desarrollo de CVD. De este modo, este capítulo se centra en evaluar si una actividad autónoma desbalanceada podría estar relacionada con el desarrollo de CVD en pacientes de SAS. Los resultados del análisis revelaron una dominancia simpática reducida en aquellos sujetos que padecían SAS y CVD, en comparación con aquellos sin CVD. Además, un análisis retrospectivo en una base de datos de sujetos con SAS que desarollarán CVD en el futuro también reveló una actividad simpática reducida, sugiriendo que un ANS desbalanceado podría constituir un factor de riesgo adicional para el desarrollo de CVD en pacientes de SAS.La tercera parte está formada por los capítulos 6 y 7, y está centrada en diferentes aplicaciones del análisis de señales cardiorespiratorias en el campo de las ciencias del deporte. El capítulo 6 aborda la estimación del volumen tidal (TV) a partir del electrocardiograma (ECG). A pesar de que una correcta monitorización de la actividad respiratoria es de gran interés en ciertas enfermedades respiratorias y en ciencias del deporte, la mayor parte de la actividad investigadora se ha centrado en la estimación de la frecuencia respiratoria, con sólo unos pocos estudios centrados en el TV, la mayoría de los cuales se basan en técnicas no relacionadas con el ECG. En este capítulo se propone un marco de trabajo para la estimación del TV en reposo y durante una prueba de esfuerzo en tapiz rodante utilizando únicamente parámetros derivados del ECG. Errores de estimación del 14% en la mayoría de los casos y del 6% en algunos sugieren que el TV puede estimarse a partir del ECG, incluso en condiciones no estacionarias.Por último, en el capítulo 7 se propone una metodología novedosa para la estimación del umbral anaeróbico (AT) a partir del análisis de las dinámicas de repolarización ventricular. El AT representa la frontera a partir de la cual el sistema cardiovascular limita la actividad física de resistencia, y aunque fue inicialmente concebido para la evaluación de la capacidad física de pacientes con CVD, también resulta de gran interés en el campo de las ciencias del deporte, permitiendo diseñar mejores rutinas de entrenamiento o para prevenir el sobre-entrenamiento. Sin embargo, la evaluación del AT requiere de técnicas invasivas o de dispositivos incómodos. En este capítulo, el AT fue estimado a partir del análisis de las variaciones de las dinámicas de repolarización ventricular durante una prueba de esfuerzo en cicloergómetro. Errores de estimación de 25 W, correspondientesa 1 minuto en este estudio, en un 63% de los sujetos (y menores que 50 W en un 74% de ellos) sugieren que el AT puede estimarse de manera no invasiva, utilizando únicamente registros de ECG.<br /

Extraction and Detection of Fetal Electrocardiograms from Abdominal Recordings

The non-invasive fetal ECG (NIFECG), derived from abdominal surface electrodes, offers novel diagnostic possibilities for prenatal medicine. Despite its straightforward applicability, NIFECG signals are usually corrupted by many interfering sources. Most significantly, by the maternal ECG (MECG), whose amplitude usually exceeds that of the fetal ECG (FECG) by multiple times. The presence of additional noise sources (e.g. muscular/uterine noise, electrode motion, etc.) further affects the signal-to-noise ratio (SNR) of the FECG. These interfering sources, which typically show a strong non-stationary behavior, render the FECG extraction and fetal QRS (FQRS) detection demanding signal processing tasks. In this thesis, several of the challenges regarding NIFECG signal analysis were addressed. In order to improve NIFECG extraction, the dynamic model of a Kalman filter approach was extended, thus, providing a more adequate representation of the mixture of FECG, MECG, and noise. In addition, aiming at the FECG signal quality assessment, novel metrics were proposed and evaluated. Further, these quality metrics were applied in improving FQRS detection and fetal heart rate estimation based on an innovative evolutionary algorithm and Kalman filtering signal fusion, respectively. The elaborated methods were characterized in depth using both simulated and clinical data, produced throughout this thesis. To stress-test extraction algorithms under ideal circumstances, a comprehensive benchmark protocol was created and contributed to an extensively improved NIFECG simulation toolbox. The developed toolbox and a large simulated dataset were released under an open-source license, allowing researchers to compare results in a reproducible manner. Furthermore, to validate the developed approaches under more realistic and challenging situations, a clinical trial was performed in collaboration with the University Hospital of Leipzig. Aside from serving as a test set for the developed algorithms, the clinical trial enabled an exploratory research. This enables a better understanding about the pathophysiological variables and measurement setup configurations that lead to changes in the abdominal signal's SNR. With such broad scope, this dissertation addresses many of the current aspects of NIFECG analysis and provides future suggestions to establish NIFECG in clinical settings.:Abstract Acknowledgment Contents List of Figures List of Tables List of Abbreviations List of Symbols (1)Introduction 1.1)Background and Motivation 1.2)Aim of this Work 1.3)Dissertation Outline 1.4)Collaborators and Conflicts of Interest (2)Clinical Background 2.1)Physiology 2.1.1)Changes in the maternal circulatory system 2.1.2)Intrauterine structures and feto-maternal connection 2.1.3)Fetal growth and presentation 2.1.4)Fetal circulatory system 2.1.5)Fetal autonomic nervous system 2.1.6)Fetal heart activity and underlying factors 2.2)Pathology 2.2.1)Premature rupture of membrane 2.2.2)Intrauterine growth restriction 2.2.3)Fetal anemia 2.3)Interpretation of Fetal Heart Activity 2.3.1)Summary of clinical studies on FHR/FHRV 2.3.2)Summary of studies on heart conduction 2.4)Chapter Summary (3)Technical State of the Art 3.1)Prenatal Diagnostic and Measuring Technique 3.1.1)Fetal heart monitoring 3.1.2)Related metrics 3.2)Non-Invasive Fetal ECG Acquisition 3.2.1)Overview 3.2.2)Commercial equipment 3.2.3)Electrode configurations 3.2.4)Available NIFECG databases 3.2.5)Validity and usability of the non-invasive fetal ECG 3.3)Non-Invasive Fetal ECG Extraction Methods 3.3.1)Overview on the non-invasive fetal ECG extraction methods 3.3.2)Kalman filtering basics 3.3.3)Nonlinear Kalman filtering 3.3.4)Extended Kalman filter for FECG estimation 3.4)Fetal QRS Detection 3.4.1)Merging multichannel fetal QRS detections 3.4.2)Detection performance 3.5)Fetal Heart Rate Estimation 3.5.1)Preprocessing the fetal heart rate 3.5.2)Fetal heart rate statistics 3.6)Fetal ECG Morphological Analysis 3.7)Problem Description 3.8)Chapter Summary (4)Novel Approaches for Fetal ECG Analysis 4.1)Preliminary Considerations 4.2)Fetal ECG Extraction by means of Kalman Filtering 4.2.1)Optimized Gaussian approximation 4.2.2)Time-varying covariance matrices 4.2.3)Extended Kalman filter with unknown inputs 4.2.4)Filter calibration 4.3)Accurate Fetal QRS and Heart Rate Detection 4.3.1)Multichannel evolutionary QRS correction 4.3.2)Multichannel fetal heart rate estimation using Kalman filters 4.4)Chapter Summary (5)Data Material 5.1)Simulated Data 5.1.1)The FECG Synthetic Generator (FECGSYN) 5.1.2)The FECG Synthetic Database (FECGSYNDB) 5.2)Clinical Data 5.2.1)Clinical NIFECG recording 5.2.2)Scope and limitations of this study 5.2.3)Data annotation: signal quality and fetal amplitude 5.2.4)Data annotation: fetal QRS annotation 5.3)Chapter Summary (6)Results for Data Analysis 6.1)Simulated Data 6.1.1)Fetal QRS detection 6.1.2)Morphological analysis 6.2)Own Clinical Data 6.2.1)FQRS correction using the evolutionary algorithm 6.2.2)FHR correction by means of Kalman filtering (7)Discussion and Prospective 7.1)Data Availability 7.1.1)New measurement protocol 7.2)Signal Quality 7.3)Extraction Methods 7.4)FQRS and FHR Correction Algorithms (8)Conclusion References (A)Appendix A - Signal Quality Annotation (B)Appendix B - Fetal QRS Annotation (C)Appendix C - Data Recording GU

Statistical Coding and Decoding of Heartbeat Intervals

The heart integrates neuroregulatory messages into specific bands of frequency, such that the overall amplitude spectrum of the cardiac output reflects the variations of the autonomic nervous system. This modulatory mechanism seems to be well adjusted to the unpredictability of the cardiac demand, maintaining a proper cardiac regulation. A longstanding theory holds that biological organisms facing an ever-changing environment are likely to evolve adaptive mechanisms to extract essential features in order to adjust their behavior. The key question, however, has been to understand how the neural circuitry self-organizes these feature detectors to select behaviorally relevant information. Previous studies in computational perception suggest that a neural population enhances information that is important for survival by minimizing the statistical redundancy of the stimuli. Herein we investigate whether the cardiac system makes use of a redundancy reduction strategy to regulate the cardiac rhythm. Based on a network of neural filters optimized to code heartbeat intervals, we learn a population code that maximizes the information across the neural ensemble. The emerging population code displays filter tuning proprieties whose characteristics explain diverse aspects of the autonomic cardiac regulation, such as the compromise between fast and slow cardiac responses. We show that the filters yield responses that are quantitatively similar to observed heart rate responses during direct sympathetic or parasympathetic nerve stimulation. Our findings suggest that the heart decodes autonomic stimuli according to information theory principles analogous to how perceptual cues are encoded by sensory systems

Modelling the interaction between induced pluripotent stem cells derived cardiomyocytes patches and the recipient hearts

Cardiovascular diseases are the main cause of death worldwide. The single biggest killer is represented by ischemic heart disease. Myocardial infarction causes the formation of non-conductive and non-contractile, scar-like tissue in the heart, which can hamper the heart's physiological function and cause pathologies ranging from arrhythmias to heart failure. The heart can not recover the tissue lost due to myocardial infarction due to the myocardium's limited ability to regenerate. The only available treatment is heart transpalant, which is limited by the number of donors and can elicit an adverse response from the recipients immune system. Recently, regenerative medicine has been proposed as an alternative approach to help post-myocardial infarction hearts recover their functionality. Among the various techniques, the application of cardiac patches of engineered heart tissue in combination with electroactive materials constitutes a promising technology. However, many challenges need to be faced in the development of this treatment. One of the main concerns is represented by the immature phenotype of the stem cells-derived cardiomyocytes used to fabricate the engineered heart tissue. Their electrophysiological differences with respect to the host myocardium may contribute to an increased arrhythmia risk. A large number of animal experiments are needed to optimize the patches' characteristics and to better understand the implications of the electrical interaction between patches and host myocardium. In this Thesis we leveraged cardiac computational modelling to simulate \emph{in silico} electrical propagation in scarred heart tissue in the presence of a patch of engineered heart tissue and conductive polymer engrafted at the epicardium. This work is composed by two studies. In the first study we designed a tissue model with simplified geometry and used machine learning and global sensitivity analysis techniques to identify engineered heart tissue patch design variables that are important for restoring physiological electrophysiology in the host myocardium. Additionally, we showed how engineered heart tissue properties could be tuned to restore physiological activation while reducing arrhythmic risk. In the second study we moved to more realistic geometries and we devised a way to manipulate ventricle meshes obtained from magnetic resonance images to apply \emph{in silico} engineered heart tissue epicardial patches. We then investigated how patches with different conduction velocity and action potential duration influence the host ventricle electrophysiology. Specifically, we showed that appropriately located patches can reduce the predisposition to anatomical isthmus mediated re-entry and that patches with a physiological action potential duration and higher conduction velocity were most effective in reducing this risk. We also demonstrated that patches with conduction velocity and action potential duration typical of immature stem cells-derived cardiomyocytes were associated with the onset of sustained functional re-entry in an ischemic cardiomyopathy model with a large transmural scar. Finally, we demonstrated that patches electrically coupled to host myocardium reduce the likelihood of propagation of focal ectopic impulses. This Thesis demonstrates how computational modelling can be successfully applied to the field of regenerative medicine and constitutes the first step towards the creation of patient-specific models for developing and testing patches for cardiac regeneration.Open Acces