Search

When pathogens enter into the body, the natural defense reaction is an immunological process triggered by the production of cytokines by macrophages. This reaction is called inflammation and it is an important and complex physiological response of living tissues to all harmful damages. As a consequence of the immune response, the most common symptoms are edema, redness, fever, and pain. All of these four symptoms are related to the physiological process of the immune defense system. For instance, edema may be caused by increased blood vessel wall permeability which leads to a swell at the infection site. So, edema is a consequence of interstitial fluid dynamics and its interactions with the immune system. In this present study, this physiological process is mathematically described by a nonlinear system of partial differential equations (PDE) based on a porous media approach. The model describes the interaction between the neutrophils and the pathogen, which in our case is non-specific bacteria. In addition, a poroelastic model is coupled with the inflammatory one resulting in an hydro-mechanical model for inflammatory edema. The new suggested model was able to reproduce key aspects of an inflammatory edema formation. [ABSTRACT FROM AUTHOR]

Arrhythmias in cardiac tissue are generally associated with irregular electrical wave propagation in the heart. Cardiac tissue is formed by a discrete cell network, which is often heterogeneous. Recently, it was shown in simulations of two-dimensional (2D) discrete models of cardiac tissue that a wave crossing a fibrotic, heterogeneous region may produce reentry and transient or persistent ectopic activity provided the fraction of conducting connections is just above the percolation threshold. Here, we investigate the occurrence of these phenomena in three-dimensions by simulations of a discrete model representing a thin slab of cardiac tissue. This is motivated (i) by the necessity to study the relevance and properties of the percolation-related mechanism for the emergence of microreentries in three dimensions and (ii) by the fact that atrial tissue is quite thin in comparison with ventricular tissue. Here, we simplify the model by neglecting details of tissue anatomy, e. g. geometries of atria or ventricles and the anisotropy in the conductivity. Hence, our modeling study is confined to the investigation of the effect of the tissue thickness as well as to the comparison of the dynamics of electrical excitation in a 2D layer with the one in a 3D slab. Our results indicate a strong and non-trivial effect of the thickness even for thin tissue slabs on the probability of microreentries and ectopic beat generation. The strong correlation of the occurrence of microreentry with the percolation threshold reported earlier in 2D layers persists in 3D slabs. Finally, a qualitative agreement of 3D simulated electrograms in the fibrotic region with the experimentally observed complex fractional atrial electrograms (CFAE) as well as strong difference between simulated electrograms in 2D and 3D were found for the cases where reentry and ectopic activity were triggered by the micro-fibrotic region. [ABSTRACT FROM AUTHOR]

This paper uses recurrence quantification analysis (RQA) combined with entropy measures and organization indices to characterize arrhythmic patterns and dynamics in computer simulations of cardiac tissue. We performed different simulations of cardiac tissues of sizes comparable to the human heart atrium. In these simulations, we observed four classic arrhythmic patterns: a spiral wave anchored to a highly fibrotic region resulting in sustained re-entry, a meandering spiral wave, fibrillation, and a spiral wave anchored to a scar region that breaks up into wavelets away from the main rotor. A detailed analysis revealed that, within the same simulation, maps of RQA metrics could differentiate regions with regular AP propagation from ones with chaotic activity. In particular, the combination of two RQA metrics, the length of the longest diagonal string of recurrence points and the mean length of diagonal lines, was able to identify the location of rotor tips, which are the active elements that maintain spiral waves and fibrillation. By proposing low-dimensional models based on the mean value and spatial correlation of metrics calculated from membrane potential time series, we identify RQA-based metrics that successfully separate the four different types of cardiac arrhythmia into distinct regions of the feature space, and thus might be used for automatic classification, in particular distinguishing between fibrillation driven by self-sustaining chaos and that created by a persistent rotor and wavebreak. We also discuss the practical applicability of such an approach. [ABSTRACT FROM AUTHOR]

