Arrhythmogenesis in the heart: Multiscale modeling of the effects of defibrillation shocks and the role of electrophysiological heterogeneity

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Effect of Heart Structure on Ventricular Fibrillation in the Rabbit: A Simulation Study

Ventricular fibrillation (VF) is a lethal condition that affects millions worldwide. The mechanism underlying VF is unstable reentrant electrical waves rotating around lines called filaments. These complex spatio-temporal patterns can be studied using both experimental and numerical methods. Computer simulations provide unique insights including high resolution dynamics throughout the heart and systematic control of quantities such as fiber orientation and cellular kinetics that are not feasible experimentally. Here we study filament dynamics using two bi-ventricular 3-D high-resolution rabbit heart geometries, one with detailed fine structure and another without fine structure. We studied filament dynamics using anisotropic and isotropic conductivities, and with four cellular action potential models with different recovery kinetics. Spiral wave dynamics observed in isotropic two-dimensional sheets were not predictive of the behavior in the whole heart. In 2-D the four cell models exhibited stable reentry, meandering spiral waves, and spiral-wave breakup. In the whole heart with fine structure, all simulation results exhibited complex dynamics reminiscent of fibrillation observed experimentally. In the whole heart without fine structure, anisotropy acted to destabilize filament dynamics although the number of filaments was reduced compared to the heart with structure. In addition, in isotropic hearts without structure the two cell models that exhibited meandering spiral waves in 2-D, stabilized into figure-of-eight surface patterns. We also studied the sensitivity of filament dynamics to computer system configuration and initial conditions. After large simulation times, different macroscopic results sometimes occurred across different system configurations, likely due to a lack of bitwise reproducibility. The study conclusions were insensitive to initial condition perturbations, however, the exact number of filaments over time and their trends were altered by these changes. In summary, we present the following new results. First, we provide a new cell model that resembles the surface patterns of VF in the rabbit heart both qualitatively and quantitatively. Second, filament dynamics in the whole heart cannot be predicted from spiral wave dynamics in 2-D and we identified anisotropy as one destabilizing factor. Third, the exact dynamics of filaments are sensitive to a variety of factors, so we suggest caution in their interpretation and their quantitative analyses.

A bilobal model of Ca2+-dependent inactivation to probe the physiology of L-type Ca2+ channels

L-type calcium channels undergo Ca²⁺-dependent inactivation (CDI) in order to precisely control the entry of Ca²⁺ into cells such as cardiomyocytes. Limpitikul et al. develop a bilobal model of CDI and use it to understand the pathogenesis of arrhythmias associated with mutations in CaM.

L-type calcium channels (LTCCs) are critical elements of normal cardiac function, playing a major role in orchestrating cardiac electrical activity and initiating downstream signaling processes. LTCCs thus use feedback mechanisms to precisely control calcium (Ca²⁺) entry into cells. Of these, Ca²⁺-dependent inactivation (CDI) is significant because it shapes cardiac action potential duration and is essential for normal cardiac rhythm. This important form of regulation is mediated by a resident Ca²⁺ sensor, calmodulin (CaM), which is comprised of two lobes that are each capable of responding to spatially distinct Ca²⁺ sources. Disruption of CaM-mediated CDI leads to severe forms of long-QT syndrome (LQTS) and life-threatening arrhythmias. Thus, a model capable of capturing the nuances of CaM-mediated CDI would facilitate increased understanding of cardiac (patho)physiology. However, one critical barrier to achieving a detailed kinetic model of CDI has been the lack of quantitative data characterizing CDI as a function of Ca²⁺. This data deficit stems from the experimental challenge of uncoupling the effect of channel gating on Ca²⁺ entry. To overcome this obstacle, we use photo-uncaging of Ca²⁺ to deliver a measurable Ca²⁺ input to CaM/LTCCs, while simultaneously recording CDI. Moreover, we use engineered CaMs with Ca²⁺ binding restricted to a single lobe, to isolate the kinetic response of each lobe. These high-resolution measurements enable us to build mathematical models for each lobe of CaM, which we use as building blocks for a full-scale bilobal model of CDI. Finally, we use this model to probe the pathogenesis of LQTS associated with mutations in CaM (calmodulinopathies). Each of these models accurately recapitulates the kinetics and steady-state properties of CDI in both physiological and pathological states, thus offering powerful new insights into the mechanistic alterations underlying cardiac arrhythmias.

