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

Evaluating pro-arrhythmogenic effects of the T634S-hERG mutation: insights from a simulation study

A mutation to serine of a conserved threonine (T634S) in the hERG K+ channel S6 pore region has been identified as a variant of uncertain significance, showing a loss-of-function effect. However, its potential consequences for ventricular excitation and arrhythmogenesis have not been reported. This study evaluated possible functional effects of the T634S-hERG mutation on ventricular excitation and arrhythmogenesis by using multi-scale computer models of the human ventricle. A Markov chain model of the rapid delayed rectifier potassium current (IKr) was reconstructed for wild-type and T634S-hERG mutant conditions and incorporated into the ten Tusscher et al. models of human ventricles at cell and tissue (1D, 2D and 3D) levels. Possible functional impacts of the T634S-hERG mutation were evaluated by its effects on action potential durations (APDs) and their rate-dependence (APDr) at the cell level; and on the QT interval of pseudo-ECGs, tissue vulnerability to unidirectional conduction block (VW), spiral wave dynamics and repolarization dispersion at the tissue level. It was found that the T634S-hERG mutation prolonged cellular APDs, steepened APDr, prolonged the QT interval, increased VW, destablized re-entry and augmented repolarization dispersion across the ventricle. Collectively, these results imply potential pro-arrhythmic effects of the T634S-hERG mutation, consistent with LQT2.

Different Ventricular Fibrillation Types in Low-Dimensional Latent Spaces

The causes of ventricular fibrillation (VF) are not yet elucidated, and it has been proposed that different mechanisms might exist. Moreover, conventional analysis methods do not seem to provide time or frequency domain features that allow for recognition of different VF patterns in electrode-recorded biopotentials. The present work aims to determine whether low-dimensional latent spaces could exhibit discriminative features for different mechanisms or conditions during VF episodes. For this purpose, manifold learning using autoencoder neural networks was analyzed based on surface ECG recordings. The recordings covered the onset of the VF episode as well as the next 6 min, and comprised an experimental database based on an animal model with five situations, including control, drug intervention (amiodarone, diltiazem, and flecainide), and autonomic nervous system blockade. The results show that latent spaces from unsupervised and supervised learning schemes yielded moderate though quite noticeable separability among the different types of VF according to their type or intervention. In particular, unsupervised schemes reached a multi-class classification accuracy of 66%, while supervised schemes improved the separability of the generated latent spaces, providing a classification accuracy of up to 74%. Thus, we conclude that manifold learning schemes can provide a valuable tool for studying different types of VF while working in low-dimensional latent spaces, as the machine-learning generated features exhibit separability among different VF types. This study confirms that latent variables are better VF descriptors than conventional time or domain features, making this technique useful in current VF research on elucidation of the underlying VF mechanisms.

On the Role of Ionic Modeling on the Signature of Cardiac Arrhythmias for Healthy and Diseased Hearts

Computational cardiology is rapidly becoming the gold standard for innovative medical treatments and device development. Despite a worldwide effort in mathematical and computational modeling research, the complexity and intrinsic multiscale nature of the heart still limit our predictability power raising the question of the optimal modeling choice for large-scale whole-heart numerical investigations. We propose an extended numerical analysis among two different electrophysiological modeling approaches: a simplified phenomenological one and a detailed biophysical one. To achieve this, we considered three-dimensional healthy and infarcted swine heart geometries. Heterogeneous electrophysiological properties, fine-tuned DT-MRI -based anisotropy features, and non-conductive ischemic regions were included in a custom-built finite element code. We provide a quantitative comparison of the electrical behaviors during steady pacing and sustained ventricular fibrillation for healthy and diseased cases analyzing cardiac arrhythmias dynamics. Action potential duration (APD) restitution distributions, vortex filament counting, and pseudo-electrocardiography (ECG) signals were numerically quantified, introducing a novel statistical description of restitution patterns and ventricular fibrillation sustainability. Computational cost and scalability associated with the two modeling choices suggests that ventricular fibrillation signatures are mainly controlled by anatomy and structural parameters, rather than by regional restitution properties. Finally, we discuss limitations and translational perspectives of the different modeling approaches in view of large-scale whole-heart in silico studies.

Prediction of Cardiac Mechanical Performance From Electrical Features During Ventricular Tachyarrhythmia Simulation Using Machine Learning Algorithms

In ventricular tachyarrhythmia, electrical instability features including action potential duration, dominant frequency, phase singularity, and filaments are associated with mechanical contractility. However, there are insufficient studies on estimated mechanical contractility based on electrical features during ventricular tachyarrhythmia using a stochastic model. In this study, we predicted cardiac mechanical performance from features of electrical instability during ventricular tachyarrhythmia simulation using machine learning algorithms, including support vector regression (SVR) and artificial neural network (ANN) models. We performed an electromechanical tachyarrhythmia simulation and extracted 12 electrical instability features and two mechanical properties, including stroke volume and the amplitude of myocardial tension (ampTens). We compared predictive performance according to kernel types of the SVR model and the number of hidden layers of the ANN model. In the SVR model, the prediction accuracies of stroke volume and ampTens were the highest when using the polynomial kernel and linear kernel, respectively. The predictive performance of the ANN model was better than that of the SVR model. The prediction accuracies were the highest when the ANN model consisted of three hidden layers. Accordingly, we propose the ANN model with three hidden layers as an optimal model for predicting cardiac mechanical contractility in ventricular tachyarrhythmia. The results of this study are expected to be used to indirectly estimate the hemodynamic response from the electrical cardiac map measured by the optical mapping system during cardiac surgery, as well as cardiac contractility under normal sinus rhythm conditions.

Integration of activation maps of epicardial veins in computational cardiac electrophysiology

In this work we address the issue of validating the monodomain equation used in combination with the Bueno-Orovio ionic model for the prediction of the activation times in cardiac electro-physiology of the left ventricle. To this aim, we consider four patients who suffered from Left Bundle Branch Block (LBBB). We use activation maps performed at the septum as input data for the model and maps at the epicardial veins for the validation. In particular, a first set (half) of the latter are used to estimate the conductivities of the patient and a second set (the remaining half) to compute the errors of the numerical simulations. We find an excellent agreement between measures and numerical results. Our validated computational tool could be used to accurately predict activation times at the epicardial veins with a short mapping, i.e. by using only a part (the most proximal) of the standard acquisition points, thus reducing the invasive procedure and exposure to radiation.

Sensitivity of a data-assimilation system for reconstructing three-dimensional cardiac electrical dynamics