The electrical activity of adult mouse and rat hearts has been analyzed extensively, often as a prerequisite for genetic engineering studies or for the development of rodent models of human diseases. Some aspects of the initiation and conduction of the cardiac action potential in rodents closely resemble those in large mammals. However, rodents have a much higher heart rate and their ventricular action potential is triangular and very short. As a consequence, an interpretation of the electrocardiogram in the mouse and rat remains difficult and controversial. In this study, optical mapping techniques have been applied to an in vitro left ventricular adult rat preparation to obtain patterns of conduction and action potential duration measurements from the epicardial surface. This information has been combined with previously published mathematical models of the rat ventricular myocyte to develop a bidomain model for action potential propagation and electrogram formation in the rat left ventricle. Important insights into the basis for the repolarization waveform in the ventricular electrogram of the adult rat have been obtained. Notably, our model demonstrated that the biphasic shape of the rat ventricular repolarization wave can be explained in terms of the transmural and apex-to-base gradients in action potential duration that exist in the rat left ventricle. [ABSTRACT FROM AUTHOR]

The bidomain equations are widely used for the simulation of electrical activity in cardiac tissue but are computationally expensive, limiting the size of the problem which can be modeled. The purpose of this study is to determine more efficient ways to solve the elliptic portion of the bidomain equations, the most computationally expensive part of the computation. Specifically, we assessed the performance of a parallel multigrid (MG) preconditioner for a conjugate gradient solver. We employed an operator splitting technique, dividing the computation in a parabolic equation, an elliptical equation, and a nonlinear system of ordinary differential equations at each time step. The elliptic equation was solved by the preconditioned conjugate gradient method, and the traditional block incomplete LU parallel preconditioner (ILU) was compared to MG. Execution time was minimized for each put- conditioner by adjusting the fill-in factor for ILU, and by choosing the optimal number of levels for MG. The parallel implementation was based on the PETSc library and we report results for up to 16 nodes on a distributed cluster, for two and three dimensional simulations. A direct solver was also available to compare results for single processor runs. MG was found to solve the system in one third of the time required by ILU but required about 40% more memory. Thus, MG offered an attractive tradeoff between memory usage and speed, since its performance lay between those of the classic iterative methods (slow and low memory consumption) and direct methods (fast and high memory consumption). Results suggest the MG preconditioner is well suited for quickly and accurately solving the bidomain equations. [ABSTRACT FROM AUTHOR]

Abstract: Computational modeling of the heart is a subject of substantial medical and scientific interest, which may contribute to increase the understanding of several phenomena associated with cardiac physiological and pathological states. Modeling the mechanics of the heart have led to considerable insights, but it still represents a complex and a demanding computational problem, especially in a strongly coupled electromechanical setting. Passive cardiac tissue is commonly modeled as hyperelastic and is characterized by quasi‐incompressible, orthotropic, and nonlinear material behavior. These factors are known to be very challenging for the numerical solution of the model. The near‐incompressibility is known to cause numerical issues such as the well‐known locking phenomenon and ill‐conditioning of the stiffness matrix. In this work, the augmented Lagrangian method is used to handle the nearly incompressible condition. This approach can potentially improve computational performance by reducing the condition number of the stiffness matrix and thereby improving the convergence of iterative solvers. We also improve the performance of iterative solvers by the use of an algebraic multigrid preconditioner. Numerical results of the augmented Lagrangian method combined with a preconditioned iterative solver for a cardiac mechanics benchmark suite are presented to show its improved performance. [ABSTRACT FROM AUTHOR]

Cardiac Purkinje networks are a fundamental part of the conduction system and are known to initiate a variety of cardiac arrhythmias. However, patient-specific modeling of Purkinje networks remains a challenge due to their high morphological complexity. This work presents a novel method based on optimization principles for the generation of Purkinje networks that combines geometric and activation accuracy in branch size, bifurcation angles, and Purkinje-ventricular-junction activation times. Three biventricular meshes with increasing levels of complexity are used to evaluate the performance of our approach. Purkinje-tissue coupled monodomain simulations are executed to evaluate the generated networks in a realistic scenario using the most recent Purkinje/ventricular human cellular models and physiological values for the Purkinje-ventricular-junction characteristic delay. The results demonstrate that the new method can generate patient-specific Purkinje networks with controlled morphological metrics and specified local activation times at the Purkinje-ventricular junctions. [ABSTRACT FROM AUTHOR]

