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 K + transport. The resulting model is able to reproduce and explain the triphasic time course of extracellular K + concentration within the ischemic zone, with values of [ K + ] o close to 14 mmol/L in the central ischemic zone after 30 min. In addition, the formation of a [ K + ] o border zone of approximately 1.2 cm 15 min after the onset of ischemia is predicted by the model. Our results indicate that the primary rising phase of [ K + ] o is mainly due to the imbalance between K + efflux, that increases slightly, and K + influx, that follows a reduction of the NaK pump activity by more than 50%. The onset of the plateau phase is caused by the appearance of electrical alternans (a novel mechanism identified by the model), which cause an abrupt reduction in the K + efflux. After the plateau, the secondary rising phase of [ K + ] o is caused by a subsequent imbalance between the K + influx, which continues to decrease slowly, and the K + efflux, which remains almost constant. Further, the study shows that the modulation of these mechanisms by the electrotonic coupling is the main responsible for the formation of the ischemic border zone in tissue, with K + transport playing only a minor role. Finally, the results of the model indicate that the injury current established between the healthy and the altered tissue is not sufficient to depolarize non-ischemic cells within the healthy tissue.

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. Moreover, the role of the His-Purkinje system (HPS) in the initiation and maintenance of arrhythmias is not completely understood. In the present work, we investigate how the three components of ischemia affect the vulnerable window (VW) for reentry using computational simulations. In addition, we analyze the role of the HPS on arrhythmogenesis. A 3D biventricular/torso human model that includes a realistic geometry of the central and border ischemic zones with one of the most electrophysiologically detailed model of ischemia to date, as well as a realistic cardiac conduction system, were used to assess the VW for reentry. Four scenarios of ischemic severity corresponding to different minutes after coronary artery occlusion were simulated. Our results suggest that ischemic severity plays an important role in the generation of reentries. Indeed, this is the first 3D simulation study to show that ventricular arrhythmias could be generated under moderate ischemic conditions, but not in mild and severe ischemia. Moreover, our results show that anoxia is the ischemic component with the most significant effect on the width of the VW. Thus, a change in the level of anoxia from moderate to severe leads to a greater increment in the VW (40 ms), in comparison with the increment of 20 ms and 35 ms produced by the individual change in the level of hyperkalemia and acidosis, respectively. Finally, the HPS was a necessary element for the generation of approximately 17% of reentries obtained. The retrograde conduction from the myocardium to HPS in the ischemic region, conduction blocks in discrete sections of the HPS, and the degree of ischemia affecting Purkinje cells, are suggested as mechanisms that favor the generation of ventricular arrhythmias.