Modelling of cardiac electrical behaviour has led to important mechanistic insights, but important challenges, including uncertainty in model formulations and parameter values, make it difficult to obtain quantitatively accurate results. An alternative approach is combining models with observations from experiments to produce a data-informed reconstruction of system states over time. Here, we extend our earlier data-assimilation studies using an ensemble Kalman filter to reconstruct a three-dimensional time series of states with complex spatio-temporal dynamics using only surface observations of voltage. We consider the effects of several algorithmic and model parameters on the accuracy of reconstructions of known scroll-wave truth states using synthetic observations. In particular, we study the algorithm's sensitivity to parameters governing different parts of the process and its robustness to several model-error conditions. We find that the algorithm can achieve an acceptable level of error in many cases, with the weakest performance occurring for model-error cases and more extreme parameter regimes with more complex dynamics. Analysis of the poorest-performing cases indicates an initial decrease in error followed by an increase when the ensemble spread is reduced. Our results suggest avenues for further improvement through increasing ensemble spread by incorporating additive inflation or using a parameter or multi-model ensemble. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Self-restoration of cardiac excitation rhythm by anti-arrhythmic ion channel gating

Homeostatic regulation protects organisms against hazardous physiological changes. However, such regulation is limited in certain organs and associated biological processes. For example, the heart fails to self-restore its normal electrical activity once disturbed, as with sustained arrhythmias. Here we present proof-of-concept of a biological self-restoring system that allows automatic detection and correction of such abnormal excitation rhythms. For the heart, its realization involves the integration of ion channels with newly designed gating properties into cardiomyocytes. This allows cardiac tissue to i) discriminate between normal rhythm and arrhythmia based on frequency-dependent gating and ii) generate an ionic current for termination of the detected arrhythmia. We show in silico, that for both human atrial and ventricular arrhythmias, activation of these channels leads to rapid and repeated restoration of normal excitation rhythm. Experimental validation is provided by injecting the designed channel current for arrhythmia termination in human atrial myocytes using dynamic clamp.

Relationship Between Electrical Instability and Pumping Performance During Ventricular Tachyarrhythmia: Computational Study

There are representative electrical parameters for understanding the mechanism of reentrant waves in studies on tachyarrhythmia, namely the action potential duration (APD), dominant frequency, phase singularity, and filament. However, there are no studies that have directly identified the correlation between these electrophysiological parameters and cardiac contractility. Therefore, we have identified individual and integrative correlations between these electrical phenomena and contractility during tachyarrhythmia by deriving regression equations and also investigated the electrophysiological parameters affecting cardiac contractility during tachyarrhythmia. We simulated ventricular tachyarrhythmia with 48 types of electrical patterns by applying four reentry generation methods and changing the electrical conductivity of the potassium channel, which has the greatest effect on ventricular tissue. The mechanical responses reflecting electrical complexity were obtained through deterministic simulations of excitation-contraction coupling. We used the stroke volume and amplitude of myocardial tension (ampTens) as the variables representing contractility. We derived stochastic models through single- and multivariable regression analyses to identify the electrical parameters affecting contractility during tachyarrhythmia. In single-variable regression analysis, the APD, dominant frequency, and filament, excluding phase singularity, have statistically significant correlations with the stroke volume and ampTens. Among them, the APD has the maximum influence on these two mechanical parameters (standard beta coefficient: 0.859 for stroke volume, 0.930 for ampTens). The stochastic model using all four electrical parameters fails to accurately predict contractility owing to the multicollinearity between the APD and dominant frequency. We have rederived the multi-variable stochastic model using three electrical parameters without the APD. The filament has the greatest effect on the stroke volume stochastically (standard beta coefficient: 0.853 and 0.752). The dominant frequency has the greatest effect on ampTens statistically (standard beta coefficient: -0.813). We conclude that among the electrical parameters, the APD has the highest individual influence on mechanical contraction, and the filament has the highest integrative influence in both statistical terms.

Phase Entrainment of Induced Ventricular Fibrillation: A Human Feasibility and Proof of Concept Study

Cardioversion and defibrillation by a single high energy shock applied by myocardial or body surface electrodes is painful, causes long term tissue damage, and is associated with worsening long term outcomes, but is almost always required for treatment of ventricular fibrillation . As a initial step towards developing methods that can terminate ventricular arrhythmias painlessly, we aim to determine if pacing stimuli at a rate of 5/s applied via an implantable cardiac defibrillator (ICD) can modify human ventricular fibrillation. In 8 patients undergoing defibrillation testing of a new/exchanged intracardiac defibrillator, five seconds of pacing at five stimuli per second was applied during the 10-20 seconds of induced ventricular fibrillation before the defibrillation shock was automatically applied, and the cardiac electrograms recorded and analyzed. The high frequency pacing did not entrain the ventricular fibrillation, but altered the dominant frequency in all 8 patients, and modulated the phase computed via the Hilbert Transform, in four of the patients. In this pilot study we demonstrate that high frequency pacing applied via ICD electrodes during VF can alter the dominant frequency and modulate the probability density of the phase of the electrogram of the ventricular fibrillation.

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.

Validation and Trustworthiness of Multiscale Models of Cardiac Electrophysiology

Computational models of cardiac electrophysiology have a long history in basic science applications and device design and evaluation, but have significant potential for clinical applications in all areas of cardiovascular medicine, including functional imaging and mapping, drug safety evaluation, disease diagnosis, patient selection, and therapy optimisation or personalisation. For all stakeholders to be confident in model-based clinical decisions, cardiac electrophysiological (CEP) models must be demonstrated to be trustworthy and reliable. Credibility, that is, the belief in the predictive capability, of a computational model is primarily established by performing validation, in which model predictions are compared to experimental or clinical data. However, there are numerous challenges to performing validation for highly complex multi-scale physiological models such as CEP models. As a result, credibility of CEP model predictions is usually founded upon a wide range of distinct factors, including various types of validation results, underlying theory, evidence supporting model assumptions, evidence from model calibration, all at a variety of scales from ion channel to cell to organ. Consequently, it is often unclear, or a matter for debate, the extent to which a CEP model can be trusted for a given application. The aim of this article is to clarify potential rationale for the trustworthiness of CEP models by reviewing evidence that has been (or could be) presented to support their credibility. We specifically address the complexity and multi-scale nature of CEP models which makes traditional model evaluation difficult. In addition, we make explicit some of the credibility justification that we believe is implicitly embedded in the CEP modeling literature. Overall, we provide a fresh perspective to CEP model credibility, and build a depiction and categorisation of the wide-ranging body of credibility evidence for CEP models. This paper also represents a step toward the extension of model evaluation methodologies that are currently being developed by the medical device community, to physiological models.

Computational Cardiac Modeling Reveals Mechanisms of Ventricular Arrhythmogenesis in Long QT Syndrome Type 8: CACNA1C R858H Mutation Linked to Ventricular Fibrillation