Introduction: The apparent location of the myocytes (M cells) with the longest action potential duration (APD) in a canine left ventricular (LV) wedge have been reported to shift after application of a sea anemone toxin, ATX-II. This toxin slows inactivation of INa and thus prolongs APD. Thus, M cells may exhibit dynamic functional states, rather than being a static, anatomically discrete, myocyte population. In this study, we attempted to further define and understand this phenomenon using a mathematical model of the human ventricular myocyte action potential incorporated into an in silico “wedge” preparation. Our simulations demonstrate that even under conditions of a fixed population and ratio of epicardial, M, and endocardial myocytes, the apparent anatomical position (transmural location) of the myocytes with the longest APD can shift following ATX-II treatment. This arises because the ATX-II effect, modeled as a small increase in the late or persistent Na+ current, and consequent prolongation of APD significantly changes the electrotonic interactions between ventricular myocytes in this LV wedge preparation. Methods and Results: This LV wedge model is based on bidomain equations. It corresponds to a rectangular tissue immersed in a passive and isotropic medium that represents the superfusion bath. In this theoretical work, the three known different and discrete populations of myocytes in the human left ventricle have been included: the epicardial, M, and endocardial cells. The effects of ATX-II on INa were simulated by altering the voltage-dependent steady-state inactivation of the parameters h (fast gate) and j (slow gate). As a result, in these ATX-II simulations a persistent late Na+ current was generated in all three types of ventricular myocytes. However, the APDs were prolonged in a heterogeneous pattern. Our simulations demonstrate that after the ATX-II effects develop, alterations in transmural electrotonic interactions can produce changes in the transmural location of myocytes with the longest APD. Conclusions: The combination of intercellular electrotonic interactions, which tend to reduce and smooth out the discrete transmural APD variations, and the heterogeneous effects of ATX-II, which preferentially prolong the APD of M cells, can shift the location of the ventricular myocytes. This shift results in significantly altered transmural patterns of action potential durations, which would be expected to change localized refractory period and excitability. These cellular changes give rise to alterations in the corresponding surface electrograms and may change the overall substrates for conduction and rhythm disturbances. [ABSTRACT FROM AUTHOR]

The development of modeling structures at the channel level that can integrate subcellular and cell models and properly reproduce different experimental data is of utmost importance in cardiac electrophysiology. In contrast to gate-based models, Markov Chain models are well suited to promote the integration of the subcellular level of the cardiomyocyte to the whole cell. In this paper, we develop Markov Chain models for the L-type Calcium current that can reproduce the electrophysiology of two established human models for the ventricular and Purkinje cells. In addition, instead of presenting a single set of parameters, we present a collection of set of parameters employing Differential Evolution algorithms that can properly reproduce very different protocol data. We show the importance of using an ensemble of a set of parameter values to obtain proper results when considering a second protocol that suppresses calcium inactivation and mimics a pathological condition. We discuss how model discrepancy, data availability, and parameter identifiability can influence the choice of the size of the collection. In summary, we have modified two cardiac models by proposing new Markov Chain models for the L-type Calcium. We keep the original whole-cell dynamics by reproducing the same characteristic action potential and calcium dynamics, whereas the Markov chain-based description of the L-type Calcium channels allows novel small spatial scale simulations of subcellular processes. Finally, the use of collections of parameters was crucial for addressing model discrepancy, identifiability issues, and avoiding fitting parameters overly precisely, i.e., overfitting. [ABSTRACT FROM AUTHOR]

In recent years, there has been an increasing interest in the mathematical and computational modeling of the human immune system (HIS). Computational models of HIS dynamics may contribute to a better understanding of the relationship between complex phenomena and immune response; in addition, computational models will support the development of new drugs and therapies for different diseases. However, modeling the HIS is an extremely difficult task that demands a huge amount of work to be performed by multidisciplinary teams. In this study, our objective is to model the spatio-temporal dynamics of representative cells and molecules of the HIS during an immune response after the injection of lipopolysaccharide (LPS) into a section of tissue. LPS constitutes the cellular wall of Gram-negative bacteria, and it is a highly immunogenic molecule, which means that it has a remarkable capacity to elicit strong immune responses. We present a descriptive, mechanistic and deterministic model that is based on partial differential equations (PDE). Therefore, this model enables the understanding of how the different complex phenomena interact with structures and elements during an immune response. In addition, the model's parameters reflect physiological features of the system, which makes the model appropriate for general use. [ABSTRACT FROM AUTHOR]