Towards a Unified Theory of Calmodulin Regulation (Calmodulation) of Voltage-Gated Calcium and Sodium Channels

Voltage-gated Na and Ca(2+) channels represent two major ion channel families that enable myriad biological functions including the generation of action potentials and the coupling of electrical and chemical signaling in cells. Calmodulin regulation (calmodulation) of these ion channels comprises a vital feedback mechanism with distinct physiological implications. Though long-sought, a shared understanding of the channel families remained elusive for two decades as the functional manifestations and the structural underpinnings of this modulation often appeared to diverge. Here, we review recent advancements in the understanding of calmodulation of Ca(2+) and Na channels that suggest a remarkable similarity in their regulatory scheme. This interrelation between the two channel families now paves the way towards a unified mechanistic framework to understand vital calmodulin-dependent feedback and offers shared principles to approach related channelopathic diseases. An exciting era of synergistic study now looms.

Computational rabbit models to investigate the initiation, perpetuation, and termination of ventricular arrhythmia

Current understanding of cardiac electrophysiology has been greatly aided by computational work performed using rabbit ventricular models. This article reviews the contributions of multiscale models of rabbit ventricles in understanding cardiac arrhythmia mechanisms. This review will provide an overview of multiscale modeling of the rabbit ventricles. It will then highlight works that provide insights into the role of the conduction system, complex geometric structures, and heterogeneous cellular electrophysiology in diseased and healthy rabbit hearts to the initiation and maintenance of ventricular arrhythmia. Finally, it will provide an overview on the contributions of rabbit ventricular modeling on understanding the mechanisms underlying shock-induced defibrillation.

How computer simulations of the human heart can improve anti-arrhythmia therapy

Over the last decade, the state-of-the-art in cardiac computational modelling has progressed rapidly. The electrophysiological function of the heart can now be simulated with a high degree of detail and accuracy, opening the doors for simulation-guided approaches to anti-arrhythmic drug development and patient-specific therapeutic interventions. In this review, we outline the basic methodology for cardiac modelling, which has been developed and validated over decades of research. In addition, we present several recent examples of how computational models of the human heart have been used to address current clinical problems in cardiac electrophysiology. We will explore the use of simulations to improve anti-arrhythmic pacing and defibrillation interventions; to predict optimal sites for clinical ablation procedures; and to aid in the understanding and selection of arrhythmia risk markers. Together, these studies illustrate how the tremendous advances in cardiac modelling are poised to revolutionize medical treatment and prevention of arrhythmia.

Filament Dynamics during Simulated Ventricular Fibrillation in a High-Resolution Rabbit Heart

The mechanisms underlying ventricular fibrillation (VF) are not well understood. The electrical activity on the heart surface during VF has been recorded extensively in the experimental setting and in some cases clinically; however, corresponding transmural activation patterns are prohibitively difficult to measure. In this paper, we use a high-resolution biventricular heart model to study three-dimensional electrical activity during fibrillation, focusing on the driving sources of VF: “filaments,” the organising centres of unstable reentrant scroll waves. We show, for the first time, specific 3D filament dynamics during simulated VF in a whole heart geometry that includes fine-scale anatomical structures. Our results suggest that transmural activity is much more complex than what would be expected from surface observations alone. We present examples of complex intramural activity, including filament breakup and reattachment, anchoring to the thin right ventricular apex; rapid transitions among various filament shapes; and filament lengths much greater than wall thickness. We also present evidence for anatomy playing a major role in VF development and coronary vessels and trabeculae influencing filament dynamics. Overall, our results indicate that intramural activity during simulated VF is extraordinarily complex and suggest that further investigation of 3D filaments is necessary to fully comprehend recorded surface patterns.

Computational modeling of cardiac optogenetics: Methodology overview & review of findings from simulations

