Effect of Global Cardiac Ischemia on Human Ventricular Fibrillation: Insights from a Multi-scale Mechanistic Model of the Human Heart

Acute regional ischemia in the heart can lead to cardiac arrhythmias such as ventricular fibrillation (VF), which in turn compromise cardiac output and result in secondary global cardiac ischemia. The secondary ischemia may influence the underlying arrhythmia mechanism. A recent clinical study documents the effect of global cardiac ischaemia on the mechanisms of VF. During 150 seconds of global ischemia the dominant frequency of activation decreased, while after reperfusion it increased rapidly. At the same time the complexity of epicardial excitation, measured as the number of epicardical phase singularity points, remained approximately constant during ischemia. Here we perform numerical studies based on these clinical data and propose explanations for the observed dynamics of the period and complexity of activation patterns. In particular, we study the effects on ischemia in pseudo-1D and 2D cardiac tissue models as well as in an anatomically accurate model of human heart ventricles. We demonstrate that the fall of dominant frequency in VF during secondary ischemia can be explained by an increase in extracellular potassium, while the increase during reperfusion is consistent with washout of potassium and continued activation of the ATP-dependent potassium channels. We also suggest that memory effects are responsible for the observed complexity dynamics. In addition, we present unpublished clinical results of individual patient recordings and propose a way of estimating extracellular potassium and activation of ATP-dependent potassium channels from these measurements

Incorporating Inductances in Tissue-Scale Models of Cardiac Electrophysiology

In standard models of cardiac electrophysiology, including the bidomain and monodomain models, local perturbations can propagate at infinite speed. We address this unrealistic property by developing a hyperbolic bidomain model that is based on a generalization of Ohm's law with a Cattaneo-type model for the fluxes. Further, we obtain a hyperbolic monodomain model in the case that the intracellular and extracellular conductivity tensors have the same anisotropy ratio. In one spatial dimension, the hyperbolic monodomain model is equivalent to a cable model that includes axial inductances, and the relaxation times of the Cattaneo fluxes are strictly related to these inductances. A purely linear analysis shows that the inductances are negligible, but models of cardiac electrophysiology are highly nonlinear, and linear predictions may not capture the fully nonlinear dynamics. In fact, contrary to the linear analysis, we show that for simple nonlinear ionic models, an increase in conduction velocity is obtained for small and moderate values of the relaxation time. A similar behavior is also demonstrated with biophysically detailed ionic models. Using the Fenton-Karma model along with a low-order finite element spatial discretization, we numerically analyze differences between the standard monodomain model and the hyperbolic monodomain model. In a simple benchmark test, we show that the propagation of the action potential is strongly influenced by the alignment of the fibers with respect to the mesh in both the parabolic and hyperbolic models when using relatively coarse spatial discretizations. Accurate predictions of the conduction velocity require computational mesh spacings on the order of a single cardiac cell. We also compare the two formulations in the case of spiral break up and atrial fibrillation in an anatomically detailed model of the left atrium, and [...].Comment: 20 pages, 12 figure

Electrophysiology Model for a Human Heart with Ischemic Scar and Realistic Purkinje Network

The role of Purkinje fibres in the onset of arrhythmias is controversial and computer simulations may shed light on possible arrhythmic mechanisms involving the Purkinje fibres. However, few computational modelling studies currently include a detailed Purkinje network as part of the model. We present a coupled Purkinje-myocardium electrophysiology model that includes an explicit model for the ischemic scar plus a detailed Purkinje network, and compare simulated activation times to those obtained by electro-anatomical mapping in vivo during sinus rhythm pacing. The results illustrate the importance of using sufficiently dense Purkinje networks in patient-specific studies to capture correctly the myocardial early activation that may be influenced by surviving Purkinje fibres in the infarct region

Fully-Coupled Electromechanical Simulations of the LV Dog Anatomy Using HPC: Model Testing and Verification

