A model of canine purkinje cell electrophysiology and Ca(2+) cycling: rate dependence, triggered activity, and comparison to ventricular myocytes

Circadian regulation of sinoatrial nodal cell pacemaking function: Dissecting the roles of autonomic control, body temperature, and local circadian rhythmicity

Strong circadian (~24h) rhythms in heart rate (HR) are critical for flexible regulation of cardiac pacemaking function throughout the day. While this circadian flexibility in HR is sustained in diverse conditions, it declines with age, accompanied by reduced maximal HR performance. The intricate regulation of circadian HR involves the orchestration of the autonomic nervous system (ANS), circadian rhythms of body temperature (CRBT), and local circadian rhythmicity (LCR), which has not been fully understood. Here, we developed a mathematical model describing ANS, CRBT, and LCR in sinoatrial nodal cells (SANC) that accurately captures distinct circadian patterns in adult and aged mice. Our model underscores how the alliance among ANS, CRBT, and LCR achieves circadian flexibility to cover a wide range of firing rates in SANC, performance to achieve maximal firing rates, while preserving robustness to generate rhythmic firing patterns irrespective of external conditions. Specifically, while ANS dominates in promoting SANC flexibility and performance, CRBT and LCR act as primary and secondary boosters, respectively, to further enhance SANC flexibility and performance. Disruption of this alliance with age results in impaired SANC flexibility and performance, but not robustness. This unexpected outcome is primarily attributed to the age-related reduction in parasympathetic activities, which maintains SANC robustness while compromising flexibility. Our work sheds light on the critical alliance of ANS, CRBT, and LCR in regulating time-of-day cardiac pacemaking function and dysfunction, offering insights into novel therapeutic targets for the prevention and treatment of cardiac arrhythmias.

Moderator Band and Ventricular Tachycardia: Structural or Functional Substrate?

The moderator band (MB) is an intracavitary structure of the right ventricle composed of muscular fibers encompassing specialized Purkinje fibers, separated each other by collagen and adipose tissue. In the last decades, premature ventricular complexes originating within the Purkinje network have been implicated in the genesis of life-threatening arrhythmias. However, right Purkinje network arrhythmias have been much less reported in the literature compared to the left counterpart. The MB has unique anatomical and electrophysiological properties, which may account for its arrhythmogenicity and may be responsible for a significant portion of idiopathic ventricular fibrillation. MB embodies autonomic nervous system cells, with important implications in arrhythmogenesis. Some idiopathic ventricular arrhythmias, defined as the absence of any identifiable structural heart disorder, can begin from this site. Due to these complex structural and functional peculiarities strictly interplayed each other, it is arduous to determine the precise mechanism underlying MB arrhythmias. MB-related arrhythmias should be differentiated from other right Purkinje fibers arrhythmias because of the opportunity for intervention and the unusual site for the ablation poorly described in the literature. In the current paper, we report the characteristics and electrical properties of the MB, their involvement in arrhythmogenesis, clinical and electrophysiological peculiarities of MB-related arrhythmias, and current treatment options.

Spatiotemporal patterns of early afterdepolarizations underlying abnormal T-wave morphologies in a tissue model of the Purkinje-ventricular system

Sudden cardiac death (SCD) is a leading cause of death worldwide, and the majority of SCDs are caused by acute ventricular arrhythmias (VAs). Early afterdepolarizations (EADs) are an important trigger of VA under pathological conditions, e.g., inherited or acquired long QT syndrome (LQTS). However, it remains unclear how EAD events at the cellular level are spatially organized at the tissue level to induce and maintain ventricular arrhythmias and whether the spatial-temporal patterns of EADs at the tissue level are associated with abnormal T-wave morphologies that are often observed in LQTS, such as broad-based, notched or bifid; late appearance; and pointed T-waves. Here, a tissue model of the Purkinje-ventricular system (PVS) was developed to quantitatively investigate the complex spatial-temporal dynamics of EADs during T-wave abnormalities. We found that (1) while major inhibition of ICaL can substantially reduce the excitability of the PVS leading to conduction failures, moderate ICaL inhibition can promote occurrences of AP alternans at short cycle lengths (CLs), and EAD events preferentially occur with a major reduction of IKr (>50%) at long CLs; (2) with a minor reduction of ICaL, spatially synchronized steady-state EAD events with inverted and biphasic T-waves can be "weakened" into beat-to-beat concurrences of spatially synchronized EADs and T-wave alternans, and as pacing CLs increase, beat-to-beat concurrences of localized EADs with late-appearing and pointed T-wave morphologies can be observed; (3) under certain conditions, localized EAD events in the midmyocardium may trigger slow uni-directional electric propagation with inverted (antegrade) or upright (retrograde) broad-based T-waves; (4) spatially discordant EADs were typically characterized by desynchronized spontaneous onsets of EAD events between two groups of PVS tissues with biphasic T-wave morphologies, and they can evolve into spatially discordant oscillating EAD patterns with sustained or self-terminated alternating EAD and electrocardiogram (ECG) patterns. Our results provide new insights into the spatiotemporal aspects of the onset and development of EADs and suggest possible mechanistic links between the complex spatial dynamics of EADs and T-wave morphologies.

Models of the cardiac L-type calcium current: A quantitative review

The L-type calcium current (I CaL ) plays a critical role in cardiac electrophysiology, and models ofI CaL are vital tools to predict arrhythmogenicity of drugs and mutations. Five decades of measuring and modelingI CaL have resulted in several competing theories (encoded in mathematical equations). However, the introduction of new models has not typically been accompanied by a data-driven critical comparison with previous work, so that it is unclear which model is best suited for any particular application. In this review, we describe and compare 73 published mammalianI CaL models and use simulated experiments to show that there is a large variability in their predictions, which is not substantially diminished when grouping by species or other categories. We provide model code for 60 models, list major data sources, and discuss experimental and modeling work that will be required to reduce this huge list of competing theories and ultimately develop a community consensus model ofI CaL . This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Molecular and Cellular Physiology.

A Novel In Silico Electromechanical Model of Human Ventricular Cardiomyocyte

Contractility has become one of the main readouts in computational and experimental studies on cardiomyocytes. Following this trend, we propose a novel mathematical model of human ventricular cardiomyocytes electromechanics, BPSLand, by coupling a recent human contractile element to the BPS2020 model of electrophysiology. BPSLand is the result of a hybrid optimization process and it reproduces all the electrophysiology experimental indices captured by its predecessor BPS2020, simultaneously enabling the simulation of realistic human active tension and its potential abnormalities. The transmural heterogeneity in both electrophysiology and contractility departments was simulated consistent with previous computational and in vitro studies. Furthermore, our model could capture delayed afterdepolarizations (DADs), early afterdepolarizations (EADs), and contraction abnormalities in terms of aftercontractions triggered by either drug action or special pacing modes. Finally, we further validated the mechanical results of the model against previous experimental and in silico studies, e.g., the contractility dependence on pacing rate. Adding a new level of applicability to the normative models of human cardiomyocytes, BPSLand represents a robust, fully-human in silico model with promising capabilities for translational cardiology.

Bifurcations and Proarrhythmic Behaviors in Cardiac Electrical Excitations