Analysis of the response of human iPSC-derived cardiomyocyte tissue to ICaL 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 (ICaL) block in hiPSC-CMs. For this, hiPSC-CM preparations were cultured and an equivalent virtual tissue was modeled. Re-entry patterns of electrical activation were induced and several biomarkers were obtained before and after ICaL block. The virtual hiPSC-CM simulations were also reproduced using a tissue composed of adult ventricular cardiomyocytes (hAdultV-CMs). The analysis of phases, currents and safety factor for propagation showed an increased size of the re-entry core when ICaL was blocked as a result of depressed cellular excitability. The bigger wavefront curvature yielded reductions of 12.2%, 6.9%, and 4.2% in the frequency of the re-entry for hiPSC-CM cultures, virtual hiPSC-CM, and hAdultV-CM tissues, respectively. Furthermore, ICaL block led to a 47.8% shortening of the vulnerable window for re-entry in the virtual hiPSC-CM tissue and to re-entry vanishment in hAdultV-CM tissue. The consistent behavior between in-vitro and in-silico hiPSC-CMs and between in-silico hiPSC-CMs and hAdultV-CMs evidences that virtual hiPSC-CM tissues are suitable for assessing cardiac efficacy, as done in the present study through the analysis of ICaL block.

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. In the last decade, intensive research in the field of computational cardiac electrophysiology (EP) has demonstrated the ability of three-dimensional (3D) cardiac computational models to perform in-silico EP studies. However, the personalization and modeling of certain key components remain challenging, particularly in the case of the infarct border zone (BZ). In this study, we used a clinical dataset from a patient with a history of infarct-related VT to build an image-based 3D ventricular model aimed at computational simulation of cardiac EP, including detailed patient-specific cardiac anatomy and infarct scar geometry. We modeled the BZ in eight different ways by combining the presence or absence of electrical remodeling with four different levels of image-based patchy fibrosis (0, 10, 20, and 30%). A 3D torso model was also constructed to compute the ECG. Patient-specific sinus activation patterns were simulated and validated against the patient's ECG. Subsequently, the pacing protocol used to induce reentrant VTs in the EP laboratory was reproduced in-silico. The clinical VT was induced with different versions of the model and from different pacing points, thus identifying the slow conducting channel responsible for such VT. Finally, the real patient's ECG recorded during VT episodes was used to validate our simulation results and to assess different strategies to model the BZ. Our study showed that reduced conduction velocities and heterogeneity in action potential duration in the BZ are the main factors in promoting reentrant activity. Either electrical remodeling or fibrosis in a degree of at least 30% in the BZ were required to initiate VT. Moreover, this proof-of-concept study confirms the feasibility of developing 3D computational models for cardiac EP able to reproduce cardiac activation in sinus rhythm and during VT, using exclusively non-invasive clinical data.

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. The computational analysis of the incidence of Ca2+ and/or repolarization alternans under HF conditions in the presence of fibroblasts could provide a better understanding of the mechanisms leading to HF arrhythmias and contractile function disorders. METHODS AND FINDINGS: The goal of the present study was to investigate in silico the mechanisms leading to the formation of Ca-alt in failing human ventricular myocytes and tissues with disperse fibroblast distributions. The contribution of ionic currents variability to alternans formation at the cellular level was analyzed and the results show that in normal ventricular tissue, altered Ca2+ dynamics lead to Ca-alt, which precede APD alternans and can be aggravated by the presence of fibroblasts. Electrophysiological remodeling of failing tissue alone is sufficient to develop alternans. The incidence of alternans is reduced when fibroblasts are present in failing tissue due to significantly depressed Ca2+ transients. The analysis of the underlying ionic mechanisms suggests that Ca-alt are driven by Ca2+-handling protein and Ca2+ cycling dysfunctions in the junctional sarcoplasmic reticulum and that their contribution to alternans occurrence depends on the cardiac remodeling conditions and on myocyte-fibroblast interactions. CONCLUSION: It can thus be concluded that fibroblasts modulate the formation of Ca-alt in human ventricular tissue subjected to heart failure-related electrophysiological remodeling. Pharmacological therapies should thus consider the extent of both the electrophysiological and structural remodeling present in the failing heart.

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. However, around one third do not respond favorably. Therefore, optimal location of pacing leads, timing delays between leads and/or choosing related biomarkers is crucial to achieve the best possible degree of ventricular synchrony during CRT application. In this study, computational modeling is used to predict the optimal location and delay of pacing leads to improve CRT response. We use a 3D electrophysiological computational model of the heart and torso to get insight into the changes in the activation patterns obtained when the heart is paced from different regions and for different atrioventricular and interventricular delays. The model represents a heart with left bundle branch block and heart failure, and allows a detailed and accurate analysis of the electrical changes observed simultaneously in the myocardium and in the QRS complex computed in the precordial leads. Computational simulations were performed using a modified version of the O'Hara et al. action potential model, the most recent mathematical model developed for human ventricular electrophysiology. The optimal location for the pacing leads was determined by QRS maximal reduction. Additionally, the influence of Purkinje system on CRT response was assessed and correlation analysis between several parameters of the QRS was made. Simulation results showed that the right ventricle (RV) upper septum near the outflow tract is an alternative location to the RV apical lead. Furthermore, LV endocardial pacing provided better results as compared to epicardial stimulation. Finally, the time to reach the 90% of the QRS area was a good predictor of the instant at which 90% of the ventricular tissue was activated. Thus, the time to reach the 90% of the QRS area is suggested as an additional index to assess CRT effectiveness to improve biventricular synchrony.

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