Cardiac optogenetics is emerging as an exciting new potential avenue to enable spatiotemporally precise control of excitable cells and tissue in the heart with low-energy optical stimuli. This approach involves the expression of exogenous light-sensitive proteins (opsins) in target heart tissue via viral gene or cell delivery. Preliminary experiments in optogenetically-modified cells, tissue, and organisms have made great strides towards demonstrating the feasibility of basic applications, including the use of light stimuli to pace or disrupt reentrant activity. However, it remains unknown whether techniques based on this intriguing technology could be scaled up and used in humans for novel clinical applications, such as pain-free optical defibrillation or dynamic modulation of action potential shape. A key step towards answering such questions is to explore potential optogenetics-based therapies using sophisticated computer simulation tools capable of realistically representing opsin delivery and light stimulation in biophysically detailed, patient-specific models of the human heart. This review provides (1) a detailed overview of the methodological developments necessary to represent optogenetics-based solutions in existing virtual heart platforms and (2) a survey of findings that have been derived from such simulations and a critical assessment of their significance with respect to the progress of the field.

Calmodulin mutations associated with long QT syndrome prevent inactivation of cardiac L-type Ca(2+) currents and promote proarrhythmic behavior in ventricular myocytes

Recent work has identified missense mutations in calmodulin (CaM) that are associated with severe early-onset long-QT syndrome (LQTS), leading to the proposition that altered CaM function may contribute to the molecular etiology of this subset of LQTS. To date, however, no experimental evidence has established these mutations as directly causative of LQTS substrates, nor have the molecular targets of CaM mutants been identified. Here, therefore, we test whether expression of CaM mutants in adult guinea-pig ventricular myocytes (aGPVM) induces action-potential prolongation, and whether affiliated alterations in the Ca(2+) regulation of L-type Ca(2+) channels (LTCC) might contribute to such prolongation. In particular, we first overexpressed CaM mutants in aGPVMs, and observed both increased action potential duration (APD) and heightened Ca(2+) transients. Next, we demonstrated that all LQTS CaM mutants have the potential to strongly suppress Ca(2+)/CaM-dependent inactivation (CDI) of LTCCs, whether channels were heterologously expressed in HEK293 cells, or present in native form within myocytes. This attenuation of CDI is predicted to promote action-potential prolongation and boost Ca(2+) influx. Finally, we demonstrated how a small fraction of LQTS CaM mutants (as in heterozygous patients) would nonetheless suffice to substantially diminish CDI, and derange electrical and Ca(2+) profiles. In all, these results highlight LTCCs as a molecular locus for understanding and treating CaM-related LQTS in this group of patients.

New insights into defibrillation of the heart from realistic simulation studies

Cardiac defibrillation, as accomplished nowadays by automatic, implantable devices, constitutes the most important means of combating sudden cardiac death. Advancing our understanding towards a full appreciation of the mechanisms by which a shock interacts with the heart, particularly under diseased conditions, is a promising approach to achieve an optimal therapy. The aim of this article is to assess the current state-of-the-art in whole-heart defibrillation modelling, focusing on major insights that have been obtained using defibrillation models, primarily those of realistic heart geometry and disease remodelling. The article showcases the contributions that modelling and simulation have made to our understanding of the defibrillation process. The review thus provides an example of biophysically based computational modelling of the heart (i.e. cardiac defibrillation) that has advanced the understanding of cardiac electrophysiological interaction at the organ level, and has the potential to contribute to the betterment of the clinical practice of defibrillation.

Sodium current reduction unmasks a structure-dependent substrate for arrhythmogenesis in the normal ventricles

Background

Organ-scale arrhythmogenic consequences of source-sink mismatch caused by impaired excitability remain unknown, hindering the understanding of pathophysiology in disease states like Brugada syndrome and ischemia.

Objective

We sought to determine whether sodium current (I_Na) reduction in the structurally normal heart unmasks a regionally heterogeneous substrate for the induction of sustained arrhythmia by premature ventricular contractions (PVCs).

Methods

We conducted simulations in rabbit ventricular computer models with 930 unique combinations of PVC location (10 sites) and coupling interval (250–400 ms), I_Na reduction (30 or 40% of normal levels), and post-PVC sinus rhythm (arrested or persistent). Geometric characteristics and source-sink mismatch were quantitatively analyzed by calculating ventricular wall thickness and a newly formulated 3D safety factor (SF), respectively.