The heart is a hierarchical dynamic system consisting of molecules, cells, and tissues, and acts as a pump for blood circulation. The pumping function depends critically on the preceding electrical activity, and disturbances in the pattern of excitation propagation lead to cardiac arrhythmia and pump failure. Excitation phenomena in cardiomyocytes have been modeled as a nonlinear dynamical system. Because of the nonlinearity of excitation phenomena, the system dynamics could be complex, and various analyses have been performed to understand the complex dynamics. Understanding the mechanisms underlying proarrhythmic responses in the heart is crucial for developing new ways to prevent and control cardiac arrhythmias and resulting contractile dysfunction. When the heart changes to a pathological state over time, the action potential (AP) in cardiomyocytes may also change to a different state in shape and duration, often undergoing a qualitative change in behavior. Such a dynamic change is called bifurcation. In this review, we first summarize the contribution of ion channels and transporters to AP formation and our knowledge of ion-transport molecules, then briefly describe bifurcation theory for nonlinear dynamical systems, and finally detail its recent progress, focusing on the research that attempts to understand the developing mechanisms of abnormal excitations in cardiomyocytes from the perspective of bifurcation phenomena.

The Purkinje-myocardial junction is the anatomic origin of ventricular arrhythmia in CPVT

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is an arrhythmia syndrome caused by gene mutations that render RYR2 Ca release channels hyperactive, provoking spontaneous Ca release and delayed afterdepolarizations (DADs). What remains unknown is the cellular source of ventricular arrhythmia triggered by DADs: Purkinje cells in the conduction system or ventricular cardiomyocytes in the working myocardium. To answer this question, we used a genetic approach in mice to knock out cardiac calsequestrin either in Purkinje cells or in ventricular cardiomyocytes. Total loss of calsequestrin in the heart causes a severe CPVT phenotype in mice and humans. We found that loss of calsequestrin only in ventricular myocytes produced a full-blown CPVT phenotype, whereas mice with loss of calsequestrin only in Purkinje cells were comparable to WT mice. Subendocardial chemical ablation or restoration of calsequestrin expression in subendocardial cardiomyocytes neighboring Purkinje cells was sufficient to protect against catecholamine-induced arrhythmias. In silico modeling demonstrated that DADs in ventricular myocardium can trigger full action potentials in the Purkinje fiber, but not vice versa. Hence, ectopic beats in CPVT are likely generated at the Purkinje-myocardial junction via a heretofore unrecognized tissue mechanism, whereby DADs in the ventricular myocardium trigger full action potentials in adjacent Purkinje cells.

Electrocardiographic Characteristics and Catheter Ablation of Ventricular Arrhythmias Originating From the Moderator Band in Children

AIMS

To investigate the electrocardiographic (ECG) characteristics and catheter ablation of ventricular arrhythmias (VAs) originating from the moderator band (MB) in children.

METHODS

A total of six children who had VAs originating from the MB-as confirmed by electrophysiological study-and who underwent catheter ablation between January 2016 and December 2020 were retrospectively reviewed. During the procedure, a three-dimensional electroanatomic mapping system was used to facilitate three-dimensional anatomical reconstruction, mapping and ablation. Patients' clinical characteristics, ECG features and procedural data were collected and analyzed.

RESULTS

The mean age was 8.4 ± 2.6 years (range: 5.3-11 years) and mean weight was 27.7 ± 11.4 kg (range: 17-47 kg). Four patients presented with frequent premature ventricular contraction (PVC), one patient presented with frequent PVC and non-sustained ventricular tachycardia, and one patient presented with sustained monomorphic ventricular tachycardia. The QRS duration averaged 126.3 ± 4.6 ms. In all patients, the VAs had left bundle branch block QRS with left superior frontal plane axes, rapid downstrokes of the QRS in the precordial leads, and late precordial transitions (>V₄). During the same period, 10 cases of VAs originated from the posterior-lateral wall of the tricuspid annulus, with a mean QRS duration of 152.8 ± 6.4 ms. Compared to that, VAs of MB origin have narrower QRS widths, downstroke slopes in the inferior lead, sharper downstroke slopes in the precordial lead, and smaller R-wave amplitudes in the V₆ lead. All patients experienced immediate ablation success with activations earlier than QRS by 26.0 ± 3.5 ms, and no procedural complications occurring. Only one case had recurrent PVC during a follow-up period ranging from 6 to 36 months.

CONCLUSION

MB VAs in children have distinctive ECG morphology and electrophysiological characteristics. Catheter ablation using a three-dimensional electroanatomic mapping system is safe and effective in these patients.

Ranolazine: An Old Drug with Emerging Potential; Lessons from Pre-Clinical and Clinical Investigations for Possible Repositioning

Ischemic heart disease is a significant public health problem with high mortality and morbidity. Extensive scientific investigations from basic sciences to clinics revealed multilevel alterations from metabolic imbalance, altered electrophysiology, and defective Ca²⁺/Na⁺ homeostasis leading to lethal arrhythmias. Despite the recent identification of numerous molecular targets with potential therapeutic interest, a pragmatic observation on the current pharmacological R&D output confirms the lack of new therapeutic offers to patients. By contrast, from recent trials, molecules initially developed for other fields of application have shown cardiovascular benefits, as illustrated with some anti-diabetic agents, regardless of the presence or absence of diabetes, emphasizing the clear advantage of “old” drug repositioning. Ranolazine is approved as an antianginal agent and has a favorable overall safety profile. This drug, developed initially as a metabolic modulator, was also identified as an inhibitor of the cardiac late Na⁺ current, although it also blocks other ionic currents, including the hERG/Ikr K⁺ current. The latter actions have been involved in this drug’s antiarrhythmic effects, both on supraventricular and ventricular arrhythmias (VA). However, despite initial enthusiasm and promising development in the cardiovascular field, ranolazine is only authorized as a second-line treatment in patients with chronic angina pectoris, notwithstanding its antiarrhythmic properties. A plausible reason for this is the apparent difficulty in linking the clinical benefits to the multiple molecular actions of this drug. Here, we review ranolazine’s experimental and clinical knowledge on cardiac metabolism and arrhythmias. We also highlight advances in understanding novel effects on neurons, the vascular system, skeletal muscles, blood sugar control, and cancer, which may open the way to reposition this “old” drug alone or in combination with other medications.

On the nature of delays allowing anatomical re-entry involving the Purkinje network: a simulation study

AIMS

Clinical observations suggest that the Purkinje network can be part of anatomical re-entry circuits in monomorphic or polymorphic ventricular arrhythmias. However, significant conduction delay is needed to support anatomical re-entry given the high conduction velocity within the Purkinje network.

METHODS AND RESULTS

We investigated, in computer models, whether damage rendering the Purkinje network as either an active lesion with slow conduction or a passive lesion with no excitable ionic channel, could explain clinical observations. Active lesions had compromised sodium current and a severe reduction in gap junction coupling, while passive lesions remained coupled by gap junctions, but modelled the membrane as a fixed resistance. Both types of tissue could provide significant delays of over 100 ms. Electrograms consistent with those obtained clinically were reproduced. However, passive tissue could not support re-entry as electrotonic coupling across the delay effectively increased the proximal refractory period to an extremely long interval. Active tissue, conversely, could robustly maintain re-entry.

CONCLUSION

Formation of anatomical re-entry using the Purkinje network is possible through highly reduced gap junctional coupling leading to slowed conduction.

Theoretical Investigation of the Mechanism by which A Gain-of-Function Mutation of the TRPM4 Channel Causes Conduction Block

