The mechanisms of potassium loss in acute myocardial ischemia: New insights from computational simulations.

Acute myocardial ischemia induces hyperkalemia (accumulation of extracellular potassium), a major perpetrator of lethal reentrant ventricular arrhythmias. Despite considerable experimental efforts to explain this pathology in the last decades, the intimate mechanisms behind hyperkalemia remain partially unknown. In order to investigate these mechanisms, we developed a novel computational model of acute myocardial ischemia which couples a) an electrophysiologically detailed human cardiomyocyte model that incorporates modifications to account for ischemia-induced changes in transmembrane currents, with b) a model of cardiac tissue and extracellular transport.

Analysis of vulnerability to reentry in acute myocardial ischemia using a realistic human heart model.

Electrophysiological alterations of the myocardium caused by acute ischemia constitute a pro-arrhythmic substrate for the generation of potentially lethal arrhythmias. Experimental evidence has shown that the main components of acute ischemia that induce these electrophysiological alterations are hyperkalemia, hypoxia (or anoxia in complete artery occlusion), and acidosis. However, the influence of each ischemic component on the likelihood of reentry is not completely established.

Analysis of the response of human iPSC-derived cardiomyocyte tissue to I block. A combined in vitro and in silico approach.

The high incidence of cardiac arrythmias underlines the need for the assessment of pharmacological therapies. In this field of drug efficacy, as in the field of drug safety highlighted by the Comprehensive in Vitro Proarrhythmia Assay initiative, new pillars for research have become crucial: firstly, the integration of in-silico experiments, and secondly the evaluation of fully integrated biological systems, such as human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). In this study, we therefore aimed to combine in-vitro experiments and in-silico simulations to evaluate the antiarrhythmic effect of L-type calcium current (I) block in hiPSC-CMs.

Mechanistic investigation of Ca2+ alternans in human heart failure and its modulation by fibroblasts.

Background: Heart failure (HF) is characterized, among other factors, by a progressive loss of contractile function and by the formation of an arrhythmogenic substrate, both aspects partially related to intracellular Ca2+ cycling disorders. In failing hearts both electrophysiological and structural remodeling, including fibroblast proliferation, contribute to changes in Ca2+ handling which promote the appearance of Ca2+ alternans (Ca-alt). Ca-alt in turn give rise to repolarization alternans, which promote dispersion of repolarization and contribute to reentrant activity.

Personalized Cardiac Computational Models: From Clinical Data to Simulation of Infarct-Related Ventricular Tachycardia.

In the chronic stage of myocardial infarction, a significant number of patients develop life-threatening ventricular tachycardias (VT) due to the arrhythmogenic nature of the remodeled myocardium. Radiofrequency ablation (RFA) is a common procedure to isolate reentry pathways across the infarct scar that are responsible for VT. Unfortunately, this strategy show relatively low success rates; up to 50% of patients experience recurrent VT after the procedure.

Optimization of Lead Placement in the Right Ventricle During Cardiac Resynchronization Therapy. A Simulation Study.

Patients suffering from heart failure and left bundle branch block show electrical ventricular dyssynchrony causing an abnormal blood pumping. Cardiac resynchronization therapy (CRT) is recommended for these patients. Patients with positive therapy response normally present QRS shortening and an increased left ventricle (LV) ejection fraction.

Ca Cycling Impairment in Heart Failure Is Exacerbated by Fibrosis: Insights Gained From Mechanistic Simulations.

Heart failure (HF) is characterized by altered Ca cycling, resulting in cardiac contractile dysfunction. Failing myocytes undergo electrophysiological remodeling, which is known to be the main cause of abnormal Ca homeostasis. However, structural remodeling, specifically proliferating fibroblasts coupled to myocytes in the failing heart, could also contribute to Ca cycling impairment.

Sensitivity analysis revealing the effect of modulating ionic mechanisms on calcium dynamics in simulated human heart failure.

Abnormal intracellular Ca2+ handling is the major contributor to the depressed cardiac contractility observed in heart failure. The electrophysiological remodeling associated with this pathology alters both the action potential and the Ca2+ dynamics, leading to a defective excitation-contraction coupling that ends in mechanical dysfunction. The importance of maintaining a correct intracellular Ca2+ concentration requires a better understanding of its regulation by ionic mechanisms.

