Simulation of ventricular rate control during atrial fibrillation using ionic channel blockers.

BACKGROUND: The atrioventricular (AV) node is the only compartment that conducts an electrical impulse between the atria and the ventricles. The main role of the AV node is to facilitate efficient pumping by conducting excitation slowly between the two chambers as well as reduce the ventricular rate during atrial fibrillation (AF). METHODS: Using computer simulations, we investigated excitation conduction from the right atrium to the bundle of His during high-rate atrial excitation with or without partial blocking of the calcium or potassium ionic current. RESULTS: Our simulations revealed differences in rate reduction and repolarization effects between calcium and potassium current blocking and high degree of potassium current blocking required to reduce the ventricular rate during AF. CONCLUSIONS: Our simulation results explain why potassium current blockers are not recommended for controlling ventricular rate during AF.

Characterization and influence of cardiac background sodium current in the atrioventricular node.

Background inward sodium current (IB,Na) that influences cardiac pacemaking has been comparatively under-investigated. The aim of this study was to determine for the first time the properties and role of IB,Na in cells from the heart's secondary pacemaker, the atrioventricular node (AVN). Myocytes were isolated from the AVN of adult male rabbits and mice using mechanical and enzymatic dispersion. Background current was measured using whole-cell patch clamp and monovalent ion substitution with major voltage- and time-dependent conductances inhibited. In the absence of a selective pharmacological inhibitor of IB,Na, computer modelling was used to assess the physiological contribution of IB,Na. Net background current during voltage ramps was linear, reversing close to 0mV. Switching between Tris- and Na(+)-containing extracellular solution in rabbit and mouse AVN cells revealed an inward IB,Na, with an increase in slope conductance in rabbit cells at -50mV from 0.54±0.03 to 0.91±0.05nS (mean±SEM; n=61 cells). IB,Na magnitude varied in proportion to [Na(+)]o. Other monovalent cations could substitute for Na(+) (Rb(+)>K(+)>Cs(+)>Na(+)>Li(+)). The single-channel conductance with Na(+) as charge carrier estimated from noise-analysis was 3.2±1.2pS (n=6). Ni(2+) (10mM), Gd(3+) (100µM), ruthenium red (100µM), or amiloride (1mM) produced modest reductions in IB,Na. Flufenamic acid was without significant effect, whilst La(3+) (100µM) or extracellular acidosis (pH6.3) inhibited the current by >60%. Under the conditions of our AVN cell simulations, removal of IB,Na arrested spontaneous activity and, in a simulated 1D-strand, reduced conduction velocity by ~20%. IB,Na is carried by distinct low conductance monovalent non-selective cation channels and can influence AVN spontaneous activity and conduction.

Splicing misregulation of SCN5A contributes to cardiac-conduction delay and heart arrhythmia in myotonic dystrophy.

Myotonic dystrophy (DM) is caused by the expression of mutant RNAs containing expanded CUG repeats that sequester muscleblind-like (MBNL) proteins, leading to alternative splicing changes. Cardiac alterations, characterized by conduction delays and arrhythmia, are the second most common cause of death in DM. Using RNA sequencing, here we identify novel splicing alterations in DM heart samples, including a switch from adult exon 6B towards fetal exon 6A in the cardiac sodium channel, SCN5A. We find that MBNL1 regulates alternative splicing of SCN5A mRNA and that the splicing variant of SCN5A produced in DM presents a reduced excitability compared with the control adult isoform. Importantly, reproducing splicing alteration of Scn5a in mice is sufficient to promote heart arrhythmia and cardiac-conduction delay, two predominant features of myotonic dystrophy. In conclusion, misregulation of the alternative splicing of SCN5A may contribute to a subset of the cardiac dysfunctions observed in myotonic dystrophy.

Three-dimensional computer model of the right atrium including the sinoatrial and atrioventricular nodes predicts classical nodal behaviours.