In the heart, TRPM4 is most abundantly distributed in the conduction system. Previously, a single mutation, 'E7K', was identified in its distal N-terminus to cause conduction disorder because of enhanced cell-surface expression. It remains, however, unclear how this expression increase leads to conduction failure rather than abnormally enhanced cardiac excitability. To address this issue theoretically, we mathematically formulated the gating kinetics of the E7K-mutant TRPM4 channel by a combined use of voltage jump analysis and ionomycin-perforated cell-attached recording technique and incorporated the resultant rate constants of opening and closing into a human Purkinje fiber single-cell action potential (AP) model (Trovato model) to perform 1D-cable simulations. The results from TRPM4 expressing HEK293 cells showed that as compared with the wild-type, the open state is much preferred in the E7K mutant with increased voltage-and Ca²⁺-sensitivities. These theoretical predictions were confirmed by power spectrum and single channel analyses of expressed wild-type and E7K-mutant TRPM4 channels. In our modified Trovato model, the facilitated opening of the E7K mutant channel markedly prolonged AP duration with concomitant depolarizing shifts of the resting membrane potential in a manner dependent on the channel density (or maximal activity). This was, however, little evident in the wild-type TRPM4 channel. Moreover, 1D-cable simulations with the modified Trovato model revealed that increasing the density of E7K (but not of wild-type) TRPM4 channels progressively reduced AP conduction velocity eventually culminating in complete conduction block. These results clearly suggest the brady-arrhythmogenicity of the E7K mutant channel which likely results from its pathologically enhanced activity.

Evolution of mathematical models of cardiomyocyte electrophysiology

Advanced computational techniques and mathematical modeling have become more and more important to the study of cardiac electrophysiology. In this review, we provide a brief history of the evolution of cardiomyocyte electrophysiology models and highlight some of the most important ones that had a major impact on our understanding of the electrical activity of the myocardium and associated transmembrane ion fluxes in normal and pathological states. We also present the use of these models in the study of various arrhythmogenesis mechanisms, particularly the integration of experimental pharmacology data into advanced humanized models for in silico proarrhythmogenic risk prediction as an essential component of the Comprehensive in vitro Proarrhythmia Assay (CiPA) drug safety paradigm.

Delayed afterdepolarization-induced triggered activity in cardiac purkinje cells mediated through cytosolic calcium diffusion waves

Cardiac Purkinje cells (PCs) are more susceptible to action potential abnormalities as compared to ventricular myocytes (VMs), which could be associated with their distinct intracellular calcium handling. We developed a detailed biophysical model of a mouse cardiac PC, which importantly reproduces the experimentally observed biphasic cytosolic calcium waves. The model includes a stochastic gating formulation for the opening and closing of ryanodine receptor (RyR) channels, simulated with a Monte Carlo method, to accurately reproduce cytosolic calcium wave propagation and the effects of spontaneous calcium release events. Simulations predict that during an action potential, smaller cytosolic calcium wavelets propagated from the sarcolemma towards the center of the cell and initiated larger magnitude cell‐wide calcium waves via a calcium‐induced‐calcium release mechanism. In the presence of RyR mutations, frequent spontaneous calcium leaks from sarcoplasmic reticulum (SR) initiated calcium waves, which upon reaching the cell periphery produced delayed afterdepolarizations (DADs) via sodium‐calcium exchanger (NCX) and T‐type calcium (I_CaT) channel activation. In the presence of isoproterenol‐mediated effects, DADs induced triggered activity by reactivation of fast sodium channels. Based on our model, we found that the activation of either L‐type calcium channels (I_CaL), I_CaT, sodium‐potassium exchanger (I_NaK) or NCX is sufficient for occurrence of triggered activity; however, a partial blockade of I_CaT or I_NaK is essential for its successful termination. Our modeling study highlights valuable insights into the mechanisms of DAD‐induced triggered activity mediated via cytosolic calcium waves in cardiac PCs and may elucidate the increased arrhythmogeneity in PCs.

Simulation of the Effects of Extracellular Calcium Changes Leads to a Novel Computational Model of Human Ventricular Action Potential With a Revised Calcium Handling

The importance of electrolyte concentrations for cardiac function is well established. Electrolyte variations can lead to arrhythmias onset, due to their important role in the action potential (AP) genesis and in maintaining cell homeostasis. However, most of the human AP computer models available in literature were developed with constant electrolyte concentrations, and fail to simulate physiological changes induced by electrolyte variations. This is especially true for Ca2+, even in the O'Hara-Rudy model (ORd), one of the most widely used models in cardiac electrophysiology. Therefore, the present work develops a new human ventricular model (BPS2020), based on ORd, able to simulate the inverse dependence of AP duration (APD) on extracellular Ca2+ concentration ([Ca2+]o), and APD rate dependence at 4 mM extracellular K+. The main changes needed with respect to ORd are: (i) an increased sensitivity of L-type Ca2+ current inactivation to [Ca2+]o; (ii) a single compartment description of the sarcoplasmic reticulum; iii) the replacement of Ca2+ release. BPS2020 is able to simulate the physiological APD-[Ca2+]o relationship, while also retaining the well-reproduced properties of ORd (APD rate dependence, restitution, accommodation and current block effects). We also used BPS2020 to generate an experimentally-calibrated population of models to investigate: (i) the occurrence of repolarization abnormalities in response to hERG current block; (ii) the rate adaptation variability; (iii) the occurrence of alternans and delayed after-depolarizations at fast pacing. Our results indicate that we successfully developed an improved version of ORd, which can be used to investigate electrophysiological changes and pro-arrhythmic abnormalities induced by electrolyte variations and current block at multiple rates and at the population level.

Human Purkinje in silico model enables mechanistic investigations into automaticity and pro-arrhythmic abnormalities

Cardiac Purkinje cells (PCs) are implicated in lethal arrhythmias caused by cardiac diseases, mutations, and drug action. However, the pro-arrhythmic mechanisms in PCs are not entirely understood, particularly in humans, as most investigations are conducted in animals. The aims of this study are to present a novel human PCs electrophysiology biophysically-detailed computational model, and to disentangle ionic mechanisms of human Purkinje-related electrophysiology, pacemaker activity and arrhythmogenicity. The new Trovato2020 model incorporates detailed Purkinje-specific ionic currents and Ca²⁺ handling, and was developed, calibrated and validated using human experimental data acquired at multiple frequencies, both in control conditions and following drug application. Multiscale investigations were performed in a Purkinje cell, in fibre and using an experimentally-calibrated population of PCs to evaluate biological variability. Simulations demonstrate the human Purkinje Trovato2020 model is the first one to yield: (i) all key AP features consistent with human Purkinje recordings; (ii) Automaticity with funny current up-regulation (iii) EADs at slow pacing and with 85% hERG block; (iv) DADs following fast pacing; (v) conduction velocity of 160 cm/s in a Purkinje fibre, as reported in human. The human in silico PCs population highlights that: (1) EADs are caused by I_CaL reactivation in PCs with large inward currents; (2) DADs and triggered APs occur in PCs experiencing Ca²⁺ accumulation, at fast pacing, caused by large L-type calcium current and small Na⁺/Ca²⁺ exchanger. The novel human Purkinje model unlocks further investigations into the role of cardiac Purkinje in ventricular arrhythmias through computer modeling and multiscale simulations.

•

A human in silico AP model was developed to investigate arrhythmia in cardiac Purkinje.

•

The new Purkinje model enables multiscale investigations from single cell to tissue.

•

Populations of human Purkinje models reproduce and explain experimental variability.

•

Ca²⁺-current reactivation triggers EADs in virtual Purkinje cells with weak repolarisation reserve.

•

Ca²⁺ accumulation caused by increased Ca²⁺ and NCX currents triggers DADs.

Transmural and rate-dependent profiling of drug-induced arrhythmogenic risks through in silico simulations of multichannel pharmacology