Results

Reducing I_Na to 30% of its normal level created a substrate for sustained arrhythmia induction by establishing large regions of critical source-sink mismatch (SF<1) for ectopic wavefronts propagating from thin to thick tissue. In the same simulations but with 40% of normal I_Na, PVCs did not induce reentry because the volume of tissue with SF<1 was >95% smaller. Likewise, when post-PVC sinus activations were persistent instead of arrested, no ectopic excitations initiated sustained reentry because sinus activation breakthroughs engulfed the excitable gap.

Conclusion

Our new SF formulation can quantify ectopic wavefront propagation robustness in geometrically complex 3D tissue with impaired excitability. This novel methodology was applied to show that I_Na reduction precipitates source-sink mismatch, creating a potent substrate for sustained arrhythmia induction by PVCs originating near regions of ventricular wall expansion, such as the RV outflow tract.

Advances in modeling ventricular arrhythmias: from mechanisms to the clinic

Modern cardiovascular research has increasingly recognized that heart models and simulation can help interpret an array of experimental data and dissect important mechanisms and interrelationships, with developments rooted in the iterative interaction between modeling and experimentation. This article reviews the progress made in simulating cardiac electrical behavior at the level of the organ and, specifically, in the development of models of ventricular arrhythmias and fibrillation, as well as their termination (defibrillation). The ability to construct multiscale models of ventricular arrhythmias, representing integrative behavior from the molecule to the entire organ, has enabled mechanistic inquiry into the dynamics of ventricular arrhythmias in the diseased myocardium, in understanding drug-induced proarrhythmia, and in the development of new modalities for defibrillation, to name a few. In this article, we also review the initial use of ventricular models of arrhythmia in personalized diagnosis, treatment planning, and prevention of sudden cardiac death. Implementing individualized cardiac simulations at the patient bedside is poised to become one of the most thrilling examples of computational science and engineering approaches in translational medicine.

Tachycardia in post-infarction hearts: insights from 3D image-based ventricular models

Ventricular tachycardia, a life-threatening regular and repetitive fast heart rhythm, frequently occurs in the setting of myocardial infarction. Recently, the peri-infarct zones surrounding the necrotic scar (termed gray zones) have been shown to correlate with ventricular tachycardia inducibility. However, it remains unknown how the latter is determined by gray zone distribution and size. The goal of this study is to examine how tachycardia circuits are maintained in the infarcted heart and to explore the relationship between the tachycardia organizing centers and the infarct gray zone size and degree of heterogeneity. To achieve the goals of the study, we employ a sophisticated high-resolution electrophysiological model of the infarcted canine ventricles reconstructed from imaging data, representing both scar and gray zone. The baseline canine ventricular model was also used to generate additional ventricular models with different gray zone sizes, as well as models in which the gray zone was represented as different heterogeneous combinations of viable tissue and necrotic scar. The results of the tachycardia induction simulations with a number of high-resolution canine ventricular models (22 altogether) demonstrated that the gray zone was the critical factor resulting in arrhythmia induction and maintenance. In all models with inducible arrhythmia, the scroll-wave filaments were contained entirely within the gray zone, regardless of its size or the level of heterogeneity of its composition. The gray zone was thus found to be the arrhythmogenic substrate that promoted wavebreak and reentry formation. We found that the scroll-wave filament locations were insensitive to the structural composition of the gray zone and were determined predominantly by the gray zone morphology and size. The findings of this study have important implications for the advancement of improved criteria for stratifying arrhythmia risk in post-infarction patients and for the development of new approaches for determining the ablation targets of infarct-related tachycardia.

Effects of pacing magnitudes and forms on bistability width in a modeled ventricular tissue

Bistability in periodically paced cardiac tissue is relevant to cardiac arrhythmias and its control. In the present paper, one-dimensional tissue of the phase I Luo-Rudy model is numerically investigated. The effects of various parameters of pacing signals on bistability width are studied. The following conclusions are obtained: (i) Pacing can be classified into two types: pulsatile and sinusoidal types. Pulsatile pacing reduces bistability width as its magnitude is increased. Sinusoidal pacing increases the width as its amplitude is increased. (ii) In a pacing period the hyperpolarizing part plays a more important role than the depolarizing part. Variations of the hyperpolarizing ratio in a period evidently change the width of bistability and its variation tendency. (iii) A dynamical mechanism is proposed to qualitatively explain the phenomena, which reveals the reason for the different effects of pulsatile and sinusoidal pacing on bistability. The methods for changing bistability width by external pacing may help control arrhythmias in cardiology.