The aim of the study was to develop a three-dimensional (3D) anatomically-detailed model of the rabbit right atrium containing the sinoatrial and atrioventricular nodes to study the electrophysiology of the nodes. A model was generated based on 3D images of a rabbit heart (atria and part of ventricles), obtained using high-resolution magnetic resonance imaging. Segmentation was carried out semi-manually. A 3D right atrium array model (∼3.16 million elements), including eighteen objects, was constructed. For description of cellular electrophysiology, the Rogers-modified FitzHugh-Nagumo model was further modified to allow control of the major characteristics of the action potential with relatively low computational resource requirements. Model parameters were chosen to simulate the action potentials in the sinoatrial node, atrial muscle, inferior nodal extension and penetrating bundle. The block zone was simulated as passive tissue. The sinoatrial node, crista terminalis, main branch and roof bundle were considered as anisotropic. We have simulated normal and abnormal electrophysiology of the two nodes. In accordance with experimental findings: (i) during sinus rhythm, conduction occurs down the interatrial septum and into the atrioventricular node via the fast pathway (conduction down the crista terminalis and into the atrioventricular node via the slow pathway is slower); (ii) during atrial fibrillation, the sinoatrial node is protected from overdrive by its long refractory period; and (iii) during atrial fibrillation, the atrioventricular node reduces the frequency of action potentials reaching the ventricles. The model is able to simulate ventricular echo beats. In summary, a 3D anatomical model of the right atrium containing the cardiac conduction system is able to simulate a wide range of classical nodal behaviours.

Biophysical characterization of the underappreciated and important relationship between heart rate variability and heart rate.

Heart rate (HR) variability (HRV; beat-to-beat changes in the R-wave to R-wave interval) has attracted considerable attention during the past 30+ years (PubMed currently lists >17 000 publications). Clinically, a decrease in HRV is correlated to higher morbidity and mortality in diverse conditions, from heart disease to fetal distress. It is usually attributed to fluctuation in cardiac autonomic nerve activity. We calculated HRV parameters from a variety of cardiac preparations (including humans, living animals, Langendorff-perfused heart, and single sinoatrial nodal cell) in diverse species, combining this with data from previously published articles. We show that regardless of conditions, there is a universal exponential decay-like relationship between HRV and HR. Using 2 biophysical models, we develop a theory for this and confirm that HRV is primarily dependent on HR and cannot be used in any simple way to assess autonomic nerve activity to the heart. We suggest that the correlation between a change in HRV and altered morbidity and mortality is substantially attributable to the concurrent change in HR. This calls for re-evaluation of the findings from many articles that have not adjusted properly or at all for HR differences when comparing HRV in multiple circumstances.

Importance of gradients in membrane properties and electrical coupling in sinoatrial node pacing.

The sinoatrial node (SAN) is heterogeneous in terms of cell size, ion channels, current densities, connexins and electrical coupling. For example, Nav1.5 (responsible for INa) and Cx43 (responsible for electrical coupling) are absent from the centre of the SAN (normally the leading pacemaker site), but present in the periphery (at SAN-atrial muscle junction). To test whether the heterogeneity is important for the functioning of the SAN, one- and two-dimensional models of the SAN and surrounding atrial muscle were created. Normal functioning of the SAN (in terms of cycle length, position of leading pacemaker site, conduction times, activation and repolarization sequences and space constants) was observed when, from the centre to the periphery, (i) cell characteristics (cell size and ionic current densities) were changed in a gradient fashion from a central-type (lacking INa) to a peripheral-type (possessing INa) and (ii) coupling conductance was increased in a gradient fashion. We conclude that the heterogeneous nature of the node is important for its normal functioning. The presence of Nav1.5 and Cx43 in the periphery may be essential for the node to be able to drive the atrial muscle: Nav1.5 provides the necessary depolarizing current and Cx43 delivers it to the atrial muscle.

Simulation study of complex action potential conduction in atrioventricular node.

The atrioventricular (AV) node, which is located between the atria and ventricles of the heart, acts as important roles in cardiac excitation conduction between the two chambers. Although there are multiple conduction pathways in the AV node, the structure of the AV node has not been clarified. In this study, we constructed a one-dimensional model of the AV node and simulated excitation conduction between the right atrium and the bundle of His via the AV node. We also investigated several characteristics of the AV node: (1) responses of the AV node to high-rate excitation in the right atrium, (2) the AV nodal reentrant beat induced by premature stimulus, and (3) ventricular rate control during atrial fibrillation with various methods. Our simulation results suggest that multiple conduction pathways act as important roles in controlling the ventricular rate. The one-dimensional model constructed in this study may be useful to analyze complex conduction patterns in the AV node.