Bacterial infections can be of two types: acute or chronic. The chronic bacterial infections are characterized by being a large bacterial infection and/or an infection where the bacteria grows rapidly. In these cases, the immune response is not capable of completely eliminating the infection which may lead to the formation of a pattern known as microabscess (or abscess). The microabscess is characterized by an area comprising fluids, bacteria, immune cells (mainly neutrophils), and many types of dead cells. This distinct pattern of formation can only be numerically reproduced and studied bymodels that capture the spatiotemporal dynamics of the human immune system (HIS). In this context, our work aims to develop and implement an initial computational model to study the process of microabscess formation during a bacterial infection. [ABSTRACT FROM AUTHOR]

The bidomain equations are considered to be one of the most complete descriptions of the electrical activity in cardiac tissue, but large scale simulations, as resulting from discretization of an entire heart, remain a computational challenge due to the elliptic portion of the problem, the part associated with solving the extracellular potential. In such cases, the use of iterative solvers and parallel computing environments are mandatory to make parameter studies feasible. The preconditioned conjugate gradient (PCG) method is a standard choice for this problem. Although robust, its efficiency greatly depends on the choice of preconditioner. On structured grids, it has been demonstrated that a geometric multigrid preconditioner performs significantly better than an incomplete LU (ILU) preconditioner. However, unstructured grids are often preferred to better represent organ boundaries and allow for coarser discretization in the bath far from cardiac surfaces. Under these circumstances, algebraic multigrid (AMG) methods are advantageous since they compute coarser levels directly from the system matrix itself, thus avoiding the complexity of explicitly generating coarser, geometric grids. In this paper, the performance of an AMG preconditioner (BoomerAMG) is compared with that of the standard ILU preconditioner and a direct solver. BoomerAMG is used in two different ways, as a preconditioner and as a standalone solver. Two 3-D simulation examples modeling the induction of arrhythmias in rabbit ventricles were used to measure performance in both sequential and parallel simulations. It is shown that the AMG preconditioner is very well suited for the solution of the bidomain equation, being clearly superior to ILU preconditioning in all regards, with speedups by factors in the range 5.9-7.7. [ABSTRACT FROM AUTHOR]

Background: An effective yellow fever (YF) vaccine has been available since 1937. Nevertheless, questions regarding its use remain poorly understood, such as the ideal dose to confer immunity against the disease, the need for a booster dose, the optimal immunisation schedule for immunocompetent, immunosuppressed, and pediatric populations, among other issues. This work aims to demonstrate that computational tools can be used to simulate different scenarios regarding YF vaccination and the immune response of individuals to this vaccine, thus assisting the response of some of these open questions. Results: This work presents the computational results obtained by a mathematical model of the human immune response to vaccination against YF. Five scenarios were simulated: primovaccination in adults and children, booster dose in adult individuals, vaccination of individuals with autoimmune diseases under immunomodulatory therapy, and the immune response to different vaccine doses. Where data were available, the model was able to quantitatively replicate the levels of antibodies obtained experimentally. In addition, for those scenarios where data were not available, it was possible to qualitatively reproduce the immune response behaviours described in the literature. Conclusions: Our simulations show that the minimum dose to confer immunity against YF is half of the reference dose. The results also suggest that immunological immaturity in children limits the induction and persistence of long-lived plasma cells are related to the antibody decay observed experimentally. Finally, the decay observed in the antibody level after ten years suggests that a booster dose is necessary to keep immunity against YF. [ABSTRACT FROM AUTHOR]

Background: Myocarditis is defined as the inflammation of the myocardium, i.e. the cardiac muscle. Among the reasons that lead to this disease, we may include infections caused by a virus, bacteria, protozoa, fungus, and others. One of the signs of the inflammation is the formation of edema, which may be a consequence of the interaction between interstitial fluid dynamics and immune response. This complex physiological process was mathematically modeled using a nonlinear system of partial differential equations (PDE) based on porous media approach. By combing a model based on Biot’s poroelasticity theory with a model for the immune response we developed a new hydro-mechanical model for inflammatory edema. To verify this new computational model, T2 parametric mapping obtained by Magnetic Resonance (MR) imaging was used to identify the region of edema in a patient diagnosed with unspecific myocarditis. Results: A patient-specific geometrical model was created using MRI images from the patient with myocarditis. With this model, edema formation was simulated using the proposed hydro-mechanical mathematical model in a two-dimensional domain. The computer simulations allowed us to correlate spatiotemporal dynamics of representative cells of the immune systems, such as leucocytes and the pathogen, with fluid accumulation and cardiac tissue deformation. Conclusions: This study demonstrates that the proposed mathematical model is a very promising tool to better understand edema formation in myocarditis. Simulations obtained from a patient-specific model reproduced important aspects related to the formation of cardiac edema, its area, position, and shape, and how these features are related to immune response. [ABSTRACT FROM AUTHOR]