Functional analysis of the L-type calcium channel has shown that the CACNA1C R858H mutation associated with severe QT interval prolongation may lead to ventricular fibrillation (VF). This study investigated multiple potential mechanisms by which the CACNA1C R858H mutation facilitates and perpetuates VF. The Ten Tusscher-Panfilov (TP06) human ventricular cell models incorporating the experimental data on the kinetic properties of L-type calcium channels were integrated into one-dimensional (1D) fiber, 2D sheet, and 3D ventricular models to investigate the pro-arrhythmic effects of CACNA1C mutations by quantifying changes in intracellular calcium handling, action potential profiles, action potential duration restitution (APDR) curves, dispersion of repolarization (DOR), QT interval and spiral wave dynamics. R858H “mutant” L-type calcium current (I_CaL) augmented sarcoplasmic reticulum calcium content, leading to the development of afterdepolarizations at the single cell level and focal activities at the tissue level. It also produced inhomogeneous APD prolongation, causing QT prolongation and repolarization dispersion amplification, rendering R858H “mutant” tissue more vulnerable to the induction of reentry compared with other conditions. In conclusion, altered I_CaL due to the CACNA1C R858H mutation increases arrhythmia risk due to afterdepolarizations and increased tissue vulnerability to unidirectional conduction block. However, the observed reentry is not due to afterdepolarizations (not present in our model), but rather to a novel blocking mechanism.

Optogenetic manipulation of anatomical re-entry by light-guided generation of a reversible local conduction block

Aims

Anatomical re-entry is an important mechanism of ventricular tachycardia, characterized by circular electrical propagation in a fixed pathway. It's current investigative and therapeutic approaches are non-biological, rather unspecific (drugs), traumatizing (electrical shocks), or irreversible (ablation). Optogenetics is a new biological technique that allows reversible modulation of electrical function with unmatched spatiotemporal precision using light-gated ion channels. We therefore investigated optogenetic manipulation of anatomical re-entry in ventricular cardiac tissue.

Methods and results

Transverse, 150-μm-thick ventricular slices, obtained from neonatal rat hearts, were genetically modified with lentiviral vectors encoding Ca2+-translocating channelrhodopsin (CatCh), a light-gated depolarizing ion channel, or enhanced yellow fluorescent protein (eYFP) as control. Stable anatomical re-entry was induced in both experimental groups. Activation of CatCh was precisely controlled by 470-nm patterned illumination, while the effects on anatomical re-entry were studied by optical voltage mapping. Regional illumination in the pathway of anatomical re-entry resulted in termination of arrhythmic activity only in CatCh-expressing slices by establishing a local and reversible, depolarization-induced conduction block in the illuminated area. Systematic adjustment of the size of the light-exposed area in the re-entrant pathway revealed that re-entry could be terminated by either wave collision or extinction, depending on the depth (transmurality) of illumination. In silico studies implicated source-sink mismatches at the site of subtransmural conduction block as an important factor in re-entry termination.

Conclusions

Anatomical re-entry in ventricular tissue can be manipulated by optogenetic induction of a local and reversible conduction block in the re-entrant pathway, allowing effective re-entry termination. These results provide distinctively new mechanistic insight into re-entry termination and a novel perspective for cardiac arrhythmia management.

Spatiotemporal Progression of Early Human Ventricular Fibrillation

Objectives

The objective of this study was to evaluate the spatio-temporal organization and progression of human ventricular fibrillation (VF) in the left (LV) and right (RV) ventricles.

Background

Studies suggest that localized sources contribute to VF maintenance, but the evolution of VF episodes has not been quantified.

Methods

Synchrony between electrograms recorded from 25 patients with induced VF is computed and used to define the Asynchronous Index (ASI), indicating regions which are out-of-step with surrounding tissue. Computer simulations show that ASI can identify the location of VF-maintaining sources, where larger values of ASI_max correlate with more stable sources.

Results

Automated synchrony analysis shows elevated values of ASI in a majority of self-terminating episodes (LV: 8/9, RV: 7/8) and sustained episodes (LV: 11/11, RV: 12/12). The locations of ASI_max in sustained episodes co-localize with rotor cores when rotational activity is simultaneously present in phase maps (LV: 8/8, RV: 5/7, p<.05). The distribution of ASI_max differentiates self-terminating from sustained episodes (mean ASI_max = 0.60±0.14 and 0.70±0.16, respectively; p=0.01). Across sustained episodes the LV exhibits an increase in ASI_max with time.

Conclusions

Quantitative analysis identifies localized asynchronous regions that correlate with sources in VF, with sustained episodes evolving to exhibit more stable activation in the LV. This successive increase in stability indicates a stabilizing agent may be responsible for perpetuating fibrillation in a "migrate-and-capture" mechanism in the LV.

Mechanisms Underlying the Emergence of Post-acidosis Arrhythmia at the Tissue Level: A Theoretical Study

Acidosis has complex electrophysiological effects, which are associated with a high recurrence of ventricular arrhythmias. Through multi-scale cardiac computer modeling, this study investigated the mechanisms underlying the emergence of post-acidosis arrhythmia at the tissue level. In simulations, ten Tusscher-Panfilov ventricular model was modified to incorporate various data on acidosis-induced alterations of cellular electrophysiology and intercellular electrical coupling. The single cell models were incorporated into multicellular one-dimensional (1D) fiber and 2D sheet tissue models. Electrophysiological effects were quantified as changes of action potential profile, sink-source interactions of fiber tissue, and the vulnerability of tissue to the genesis of unidirectional conduction that led to initiation of re-entry. It was shown that acidosis-induced sarcoplasmic reticulum (SR) calcium load contributed to delayed afterdepolarizations (DADs) in single cells. These DADs may be synchronized to overcome the source-sink mismatch arising from intercellular electrotonic coupling, and produce a premature ventricular complex (PVC) at the tissue level. The PVC conduction can be unidirectionally blocked in the transmural ventricular wall with altered electrical heterogeneity, resulting in the genesis of re-entry. In conclusion, altered source-sink interactions and electrical heterogeneity due to acidosis-induced cellular electrophysiological alterations may increase susceptibility to post-acidosis ventricular arrhythmias.

Spiral-wave dynamics in a mathematical model of human ventricular tissue with myocytes and Purkinje fibers

We present systematic numerical studies of the possible effects of the coupling of human endocardial and Purkinje cells at cellular and two-dimensional tissue levels. We find that the autorhythmic-activity frequency of the Purkinje cell in a composite decreases with an increase in the coupling strength; this can even eliminate the autorhythmicity. We observe a delay between the beginning of the action potentials of endocardial and Purkinje cells in a composite; such a delay increases as we decrease the diffusive coupling, and eventually a failure of transmission occurs. An increase in the diffusive coupling decreases the slope of the action-potential-duration-restitution curve of an endocardial cell in a composite. By using a minimal model for the Purkinje network, in which we have a two-dimensional, bilayer tissue, with a layer of Purkinje cells on top of a layer of endocardial cells, we can stabilize spiral-wave turbulence; however, for a sparse distribution of Purkinje-ventricular junctions, at which these two layers are coupled, we can also obtain additional focal activity and many complex transient regimes. We also present additional effects resulting from the coupling of Purkinje and endocardial layers and discuss the relation of our results to the studies performed in anatomically accurate models of the Purkinje network.