Placement of implantable cardioverter-defibrillators in paediatric and congenital heart defect patients: a pipeline for model generation and simulation prediction of optimal configurations

There is currently no reliable way of predicting the optimal implantable cardioverter-defibrillator (ICD) placement in paediatric and congenital heart defect (CHD) patients. This study aimed to: (1) develop a new image processing pipeline for constructing patient-specific heart-torso models from clinical magnetic resonance images (MRIs); (2) use the pipeline to determine the optimal ICD configuration in a paediatric tricuspid valve atresia patient; (3) establish whether the widely used criterion of shock-induced extracellular potential (Φe) gradients ≥5 V cm(-1) in ≥95% of ventricular volume predicts defibrillation success. A biophysically detailed heart-torso model was generated from patient MRIs. Because transvenous access was impossible, three subcutaneous and three epicardial lead placement sites were identified along with five ICD scan locations. Ventricular fibrillation was induced, and defibrillation shocks were applied from 11 ICD configurations to determine defibrillation thresholds (DFTs). Two configurations with epicardial leads resulted in the lowest DFTs overall and were thus considered optimal. Three configurations shared the lowest DFT among subcutaneous lead ICDs. The Φe gradient criterion was an inadequate predictor of defibrillation success, as defibrillation failed in numerous instances even when 100% of the myocardium experienced such gradients. In conclusion, we have developed a new image processing pipeline and applied it to a CHD patient to construct the first active heart-torso model from clinical MRIs.

Computational cardiology: the heart of the matter

This paper reviews the newest developments in computational cardiology. It focuses on the contribution of cardiac modeling to the development of new therapies as well as the advancement of existing ones for cardiac arrhythmias and pump dysfunction. Reviewed are cardiac modeling efforts aimed at advancing and optimizing existent therapies for cardiac disease (defibrillation, ablation of ventricular tachycardia, and cardiac resynchronization therapy) and at suggesting novel treatments, including novel molecular targets, as well as efforts to use cardiac models in stratification of patients likely to benefit from a given therapy, and the use of models in diagnostic procedures.

Modeling atrial arrhythmias: impact on clinical diagnosis and therapies

Atrial arrhythmias are the most frequent sustained rhythm disorders in humans and often lead to severe complications such as heart failure and stroke. Despite the important insights provided by animal models into the mechanisms of atrial arrhythmias, direct translation of experimental findings to new therapies in patients has not been straightforward. With the advances in computer technology, large-scale electroanatomical computer models of the atria that integrate information from the molecular to organ scale have reached a level of sophistication that they can be used to interpret the outcome of experimental and clinical studies and aid in the rational design of therapies. This paper reviews the state-of-the-art of computer models of the electrical dynamics of the atria and discusses the evolving role of simulation in assisting the clinical diagnosis and treatment of atrial arrhythmias.

Coupled personalization of cardiac electrophysiology models for prediction of ischaemic ventricular tachycardia

In order to translate the important progress in cardiac electrophysiology modelling of the last decades into clinical applications, there is a requirement to make macroscopic models that can be used for the planning and performance of the clinical procedures. This requires model personalization, i.e. estimation of patient-specific model parameters and computations compatible with clinical constraints. Simplified macroscopic models can allow a rapid estimation of the tissue conductivity, but are often unreliable to predict arrhythmias. Conversely, complex biophysical models are more complete and have mechanisms of arrhythmogenesis and arrhythmia sustainibility, but are computationally expensive and their predictions at the organ scale still have to be validated. We present a coupled personalization framework that combines the power of the two kinds of models while keeping the computational complexity tractable. A simple eikonal model is used to estimate the conductivity parameters, which are then used to set the parameters of a biophysical model, the Mitchell-Schaeffer (MS) model. Additional parameters related to action potential duration restitution curves for the tissue are further estimated for the MS model. This framework is applied to a clinical dataset derived from a hybrid X-ray/magnetic resonance imaging and non-contact mapping procedure on a patient with heart failure. This personalized MS model is then used to perform an in silico simulation of a ventricular tachycardia (VT) stimulation protocol to predict the induction of VT. This proof of concept opens up possibilities of using VT induction modelling in order to both assess the risk of VT for a given patient and also to plan a potential subsequent radio-frequency ablation strategy to treat VT.