Structure, function and clinical relevance of the cardiac conduction system, including the atrioventricular ring and outflow tract tissues.

It is now over 100years since the discovery of the cardiac conduction system, consisting of three main parts, the sinus node, the atrioventricular node and the His-Purkinje system. The system is vital for the initiation and coordination of the heartbeat. Over the last decade, immense strides have been made in our understanding of the cardiac conduction system and these recent developments are reviewed here. It has been shown that the system has a unique embryological origin, distinct from that of the working myocardium, and is more extensive than originally thought with additional structures: atrioventricular rings, a third node (so called retroaortic node) and pulmonary and aortic sleeves. It has been shown that the expression of ion channels, intracellular Ca(2+)-handling proteins and gap junction channels in the system is specialised (different from that in the ordinary working myocardium), but appropriate to explain the functioning of the system, although there is continued debate concerning the ionic basis of pacemaking. We are beginning to understand the mechanisms (fibrosis and remodelling of ion channels and related proteins) responsible for dysfunction of the system (bradycardia, heart block and bundle branch block) associated with atrial fibrillation and heart failure and even athletic training. Equally, we are beginning to appreciate how naturally occurring mutations in ion channels cause congenital cardiac conduction system dysfunction. Finally, current therapies, the status of a new therapeutic strategy (use of a specific heart rate lowering drug) and a potential new therapeutic strategy (biopacemaking) are reviewed.

Connexins and the atrioventricular node.

The structure and functioning of the atrioventricular (AV) node has remained mysterious owing to its high degree of complexity. In this review article, we integrate advances in knowledge regarding connexin expression in the AV node. Complex patterning of 4 different connexin isoforms with single channel conductances ranging from ultralow to high explains the dual pathway electrophysiology of the AV node, the presence of 2 nodal extensions, longitudinal dissociation in the penetrating bundle, and, most importantly, how the AV node maintains slow conduction between the atria and the ventricles. It is shown that the complex patterning of connexins is the consequence of the embryonic development of the cardiac conduction system. Finally, it is argued that connexin dysregulation may be responsible for AV node dysfunction.

Pharmacological blockade of IKs destabilizes spiral-wave reentry under ß-adrenergic stimulation in favor of its early termination.

We tested a hypothesis that an enhancement of I(Ks) may play a pivotal role in ventricular proarrhythmia under high sympathetic activity. A 2-dimensional ventricular muscle layer was prepared in rabbit hearts, and action potential signals were analyzed by optical mapping. During constant stimulation, isoproterenol (ISP, 0.1 µM) significantly shortened action potential duration (APD); chromanol 293B (30 µM), a selective I(Ks)-blocker, reversed the APD shortening. VTs induced in the presence of ISP lasted longer than in the control, and this was reversed by 293B. E-4031 (0.1 µM), a selective I(Kr)-blocker, did not cause such reversal. Spiral-wave (SW) reentry with ISP was characterized by more stable rotation around a shorter functional block line (FBL) than in the control. After application of 293B, SW reentry was destabilized, and rotation around a longer FBL with prominent drift reappeared. The APD abbreviation by ISP close to the rotation center was more pronounced than in the periphery, leading to an opposite APD gradient (center < periphery) compared with controls. This effect was also reversed by 293B. In conclusion, ß-adrenergic stimulation stabilizes SW reentry most likely though an enhancement of I(Ks). Blockade of I(Ks) may be a promising therapeutic modality in prevention of ventricular tachyarrhythmias under high sympathetic activity.

Ageing-dependent remodelling of ion channel and Ca2+ clock genes underlying sino-atrial node pacemaking.