Pattern formation has been widely observed in extended chemical and biological processes. Although the biochemical systems are highly heterogeneous, homogenized continuum approaches formed by partial differential equations have been employed frequently. Such approaches are usually justified by the difference of scales between the heterogeneities and the characteristic spatial size of the patterns. Under different conditions, for example, under weak coupling, discrete models are more adequate. However, discrete models may be less manageable, for instance, in terms of numerical implementation and mesh generation, than the associated continuum models. Here we study a model to approach discreteness which permits the computer implementation on general unstructured meshes. The model is cast as a partial differential equation but with a parameter that depends not only on heterogeneities sizes, as in the case of quasicontinuum models, but also on the discretization mesh. Therefore, we refer to it as a discretization-dependent model. We validate the approach in a generic excitable media that simulates three different phenomena: the propagation of action membrane potential in cardiac tissue, in myelinated axons of neurons, and concentration waves in chemical microemulsions. [ABSTRACT FROM AUTHOR]

Cardiac electrical alternans is a period-2 dynamical behavior with alternating long and short action potential durations (APD) that often precedes dangerous arrhythmias associated with cardiac arrest. Despite the importance of alternans, many current ordinary differential equations models of cardiac electrophysiology do not produce alternans, thereby limiting the use of these models for studying the mechanisms that underlie this condition. Because delay differential equations (DDEs) commonly induce complex dynamics in other biological systems, we investigate whether incorporating DDEs can lead to alternans development in cardiac models by studying the Fox et al. canine ventricular action potential model. After suppressing the alternans in the original model, we show that alternans can be obtained by introducing DDEs in the model gating variables, and we quantitatively compare the DDE-induced alternans with the alternans present in the original model. We analyze the behavior of the voltage, currents, and gating variables of the model to study the effects of the delays and to determine how alternans develops in that setting, and we discuss the mathematical and physiological implications of our findings. In future work, we aim to apply our approach to induce alternans in models that do not naturally exhibit such dynamics. [ABSTRACT FROM AUTHOR]

Foam injection in porous media is often used to control the gas fingering in multi-phase flow. Mathematical models of foam dynamics involve non-Newtonian formulations. To numerically simulate these complex phenomena, experimental data is gathered and used to estimate the parameter values of models via optimization techniques. The present work improves this procedure by introducing a new objective function based on the mobility reduction factor and does not require further experimental observations other than those usually obtained in core-flooding experiments. A series of numerical experiments were carried out to show the features and robustness of the proposed approach. In particular, the identifiability analysis results show that key parameters of the models are practically non-identifiable when using traditional objective functions that rely only on apparent viscosity. When the proposed objective function is used for parameter estimation, the identifiability issues are solved. In addition, the uncertainty quantification analysis revealed that lower uncertainty is achieved when using the new objective function when compared to the one that uses only apparent viscosity. In summary, we show that the new objective function generates better-calibrated models with high fidelity and low uncertainties and alleviates parameter non-identifiability issues. Therefore, the results suggest the new objective function is better suited for calibrating foam displacement models for enhanced oil recovery. • A new objective function is proposed for the estimation of parameters for models of non-Newtonian foam dynamics. • Inverse and forward uncertainty quantification, sensitivity analysis, and identifiability methods assessed the new objective function. • The minimization of the new objective function generates better-calibrated models with high fidelity and low uncertainties. [ABSTRACT FROM AUTHOR]

Advanced medical imaging technologies provide a wealth of information on cardiac anatomy and structure at a paracellular resolution, allowing to identify microstructural discontinuities which disrupt the intracellular matrix. Current state-of-the-art computer models built upon such datasets account for increasingly finer anatomical details, however, structural discontinuities at the paracellular level are typically discarded in the model generation process, owing to the significant costs which incur when using high resolutions for explicit representation. In this study, a novel discontinuous finite element (dFE) approach for discretizing the bidomain equations is presented, which accounts for fine-scale structures in a computer model without the need to increase spatial resolution. In the dFE method, this is achieved by imposing infinitely thin lines of electrical insulation along edges of finite elements which approximate the geometry of discontinuities in the intracellular matrix. Simulation results demonstrate that the dFE approach accounts for effects induced by microscopic size scale discontinuities, such as the formation of microscopic virtual electrodes, with vast computational savings as compared to high resolution continuous finite element models. Moreover, the method can be implemented in any standard continuous finite element code with minor effort. [ABSTRACT FROM PUBLISHER]