Comparison between Hodgkin-Huxley and Markov formulations of cardiac ion channels.

When simulating the macroscopic current flowing through cardiac ion channels, two mathematical formalisms can be adopted: the Hodgkin-Huxley model (HHM) formulation, which describes openings and closings of channel 'gates', or the Markov model (MM) formulation, based on channel 'state' transitions. The latter was first used in 1995 to simulate the effects of mutations in ionic currents and, since then, its use has been extended to wild-type channels also. While the MMs better describe the actual behavior of ion channels, they are mathematically more complex than HHMs in terms of parameter estimation and identifiability and are computationally much more demanding, which can dramatically increase computational time in large-scale (e.

Three-dimensional cardiac computational modelling: methods, features and applications.

The combination of computational models and biophysical simulations can help to interpret an array of experimental data and contribute to the understanding, diagnosis and treatment of complex diseases such as cardiac arrhythmias. For this reason, three-dimensional (3D) cardiac computational modelling is currently a rising field of research. The advance of medical imaging technology over the last decades has allowed the evolution from generic to patient-specific 3D cardiac models that faithfully represent the anatomy and different cardiac features of a given alive subject.

Electrophysiological and structural remodeling in heart failure modulate arrhythmogenesis. 1D simulation study.

Background: Heart failure is a final common pathway or descriptor for various cardiac pathologies. It is associated with sudden cardiac death, which is frequently caused by ventricular arrhythmias. Electrophysiological remodeling, intercellular uncoupling, fibrosis and autonomic imbalance have been identified as major arrhythmogenic factors in heart failure etiology and progression.

In silico ischaemia-induced reentry at the Purkinje-ventricle interface.

Aims: This computational modelling work illustrates the influence of hyperkalaemia and electrical uncoupling induced by defined ischaemia on action potential (AP) propagation and the incidence of reentry at the Purkinje-ventricle interface in mammalian hearts.

Methods And Results: Unidimensional and bidimensional models of the Purkinje-ventricle subsystem, including ischaemic conditions (defined as phase 1B) in the ventricle and an ischaemic border zone, were developed by altering several important electrophysiological parameters of the Luo-Rudy AP model of the ventricular myocyte. Purkinje electrical activity was modelled using the equations of DiFrancesco and Noble.

Multiscale computational analysis of the bioelectric consequences of myocardial ischaemia and infarction.

Ischaemic heart disease is considered as the single most frequent cause of death, provoking more than 7 000 000 deaths every year worldwide. A high percentage of patients experience sudden cardiac death, caused in most cases by tachyarrhythmic mechanisms associated to myocardial ischaemia and infarction. These diseases are difficult to study using solely experimental means due to their complex dynamics and unstable nature.

Non-uniform dispersion of the source-sink relationship alters wavefront curvature.

The distribution of cellular source-sink relationships plays an important role in cardiac propagation. It can lead to conduction slowing and block as well as wave fractionation. It is of great interest to unravel the mechanisms underlying evolution in wavefront geometry.

Modeling the different sections of the cardiac conduction system to obtain realistic electrocardiograms.

The cardiac conduction system is divided in different sections that play an important role in the cardiac depolarization sequence and define the morphology of the electrocardiogram. In this study we have built several configurations for each section based on anatomical descriptions. The effect of the morphology of the bundle branches, and the density of both Purkinje branches and Purkinje-myocardial junctions (PMJ) has been studied by comparing the pseudo-ECGs obtained with the standard precordial leads of the electrocardiogram.

Dominant frequency and organization index maps in a realistic three-dimensional computational model of atrial fibrillation.

Aims: To study, using simulation, the spectral characteristics of different patterns of atrial fibrillation (AF) at high spatial resolution. Dominant frequency (DF) and organization index (OI) maps have been used to approximate the location of the focal source of high frequency during AF events.

Methods And Results: A realistic three-dimensional model of the human atria that includes fibre orientation, electrophysiological heterogeneity, and anisotropy was implemented.

Simulation and mechanistic investigation of the arrhythmogenic role of the late sodium current in human heart failure.

Heart failure constitutes a major public health problem worldwide. The electrophysiological remodeling of failing hearts sets the stage for malignant arrhythmias, in which the role of the late Na(+) current (I(NaL)) is relevant and is currently under investigation. In this study we examined the role of I(NaL) in the electrophysiological phenotype of ventricular myocytes, and its proarrhythmic effects in the failing heart.

Interaction of specialized cardiac conduction system with antiarrhythmic drugs: a simulation study.

The use of antiarrhythmic drugs is common to treat heart rhythm disorders. Computational modeling and simulation are promising tools that could be used to investigate the effects of specific drugs on cardiac electrophysiology. In this paper, we study the multiscale effects of dofetilide, a drug that blocks IKr, from cellular to organ level paying special attention to its effect on heart structures, in particular the specialized cardiac conduction system (CCS).

Systematic characterization of the ionic basis of rabbit cellular electrophysiology using two ventricular models.

Several mathematical models of rabbit ventricular action potential (AP) have been proposed to investigate mechanisms of arrhythmias and excitation-contraction coupling. Our study aims at systematically characterizing how ionic current properties modulate the main cellular biomarkers of arrhythmic risk using two widely-used rabbit ventricular models, and comparing simulation results using the two models with experimental data available for rabbit. A sensitivity analysis of AP properties, Ca²⁺ and Na⁺ dynamics, and their rate dependence to variations (±15% and ±30%) in the main transmembrane current conductances and kinetics was performed using the Shannon et al.

Sex and age related differences in drug induced QT prolongation by dofetilide under reduced repolarization reserve in simulated ventricular cells.

The aim of this study was to investigate sex and age related differences in drug induced QT prolongation by dofetilide under reduced repolarization reserve in simulated ventricular cells. Left ventricular endocardial action potentials were simulated using a modified Luo Rudy model. Sex, age and regional differences in currents I(CaL), IK(r), IK(s), and I(to) were incorporated into the model by modifying the equations representing them.

Effects of the antiarrhythmic drug dofetilide on transmural dispersion of repolarization in ventriculum. A computer modeling study.

Dofetilide is a class-III drug that inhibits the rapid component of the delayed potassium current ( I(Kr)). Experimental studies have shown that the different layers of ventricular muscle present differences in action potential duration (APD) and different responses to class III agents. It has been suggested that it contributes to APD heterogeneity in the ventricles.

Exploring the role of pH in modulating the effects of lidocaine in virtual ischemic tissue.

Lidocaine is a class I antiarrhytmic drug that blocks Na(+) channels and exists in both neutral and charged forms at a physiological pH. In this work, a mathematical model of pH and the frequency-modulated effects of lidocaine has been developed and incorporated into the Luo-Rudy model of the ventricular action potential. We studied the effects of lidocaine on Na(+) current, maximum upstroke velocity, and conduction velocity and demonstrated that a decrease of these parameters was dependent on pH, frequency, and concentration.

Adaptive macro finite elements for the numerical solution of monodomain equations in cardiac electrophysiology.

Many problems in biology and engineering are governed by anisotropic reaction-diffusion equations with a very rapidly varying reaction term. This usually implies the use of very fine meshes and small time steps in order to accurately capture the propagating wave while avoiding the appearance of spurious oscillations in the wave front. This work develops a family of macro finite elements amenable for solving anisotropic reaction-diffusion equations with stiff reactive terms.

The relative role of refractoriness and source-sink relationship in reentry generation during simulated acute ischemia.

During acute myocardial ischemia, reentrant episodes may lead to ventricular fibrillation (VF), giving rise to potentially mortal arrhythmias. VF has been traditionally related to dispersion of refractoriness and more recently to the source-sink relationship. Our goal is to theoretically investigate the relative role of dispersion of refractoriness and source-sink mismatch in vulnerability to reentry in the specific situation of regional myocardial acute ischemia.

A grid computing-based approach for the acceleration of simulations in cardiology.

This paper combines high-performance computing and grid computing technologies to accelerate multiple executions of a biomedical application that simulates the action potential propagation on cardiac tissues. First, a parallelization strategy was employed to accelerate the execution of simulations on a cluster of personal computers (PCs). Then, grid computing was employed to concurrently perform the multiple simulations that compose the cardiac case studies on the resources of a grid deployment, by means of a service-oriented approach.