Effect of myocyte-fibroblast coupling on the onset of pathological dynamics in a model of ventricular tissue

Managing lethal cardiac arrhythmias is one of the biggest challenges in modern cardiology, and hence it is very important to understand the factors underlying such arrhythmias. While early afterdepolarizations (EAD) of cardiac cells is known to be one such arrhythmogenic factor, the mechanisms underlying the emergence of tissue level arrhythmias from cellular level EADs is not fully understood. Another known arrhythmogenic condition is fibrosis of cardiac tissue that occurs both due to aging and in many types of heart diseases. In this paper we describe the results of a systematic in-silico study, using the TNNP model of human cardiac cells and MacCannell model for (myo)fibroblasts, on the possible effects of diffuse fibrosis on arrhythmias occurring via EADs. We find that depending on the resting potential of fibroblasts (VFR), M-F coupling can either increase or decrease the region of parameters showing EADs. Fibrosis increases the probability of occurrence of arrhythmias after a single focal stimulation and this effect increases with the strength of the M-F coupling. While in our simulations, arrhythmias occur due to fibrosis induced ectopic activity, we do not observe any specific fibrotic pattern that promotes the occurrence of these ectopic sources.

A novel algorithm for ventricular arrhythmia classification using a fuzzy logic approach

In the present study, it has been shown that an unnecessary implantable cardioverter-defibrillator (ICD) shock is often delivered to patients with an ambiguous ECG rhythm in the overlap zone between ventricular tachycardia (VT) and ventricular fibrillation (VF); these shocks significantly increase mortality. Therefore, accurate classification of the arrhythmia into VT, organized VF (OVF) or disorganized VF (DVF) is crucial to assist ICDs to deliver appropriate therapy. A classification algorithm using a fuzzy logic classifier was developed for accurately classifying the arrhythmias into VT, OVF or DVF. Compared with other studies, our method aims to combine ten ECG detectors that are calculated in the time domain and the frequency domain in addition to different levels of complexity for detecting subtle structure differences between VT, OVF and DVF. The classification in the overlap zone between VT and VF is refined by this study to avoid ambiguous identification. The present method was trained and tested using public ECG signal databases. A two-level classification was performed to first detect VT with an accuracy of 92.6 %, and then the discrimination between OVF and DVF was detected with an accuracy of 84.5 %. The validation results indicate that the proposed method has superior performance in identifying the organization level between the three types of arrhythmias (VT, OVF and DVF) and is promising for improving the appropriate therapy choice and decreasing the possibility of sudden cardiac death.

Nonlinear physics of electrical wave propagation in the heart: a review

The beating of the heart is a synchronized contraction of muscle cells (myocytes) that is triggered by a periodic sequence of electrical waves (action potentials) originating in the sino-atrial node and propagating over the atria and the ventricles. Cardiac arrhythmias like atrial and ventricular fibrillation (AF,VF) or ventricular tachycardia (VT) are caused by disruptions and instabilities of these electrical excitations, that lead to the emergence of rotating waves (VT) and turbulent wave patterns (AF,VF). Numerous simulation and experimental studies during the last 20 years have addressed these topics. In this review we focus on the nonlinear dynamics of wave propagation in the heart with an emphasis on the theory of pulses, spirals and scroll waves and their instabilities in excitable media with applications to cardiac modeling. After an introduction into electrophysiological models for action potential propagation, the modeling and analysis of spatiotemporal alternans, spiral and scroll meandering, spiral breakup and scroll wave instabilities like negative line tension and sproing are reviewed in depth and discussed with emphasis on their impact for cardiac arrhythmias.

Pro-arrhythmogenic effects of CACNA1C G1911R mutation in human ventricular tachycardia: insights from cardiac multi-scale models

Mutations in the CACNA1C gene are associated with ventricular tachycardia (VT). Although the CACNA1C mutations were well identified in patients with cardiac arrhythmias, mechanisms by which cardiac arrhythmias are generated in such genetic mutation conditions remain unclear. In this study, we identified a novel mechanism of VT resulted from enhanced repolarization dispersion which is a key factor for arrhythmias in the CACNA1C G1911R mutation using multi-scale computational models of the human ventricle. The increased calcium influx in the mutation prolonged action potential duration (APD), produced steepened action potential duration restitution (APDR) curves as well as augmented membrane potential differences among different cell types during repolarization, increasing transmural dispersion of repolarization (DOR) and the spatial and temporal heterogeneity of cardiac electrical activities. Consequentially, the vulnerability to unidirectional conduction block in response to a premature stimulus increased at tissue level in the G1911R mutation. The increased functional repolarization dispersion anchored reentrant excitation waves in tissue and organ models, facilitating the initiation and maintenance of VT due to less meandering rotor tip. Thus, the increased repolarization dispersion caused by the G1911R mutation is a primary factor that may primarily contribute to the genesis of cardiac arrhythmias in Timothy Syndrome.

Investigating a Novel Activation-Repolarisation Time Metric to Predict Localised Vulnerability to Reentry Using Computational Modelling

Exit sites associated with scar-related reentrant arrhythmias represent important targets for catheter ablation therapy. However, their accurate location in a safe and robust manner remains a significant clinical challenge. We recently proposed a novel quantitative metric (termed the Reentry Vulnerability Index, RVI) to determine the difference between activation and repolarisation intervals measured from pairs of spatial locations during premature stimulation to accurately locate the critical site of reentry formation. In the clinic, the method showed potential to identify regions of low RVI corresponding to areas vulnerable to reentry, subsequently identified as ventricular tachycardia (VT) circuit exit sites. Here, we perform an in silico investigation of the RVI metric in order to aid the acquisition and interpretation of RVI maps and optimise its future usage within the clinic. Within idealised 2D sheet models we show that the RVI produces lower values under correspondingly more arrhythmogenic conditions, with even low resolution (8 mm electrode separation) recordings still able to locate vulnerable regions. When applied to models of infarct scars, the surface RVI maps successfully identified exit sites of the reentrant circuit, even in scenarios where the scar was wholly intramural. Within highly complex infarct scar anatomies with multiple reentrant pathways, the identified exit sites were dependent upon the specific pacing location used to compute the endocardial RVI maps. However, simulated ablation of these sites successfully prevented the reentry re-initiation. We conclude that endocardial surface RVI maps are able to successfully locate regions vulnerable to reentry corresponding to critical exit sites during sustained scar-related VT. The method is robust against highly complex and intramural scar anatomies and low resolution clinical data acquisition. Optimal location of all relevant sites requires RVI maps to be computed from multiple pacing locations.

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.