Verification of electro-mechanic models of the heart require a good amount of reliable, high resolution, thorough in-vivo measurements. The detail of the mathematical models used to create simulations of the heart beat vary greatly. Generally, the objective of the simulation determines the modeling approach. However, it is important to exactly quantify the amount of error between the various approaches that can be used to simulate a heart beat by comparing them to ground truth data. The more detailed the model is, the more computing power it requires, we therefore employ a high-performance computing solver throughout this study. We aim to compare models to data measured experimentally to identify the effect of using a mathematical model of fibre orientation versus the measured fibre orientations using DT-MRI. We also use simultaneous endocardial stimuli vs an instantaneous myocardial stimulation to trigger the mechanic contraction. Our results show that synchronisation of the electrical and mechanical events in the heart beat are necessary to create a physiological timing of hemodynamic events. Synchronous activation of all of the myocardium provides an unrealistic timing of hemodynamic events in the cardiac cycle. Results also show the need of establishing a protocol to quantify the zero-pressure configuration of the left ventricular geometry to initiate the simulation protocol; however, the predicted zero-pressure configuration of the same geometry was different, depending on the origin of the fibre field employed.This work has been done with the support of the grant SEV-2011-00067 of Severo Ochoa Program, awarded by the Spanish Government to the Barcelona Supercomputing Center. Part of the research leading to these results has received funding from the Seventh Framework Programme (FP7/2007-2013) under grant agreement n 611823. It has also been partially funded from the by the Spanish Ministry of Economy and Competitiveness (TIN2011-28067).Peer ReviewedPostprint (author's final draft

Spiral-Wave Turbulence and Its Control in the Presence of Inhomogeneities in Four Mathematical Models of Cardiac Tissue

Regular electrical activation waves in cardiac tissue lead to the rhythmic contraction and expansion of the heart that ensures blood supply to the whole body. Irregularities in the propagation of these activation waves can result in cardiac arrhythmias, like ventricular tachycardia (VT) and ventricular fibrillation (VF), which are major causes of death in the industrialised world. Indeed there is growing consensus that spiral or scroll waves of electrical activation in cardiac tissue are associated with VT, whereas, when these waves break to yield spiral- or scroll-wave turbulence, VT develops into life-threatening VF: in the absence of medical intervention, this makes the heart incapable of pumping blood and a patient dies in roughly two-and-a-half minutes after the initiation of VF. Thus studies of spiral- and scroll-wave dynamics in cardiac tissue pose important challenges for in vivo and in vitro experimental studies and for in silico numerical studies of mathematical models for cardiac tissue. A major goal here is to develop low-amplitude defibrillation schemes for the elimination of VT and VF, especially in the presence of inhomogeneities that occur commonly in cardiac tissue. We present a detailed and systematic study of spiral- and scroll-wave turbulence and spatiotemporal chaos in four mathematical models for cardiac tissue, namely, the Panfilov, Luo-Rudy phase 1 (LRI), reduced Priebe-Beuckelmann (RPB) models, and the model of ten Tusscher, Noble, Noble, and Panfilov (TNNP). In particular, we use extensive numerical simulations to elucidate the interaction of spiral and scroll waves in these models with conduction and ionic inhomogeneities; we also examine the suppression of spiral- and scroll-wave turbulence by low-amplitude control pulses. Our central qualitative result is that, in all these models, the dynamics of such spiral waves depends very sensitively on such inhomogeneities. We also study two types of control schemes that have been suggested for the control of spiral turbulence, via low amplitude current pulses, in such mathematical models for cardiac tissue; our investigations here are designed to examine the efficacy of such control schemes in the presence of inhomogeneities. We find that a local pulsing scheme does not suppress spiral turbulence in the presence of inhomogeneities; but a scheme that uses control pulses on a spatially extended mesh is more successful in the elimination of spiral turbulence. We discuss the theoretical and experimental implications of our study that have a direct bearing on defibrillation, the control of life-threatening cardiac arrhythmias such as ventricular fibrillation

Share and enjoy: anatomical models database-generating and sharing cardiovascular model data using web services

Sharing data between scientists and with clinicians in cardiac research has been facilitated significantly by the use of web technologies. The potential of this technology has meant that information sharing has been routinely promoted through databases that have encouraged stakeholder participation in communities around these services. In this paper we discuss the Anatomical Model Database (AMDB) (Gianni et al. Functional imaging and modeling of the heart. Springer, Heidelberg, 2009; Gianni et al. Phil Trans Ser A Math Phys Eng Sci 368:3039–3056, 2010) which both facilitate a database-centric approach to collaboration, and also extends this framework with new capabilities for creating new mesh data. AMDB currently stores cardiac geometric models described in Gianni et al. (Functional imaging and modelling of the heart. Springer, Heidelberg, 2009), a number of additional cardiac models describing geometry and functional properties, and most recently models generated using a web service. The functional models represent data from simulations in geometric form, such as electrophysiology or mechanics, many of which are present in AMDB as part of a benchmark study. Finally, the heartgen service has been added for producing left or bi-ventricle models derived from binary image data using the methods described in Lamata et al. (Med Image Anal 15:801–813, 2011). The results can optionally be hosted on AMDB alongside other community-provided anatomical models. AMDB is, therefore, a unique database storing geometric data (rather than abstract models or image data) combined with a powerful web service for generating new geometric models

A Novel and Lethal De Novo LQT-3 Mutation in a Newborn with Distinct Molecular Pharmacology and Therapeutic Response

SCN5A encodes the alpha-subunit (Na(v)1.5) of the principle Na(+) channel in the human heart. Genetic lesions in SCN5A can cause congenital long QT syndrome (LQTS) variant 3 (LQT-3) in adults by disrupting inactivation of the Na(v)1.5 channel. Pharmacological targeting of mutation-altered Na(+) channels has proven promising in developing a gene-specific therapeutic strategy to manage specifically this LQTS variant. SCN5A mutations that cause similar channel dysfunction may also contribute to sudden infant death syndrome (SIDS) and other arrhythmias in newborns, but the prevalence, impact, and therapeutic management of SCN5A mutations may be distinct in infants compared with adults.Here, in a multidisciplinary approach, we report a de novo SCN5A mutation (F1473C) discovered in a newborn presenting with extreme QT prolongation and differential responses to the Na(+) channel blockers flecainide and mexiletine. Our goal was to determine the Na(+) channel phenotype caused by this severe mutation and to determine whether distinct effects of different Na(+) channel blockers on mutant channel activity provide a mechanistic understanding of the distinct therapeutic responsiveness of the mutation carrier. Sequence analysis of the proband revealed the novel missense SCN5A mutation (F1473C) and a common variant in KCNH2 (K897T). Patch clamp analysis of HEK 293 cells transiently transfected with wild-type or mutant Na(+) channels revealed significant changes in channel biophysics, all contributing to the proband's phenotype as predicted by in silico modeling. Furthermore, subtle differences in drug action were detected in correcting mutant channel activity that, together with both the known genetic background and age of the patient, contribute to the distinct therapeutic responses observed clinically.The results of our study provide further evidence of the grave vulnerability of newborns to Na(+) channel defects and suggest that both genetic background and age are particularly important in developing a mutation-specific therapeutic personalized approach to manage disorders in the young

Mathematical Modeling and Simulation of Ventricular Activation Sequences: Implications for Cardiac Resynchronization Therapy

Next to clinical and experimental research, mathematical modeling plays a crucial role in medicine. Biomedical research takes place on many different levels, from molecules to the whole organism. Due to the complexity of biological systems, the interactions between components are often difficult or impossible to understand without the help of mathematical models. Mathematical models of cardiac electrophysiology have made a tremendous progress since the first numerical ECG simulations in the 1960s. This paper briefly reviews the development of this field and discusses some example cases where models have helped us forward, emphasizing applications that are relevant for the study of heart failure and cardiac resynchronization therapy

Identification of IKr Kinetics and Drug Binding in Native Myocytes

Determining the effect of a compound on IKr is a standard screen for drug safety. Often the effect is described using a single IC50 value, which is unable to capture complex effects of a drug. Using verapamil as an example, we present a method for using recordings from native myocytes at several drug doses along with qualitative features of IKr from published studies of HERG current to estimate parameters in a mathematical model of the drug effect on IKr. IKr was recorded from canine left ventricular myocytes using ruptured patch techniques. A voltage command protocol was used to record tail currents at voltages from −70 to −20 mV, following activating pulses over a wide range of voltages and pulse durations. Model equations were taken from a published IKr Markov model and the drug was modeled as binding to the open state. Parameters were estimated using a combined global and local optimization algorithm based on collected data with two additional constraints on IKrI–V relation and IKr inactivation. The method produced models that quantitatively reproduce both the control IKr kinetics and dose dependent changes in the current. In addition, the model exhibited use and rate dependence. The results suggest that: (1) the technique proposed here has the practical potential to develop data-driven models that quantitatively reproduce channel behavior in native myocytes; (2) the method can capture important drug effects that cannot be reproduced by the IC50 method. Although the method was developed for IKr, the same strategy can be applied to other ion channels, once appropriate channel-specific voltage protocols and qualitative features are identified