The function of the sino-atrial node (SAN), the pacemaker of the heart, is known to decline with age, resulting in pacemaker disease in the elderly. The aim of the study was to investigate the effects of ageing on the SAN by characterizing electrophysiological changes and determining whether changes in gene expression are involved. In young and old rats, SAN function was characterized in the anaesthetized animal, isolated heart and isolated right atrium using ECG and action potential recordings; gene expression was characterized using quantitative PCR. The SAN function declined with age as follows: the intrinsic heart rate declined by 18 ± 3%; the corrected SAN recovery time increased by 43 ± 13%; and the SAN action potential duration increased by 11 ± 3% (at 75% repolarization). Gene expression in the SAN changed considerably with age, e.g. there was an age-dependent decrease in the Ca(2+) clock gene, RYR2, and changes in many ion channels (e.g. increases in Na(v)1.5, Na(v)ß1 and Ca(v)1.2 and decreases in K(v)1.5 and HCN1). In conclusion, with age, there are changes in the expression of ion channel and Ca(2+) clock genes in the SAN, and the changes may provide a partial explanation for the age-dependent decline in pacemaker function.

Anatomical and molecular mapping of the left and right ventricular His-Purkinje conduction networks.

Functioning of the cardiac conduction system depends critically on its structure and its complement of ion channels. Therefore, the aim of this study was to document both the structure and ion channel expression of the left and right ventricular His-Purkinje networks, as we have previously done for the sinoatrial and atrioventricular nodes. A three-dimensional (3D) anatomical computer model of the His-Purkinje network of the rabbit heart was constructed after staining the network by immunoenzyme labelling of a marker protein, middle neurofilament. The bundle of His is a ribbon-like structure and the architecture of the His-Purkinje network differs between the left and right ventricles. The 3D model is able to explain the breakthrough points of the action potential on the ventricular epicardium during sinus rhythm. Using quantitative PCR, the expression levels of the major ion channels were measured in the free running left and right Purkinje fibres of the rabbit heart. Expression of ion channels differs from that of the working myocardium and can explain the specialised electrical activity of the Purkinje fibres as suggested by computer simulations; the expression profile of the left Purkinje fibres is more specialised than that of the right Purkinje fibres. The structure and ion channel expression of the Purkinje fibres are highly specialised and tailored to the functioning of the system. The His-Purkinje network in the left ventricle is more developed than that in the right ventricle and this may explain its greater clinical importance.

Computer three-dimensional anatomical reconstruction of the human sinus node and a novel paranodal area.

We have previously shown in rabbit that the pacemaker of the heart (the sinus node) is widespread and matches the wide distribution of the leading pacemaker site within the right atrium. There is, however, uncertainty about the precise location of the pacemaker in human heart, and its spatial relationships with the surrounding right atrial muscle. Therefore, the aim of the current study was to investigate the distribution of the sinus node tissue in a series of healthy human hearts and, for one of the hearts to construct a computer three-dimensional anatomical model of the sinus node, including the likely orientation of myocytes. A combination of experimental techniques--diffusion tensor magnetic resonance imaging (DT-MRI), histology, immunohistochemistry, image processing and computer modelling--was used. Our data show that the sinus node was larger than in previous studies and, most importantly, we identified a previously unknown area running alongside the sinus node (the "paranodal area"), which is even more extensive than the sinus node. This area possesses properties of both nodal and atrial tissues and may have a role in pacemaking. For example, it could explain the wide spread distribution of the leading pacemaker site in human right atrium, a phenomenon known as the wandering pacemaker observed in clinics. In summary, a novel 3D anatomical reconstruction presents a new picture of the distribution of nodal cells within the human right atrium.

Mechanistic links between Na+ channel (SCN5A) mutations and impaired cardiac pacemaking in sick sinus syndrome.

RATIONALE: Familial sick sinus syndrome (SSS) has been linked to loss-of-function mutations of the SCN5A gene, which result in decreased inward Na(+) current, I(Na). However, the functional role of I(Na) in cardiac pacemaking is controversial, and mechanistic links between mutations and sinus node dysfunction in SSS are unclear. OBJECTIVE: To determine mechanisms by which the SCN5A mutations impair cardiac pacemaking. METHODS AND RESULTS: Action potential (AP) models for rabbit sinoatrial node (SAN) cells were modified to incorporate experimentally reported I(Na) changes induced by 2 groups of SCN5A gene mutations (affecting the activation and inactivation of I(Na), respectively). The cell models were incorporated into an anatomically detailed 2D model of the intact SAN-atrium. Effects of the mutations and vagal nerve activity on cardiac pacemaking at the single-cell and tissue levels were studied. Multielectrode extracellular potential recordings of activation pattern from intact SAN-atrium preparations were performed to test predictions of the models. At the single-cell level, the mutations slowed down pacemaking rates in peripheral, but not in central SAN cells that control the heart rhythm. However, in tissue simulations, the mutations not only slowed down pacemaking, but also compromised AP conduction across the SAN-atrium, leading to a possible SAN exit block or sinus arrest, the major features of SSS. Simulated vagal nerve activity amplified the bradycardiac effects of the mutations. Two groups of SCN5A mutations showed subtle differences in impairing the ability of the SAN to drive the surrounding atrium, primarily attributable to their differential effects on atrial excitability and conduction safety. Experimental data with tetrodotoxin and carbachol confirmed the simulation outcomes. CONCLUSIONS: Our study substantiates the causative link between SCN5A gene mutations and SSS and illustrates mechanisms by which the mutations impair the driving ability of the SAN.