Models of ventricular arrhythmia mechanisms

The mechanisms that initiate and sustain ventricular arrhythmias in the human heart are clinically important, but hard to study experimentally. In this study, a monodomain model of electrical activation was used to examine how dynamics of electrophysiology at the cell scale influence the surface activation patterns of VF at the tissue scale. Cellular electrophysiology was described with two variants of a phenomenological model of the human ventricular epicardial action potential. The tissue geometry was an 8.0 × 8.0 × 1.2 cm 3D tissue slab with axially symmetric anisotropy. In both cases an initial re-entrant wave fragmented into multiple wavelets of activation. The model variant with steep action potential duration restitution produced much more complex activation, with a greater average number of filaments (13.79) than the variant with less steep restitution (3.08). More complex activation was associated with proportionally fewer transmural filaments, and so the average number of epicardial wavefronts and phase singularities per filament was lower. The average number of epicardial phase singularities and wavefronts for the model variant with less steep restitution were consistent with experimental observations in the human heart. This study shows that small changes in cell scale dynamics can have a large influence on the complexity of re-entrant activation in simulated 3D tissue, as well as on the features observed on the epicardial surface.

Ventricular Tachycardia and Early Fibrillation in Patients With Brugada Syndrome and Ischemic Cardiomyopathy Show Predictable Frequency-Phase Properties on the Precordial ECG Consistent With the Respective Arrhythmogenic Substrate

Supplemental Digital Content is available in the text.

Ventricular fibrillation (VF) has been proposed to be maintained by localized high-frequency sources. We tested whether spectral-phase analysis of the precordial ECG enabled identification of periodic activation patterns generated by such sources.

Analysis of cardiac fibrillation using phase mapping

The sequence of myocardial electrical activation during fibrillation is complex and changes with each cycle. Phase analysis represents the electrical activation-recovery process as an angle. Lines of equal phase converge at a phase singularity at the center of rotation of a reentrant wave, and the identification of reentry and tracking of reentrant wavefronts can be automated. We examine the basic ideas behind phase analysis. With the exciting prospect of using phase analysis of atrial electrograms to guide ablation in the human heart, we highlight several recent developments in preprocessing electrograms so that phase can be estimated reliably.

Influence of the medium's dimensionality on defect-mediated turbulence

Spatiotemporal chaos in oscillatory and excitable media is often characterized by the presence of phase singularities called defects. Understanding such defect-mediated turbulence and its dependence on the dimensionality of a given system is an important challenge in nonlinear dynamics. This is especially true in the context of ventricular fibrillation in the heart, where the importance of the thickness of the ventricular wall is contentious. Here, we study defect-mediated turbulence arising in two different regimes in a conceptual model of excitable media and investigate how the statistical character of the turbulence changes if the thickness of the medium is changed from (quasi-) two- dimensional to three dimensional. We find that the thickness of the medium does not have a significant influence in, far from onset, fully developed turbulence while there is a clear transition if the system is close to a spiral instability. We provide clear evidence that the observed transition and change in the mechanism that drives the turbulent behavior is purely a consequence of the dimensionality of the medium. Using filament tracking, we further show that the statistical properties in the three-dimensional medium are different from those in turbulent regimes arising from filament instabilities like the negative line tension instability. Simulations also show that the presence of this unique three-dimensional turbulent dynamics is not model specific.

Quantifying inter-species differences in contractile function through biophysical modelling

Animal models and measurements are frequently used to guide and evaluate clinical interventions. In this context, knowledge of inter-species differences in physiology is crucial for understanding the limitations and relevance of animal experimental assays for informing clinical applications. Extensive effort has been put into studying the structure and function of cardiac contractile proteins and how differences in these translate into the functional properties of muscles. However, integrating this knowledge into a quantitative description, formalising and highlighting inter-species differences both in the kinetics and in the regulation of physiological mechanisms, remains challenging. In this study we propose and apply a novel approach for the quantification of inter-species differences between mouse, rat and human. Assuming conservation of the fundamental physiological mechanisms underpinning contraction, biophysically based computational models are fitted to simulate experimentally recorded phenotypes from multiple species. The phenotypic differences between species are then succinctly quantified as differences in the biophysical model parameter values. This provides the potential of quantitatively establishing the human relevance of both animal-based experimental and computational models for application in a clinical context. Our results indicate that the parameters related to the sensitivity and cooperativity of calcium binding to troponin C and the activation and relaxation rates of tropomyosin/crossbridge binding kinetics differ most significantly between mouse, rat and human, while for example the reference tension, as expected, shows only minor differences between the species. Hence, while confirming expected inter-species differences in calcium sensitivity due to large differences in the observed calcium transients, our results also indicate more unexpected differences in the cooperativity mechanism. Specifically, the decrease in the unbinding rate of calcium to troponin C with increasing active tension was much lower for mouse than for rat and human. Our results also predicted crossbridge binding to be slowest in human and fastest in mouse.

Optical Mapping of Ventricular Fibrillation Dynamics

There is very limited information regarding the dynamic patterns of the electrical activity during ventricular fibrillation (VF) in humans. Most of the data used to generate and test hypotheses regarding the mechanisms of VF come from animal models and computer simulations and the quantification of VF patterns is non-trivial. Many of the experimental recordings of the dynamic spatial patterns of VF have been obtained from mammals using "optical mapping" or "video imaging" technology in which "phase maps" are derived from high-resolution transmembrane recordings from the heart surface. The surface manifestation of the unstable reentrant waves sustaining VF can be identified as "phase singularities" and their number and location provide one measure of VF complexity. After providing a brief history of optical mapping of VF, we compare and contrast a quantitative analysis of VF patterns from the heart surface for four different animal models, hence providing physiological insight into the variety of VF dynamics among species. We found that in all four animal models the action potential duration restitution slope was actually negative during VF and that the spatial dispersion of electrophysiological parameters were not different during the first second of VF compared to pacing immediately before VF initiation. Surprisingly, our results suggest that APD restitution and spatial dispersion may not be essential causes of VF dynamics. Analyses of electrophysiological quantities in the four animal models are consistent with the idea that VF is essentially a two-dimensional phenomenon in small rabbit hearts whose size are near the boundary of the "critical mass" required to sustain VF, while VF in large pig hearts is three-dimensional and exhibits the maximal theoretical phase singularity density, and thus will not terminate spontaneously.

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.

Images as drivers of progress in cardiac computational modelling