In vitro human ether-à-go-go related gene (hERG) inhibition assay alone might provide insufficient information to discriminate "safe" from "dangerous" drugs. Here, effects of multichannel inhibition on cardiac electrophysiology were investigated using a family of cardiac cell models (Purkinje (P), endocardial (Endo), mid-myocardial (M) and epicardial (Epi)). We found that: (1) QT prolongation alone might not necessarily lead to early afterdepolarization (EAD) events, and it might be insufficient to predict arrhythmogenic liability; (2) the occurrence and onset of EAD events could be a candidate biomarker of drug-induced arrhythmogenicity; (3) M cells are more vulnerable to drug-induced arrhythmias, and can develop early afterdepolarization (EAD) at slower pacing rates; (4) the application of quinidine can cause EADs in all cell types, while I_NaL is the major depolarizing current during the generation of drug-induced EAD in P cells, I_CaL is mostly responsible in other cell types; (5) drug-induced action potential (AP) alternans with beat-to-beat variations occur at high pacing rates in P cells. These results suggested that quantitative profiling of transmural and rate-dependent properties can be essential to evaluate drug-induced arrhythmogenic risks, and may provide mechanistic insights into drug-induced arrhythmias.

Inward Rectifier Potassium Channels (Kir2.x) and Caveolin-3 Domain-Specific Interaction: Implications for Purkinje Cell-Dependent Ventricular Arrhythmias

BACKGROUND

In human cardiac ventricle, I_K1 is mainly comprised Kir2.1, but Kir2.2 and Kir2.3 heterotetramers occur and modulate I_K1. Long-QT syndrome-9-associated CAV3 mutations cause decreased Kir2.1 current density, but Kir2.x heterotetramers have not been studied. Here, we determine the effect of long-QT syndrome-9-CAV3 mutation F97C on Kir2.x homo- and heterotetramers and model-associated arrhythmia mechanisms.

METHODS AND RESULTS

Super-resolution microscopy, co-immunoprecipitation, cellular electrophysiology, on-cell Western blotting, and simulation of Purkinje and ventricular myocyte mathematical models were used. Kir2.x isoforms have unique subcellular colocalization in human cardiomyocytes and coimmunoprecipitate with Cav3. F97C-Cav3 decreased peak inward Kir2.2 current density by 50% (-120 mV; P=0.019) and peak outward by 75% (-40 mV; P<0.05) but did not affect Kir2.3 current density. FRET (Förster resonance energy transfer) efficiency for Kir2.2 with Cav3 is high, and on-cell Western blotting demonstrates decreased Kir2.2 membrane expression with F97C-Cav3. Cav3-F97C reduced peak inward and outward current density of Kir2.2/Kir2.1 or Kir2.2/Kir2.3 heterotetramers (P<0.05). Only Cav3 scaffolding and membrane domains co-immunoprecipitation with Kir2.1 and Kir2.2 and Kir2.x-N-terminal Cav3 binding motifs are required for interaction. Mathematical Purkinje, but not ventricular, myocyte model incorporating simulated current reductions, predicts spontaneous delayed after-depolarization-mediated triggered activity.

CONCLUSIONS

Kir2.x isoforms have a unique intracellular pattern of distribution in association with specific Cav3 domains and that critically depends on interaction with N-terminal Kir2.x Cav3-binding motifs. Long-QT syndrome-9-CAV3 mutation differentially regulates current density and cell surface expression of Kir2.x homomeric and heteromeric channels. Mathematical Purkinje cell model incorporating experimental findings suggests delayed after-depolarization-type triggered activity as a possible arrhythmia mechanism.

Propagation Failure by TRPM4 Overexpression

Transient receptor potential melastatin member 4 (TRPM4) channels are nonselective monovalent cationic channels found in human atria and conduction system. Overexpression of TRPM4 channels has been found in families suffering from inherited cardiac arrhythmias, notably heart block. In this study, we integrate a mathematical formulation of the TRPM4 channel into a Purkinje cell model (Pan-Rudy model). Instead of simply adding the channel to the model, a combination of existing currents equivalent to the TRPM4 current was constructed, based on TRPM4 current dynamics. The equivalent current was then replaced by the TRPM4 current to preserve the model action potential. Single-cell behavior showed early afterdepolarizations for increases in TRPM4 channel expression above twofold. In a homogeneous strand of tissue, propagation conducted faithfully for lower expression levels but failed completely for more than a doubling of TRPM4 channel expression. Only with a heterogeneous distribution of channel expression was intermittent heart block seen. This study suggests that in Purkinje fibers, TRPM4 channels may account for sodium background current (I_Nab), and that a heterogeneous expression of TRPM4 channels in the His/Purkinje system is required for type II heart block, as seen clinically.

Caveolin-3 Microdomain: Arrhythmia Implications for Potassium Inward Rectifier and Cardiac Sodium Channel

In human cardiac ventricular myocytes, caveolin-3 functions as a scaffolding and regulatory protein for signaling molecules and compartmentalizes ion channels. Our lab has recently explored this sub-cellular microdomain and found that potassium inward rectifier Kir2.x is found in association with caveolin-3. The three cardiac Kir2.x isoforms (Kir2.1, Kir2.2, and Kir2.3) are the molecular correlates of I_K1 in the heart, of which Kir2.1 is the dominant isoform in the ventricle. Kir2.1 channels assemble with Kir2.2 and Kir2.3 forming hetero-tetramers that modulate I_K1. I_K1 sets the resting membrane potential and assists with terminal phase 3 ventricular repolarization. In our studies using native human ventricular tissue, Kir2.x co-localizes with caveolin-3 and significance of the association between Kir2.x and caveolin-3 is emphasized in relation to mutations in the gene which encodes caveolin-3, CAV3, associated with Long QT Syndrome 9 (LQT9). LQT9-associated CAV3 mutations cause decreased current density in Kir2.1 and Kir2.2 as homomeric and heteromeric channels, which affects repolarization and membrane potential stability. A portion of Kir2.1 cardiac localization parallels that of the cardiac sodium channel (Nav1.5). This may have implications for Long QT9 in which CAV3 mutations cause an increase in the late current of Nav1.5 (I_Na-L) via nNOS mediated nitrosylation of Nav1.5. In iPS-CMs, expression of LQT9 CAV3 mutations resulted in action potential duration (APD) prolongation and early-after depolarizations (EADs), supporting the arrhythmogenicity of LQT9. To evaluate the combined effect of the CAV3 mutants on I_Na-L and I_K1, we studied both ventricular and Purkinje myocyte mathematical modeling. Interestingly, mathematical ventricular myocytes, similar to iPS-CMs, demonstrated EADs but no sustained arrhythmia. In contrast, Purkinje modeling demonstrated delayed-after depolarizations (DADs) driven mechanism for sustained arrhythmia, dependent on the combined loss of I_K1 and gain of I_Na-L. This finding changes the overall assumed arrhythmia phenotype for LQT9. In future studies, we are exploring caveolar micro-domain disruption in heart failure and how this effects Kir2.x and Nav1.5. Here we review the caveolae cardiac microdomain of Kir2.x and Nav1.5 and explore some of the downstream effects of caveolin-3 and caveolae disruption in specific clinical scenarios.

Alternans promotion in cardiac electrophysiology models by delay differential equations

Cardiac electrical alternans is a state of alternation between long and short action potentials and is frequently associated with harmful cardiac conditions. Different dynamic mechanisms can give rise to alternans; however, many cardiac models based on ordinary differential equations are not able to reproduce this phenomenon. A previous study showed that alternans can be induced by the introduction of delay differential equations (DDEs) in the formulations of the ion channel gating variables of a canine myocyte model. The present work demonstrates that this technique is not model-specific by successfully promoting alternans using DDEs for five cardiac electrophysiology models that describe different types of myocytes, with varying degrees of complexity. By analyzing results across the different models, we observe two potential requirements for alternans promotion via DDEs for ionic gates: (i) the gate must have a significant influence on the action potential duration and (ii) a delay must significantly impair the gate's recovery between consecutive action potentials.