Molecular architecture of the human sinus node: insights into the function of the cardiac pacemaker.

BACKGROUND: Although we know much about the molecular makeup of the sinus node (SN) in small mammals, little is known about it in humans. The aims of the present study were to investigate the expression of ion channels in the human SN and to use the data to predict electrical activity. METHODS AND RESULTS: Quantitative polymerase chain reaction, in situ hybridization, and immunofluorescence were used to analyze 6 human tissue samples. Messenger RNA (mRNA) for 120 ion channels (and some related proteins) was measured in the SN, a novel paranodal area, and the right atrium (RA). The results showed, for example, that in the SN compared with the RA, there was a lower expression of Na(v)1.5, K(v)4.3, K(v)1.5, ERG, K(ir)2.1, K(ir)6.2, RyR2, SERCA2a, Cx40, and Cx43 mRNAs but a higher expression of Ca(v)1.3, Ca(v)3.1, HCN1, and HCN4 mRNAs. The expression pattern of many ion channels in the paranodal area was intermediate between that of the SN and RA; however, compared with the SN and RA, the paranodal area showed greater expression of K(v)4.2, K(ir)6.1, TASK1, SK2, and MiRP2. Expression of ion channel proteins was in agreement with expression of the corresponding mRNAs. The levels of mRNA in the SN, as a percentage of those in the RA, were used to estimate conductances of key ionic currents as a percentage of those in a mathematical model of human atrial action potential. The resulting SN model successfully produced pacemaking. CONCLUSIONS: Ion channels show a complex and heterogeneous pattern of expression in the SN, paranodal area, and RA in humans, and the expression pattern is appropriate to explain pacemaking.

Atrial proarrhythmia due to increased inward rectifier current (I(K1)) arising from KCNJ2 mutation--a simulation study.

Atrial fibrillation (AF) has been linked to increased inward rectifier potassium current, I(K1), either due to AF-induced electrical remodelling, or from functional changes due to the Kir2.1 V93I mutation. The aim of this simulation study was to identify at cell and tissue levels' mechanisms by which increased I(K1) facilitates and perpetuates AF. The Courtemanche et al. human atrial cell action potential (AP) model was modified to incorporate reported changes in I(K1) induced by the Kir2.1 V93I mutation in both heterozygous (Het) and homozygous (Hom) mutant forms. The modified models for wild type (WT), Het and Hom conditions were incorporated into homogeneous 1D, 2D and 3D tissue models. Restitution curves of AP duration (APD), effective refractory period (ERP) and conduction velocity (CV) were computed and both the temporal and the spatial vulnerability of atrial tissue to re-entry were measured. The lifespan and tip meandering pattern of re-entry were also characterised. For comparison, parallel simulations were performed by incorporating into the Courtmanche et al. model a linear increase in maximal I(K1) conductance. It was found that the gain-in-function of V93I 'mutant'I(K1) led to abbreviated atrial APs and flattened APD, ERP and CV restitution curves. It also hyperpolarised atrial resting membrane potential and slowed down intra-atrial conduction. V93I 'mutant'I(K1) reduced the tissue's temporal vulnerability but increased spatial vulnerability to initiate and sustain re-entry, resulting in an increased overall susceptibility of atrial tissue to arrhythmogenesis. In the 2D model, spiral waves self-terminated for WT (lifespan < 3.3 s) tissue, but persisted in Het and Hom tissues for the whole simulation period (lifespan > 10 s). The tip of the spiral wave meandered more in WT tissue than in Het and Hom tissues. Increased I(K1) due to augmented maximal conductance produced similar results to those of Het and Hom Kir2.1 V93I mutant conditions. In the 3D model the dynamic behaviour of scroll waves was stabilized by increased I(K1). In conclusion, increased I(K1) current, either by the Kir2.1 V93I mutation or by augmented maximal conductance, increases atrial susceptibility to arrhythmia by increasing the lifespan of re-entrant spiral waves and the stability of scroll waves in 3D tissue, thereby facilitating initiation and maintenance of re-entrant circuits.