Computational models have become a fundamental tool in cardiac research. Models are evolving to cover multiple scales and physical mechanisms. They are moving towards mechanistic descriptions of personalised structure and function, including effects of natural variability. These developments are underpinned to a large extent by advances in imaging technologies. This article reviews how novel imaging technologies, or the innovative use and extension of established ones, integrate with computational models and drive novel insights into cardiac biophysics. In terms of structural characterization, we discuss how imaging is allowing a wide range of scales to be considered, from cellular levels to whole organs. We analyse how the evolution from structural to functional imaging is opening new avenues for computational models, and in this respect we review methods for measurement of electrical activity, mechanics and flow. Finally, we consider ways in which combined imaging and modelling research is likely to continue advancing cardiac research, and identify some of the main challenges that remain to be solved.

Developing a novel comprehensive framework for the investigation of cellular and whole heart electrophysiology in the in situ human heart: historical perspectives, current progress and future prospects

Understanding the mechanisms of fatal ventricular arrhythmias is of great importance. In view of the many electrophysiological differences that exist between animal species and humans, the acquisition of basic electrophysiological data in the intact human heart is essential to drive and complement experimental work in animal and in-silico models. Over the years techniques have been developed to obtain basic electrophysiological signals directly from the patients by incorporating these measurements into routine clinical procedures which access the heart such as cardiac catheterisation and cardiac surgery. Early recordings with monophasic action potentials provided valuable information including normal values for the in vivo human heart, cycle length dependent properties, the effect of ischaemia, autonomic nervous system activity, and mechano-electric interaction. Transmural recordings addressed the controversial issue of the mid myocardial "M" cell. More recently, the technique of multielectrode mapping (256 electrodes) developed in animal models has been extended to humans, enabling mapping of activation and repolarisation on the entire left and right ventricular epicardium in patients during cardiac surgery. Studies have examined the issue of whether ventricular fibrillation was driven by a "mother" rotor with inhomogeneous and fragmented conduction as in some animal models, or by multiple wavelets as in other animal studies; results showed that both mechanisms are operative in humans. The simpler spatial organisation of human VF has important implications for treatment and prevention. To link in-vivo human electrophysiological mapping with cellular biophysics, multielectrode mapping is now being combined with myocardial biopsies. This technique enables region-specific electrophysiology changes to be related to underlying cellular biology, for example: APD alternans, which is a precursor of VF and sudden death. The mechanism is incompletely understood but related to calcium cycling and APD restitution. Multielectrode sock mapping during incremental pacing enables epicardial sites to be identified which exhibit marked APD alternans and sites where APD alternans is absent. Whole heart electrophysiology is assessed by activation repolarisation mapping and analysis is performed immediately on-site in order to guide biopsies to specific myocardial sites. Samples are analysed for ion channel expression, Ca(2+)-handling proteins, gap junctions and extracellular matrix. This new comprehensive approach to bridge cellular and whole heart electrophysiology allowed to identify 20 significant changes in mRNA for ion channels Ca(2+)-handling proteins, a gap junction channel, a Na(+)-K(+) pump subunit and receptors (particularly Kir 2.1) between the positive and negative alternans sites.

Mechanism of reentry induction by a 9-V battery in rabbit ventricles

Although the application of a 9-V battery to the epicardial surface is a simple method of ventricular fibrillation induction, the fundamental mechanisms underlying this process remain unstudied. We used a combined experimental and modelling approach to understand how the interaction of direct current (DC) from a battery may induce reentrant activity within rabbit ventricles and its dependence on battery application timing and duration. A rabbit ventricular computational model was used to simulate 9-V battery stimulation for different durations at varying onset times during sinus rhythm. Corresponding high-resolution optical mapping measurements were conducted on rabbit hearts with DC stimuli applied via a relay system. DC application to diastolic tissue induced anodal and cathodal make excitations in both simulations and experiments. Subsequently, similar static epicardial virtual electrode patterns were formed that interacted with sinus beats but did not induce reentry. Upon battery release during diastole, break excitations caused single ectopics, similar to application, before sinus rhythm resumed. Reentry induction was possible for short battery applications when break excitations were slowed and forced to take convoluted pathways upon interaction with refractory tissue from prior make excitations or sinus beats. Short-lived reentrant activity could be induced for battery release shortly after a sinus beat for longer battery applications. In conclusion, the application of a 9-V battery to the epicardial surface induces reentry through a complex interaction of break excitations after battery release with prior induced make excitations or sinus beats.

Bridging the gap between computation and clinical biology: validation of cable theory in humans

Introduction: Computerized simulations of cardiac activity have significantly contributed to our understanding of cardiac electrophysiology, but techniques of simulations based on patient-acquired data remain in their infancy. We sought to integrate data acquired from human electrophysiological studies into patient-specific models, and validated this approach by testing whether electrophysiological responses to sequential premature stimuli could be predicted in a quantitatively accurate manner.

Methods: Eleven patients with structurally normal hearts underwent electrophysiological studies. Semi-automated analysis was used to reconstruct activation and repolarization dynamics for each electrode. This S₂ extrastimuli data was used to inform individualized models of cardiac conduction, including a novel derivation of conduction velocity restitution. Activation dynamics of multiple premature extrastimuli were then predicted from this model and compared against measured patient data as well as data derived from the ten-Tusscher cell-ionic model.