Stochastic spontaneous calcium release events and sodium channelopathies promote ventricular arrhythmias

Premature ventricular complexes (PVCs), the first initiating beats of a variety of cardiac arrhythmias, have been associated with spontaneous calcium release (SCR) events at the cell level. However, the mechanisms underlying the degeneration of such PVCs into arrhythmias are not fully understood. The objective of this study was to investigate the conditions under which SCR-mediated PVCs can lead to ventricular arrhythmias. In particular, we sought to determine whether sodium (Na⁺) current loss-of-function in the structurally normal ventricles provides a substrate for unidirectional conduction block and reentry initiated by SCR-mediated PVCs. To achieve this goal, a stochastic model of SCR was incorporated into an anatomically accurate compute model of the rabbit ventricles with the His-Purkinje system (HPS). Simulations with reduced Na⁺ current due to a negative-shift in the steady-state channel inactivation showed that SCR-mediated delayed afterdepolarizations led to PVC formation in the HPS, where the electrotonic load was lower, conduction block, and reentry in the 3D myocardium. Moreover, arrhythmia initiation was only possible when intrinsic electrophysiological heterogeneity in action potential within the ventricles was present. In conclusion, while benign in healthy individuals SCR-mediated PVCs can lead to life-threatening ventricular arrhythmias when combined with Na⁺ channelopathies.

Premature Ventricular Complexes in Apparently Normal Hearts

Premature ventricular complexes (PVCs) are consistently associated with worse prognosis and higher morbidity and mortality. This article reviews PVCs and their presentation in patients with an apparently normal heart. Patients with PVCs may be completely asymptomatic, whereas others may note severely disabling symptoms. Cardiomyopathy may occur with frequent PVCs. Diagnostic work-up is directed at obtaining 12-lead ECG to characterize QRS morphology, Holter monitor to assess frequency, and echo and advanced imaging to assess for early cardiomyopathy and exclude structural heart disease. Options for management include watchful waiting, medical therapy, or catheter ablation. Malignant variants of PVCs may induce ventricular fibrillation even in a normal heart.

Role of Cytosolic Calcium Diffusion in Murine Cardiac Purkinje Cells

Cardiac Purkinje cells (PCs) are morphologically and electrophysiologically different from ventricular myocytes and, importantly, exhibit distinct calcium (Ca²⁺) homeostasis. Recent studies suggest that PCs are more susceptible to action potential (AP) abnormalities than ventricular myocytes; however, the exact mechanisms are poorly understood. In this study, we utilized a detailed biophysical mathematical model of a murine PC to systematically examine the role of cytosolic Ca²⁺ diffusion in shaping the AP in PCs. A biphasic spatiotemporal Ca²⁺ diffusion process, as recorded experimentally, was implemented in the model. In this study, we investigated the role of cytosolic Ca²⁺ dynamics on AP and ionic current properties by varying the effective Ca²⁺ diffusion rate. It was observed that AP morphology, specifically the plateau, was affected due to changes in the intracellular Ca²⁺ dynamics. Elevated Ca²⁺ concentration in the sarcolemmal region activated inward sodium–Ca²⁺ exchanger (NCX) current, resulting in a prolongation of the AP plateau at faster diffusion rates. Artificially clamping the NCX current to control values completely reversed the alterations in the AP plateau, thus confirming the role of NCX in modifying the AP morphology. Our results demonstrate that cytosolic Ca²⁺ diffusion waves play a significant role in shaping APs of PCs and could provide mechanistic insights in the increased arrhythmogeneity of PCs.

Ranolazine is an Effective and Safe Treatment of Adults with Symptomatic Premature Ventricular Contractions due to Triggered Ectopy

Early and delayed afterdepolarizations (EAD/DAD) cause triggered ventricular ectopy. Because ranolazine (RAN) suppresses EAD/DAD, we postulated that RAN might be effective in reducing premature ventricular contractions (PVCs). To assess the effect of RAN in patients with symptomatic PVCs due to triggered ectopy and its safety and tolerability. A total of 59 patients with symptomatic PVCs were identified from full-disclosure Holters. Doses of 500 and 1,000 mg offlabel RAN, daily, were given to 34 and 66% patients, respectively, and repeat Holters were performed prospectively during mean followup of 3.1 months. The two Holters were retrospectively compared. Congestive heart failure (CHF) was defined as symptoms including dyspnea, orthopnea, paroxysmal nocturnal dyspnea, and fatigue, with a brain natriuretic peptide > 400. Systolic (heart failure with reduced ejection fraction) versus diastolic (heart failure with preserved ejection fraction, HFpEF) CHF depended upon an echocardiographic left ventricular ejection fraction (LVEF) at least 50% by apical two- and four-chamber Simpson's method (HFpEF). The mean age of the patients was 63 years, 60% were males, mean left ventricular ejection fraction was 60%, with 34% having coronary artery disease, 73% were hypertensive, 24% had type 2 diabetic, and 34% were on beta blockers. Upon repeat Holters at a mean of 3.1 months after initiating RAN, 95% (56/59) of the patients had their PVC count reduced as follows: 24% (14/59) had more than 90% decrease, 34% (20/59) had 71 to 90% decrease, and 17% (10/59) had 50 to 70% decrease. In the entire group, RAN reduced PVCs by 71% (mean: 13,329 to 3,837; p < 0.001). Ventricular bigeminy was reduced by 80% (4,168 to 851; p < 0.001), ventricular coupletswere reduced by 78% (374 to 81; p < 0.001), and ventricular tachycardiawas reduced by 91% (56 to 5; p < 0.001). The PVC reduction was dose dependent. Off-label RAN offers an effective and safe pharmacologic treatment for symptomatic triggered PVCs. A large, prospective randomized study is needed.

Constitutive Intracellular Na+ Excess in Purkinje Cells Promotes Arrhythmogenesis at Lower Levels of Stress Than Ventricular Myocytes From Mice With Catecholaminergic Polymorphic Ventricular Tachycardia

Supplemental Digital Content is available in the text.

Background—

In catecholaminergic polymorphic ventricular tachycardia (CPVT), cardiac Purkinje cells (PCs) appear more susceptible to Ca²⁺ dysfunction than ventricular myocytes (VMs). The underlying mechanisms remain unknown. Using a CPVT mouse (RyR2^{R4496C+/Cx40eGFP}), we tested whether PC intracellular Ca²⁺ ([Ca²⁺]_i) dysregulation results from a constitutive [Na⁺]_i surplus relative to VMs.

Methods and Results—

Simultaneous optical mapping of voltage and [Ca²⁺]_i in CPVT hearts showed that spontaneous Ca²⁺ release preceded pacing-induced triggered activity at subendocardial PCs. On simultaneous current-clamp and Ca²⁺ imaging, early and delayed afterdepolarizations trailed spontaneous Ca²⁺ release and were more frequent in CPVT PCs than CPVT VMs. As a result of increased activity of mutant ryanodine receptor type 2 channels, sarcoplasmic reticulum Ca²⁺ load, measured by caffeine-induced Ca²⁺ transients, was lower in CPVT VMs and PCs than respective controls, and sarcoplasmic reticulum fractional release was greater in both CPVT PCs and VMs than respective controls. [Na⁺]_i was higher in both control and CPVT PCs than VMs, whereas the density of the Na⁺/Ca²⁺ exchanger current was not different between PCs and VMs. Computer simulations using a PC model predicted that the elevated [Na⁺]_i of PCs promoted delayed afterdepolarizations, which were always preceded by spontaneous Ca²⁺ release events from hyperactive ryanodine receptor type 2 channels. Increasing [Na⁺]_i monotonically increased delayed afterdepolarization frequency. Confocal imaging experiments showed that postpacing Ca²⁺ spark frequency was highest in intact CPVT PCs, but such differences were reversed on saponin-induced membrane permeabilization, indicating that differences in [Na⁺]_i played a central role.