Computer three-dimensional reconstruction of the atrioventricular node.

Because of its complexity, the atrioventricular node (AVN), remains 1 of the least understood regions of the heart. The aim of the study was to construct a detailed anatomic model of the AVN and relate it to AVN function. The electric activity of a rabbit AVN preparation was imaged using voltage-dependent dye. The preparation was then fixed and sectioned. Sixty-five sections at 60- to 340-microm intervals were stained for histology and immunolabeled for neurofilament (marker of nodal tissue) and connexin43 (gap junction protein). This revealed multiple structures within and around the AVN, including transitional tissue, inferior nodal extension, penetrating bundle, His bundle, atrial and ventricular muscle, central fibrous body, tendon of Todaro, and valves. A 3D anatomically detailed mathematical model (approximately 13 million element array) of the AVN and surrounding atrium and ventricle, incorporating all cell types, was constructed. Comparison of the model with electric activity recorded in experiments suggests that the inferior nodal extension forms the slow pathway, whereas the transitional tissue forms the fast pathway into the AVN. In addition, it suggests the pacemaker activity of the atrioventricular junction originates in the inferior nodal extension. Computer simulation of the propagation of the action potential through the anatomic model shows how, because of the complex structure of the AVN, reentry (slow-fast and fast-slow) can occur. In summary, a mathematical model of the anatomy of the AVN has been generated that allows AVN conduction to be explored.

Imaging the heart: computer 3-dimensional anatomic models of the heart.

Since the 1960s, models of the action potential in various cardiac cell types have been developed, and since the 1990s, 3-dimensional anatomic (or geometric) models of various cardiac structures have been developed. We are approaching the time when, for one species, we should have a complete set of action potential and anatomic models for the various cardiac tissues and then we will have realized the aim of constructing a "virtual heart" with accurate anatomy and electrophysiology. However, already the two types of model are beginning to be used in tandem to reconstruct the activation sequence of the heart both during sinus rhythm and arrhythmias.

Pacing-induced spontaneous activity in myocardial sleeves of pulmonary veins after treatment with ryanodine.

BACKGROUND: Recent clinical electrophysiology studies and successful results of radiofrequency catheter ablation therapy suggest that high-frequency focal activity in the pulmonary veins (PVs) plays important roles in the initiation and perpetuation of atrial fibrillation, but the mechanisms underlying the focal arrhythmogenic activity are not understood. METHODS AND RESULTS: Extracellular potential mapping of rabbit right atrial preparations showed that ryanodine (2 micromol/L) caused a shift of the leading pacemaker from the sinoatrial node to an ectopic focus near the right PV-atrium junction. The transmembrane potential recorded from the isolated myocardial sleeve of the right PV showed typical atrial-type action potentials with a stable resting potential under control conditions. Treatment with ryanodine (0.5 to 2 micromol/L) resulted in a depolarization of the resting potential and a development of pacemaker depolarization. These changes were enhanced transiently after an increase in the pacing rate: a self-terminating burst of spontaneous action potentials (duration, 33.6+/-5.0 s; n=32) was induced by a train of rapid stimuli (3.3 Hz) applied after a brief rest period. The pacing-induced activity was attenuated by either depletion of the sarcoplasmic reticulum of Ca2+ or blockade of the sarcolemmal Na+-Ca2+ exchanger or Cl- channels and potentiated by beta-adrenergic stimulation. CONCLUSIONS: PV myocardial sleeves have the potential to generate spontaneous activity, and such arrhythmogenic activity is uncovered by modulation of intracellular Ca2+ dynamics.