Results: Activation dynamics following a premature S₃ were significantly different from those after an S₂. Patient specific models demonstrated accurate prediction of the S₃ activation wave, (Pearson's R² = 0.90, median error 4%). Examination of the modeled conduction dynamics allowed inferences into the spatial dispersion of activation delay. Further validation was performed against data from the ten-Tusscher cell-ionic model, with our model accurately recapitulating predictions of repolarization times (R² = 0.99).

Conclusions: Simulations based on clinically acquired data can be used to successfully predict complex activation patterns following sequential extrastimuli. Such modeling techniques may be useful as a method of incorporation of clinical data into predictive models.

Rotors and the dynamics of cardiac fibrillation

The objective of this article is to present a broad review of the role of cardiac electric rotors and their accompanying spiral waves in the mechanism of cardiac fibrillation. At the outset, we present a brief historical overview regarding reentry and then discuss the basic concepts and terminologies pertaining to rotors and their initiation. Thereafter, the intrinsic properties of rotors and spiral waves, including phase singularities, wavefront curvature, and dominant frequency maps, are discussed. The implications of rotor dynamics for the spatiotemporal organization of fibrillation, independent of the species being studied, are described next. The knowledge gained regarding the role of cardiac structure in the initiation or maintenance of rotors and the ionic bases of spiral waves in the past 2 decades, as well as the significance for drug therapy, is reviewed subsequently. We conclude by examining recent evidence suggesting that rotors are critical in sustaining both atrial and ventricular fibrillation in the human heart and its implications for treatment with radiofrequency ablation.

Novel nonsurgical left ventricular assist device and system

Treatment options for advanced stages of congestive heart failure remain limited. Left ventricular assist devices (LVADs) have emerged as a means to support failing circulation. However, these devices are not without significant risk such as major open chest surgery. We utilized a novel approach for device placement at the aorto-left atria continuity as a site to create a conduit capable of accommodating a percutaneous LVAD system. We designed and developed an expandable nitinol based device for placement at this site to create a shunt between the LA and aorta. Our experiments support this anatomic location as an accessible and feasible site for accommodation of an entirely percutaneous LVAD. The novelty of this approach would bypass the left ventricle, and thereby minimize complications and morbidities associated with current LVAD placement.

The functional role of electrophysiological heterogeneity in the rabbit ventricle during rapid pacing and arrhythmias

Electrophysiological heterogeneity in action potential recordings from healthy intact hearts remains highly variable and, where present, is almost entirely abolished at fast pacing rates. Consequently, the functional importance of intrinsic action potential duration (APD) heterogeneity in healthy ventricles, and particularly its role during rapidly activating reentrant arrhythmias, remain poorly understood. By incorporating both transmural and apicobasal APD heterogeneity within a biventricular rabbit computational model and comparing with an equivalent homogeneous model, we directly investigated the functional importance of intrinsic APD heterogeneity under fast pacing and arrhythmogenic protocols. Although differences in APD were significantly modulated at the tissue level during pacing and further reduced as pacing frequency increased, small differences were still noticeable. Such differences were further marginally accentuated/attenuated via electrotonic effects relative to wavefront propagation directions. The remaining small levels of APD heterogeneity under the fastest pacing frequencies resulted in arrhythmia initiation via heterogeneous conduction block, in contrast to complete block in the homogeneous model. Such induction mechanisms were more evident during premature stimuli at slower paced rhythms where intrinsic heterogeneity remained to a greater degree. During sustained arrhythmias, however, intrinsic heterogeneity made little difference to overall reentrant behavior, either visually, or in terms of duration, metrics quantifying filament/phase singularity dynamics, and global electrocardiogram characteristics. These findings suggest that, despite being important during arrhythmia initiation, intrinsic electrophysiological heterogeneity plays little functional role during rapid pacing and sustained arrhythmia dynamics in the healthy ventricle and thus questions the need to incorporate such detail in computational models when simulating rapid arrhythmias.

Automated signal pattern detection in ECG during human ventricular arrhythmias

Ventricular arrhythmias seriously affects cardiac function. Of these arrhythmias, Ventricular fibrillation is considered as a lethal cardiac condition. Recent studies have reported that ventricular arrhythmias are not completely random and may exhibit regional spatio-temporal organizations. These organizations could be indicative of reoccurring signal patterns and might be embedded within the surface electrocardiograms (ECGs) during ventricular arrhythmias. In this work, we aim to identify such reoccurring ECG signal patterns during ventricular arrhythmias. The detection of such signal patterns and their distribution could be of help in sub-classifying the affected population for better targeted diagnosis and treatment. Our analysis on 14 ECG segments (on average 3.24 minutes per segment) obtained from the MIT-BIH ventricular arrhythmia database identified three reoccurring signal patterns. A wavelet based technique was developed for automating the pattern identification process using ECGs. The proposed method achieved automated detection accuracies of 73.3%, 75.0% and 86.6% for the proposed signal patterns.

Ventricular fibrillation conduction through an isthmus of preserved myocardium between radiofrequency lesions

BACKGROUND

Selective local acceleration of myocardial activation during ventricular fibrillation (VF) contributes information on the interactions between neighboring zones during the arrhythmia. This study analyzes these interactions, centering the observations on an isthmus of myocardium between two radiofrequency (RF) lesions.

METHODS

In nine isolated rabbit hearts, a gap of preserved myocardium was established between two RF lesions in the anterolateral left ventricle (LV) wall. Before, during, and after increasing the spatial heterogeneity of VF by local myocardial stretching, VF epicardial recordings were obtained.

RESULTS

Local stretch in the anterior LV wall decreased the excitable window (17 ± 7 ms vs 26 ± 7 ms; P < 0.05) and increased the dominant frequency (DFr; 18.9 ± 5.0 Hz vs 15.2 ± 3.6 Hz; P < 0.05) in this zone, without changes in the non-stretched posterolateral zone (25 ± 4 ms vs 27 ± 6 ms, ns and 14.1 ± 2.7 Hz vs 14.3 ± 3.0 Hz, ns). The DFr ratio at both sides of the gap was inversely correlated to the excitable window ratio (R = -0.57; P = 0.002). Before (31% vs 26%), during (29% vs 22%), and after stretch suppression (35% vs 25%), the wavefronts passing through the gap from the posterolateral to the anterior LV wall were seen to predominate. The number of wavefronts that passed from the anterior to the posterolateral LV wall was related to the excitable window in this zone (R = 0.41; P = 0.03).

CONCLUSIONS

The VF acceleration induced in the stretched zone does not increase the flow of wavefronts toward the non-stretched zone in the adjacent gap of preserved myocardium. The absence of significant changes in the electrophysiological parameters of the non-stretched myocardium limits the arrival of wavefronts in this zone.

Evaluation of the complexity of myocardial activation during ventricular fibrillation. An experimental study

INTRODUCTION AND OBJECTIVES

An experimental model is used to analyze the characteristics of ventricular fibrillation in situations of variable complexity, establishing relationships among the data produced by different methods for analyzing the arrhythmia.

METHODS

In 27 isolated rabbit heart preparations studied under the action of drugs (propranolol and KB-R7943) or physical procedures (stretching) that produce different degrees of change in the complexity of myocardial activation during ventricular fibrillation, use was made of spectral, morphological, and mapping techniques to process the recordings obtained with epicardial multielectrodes.

RESULTS

The complexity of ventricular fibrillation assessed by mapping techniques was related to the dominant frequency, normalized spectral energy, signal regularity index, and their corresponding coefficients of variation, as well as the area of the regions of interest identified on the basis of these parameters. In the multivariate analysis, we used as independent variables the area of the regions of interest related to the spectral energy and the coefficient of variation of the energy (complexity index=-0.005×area of the spectral energy regions -2.234×coefficient of variation of the energy+1.578; P=.0001; r=0.68).

CONCLUSIONS

The spectral and morphological indicators and, independently, those derived from the analysis of normalized energy regions of interest provide a reliable approach to the evaluation of the complexity of ventricular fibrillation as an alternative to complex mapping techniques.

Prolongation of minimal action potential duration in sustained fibrillation decreases complexity by transient destabilization

AIMS

Sustained ventricular fibrillation (VF) is maintained by multiple stable rotors. Destabilization of sustained VF could be beneficial by affecting VF complexity (defined by the number of rotors). However, underlying mechanisms affecting VF stability are poorly understood. Therefore, the aim of this study was to correlate changes in arrhythmia complexity with changes in specific electrophysiological parameters, allowing a search for novel factors and underlying mechanisms affecting stability of sustained VF.

METHODS AND RESULTS

Neonatal rat ventricular cardiomyocyte monolayers and Langendorff-perfused adult rat hearts were exposed to increasing dosages of the gap junctional uncoupler 2-aminoethoxydiphenyl borate (2-APB) to induce arrhythmias. Ion channel blockers/openers were added to study effects on VF stability. Electrophysiological parameters were assessed by optical mapping and patch-clamp techniques. Arrhythmia complexity in cardiomyocyte cultures increased with increasing dosages of 2-APB (n > 38), leading to sustained VF: 0.0 ± 0.1 phase singularities/cm(2) in controls vs. 0.0 ± 0.1, 1.0 ± 0.9, 3.3 ± 3.2, 11.0 ± 10.1, and 54.3 ± 21.7 in 5, 10, 15, 20, and 25 µmol/L 2-APB, respectively. Arrhythmia complexity inversely correlated with wavelength. Lengthening of wavelength during fibrillation could only be induced by agents (BaCl(2)/BayK8644) increasing the action potential duration (APD) at maximal activation frequencies (minimal APD); 123 ± 32%/117 ± 24% of control. Minimal APD prolongation led to transient VF destabilization, shown by critical wavefront collision leading to rotor termination, followed by significant decreases in VF complexity and activation frequency (52%/37%). These key findings were reproduced ex vivo in rat hearts (n = 6 per group).

CONCLUSION

These results show that stability of sustained fibrillation is regulated by minimal APD. Minimal APD prolongation leads to transient destabilization of fibrillation, ultimately decreasing VF complexity, thereby providing novel insights into anti-fibrillatory mechanisms.

The role of fine-scale anatomical structure in the dynamics of reentry in computational models of the rabbit ventricles

Fine-scale anatomical structures in the heart may play an important role in sustaining cardiac arrhythmias. However, the extent of this role and how it may differ between species are not fully understood. In this study we used computational modelling to assess the impact of anatomy upon arrhythmia maintenance in the rabbit ventricles. Specifically, we quantified the dynamics of excitation wavefronts during episodes of simulated tachyarrhythmias and fibrillatory arrhythmias, defined as being respectively characterised by relatively low and high spatio-temporal disorganisation.Two computational models were used: a highly anatomically detailed MR-derived rabbit ventricular model (representing vasculature, endocardial structures) and a simplified equivalent model, constructed from the same MR-data but lacking such fine-scale anatomical features. During tachyarrhythmias, anatomically complex and simplified models showed very similar dynamics; however, during fibrillatory arrhythmias, as activation wavelength decreased, the presence of fine-scale anatomical details appeared to marginally increase disorganisation of wavefronts during arrhythmias in the complex model. Although a small amount of clustering of reentrant rotor centres (filaments) around endocardial structures was witnessed in follow-up analysis (which slightly increased during fibrillation as rotor size decreased), this was significantly less than previously reported in large animals. Importantly, no anchoring of reentrant rotors was visibly identifiable in arrhythmia movies. These differences between tachy- and fibrillatory arrhythmias suggest that the relative size of reentrant rotors with respect to anatomical obstacles governs the influence of fine-scale anatomy in the maintenance of ventricular arrhythmias in the rabbit. In conclusion, our simulations suggest that fine-scale anatomical features play little apparent role in the maintenance of tachyarrhythmias in the rabbit ventricles and, contrary to experimental reports in larger animals, appear to play only a minor role in the maintenance of fibrillatory arrhythmias. These findings also have important implications in optimising the level of detail required in anatomical computational meshes frequently used in arrhythmia investigations.

Bridging experiments, models and simulations: an integrative approach to validation in computational cardiac electrophysiology

Computational models in physiology often integrate functional and structural information from a large range of spatiotemporal scales from the ionic to the whole organ level. Their sophistication raises both expectations and skepticism concerning how computational methods can improve our understanding of living organisms and also how they can reduce, replace, and refine animal experiments. A fundamental requirement to fulfill these expectations and achieve the full potential of computational physiology is a clear understanding of what models represent and how they can be validated. The present study aims at informing strategies for validation by elucidating the complex interrelations among experiments, models, and simulations in cardiac electrophysiology. We describe the processes, data, and knowledge involved in the construction of whole ventricular multiscale models of cardiac electrophysiology. Our analysis reveals that models, simulations, and experiments are intertwined, in an assemblage that is a system itself, namely the model-simulation-experiment (MSE) system. We argue that validation is part of the whole MSE system and is contingent upon 1) understanding and coping with sources of biovariability; 2) testing and developing robust techniques and tools as a prerequisite to conducting physiological investigations; 3) defining and adopting standards to facilitate the interoperability of experiments, models, and simulations; 4) and understanding physiological validation as an iterative process that contributes to defining the specific aspects of cardiac electrophysiology the MSE system targets, rather than being only an external test, and that this is driven by advances in experimental and computational methods and the combination of both.