Conclusions—

In CPVT mice, the constitutive [Na⁺]_i excess of PCs promotes triggered activity and arrhythmogenesis at lower levels of stress than VMs.

Modeling our understanding of the His-Purkinje system

The His-Purkinje System (HPS) is responsible for the rapid electric conduction in the ventricles. It relays electrical impulses from the atrioventricular node to the muscle cells and, thus, coordinates the contraction of ventricles in order to ensure proper cardiac pump function. The HPS has been implicated in the genesis of ventricular tachycardia and fibrillation as a source of ectopic beats, as well as forming distinct portions of reentry circuitry. Despite its importance, it remains much less well characterized, structurally and functionally, than the myocardium. Notably, important differences exist with regard to cell structure and electrophysiology, including ion channels, intracellular calcium handling, and gap junctions. Very few computational models address the HPS, and the majority of organ level modeling studies omit it. This review will provide an overview of our current knowledge of structure and function (including electrophysiology) of the HPS. We will review the most recent advances in modeling of the system from the single cell to the organ level, with considerations for relevant interspecies distinctions.

Myokit: A simple interface to cardiac cellular electrophysiology

Myokit is a new powerful and versatile software tool for modeling and simulation of cardiac cellular electrophysiology. Myokit consists of an easy-to-read modeling language, a graphical user interface, single and multi-cell simulation engines and a library of advanced analysis tools accessible through a Python interface. Models can be loaded from Myokit's native file format or imported from CellML. Model export is provided to C, MATLAB, CellML, CUDA and OpenCL. Patch-clamp data can be imported and used to estimate model parameters. In this paper, we review existing tools to simulate the cardiac cellular action potential to find that current tools do not cater specifically to model development and that there is a gap between easy-to-use but limited software and powerful tools that require strong programming skills from their users. We then describe Myokit's capabilities, focusing on its model description language, simulation engines and import/export facilities in detail. Using three examples, we show how Myokit can be used for clinically relevant investigations, multi-model testing and parameter estimation in Markov models, all with minimal programming effort from the user. This way, Myokit bridges a gap between performance, versatility and user-friendliness.

Optogenetic determination of the myocardial requirements for extrasystoles by cell type-specific targeting of ChannelRhodopsin-2

Extrasystoles lead to several consequences, ranging from uneventful palpitations to lethal ventricular arrhythmias, in the presence of pathologies, such as myocardial ischemia. The role of working versus conducting cardiomyocytes, as well as the tissue requirements (minimal cell number) for the generation of extrasystoles, and the properties leading ectopies to become arrhythmia triggers (topology), in the normal and diseased heart, have not been determined directly in vivo. Here, we used optogenetics in transgenic mice expressing ChannelRhodopsin-2 selectively in either cardiomyocytes or the conduction system to achieve cell type-specific, noninvasive control of heart activity with high spatial and temporal resolution. By combining measurement of optogenetic tissue activation in vivo and epicardial voltage mapping in Langendorff-perfused hearts, we demonstrated that focal ectopies require, in the normal mouse heart, the simultaneous depolarization of at least 1,300-1,800 working cardiomyocytes or 90-160 Purkinje fibers. The optogenetic assay identified specific areas in the heart that were highly susceptible to forming extrasystolic foci, and such properties were correlated to the local organization of the Purkinje fiber network, which was imaged in three dimensions using optical projection tomography. Interestingly, during the acute phase of myocardial ischemia, focal ectopies arising from this location, and including both Purkinje fibers and the surrounding working cardiomyocytes, have the highest propensity to trigger sustained arrhythmias. In conclusion, we used cell-specific optogenetics to determine with high spatial resolution and cell type specificity the requirements for the generation of extrasystoles and the factors causing ectopies to be arrhythmia triggers during myocardial ischemia.

Selective Na(+) /Ca(2+) exchanger inhibition prevents Ca(2+) overload-induced triggered arrhythmias

BACKGROUND AND PURPOSE

Augmented Na(+) /Ca(2+) exchanger (NCX) activity may play a crucial role in cardiac arrhythmogenesis; however, data regarding the anti-arrhythmic efficacy of NCX inhibition are debatable. Feasible explanations could be the unsatisfactory selectivity of NCX inhibitors and/or the dependence of the experimental model on the degree of Ca(2+) i overload. Hence, we used NCX inhibitors SEA0400 and the more selective ORM10103 to evaluate the efficacy of NCX inhibition against arrhythmogenic Ca(2+) i rise in conditions when [Ca(2+) ]i was augmented via activation of the late sodium current (INaL ) or inhibition of the Na(+) /K(+) pump.

EXPERIMENTAL APPROACH

Action potentials (APs) were recorded from canine papillary muscles and Purkinje fibres by microelectrodes. NCX current (INCX ) was determined in ventricular cardiomyocytes utilizing the whole-cell patch clamp technique. Ca(2+) i transients (CaTs) were monitored with a Ca(2+) -sensitive fluorescent dye, Fluo-4.

KEY RESULTS

Enhanced INaL increased the Ca(2+) load and AP duration (APD). SEA0400 and ORM10103 suppressed INCX and prevented/reversed the anemone toxin II (ATX-II)-induced [Ca(2+) ]i rise without influencing APD, CaT or cell shortening, or affecting the ATX-II-induced increased APD. ORM10103 significantly decreased the number of strophanthidin-induced spontaneous diastolic Ca(2+) release events; however, SEA0400 failed to restrict the veratridine-induced augmentation in Purkinje-ventricle APD dispersion.

CONCLUSIONS AND IMPLICATIONS

Selective NCX inhibition - presumably by blocking rev INCX (reverse mode NCX current) - is effective against arrhythmogenesis caused by [Na(+) ]i -induced [Ca(2+) ]i elevation, without influencing the AP waveform. Therefore, selective INCX inhibition, by significantly reducing the arrhythmogenic trigger activity caused by the perturbed Ca(2+) i handling, should be considered as a promising anti-arrhythmic therapeutic strategy.

Stochastic spontaneous calcium release events trigger premature ventricular complexes by overcoming electrotonic load

AIMS

Premature ventricular complexes (PVCs) due to spontaneous calcium (Ca) release (SCR) events at the cell level can precipitate ventricular arrhythmias. However, the mechanistic link between SCRs and PVC formation remains incompletely understood. The aim of this study was to investigate the conditions under which delayed afterdepolarizations resulting from stochastic subcellular SCR events can overcome electrotonic source-sink mismatch, leading to PVC initiation.

METHODS AND RESULTS

A stochastic subcellular-scale mathematical model of SCR was incorporated in a realistic model of the rabbit ventricles and Purkinje system (PS). Elevated levels of diastolic sarcoplasmic reticulum Ca(2+) (CaSR) were imposed until triggered activity was observed, allowing us to compile statistics on probability, timing, and location of PVCs. At CaSR≥ 1500 µmol/L PVCs originated in the PS. When SCR was incapacitated in the PS, PVCs also emerged in the ventricles, but at a higher CaSR (≥1550 µmol/L) and with longer waiting times. For each model configuration tested, the probability of PVC occurrence increased from 0 to 100% within a well-defined critical CaSR range; this transition was much more abrupt in organ-scale models (∼50 µmol/L CaSR range) than in the tissue strand (∼100 µmol/L) or single-cell (∼450 µmol/L) models. Among PVCs originating in the PS, ∼68% were located near Purkinje-ventricular junctions (<1 mm).