Fibrosis is thought to play an important role in the formation and maintenance of atrial fibrillation (AF). The propensity of fibrosis to increase AF vulnerability depends not only on its amount, its texture plays a crucial role as well. While the detection of fibrotic tissue patches in the atria with extracellular recordings is feasible based on the analysis of electrogram fractionation, as used in clinical practice to identify ablation targets, the classification of fibrotic texture is a more challenging problem. This study seeks to establish a method for the electroanatomical characterization of the fibrotic textures based on the analysis of electrogram fractionation. The proposed method exploits the dependence of fractionation patterns on the incidence direction of wavefronts which differs significantly as a function of texture. A histologically detailed computer model of the right atrial isthmus was developed for testing the method. A stimulation protocol was conceived which generated various incidence directions for any given recording site where electrograms were computed. A classification method is derived then for discriminating three types of fibrosis, no fibrosis (control), diffuse, and patchy fibrosis. Simulation results showed that electrogram fractionation and amplitudes and their dependence upon incidence direction allow a robust discrimination between different classes of fibrosis. Finally, to minimize the technical effort, sensitivity analysis was performed to identify a minimum number of incidence directions required for robust classification. [ABSTRACT FROM PUBLISHER]

Key aspects of cardiac electrophysiology, such as slow conduction, conduction block, and saltatory effects have been the research topic of many studies since they are strongly related to cardiac arrhythmia, reentry, fibrillation, or defibrillation. However, to reproduce these phenomena the numerical models need to use subcellular discretization for the solution of the PDEs and nonuniform, heterogeneous tissue electric conductivity. Due to the high computational costs of simulations that reproduce the fine microstructure of cardiac tissue, previous studies have considered tissue experiments of small or moderate sizes and used simple cardiac cell models. In this paper, we develop a cardiac electrophysiology model that captures the microstructure of cardiac tissue by using a very fine spatial discretization (8 µm) and uses a very modern and complex cell model based on Markov chains for the characterization of ion channel's structure and dynamics. To cope with the computational challenges, the model was parallelized using a hybrid approach: cluster computing and GPGPUs (general-purpose computing on graphics processing units). Our parallel implementation of this model using a multi-GPU platform was able to reduce the execution times of the simulations from more than 6 days (on a single processor) to 21 minutes (on a small 8-node cluster equipped with 16 GPUs, i.e., 2 GPUs per node). [ABSTRACT FROM AUTHOR]

Ectopic beats are known to be involved in the initiation of a variety of cardiac arrhythmias. Although their location may vary, ectopic excitations have been found to originate from infarct areas, regions of micro-fibrosis and other heterogeneous tissues. However, the underlying mechanisms that link ectopic foci to heterogeneous tissues have yet to be fully understood. In this work, we investigate the mechanism of micro-reentry that leads to the generation of ectopic beats near infarct areas using a patient-specific heart model. The patient-specific geometrical model of the heart, including scar and peri-infarct zones, is obtained through magnetic resonance imaging (MRI). The infarct region is composed of ischemic myocytes and non-conducting cells (fibrosis, for instance). Electrophysiology is captured using an established cardiac myocyte model of the human ventricle modified to describe ischemia. The simulation results clearly reveal that ectopic beats emerge from micro-reentries that are sustained by the heterogeneous structure of the infarct regions. Because microscopic information about the heterogeneous structure of the infarct regions is not available, Monte-Carlo simulations are used to identify the probabilities of an infarct region to behave as an ectopic focus for different levels of ischemia and different percentages of non-conducting cells. From the proposed model, it is observed that ectopic beats are generated when a percentage of non-conducting cells is near a topological metric known as the percolation threshold. Although the mechanism for micro-reentries was proposed half a century ago to be a source of ectopic beats or premature ventricular contractions during myocardial infarction, the present study is the first to reproduce this mechanism in-silico using patient-specific data. [ABSTRACT FROM AUTHOR]