A novel, minimally invasive, segmental myocardial infarction with a clear healed infarct borderzone in rabbits

Ventricular arrhythmias in the setting of a healed myocardial infarction have been studied to a much lesser degree than acute and subacute infarction, due to the pericardial scarring, which results from the traditional open-chest techniques used for myocardial infarction (MI) induction. We sought to develop a segmental MI with low perioperative mortality in the rabbit that allows optimal visualization and therefore improved study of the infarction borderzone. Rabbits underwent MI using endovascular coil occlusion of the first obtuse marginal artery. Three weeks postprocedure, we evaluated our model by echocardiography and electrophysiology studies, optical mapping of isolated hearts, and histological studies. Seventeen rabbits underwent the protocol (12 MI and 5 sham) with a 92% survival to completion of the study (11 MI and 5 sham). MI rabbits demonstrated wall motion abnormalities on echocardiography while shams did not. At electrophysiological study, two MI rabbits had inducible ventricular tachycardia and one had inducible ventricular fibrillation. Isolated hearts demonstrated no pericardial scarring with a smooth, easily identifiable infarct borderzone. Optical mapping of the borderzone region showed successful mapping of peri-infarct reentry formation, with ventricular fibrillation inducible in 11 of 11 MI hearts and 1 of 5 sham hearts. We demonstrate successful high resolution mapping in the borderzone, showing delayed conduction in this region corresponding to late deflections in the QRS on ECG. We report the successful development of a minimally invasive MI via targeted coil delivery to the obtuse marginal artery with an exceptionally high rate of procedural survival and an arrhythmogenic phenotype. This model mimics human post-MI on echocardiography, gross pathology, histology, and electrophysiology.