Integrative computational models of cardiac arrhythmias -- simulating the structurally realistic heart

Simulation of cardiac electrical function, and specifically, simulation aimed at understanding the mechanisms of cardiac rhythm disorders, represents an example of a successful integrative multiscale modeling approach, uncovering emergent behavior at the successive scales in the hierarchy of structural complexity. The goal of this article is to present a review of the integrative multiscale models of realistic ventricular structure used in the quest to understand and treat ventricular arrhythmias. It concludes with the new advances in image-based modeling of the heart and the promise it holds for the development of individualized models of ventricular function in health and disease.

Numerical solution of the bidomain equations

Knowledge of cardiac electrophysiology is efficiently formulated in terms of mathematical models. However, most of these models are very complex and thus defeat direct mathematical reasoning founded on classical and analytical considerations. This is particularly so for the celebrated bidomain model that was developed almost 40 years ago for the concurrent analysis of extra- and intracellular electrical activity. Numerical simulations based on this model represent an indispensable tool for studying electrophysiology. However, complex mathematical models, steep gradients in the solutions and complicated geometries lead to extremely challenging computational problems. The greatest achievement in scientific computing over the past 50 years has been to enable the solving of linear systems of algebraic equations that arise from discretizations of partial differential equations in an optimal manner, i.e. such that the central processing unit (CPU) effort increases linearly with the number of computational nodes. Over the past decade, such optimal methods have been introduced in the simulation of electrophysiology. This development, together with the development of affordable parallel computers, has enabled the solution of the bidomain model combined with accurate cellular models, on geometries resembling a human heart. However, in spite of recent progress, the full potential of modern computational methods has yet to be exploited for the solution of the bidomain model. This paper reviews the development of numerical methods for solving the bidomain model. However, the field is huge and we thus restrict our focus to developments that have been made since the year 2000.

Towards predictive modelling of the electrophysiology of the heart

The simulation of cardiac electrical function is an example of a successful integrative multiscale modelling approach that is directly relevant to human disease. Today we stand at the threshold of a new era, in which anatomically detailed, tomographically reconstructed models are being developed that integrate from the ion channel to the electromechanical interactions in the intact heart. Such models hold high promise for interpretation of clinical and physiological measurements, for improving the basic understanding of the mechanisms of dysfunction in disease, such as arrhythmias, myocardial ischaemia and heart failure, and for the development and performance optimization of medical devices. The goal of this article is to present an overview of current state-of-art advances towards predictive computational modelling of the heart as developed recently by the authors of this article. We first outline the methodology for constructing electrophysiological models of the heart. We then provide three examples that demonstrate the use of these models, focusing specifically on the mechanisms for arrhythmogenesis and defibrillation in the heart. These include: (1) uncovering the role of ventricular structure in defibrillation; (2) examining the contribution of Purkinje fibres to the failure of the shock; and (3) using magnetic resonance imaging reconstructed heart models to investigate the re-entrant circuits formed in the presence of an infarct scar.

Organization of ventricular fibrillation in the human heart: experiments and models

Sudden cardiac death is a major health problem in the industrialized world. The lethal event is typically ventricular fibrillation (VF), during which the co-ordinated regular contraction of the heart is overthrown by a state of mechanical and electrical anarchy. Understanding the excitation patterns that sustain VF is important in order to identify potential therapeutic targets. In this paper, we studied the organization of human VF by combining clinical recordings of electrical excitation patterns on the epicardial surface during in vivo human VF with simulations of VF in an anatomically and electrophysiologically detailed computational model of the human ventricles. We find both in the computational studies and in the clinical recordings that epicardial surface excitation patterns during VF contain around six rotors. Based on results from the simulated three-dimensional excitation patterns during VF, which show that the total number of electrical sources is 1.4 +/- 0.12 times greater than the number of epicardial rotors, we estimate that the total number of sources present during clinically recorded VF is 9.0 +/- 2.6. This number is approximately fivefold fewer compared with that observed during VF in dog and pig hearts, which are of comparable size to human hearts. We explain this difference by considering differences in action potential duration dynamics across these species. The simpler spatial organization of human VF has important implications for treatment and prevention of this dangerous arrhythmia. Moreover, our findings underline the need for integrated research, in which human-based clinical and computational studies complement animal research.