Heart failure (HF) is characterized by altered Ca2+ cycling, resulting in cardiac contractile dysfunction. Failing myocytes undergo electrophysiological remodeling, which is known to be the main cause of abnormal Ca2+ homeostasis. However, structural remodeling, specifically proliferating fibroblasts coupled to myocytes in the failing heart, could also contribute to Ca2+ cycling impairment. The goal of the present study was to systematically analyze the mechanisms by which myocyte-fibroblast coupling could affect Ca2+ dynamics in normal conditions and in HF. Simulations of healthy and failing human myocytes were performed using established mathematical models, and cells were either isolated or coupled to fibroblasts. Univariate and multivariate sensitivity analyses were performed to quantify effects of ion transport pathways on biomarkers computed from intracellular [Ca2+] waveforms. Variability in ion channels and pumps was imposed and populations of models were analyzed to determine effects on Ca2+ dynamics. Our results suggest that both univariate and multivariate sensitivity analyses are valuable methodologies to shed light into the ionic mechanisms underlying Ca2+ impairment in HF, although differences between the two methodologies are observed at high parameter variability. These can result from either the fact that multivariate analyses take into account ion channels or non-linear effects of ion transport pathways on Ca2+ dynamics. Coupling either healthy or failing myocytes to fibroblasts decreased Ca2+ transients due to an indirect sink effect on action potential (AP) and thus on Ca2+ related currents. Simulations that investigated restoration of normal physiology in failing myocytes showed that Ca2+ cycling can be normalized by increasing SERCA and L-type Ca2+ current activity while decreasing Na+-Ca2+ exchange and SR Ca2+ leak. Changes required to normalize APs in failing myocytes depended on whether myocytes were coupled to fibroblasts. In conclusion, univariate and multivariate sensitivity analyses are helpful tools to understand how Ca2+ cycling is impaired in HF and how this can be exacerbated by coupling of myocytes to fibroblasts. The design of pharmacological actions to restore normal activity should take into account the degree of fibrosis in the failing heart.

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. To study the electrical activity and ionic homeostasis of failing myocytes, a modified version of the O'Hara et al. human action potential model was used, including electrophysiological remodeling. The impact of the main ionic transport mechanisms was analyzed using single-parameter sensitivity analyses, the first of which explored the modulation of electrophysiological characteristics related to Ca2+ exerted by the remodeled parameters. The second sensitivity analysis compared the potential consequences of modulating individual channel conductivities, as one of the main effects of potential drugs, on Ca2+ dynamic properties under both normal conditions and in heart failure. The first analysis revealed the important contribution of the sarcoplasmic reticulum Ca2+-ATPase (SERCA) dysfunction to the altered Ca2+ homeostasis, with the Na+/Ca2+ exchanger (NCX) and other Ca2+ cycling proteins also playing a significant role. Our results highlight the importance of improving the SR uptake function to increase Ca2+ content and restore Ca2+ homeostasis and contractility. The second sensitivity analysis highlights the different response of the failing myocyte versus the healthy myocyte to potential pharmacological actions on single channels. The result of modifying the conductances of the remodeled proteins such as SERCA and NCX in heart failure has less impact on Ca2+ modulation. These differences should be taken into account when designing drug therapies.

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.g. whole heart) simulations. We hypothesize that a HHM formulation obtained from classical patch-clamp protocols in wild-type and mutant ion channels can be used to correctly simulate cardiac action potentials and their static and dynamic properties. To validate our hypothesis, we selected two pivotal cardiac ionic currents (the rapid delayed rectifier K(+) current, IKr, and the inward Na(+) current, INa) and formulated HHMs for both wild-type and mutant channels (LQT2-linked T474I mutation for IKr and LQT3-linked ΔKPQ mutation for INa). Action potentials were then simulated using the MM and HHM versions of the currents, and the action potential waveforms, biomarkers and action potential duration rate dependence properties were compared in control conditions and in the presence of physiological variability. While small differences between ionic currents were found between the two models (correlation coefficient ρ>0.92), the simulations yielded almost identical action potentials (ρ>0.99), suggesting that HHMs may also be valid to simulate the effects of mutations affecting IKr and INa on the action potential.

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. Here we analyse sixty representative 3D cardiac computational models developed and published during the last fifty years, describing their information sources, features, development methods and online availability. This paper also reviews the necessary components to build a 3D computational model of the heart aimed at biophysical simulation, paying especial attention to cardiac electrophysiology (EP), and the existing approaches to incorporate those components. We assess the challenges associated to the different steps of the building process, from the processing of raw clinical or biological data to the final application, including image segmentation, inclusion of substructures and meshing among others. We briefly outline the personalisation approaches that are currently available in 3D cardiac computational modelling. Finally, we present examples of several specific applications, mainly related to cardiac EP simulation and model-based image analysis, showing the potential usefulness of 3D cardiac computational modelling into clinical environments as a tool to aid in the prevention, diagnosis and treatment of cardiac diseases.

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. OBJECTIVE: In this study we investigate in silico the role of electrophysiological and structural heart failure remodeling on the modulation of key elements of the arrhythmogenic substrate, i.e., electrophysiological gradients and abnormal impulse propagation. METHODS: Two different mathematical models of the human ventricular action potential were used to formulate models of the failing ventricular myocyte. This provided the basis for simulations of the electrical activity within a transmural ventricular strand. Our main goal was to elucidate the roles of electrophysiological and structural remodeling in setting the stage for malignant life-threatening arrhythmias. RESULTS: Simulation results illustrate how the presence of M cells and heterogeneous electrophysiological remodeling in the human failing ventricle modulate the dispersion of action potential duration and repolarization time. Specifically, selective heterogeneous remodeling of expression levels for the Na+/Ca2+ exchanger and SERCA pump decrease these heterogeneities. In contrast, fibroblast proliferation and cellular uncoupling both strongly increase repolarization heterogeneities. Conduction velocity and the safety factor for conduction are also reduced by the progressive structural remodeling during heart failure. CONCLUSION: An extensive literature now establishes that in human ventricle, as heart failure progresses, gradients for repolarization are changed significantly by protein specific electrophysiological remodeling (either homogeneous or heterogeneous). Our simulations illustrate and provide new insights into this. Furthermore, enhanced fibrosis in failing hearts, as well as reduced intercellular coupling, combine to increase electrophysiological gradients and reduce electrical propagation. In combination these changes set the stage for arrhythmias.

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. In the past decades, integrative computational simulation techniques have become a powerful tool to complement experimental and clinical research when trying to elucidate the intimate mechanisms of ischaemic electrophysiological processes and to aid the clinician in the improvement and optimization of therapeutic procedures. The purpose of this paper is to briefly review some of the multiscale computational models of myocardial ischaemia and infarction developed in the past 20 years, ranging from the cellular level to whole-heart simulations.

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. Our study suggests that an extracellular potassium concentration [K(+)]o >14 mM and a slight decrease in intercellular coupling induced by ischaemia in ventricle can cause conduction block from Purkinje to ventricle. Under these conditions, propagation from ventricle to Purkinje is possible. Thus, unidirectional block (UDB) and reentry can result. When conditions of UDB are met, retrograde propagation with a long delay (320 ms) may re-excite Purkinje cells, and give rise to a reentrant pathway. This induced reentry may be the origin of arrhythmias observed in phase 1B ischaemia. CONCLUSION: In a defined setting of ischaemia (phase 1B), a small amount of uncoupling between ventricular cells, as well as between Purkinje and ventricular tissue, may induce UDBs and reentry. Hyperkalaemia is also confirmed to be an important factor in the genesis of reentrant rhythms, since it regulates the range of coupling in which UDBs may be induced.

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. Our goal is to investigate the role of the source-sink relationship on wavefront geometry using computer simulations. We analyzed the role of variability in the microscopic source-sink relationship in driving changes in wavefront geometry. The electrophysiological activity of a homogeneous isotropic tissue was simulated using the ten Tusscher and Panfilov 2006 action potential model and the source-sink relationship was characterized using an improved version of the Romero et al. safety factor formulation (SFm2). Our simulations reveal that non-uniform dispersion of the cellular source-sink relationship (dispersion along the wavefront) leads to alterations in curvature. To better understand the role of the source-sink relationship in the process of wave formation, the electrophysiological activity at the initiation of excitation waves in a 1D strand was examined and the source-sink relationship was characterized using the two recently updated safety factor formulations: the SFm2 and the Boyle-Vigmond (SFVB) definitions. The electrophysiological activity at the initiation of excitation waves was intimately related to the SFm2 profiles, while the SFVB led to several counterintuitive observations. Importantly, with the SFm2 characterization, a critical source-sink relationship for initiation of excitation waves was identified, which was independent of the size of the electrode of excitation, membrane excitability, or tissue conductivity. In conclusion, our work suggests that non-uniform dispersion of the source-sink relationship alters wavefront curvature and a critical source-sink relationship profile separates wave expansion from collapse. Our study reinforces the idea that the safety factor represents a powerful tool to study the mechanisms of cardiac propagation in silico, providing a better understanding of cardiac arrhythmias and their therapy.

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. A functional model for the PMJs based on the existence of a conduction adaptation layer is also presented. Simulation results showed a large influence of the His bundle and bundle branches in the pseudo-ECG and helped to elucidate the most appropriate morphology. The functional PMJ model allowed bidirectional communication between the conduction system and the myocardium with realistic transmission delays between both mediums. These results can help to improve current conduction system models and improve depolarization sequences of activation in the ventricles.

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. A model for cellular heart failure was proposed using a modified version of Grandi et al. model for human ventricular action potential that incorporates the formulation of I(NaL). A sensitivity analysis of the model was performed and simulations of the pathological electrical activity of the cell were conducted. The proposed model for the human I(NaL) and the electrophysiological remodeling of myocytes from failing hearts accurately reproduce experimental observations. The sensitivity analysis of the modulation of electrophysiological parameters of myocytes from failing hearts due to ion channels remodeling, revealed a role for I(NaL) in the prolongation of action potential duration (APD), triangulation of the shape of the AP, and changes in Ca(2+) transient. A mechanistic investigation of intracellular Na(+) accumulation and APD shortening with increasing frequency of stimulation of failing myocytes revealed a role for the Na(+)/K(+) pump, the Na(+)/Ca(2+) exchanger and I(NaL). The results of the simulations also showed that in failing myocytes, the enhancement of I(NaL) increased the reverse rate-dependent APD prolongation and the probability of initiating early afterdepolarizations. The electrophysiological remodeling of failing hearts and especially the enhancement of the I(NaL) prolong APD and alter Ca(2+) transient facilitating the development of early afterdepolarizations. An enhanced I(NaL) appears to be an important contributor to the electrophysiological phenotype and to the dysregulation of [Ca(2+)](i) homeostasis of failing myocytes.

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. (2004) and the Mahajan et al. (2008a,b) AP rabbit models. The effects of severe transmembrane current blocks (up to 100%) on steady-state AP and calcium transients, and AP duration (APD) restitution curves were also simulated using both models. Our simulations show that, in both virtual rabbit cardiomyocytes, APD is significantly modified by most repolarization currents, AP triangulation is regulated mostly by the inward rectifier K⁺ current (I(K1)) whereas APD rate adaptation as well as [Na⁺](i) rate dependence is influenced by the Na⁺/K⁺ pump current (I(NaK)). In addition, steady-state [Ca²âº](i) levels, APD restitution properties and [Ca²âº](i) rate dependence are strongly dependent on I(NaK), the L-Type Ca²âº current (I(CaL)) and the Na⁺/Ca²âº exchanger current (I(NaCa)), although the relative role of these currents is markedly model dependent. Furthermore, our results show that simulations using both models agree with many experimentally-reported electrophysiological characteristics. However, our study shows that the Shannon et al. model mimics rabbit electrophysiology more accurately at normal pacing rates, whereas Mahajan et al. model behaves more appropriately at faster rates. Our results reinforce the usefulness of sensitivity analysis for further understanding of cellular electrophysiology and validation of cardiac AP models.

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. However, in vivo studies suggest that the strong cellular coupling reduces APD dispersion in intact heart. The aim of this paper is to study the effect of dofetilide on the action potentials (APs) in isolated ventricular cells and on APD dispersion in a strand of ventricular tissue. A mathematical model of dofetilide effects on I(Kr) has been developed and incorporated into the Luo--Rudy dynamic model of ventricular AP. Our results show that dofetilide induces in midmyocardium cells a faster time-course inhibition of I(Kr) than in endocardial or epicardial cells, and periods of instability with beat-to-beat APs variability. This behavior could favor temporal dispersion of repolarization between the different cells. The results also indicate that although dofetilide increases, the transmural gradient of APD in the ventricular wall, early afterdepolarizations (EADs) did not appear even under strong uncoupling conditions. However, reduced repolarization reserve favors the induction of EADs, even under normal coupling conditions.

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 developed elements are incorporated on a semi-implicit algorithm based on operator splitting that includes adaptive time stepping for handling the stiff reactive term. A linear system is solved on each time step to update the transmembrane potential, whereas the remaining ordinary differential equations are solved uncoupled. The method allows solving the linear system on a coarser mesh thanks to the static condensation of the internal degrees of freedom (DOF) of the macroelements while maintaining the accuracy of the finer mesh. The method and algorithm have been implemented in parallel. The accuracy of the method has been tested on two- and three-dimensional examples demonstrating excellent behavior when compared to standard linear elements. The better performance and scalability of different macro finite elements against standard finite elements have been demonstrated in the simulation of a human heart and a heterogeneous two-dimensional problem with reentrant activity. Results have shown a reduction of up to four times in computational cost for the macro finite elements with respect to equivalent (same number of DOF) standard linear finite elements as well as good scalability properties.

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. The electrical activity of a regionally ischemic tissue was simulated using a modified version of the Luo-Rudy dynamic model. Ischemic conditions were varied to simulate the time-course of acute ischemia. Our results showed that dispersion of refractoriness increased with the severity of ischemia. However, no correlation between dispersion of refractoriness and the width of the vulnerable window was found. Additionally, in approximately 50% of the reentries, unidirectional block (UDB) took place in cells completely recovered from refractoriness. We examined patterns of activation after premature stimulation and they were intimately related to the source-sink relationship, quantified by the safety factor (SF). Moreover, the isoline where the SF dropped below unity matched the area where propagation failed. It was concluded that the mismatch of the source-sink relationship, rather than solely refractoriness, was the ultimate cause of the UDB leading to reentry. The SF represents a very powerful tool to study the mechanisms responsible for reentry.