CONCLUSION

SCR events overcome source-sink mismatch to trigger PVCs at a critical CaSR threshold. Above this threshold, PVCs emerge due to increased probability and reduced variability in timing of SCR events, leading to significant diastolic depolarization. Sites of lower electronic load, such as the PS, are preferential locations for triggering.

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

The combination of computational models and biophysical simulations can help to interpret an array of experimental data and contribute to the understanding, diagnosis and treatment of complex diseases such as cardiac arrhythmias. For this reason, three-dimensional (3D) cardiac computational modelling is currently a rising field of research. The advance of medical imaging technology over the last decades has allowed the evolution from generic to patient-specific 3D cardiac models that faithfully represent the anatomy and different cardiac features of a given alive subject. Here we analyse sixty representative 3D cardiac computational models developed and published during the last fifty years, describing their information sources, features, development methods and online availability. This paper also reviews the necessary components to build a 3D computational model of the heart aimed at biophysical simulation, paying especial attention to cardiac electrophysiology (EP), and the existing approaches to incorporate those components. We assess the challenges associated to the different steps of the building process, from the processing of raw clinical or biological data to the final application, including image segmentation, inclusion of substructures and meshing among others. We briefly outline the personalisation approaches that are currently available in 3D cardiac computational modelling. Finally, we present examples of several specific applications, mainly related to cardiac EP simulation and model-based image analysis, showing the potential usefulness of 3D cardiac computational modelling into clinical environments as a tool to aid in the prevention, diagnosis and treatment of cardiac diseases.

In Silico Organ Modelling in Predicting Efficacy and Safety of New Medicines

The development of new medicines faces important challenges due to difficulties in the assessment of their efficacy and their safety in the targeted human population. In silico approaches through the use of mathematical modelling and computer simulations are increasingly being used to overcome some of the limitations of current experimental methods used in the development of new medicines. This chapter describes state-of-the-art in silico approaches for the evaluation of the safety and efficacy of medicines targeting important causes of mortality such as cardiovascular disease. Firstly, we describe the in silico multi-scale mathematical models and simulation techniques required to describe drug-induced effects on physiological systems such as the heart from the subcellular to the whole organ level. Then we illustrate the power of in silico approaches used to augment experimental and clinical investigations, by providing the framework to unravel multi-scale mechanisms underlying variability in the response to medicines and to focus on effects in human rather than animal models. We devote the last part of the chapter to discussing the process of validation of in silico models and simulations, which is key in building up their credibility.

Exploring susceptibility to atrial and ventricular arrhythmias resulting from remodeling of the passive electrical properties in the heart: a simulation approach

Under diseased conditions, remodeling of the cardiac tissue properties (“passive properties”) takes place; these are aspects of electrophysiological behavior that are not associated with active ion transport across cell membranes. Remodeling of the passive electrophysiological properties most often results from structural remodeling, such as gap junction down-regulation and lateralization, fibrotic growth infiltrating the myocardium, or the development of an infarct scar. Such structural remodeling renders atrial or ventricular tissue as a major substrate for arrhythmias. The current review focuses on these aspects of cardiac arrhythmogenesis. Due to the inherent complexity of cardiac arrhythmias, computer simulations have provided means to elucidate interactions pertinent to this spatial scale. Here we review the current state-of-the-art in modeling atrial and ventricular arrhythmogenesis as arising from the disease-induced changes in the passive tissue properties, as well as the contributions these modeling studies have made to our understanding of the mechanisms of arrhythmias in the heart. Because of the rapid advance of structural imaging methodologies in cardiac electrophysiology, we chose to present studies that have used such imaging methodologies to construct geometrically realistic models of cardiac tissue, or the organ itself, where the regional remodeling properties of the myocardium can be represented in a realistic way. We emphasize how the acquired knowledge can be used to pave the way for clinical applications of cardiac organ modeling under the conditions of structural remodeling.

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

AIMS

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

METHODS AND RESULTS

Unidimensional and bidimensional models of the Purkinje-ventricle subsystem, including ischaemic conditions (defined as phase 1B) in the ventricle and an ischaemic border zone, were developed by altering several important electrophysiological parameters of the Luo-Rudy AP model of the ventricular myocyte. Purkinje electrical activity was modelled using the equations of DiFrancesco and Noble. Our study suggests that an extracellular potassium concentration [K(+)]o >14 mM and a slight decrease in intercellular coupling induced by ischaemia in ventricle can cause conduction block from Purkinje to ventricle. Under these conditions, propagation from ventricle to Purkinje is possible. Thus, unidirectional block (UDB) and reentry can result. When conditions of UDB are met, retrograde propagation with a long delay (320 ms) may re-excite Purkinje cells, and give rise to a reentrant pathway. This induced reentry may be the origin of arrhythmias observed in phase 1B ischaemia.

CONCLUSION

In a defined setting of ischaemia (phase 1B), a small amount of uncoupling between ventricular cells, as well as between Purkinje and ventricular tissue, may induce UDBs and reentry. Hyperkalaemia is also confirmed to be an important factor in the genesis of reentrant rhythms, since it regulates the range of coupling in which UDBs may be induced.

Computational modeling of inhibition of voltage-gated Ca channels: identification of different effects on uterine and cardiac action potentials

The uterus and heart share the important physiological feature whereby contractile activation of the muscle tissue is regulated by the generation of periodic, spontaneous electrical action potentials (APs). Preterm birth arising from premature uterine contractions is a major complication of pregnancy and there remains a need to pursue avenues of research that facilitate the use of drugs, tocolytics, to limit these inappropriate contractions without deleterious actions on cardiac electrical excitation. A novel approach is to make use of mathematical models of uterine and cardiac APs, which incorporate many ionic currents contributing to the AP forms, and test the cell-specific responses to interventions. We have used three such models-of uterine smooth muscle cells (USMC), cardiac sinoatrial node cells (SAN), and ventricular cells-to investigate the relative effects of reducing two important voltage-gated Ca currents-the L-type (ICaL) and T-type (ICaT) Ca currents. Reduction of ICaL (10%) alone, or ICaT (40%) alone, blunted USMC APs with little effect on ventricular APs and only mild effects on SAN activity. Larger reductions in either current further attenuated the USMC APs but with also greater effects on SAN APs. Encouragingly, a combination of ICaL and ICaT reduction did blunt USMC APs as intended with little detriment to APs of either cardiac cell type. Subsequent overlapping maps of ICaL and ICaT inhibition profiles from each model revealed a range of combined reductions of ICaL and ICaT over which an appreciable diminution of USMC APs could be achieved with no deleterious action on cardiac SAN or ventricular APs. This novel approach illustrates the potential for computational biology to inform us of possible uterine and cardiac cell-specific mechanisms. Incorporating such computational approaches in future studies directed at designing new, or repurposing existing, tocolytics will be beneficial for establishing a desired uterine specificity of action.

PCP4 regulates Purkinje cell excitability and cardiac rhythmicity

Cardiac Purkinje cells are important triggers of ventricular arrhythmias associated with heritable and acquired syndromes; however, the mechanisms responsible for this proarrhythmic behavior are incompletely understood. Here, through transcriptional profiling of genetically labeled cardiomyocytes, we identified expression of Purkinje cell protein-4 (Pcp4), a putative regulator of calmodulin and Ca2+/calmodulin-dependent kinase II (CaMKII) signaling, exclusively within the His-Purkinje network. Using Pcp4-null mice and acquired cardiomyopathy models, we determined that reduced expression of PCP4 is associated with CaMKII activation, abnormal electrophysiology, dysregulated intracellular calcium handling, and proarrhythmic behavior in isolated Purkinje cells. Pcp4-null mice also displayed profound autonomic dysregulation and arrhythmic behavior in vivo. Together, these results demonstrate that PCP4 regulates cardiac excitability through both Purkinje cell-autonomous and central mechanisms and identify this modulator of CaMKII signaling as a potential arrhythmia-susceptibility candidate.