Intrinsic inhomogeneities and the coexistence of spirals with different periods of rotation

We propose a mechanism by which wave fronts emanating from a spiral may break far from the spiral core due to intrinsic spatial inhomogeneities. A series of computer simulations are presented to demonstrate how coupling domains, which on their own would not cause breakup, may cause a single spiral to break into many spirals.

Systems approach to understanding electromechanical activity in the human heart: a national heart, lung, and blood institute workshop summary

The National Heart, Lung, and Blood Institute (NHLBI) convened a workshop of cardiologists, cardiac electrophysiologists, cell biophysicists, and computational modelers on August 20 and 21, 2007, in Washington, DC, to advise the NHLBI on new research directions needed to develop integrative approaches to elucidate human cardiac function. The workshop strove to identify limitations in the use of data from nonhuman animal species for elucidation of human electromechanical function/activity and to identify what specific information on ion channel kinetics, calcium handling, and dynamic changes in the intracellular/extracellular milieu is needed from human cardiac tissues to develop more robust computational models of human cardiac electromechanical activity. This article summarizes the workshop discussions and recommendations on the following topics: (1) limitations of animal models and differences from human electrophysiology, (2) modeling ion channel structure/function in the context of whole-cell electrophysiology, (3) excitation-contraction coupling and regulatory pathways, (4) whole-heart simulations of human electromechanical activity, and (5) what human data are currently needed and how to obtain them. The recommendations can be found on the NHLBI Web site at http://www.nhlbi.nih.gov/meetings/workshops/electro.htm.

Effect of bundle branch block on cardiac output: a whole heart simulation study

The heart is an electrically controlled fluid pump which operates by mechanical contraction. Whole heart modelling is a computationally daunting task which must incorporate several subsystems: mechanical, electrical, and fluidic. Numerous feedback mechanisms on many levels, and operating at different scales, exist to finely control behaviour. Understanding these interactions is necessary to understand heart operation, as well as pathologies and therapies. A review of the components in such a model is given. The authors then present a framework for their electro-mechano-fluidic whole heart model based on cable methods. The model incorporates atria and ventricles, and has functioning valves with papillary muscles. The effect of altered propagation due to left and right bundle branch block on cardiac output is examined using the cable-based model. Results are compared to clinically observed phenomena. Good agreement was obtained, but tighter coupling of mechanical and electrical events is needed to fully account for behaviour. Cable-based models offer an alternative to continuum models.

Introduction: Cardiovascular physics

The number of patients suffering from cardiovascular diseases increases unproportionally high with the increase of the human population and aging, leading to very high expenses in the public health system. Therefore, the challenge of cardiovascular physics is to develop high-sophisticated methods which are able to, on the one hand, supplement and replace expensive medical devices and, on the other hand, improve the medical diagnostics with decreasing the patient's risk. Cardiovascular physics-which interconnects medicine, physics, biology, engineering, and mathematics-is based on interdisciplinary collaboration of specialists from the above scientific fields and attempts to gain deeper insights into pathophysiology and treatment options. This paper summarizes advances in cardiovascular physics with emphasis on a workshop held in Bad Honnef, Germany, in May 2005. The meeting attracted an interdisciplinary audience and led to a number of papers covering the main research fields of cardiovascular physics, including data analysis, modeling, and medical application. The variety of problems addressed by this issue underlines the complexity of the cardiovascular system. It could be demonstrated in this Focus Issue, that data analyses and modeling methods from cardiovascular physics have the ability to lead to significant improvements in different medical fields. Consequently, this Focus Issue of Chaos is a status report that may invite all interested readers to join the community and find competent discussion and cooperation partners.