Idiopathic ventricular arrhythmias originating from the moderator band: Electrocardiographic characteristics and treatment by catheter ablation

BACKGROUND

The moderator band (MB) can be a source of premature ventricular contractions (PVCs), monomorphic ventricular tachycardia (VT), and idiopathic ventricular fibrillation (IVF).

OBJECTIVE

The purpose of this study was to define the electrocardiographic (ECG) characteristics and procedural techniques to successfully identify and ablate MB PVCs/VT.

METHODS

In 10 patients with left bundle branch block morphology PVCs/VT, electroanatomic mapping in conjunction with intracardiac echocardiography (ICE) localized the site of origin of the PVCs to the MB. Clinical characteristics of the patients, ECG features, and procedural data were collected and analyzed.

RESULTS

Seven patients presented with IVF and 3 presented with monomorphic VT. In all patients, the ventricular arrhythmias (VAs) had a left bundle branch block QRS with a late precordial transition (>V4), a rapid downstroke of the QRS in the precordial leads, and a left superior frontal plane axis. Mean QRS duration was 152.7 ± 15.2 ms. Six patients required a repeat procedure. After mean follow-up of 21.5 ± 11.6 months, all patients were free of sustained VAs, with only 1 patient requiring antiarrhythmic drug therapy and 1 patient having isolated PVCs no longer inducing VF. There were no procedural complications.

CONCLUSION

VAs originating from the MB have a distinctive morphology and often are associated with PVC-induced ventricular fibrillation. Catheter ablation can be safely performed and is facilitated by ICE imaging.

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.

Effects of maximal sodium and potassium conductance on the stability of Hodgkin-Huxley model

Hodgkin-Huxley (HH) equation is the first cell computing model in the world and pioneered the use of model to study electrophysiological problems. The model consists of four differential equations which are based on the experimental data of ion channels. Maximal conductance is an important characteristic of different channels. In this study, mathematical method is used to investigate the importance of maximal sodium conductance g-Na and maximal potassium conductance g-K. Applying stability theory, and taking g-Na and g-K as variables, we analyze the stability and bifurcations of the model. Bifurcations are found when the variables change, and bifurcation points and boundary are also calculated. There is only one bifurcation point when g-Na is the variable, while there are two points when g-K is variable. The (g-Na,  g-K) plane is partitioned into two regions and the upper bifurcation boundary is similar to a line when both g-Na and g-K are variables. Numerical simulations illustrate the validity of the analysis. The results obtained could be helpful in studying relevant diseases caused by maximal conductance anomaly.

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.

Evoked centripetal Ca(2+) mobilization in cardiac Purkinje cells: insight from a model of three Ca(2+) release regions

Despite strong suspicion that abnormal Ca(2+) handling in Purkinje cells (P-cells) is implicated in life-threatening forms of ventricular tachycardias, the mechanism underlying the Ca(2+) cycling of these cells under normal conditions is still unclear. There is mounting evidence that P-cells have a unique Ca(2+) handling system. Notably complex spontaneous Ca(2+) activity was previously recorded in canine P-cells and was explained by a mechanistic hypothesis involving a triple layered system of Ca(2+) release channels. Here we examined the validity of this hypothesis for the electrically evoked Ca(2+) transient which was shown, in the dog and rabbit, to occur progressively from the periphery to the interior of the cell. To do so, the hypothesis was incorporated in a model of intracellular Ca(2+) dynamics which was then used to reproduce numerically the Ca(2+) activity of P-cells under stimulated conditions. The modelling was thus performed through a 2D computational array that encompassed three distinct Ca(2+) release nodes arranged, respectively, into three consecutive adjacent regions. A system of partial differential equations (PDEs) expressed numerically the principal cellular functions that modulate the local cytosolic Ca(2+) concentration (Cai). The apparent node-to-node progression of elevated Cai was obtained by combining Ca(2+) diffusion and 'Ca(2+)-induced Ca(2+) release'. To provide the modelling with a reliable experimental reference, we first re-examined the Ca(2+) mobilization in swine stimulated P-cells by 2D confocal microscopy. As reported earlier for the dog and rabbit, a centripetal Ca(2+) transient was readily visible in 22 stimulated P-cells from six adult Yucatan swine hearts (pacing rate: 0.1 Hz; pulse duration: 25 ms, pulse amplitude: 10% above threshold; 1 mm Ca(2+); 35°C; pH 7.3). An accurate replication of the observed centripetal Ca(2+) propagation was generated by the model for four representative cell examples and confirmed by statistical comparisons of simulations against cell data. Selective inactivation of Ca(2+) release regions of the computational array showed that an intermediate layer of Ca(2+) release nodes with an ~30-40% lower Ca(2+) activation threshold was required to reproduce the phenomenon. Our computational analysis was therefore fully consistent with the activation of a triple layered system of Ca(2+) release channels as a mechanism of centripetal Ca(2+) signalling in P-cells. Moreover, the model clearly indicated that the intermediate Ca(2+) release layer with increased sensitivity for Ca(2+) plays an important role in the specific intracellular Ca(2+) mobilization of Purkinje fibres and could therefore be a relevant determinant of cardiac conduction.

Three-dimensional impulse propagation in myocardium: arrhythmogenic mechanisms at the tissue level

Impulse propagation in the heart depends on the excitability of individual cardiomyocytes, impulse transmission between adjacent myocytes, and the 3-dimensional arrangement of those cells. Here, we review the role of each of these factors in normal and aberrant cardiac electric activation, with particular emphasis on the effects of 3-dimensional myocyte architecture at the tissue scale. The analysis draws on findings from in vivo and in vitro experiments, as well as biophysically based computer models that have been used to integrate and interpret these experimental data. It indicates that discontinuous arrangement of myocytes and extracellular connective tissue at the tissue scale can give rise to current source-to-sink mismatch, spatiotemporal distribution of refractoriness, and rate-sensitive electric instability, which contribute to the initiation and maintenance of reentrant cardiac arrhythmia. This exacerbates the risk of rhythm disturbance associated with heart disease. We conclude that structure-based, multiscale computer models that incorporate accurate information about local cellular electric activity provide a powerful platform for investigating the basis of reentrant cardiac arrhythmia. However, it is important that these models capture key features of structure and related electric function at the tissue scale.

The arrhythmogenic consequences of increasing late INa in the cardiomyocyte

This review presents the roles of cardiac sodium channel NaV1.5 late current (late INa) in generation of arrhythmic activity. The assumption of the authors is that proper Na(+) channel function is necessary to the maintenance of the transmembrane electrochemical gradient of Na(+) and regulation of cardiac electrical activity. Myocyte Na(+) channels' openings during the brief action potential upstroke contribute to peak INa and initiate excitation-contraction coupling. Openings of Na(+) channels outside the upstroke contribute to late INa, a depolarizing current that persists throughout the action potential plateau. The small, physiological late INa does not appear to be critical for normal electrical or contractile function in the heart. Late INa does, however, reduce the net repolarizing current, prolongs action potential duration, and increases cellular Na(+) loading. An increase of late INa, due to acquired conditions (e.g. heart failure) or inherited Na(+) channelopathies, facilitates the formation of early and delayed afterpolarizations and triggered arrhythmias, spontaneous diastolic depolarization, and cellular Ca(2+) loading. These in turn increase the spatial and temporal dispersion of repolarization time and may lead to reentrant arrhythmias.