Computational modelling of mouse atrio ventricular node action potential and automaticity.

The atrioventricular node (AVN) is a crucial component of the cardiac conduction system. Despite its pivotal role in regulating the transmission of electrical signals between atria and ventricles, a comprehensive understanding of the cellular electrophysiological mechanisms governing AVN function has remained elusive. This paper presents a detailed computational model of mouse AVN cell action potential (AP). Our model builds upon previous work and introduces several key refinements, including accurate representation of membrane currents and exchangers, calcium handling, cellular compartmentalization, dynamic update of intracellular ion concentrations, and calcium buffering. We recalibrated and validated the model against existing and unpublished experimental data. In control conditions, our model reproduces the AVN AP experimental features, (e.g. rate = 175 bpm, experimental range [121, 191] bpm). Notably, our study sheds light on the contribution of L-type calcium currents, through both Cav1.2 and Cav1.3 channels, in AVN cells. The model replicates several experimental observations, including the cessation of firing upon block of Cav1.3 or INa,r current. If block induces a reduction in beating rate of 11%. In summary, this work presents a comprehensive computational model of mouse AVN cell AP, offering a valuable tool for investigating pacemaking mechanisms and simulating the impact of ionic current blockades. By integrating calcium handling and refining formulation of ionic currents, our model advances understanding of this critical component of the cardiac conduction system, providing a platform for future developments in cardiac electrophysiology. KEY POINTS: This paper introduces a comprehensive computational model of mouse atrioventricular node (AVN) cell action potentials (APs). Our model is based on the electrophysiological data from isolated mouse AVN cells and exhibits an action potential and calcium transient that closely match the experimental records. By simulating the effects of blocking specific ionic currents, the model effectively predicts the roles of L-type Cav1.2 and Cav1.3 channels, T-type calcium channels, sodium currents (TTX-sensitive and TTX-resistant), and the funny current (If) in AVN pacemaking. The study also emphasizes the significance of other ionic currents, including IKr, Ito, IKur, in regulating AP characteristics and cycle length in AVN cells. The model faithfully reproduces the rate dependence of action potentials under pacing, opening the possibility of use in impulse propagation models. The population-of-models approach showed the robustness of this new AP model in simulating a wide spectrum of cellular pacemaking in AVN.

A GABAergic system in atrioventricular node pacemaker cells controls electrical conduction between the atria and ventricles.

Physiologically, the atria contract first, followed by the ventricles, which is the prerequisite for normal blood circulation. The above phenomenon of atrioventricular sequential contraction results from the characteristically slow conduction of electrical excitation of the atrioventricular node (AVN) between the atria and the ventricles. However, it is not clear what controls the conduction of electrical excitation within AVNs. Here, we find that AVN pacemaker cells (AVNPCs) possess an intact intrinsic GABAergic system, which plays a key role in electrical conduction from the atria to the ventricles. First, along with the discovery of abundant GABA-containing vesicles under the surface membranes of AVNPCs, key elements of the GABAergic system, including GABA metabolic enzymes, GABA receptors, and GABA transporters, were identified in AVNPCs. Second, GABA synchronously elicited GABA-gated currents in AVNPCs, which significantly weakened the excitability of AVNPCs. Third, the key molecular elements of the GABAergic system markedly modulated the conductivity of electrical excitation in the AVN. Fourth, GABAA receptor deficiency in AVNPCs accelerated atrioventricular conduction, which impaired the AVN's protective potential against rapid ventricular frequency responses, increased susceptibility to lethal ventricular arrhythmias, and decreased the cardiac contractile function. Finally, interventions targeting the GABAergic system effectively prevented the occurrence and development of atrioventricular block. In summary, the endogenous GABAergic system in AVNPCs determines the slow conduction of electrical excitation within AVNs, thereby ensuring sequential atrioventricular contraction. The endogenous GABAergic system shows promise as a novel intervention target for cardiac arrhythmias.

Development of the Cardiac Conduction System.

The electrical impulses that coordinate the sequential, rhythmic contractions of the atria and ventricles are initiated and tightly regulated by the specialized tissues of the cardiac conduction system. In the mature heart, these impulses are generated by the pacemaker cardiomyocytes of the sinoatrial node, propagated through the atria to the atrioventricular node where they are delayed and then rapidly propagated to the atrioventricular bundle, right and left bundle branches, and finally, the peripheral ventricular conduction system. Each of these specialized components arise by complex patterning events during embryonic development. This chapter addresses the origins and transcriptional networks and signaling pathways that drive the development and maintain the function of the cardiac conduction system.

Atrioventricular and Ventricular Functional Interdependence in Individuals Without Overt Cardiac Disease.

Background Left atrial (LA) and right ventricular (RV) performance play an integral role in the pathophysiology and prognosis of heart failure. We hypothesized that subclinical left ventricular dysfunction adversely affects LA/RV geometry and function even in a preclinical setting. This study aimed to investigate the atrioventricular and ventricular functional interdependence in a community-based cohort without overt cardiovascular disease. Methods and Results Left ventricular global longitudinal strain (LVGLS), RV free-wall longitudinal strain and LA phasic strain were assessed by speckle-tracking echocardiography in 1080 participants (600 men; 62±12 years) between 2014 and 2018. One hundred and forty-three participants (13.2%) had an abnormal LVGLS (>-18.6%). LA reservoir strain, conduit strain, and RV free-wall longitudinal strain were significantly decreased in abnormal LVGLS group compared with normal LVGLS group (all P<0.001). LA and RV dysfunction (LA reservoir strain<31.4% and RVLS>-19.2%) were present in 18.9% and 19.6% of participants with abnormal LVGLS. Decreased LVGLS was associated with worse LA reservoir strain, conduit strain and RV free-wall longitudinal strain (standardized ß=-0.20, -0.19 and 0.11 respectively, all P<0.01) independent of cardiovascular risk factors. LA and/or RV dysfunction concomitant with abnormal LVGLS carried significantly increased risk of elevated B-type natriuretic peptide levels (>28.6 pg/mL for men and >44.4 pg/mL for women) compared with normal LVGLS (odds ratio, 2.01; P=0.030). Conclusions LA/RV dysfunction was present in 20% individuals with abnormal LVGLS and multi-chamber impairment was associated with elevated B-type natriuretic peptide level, which may provide valuable insights for a better understanding of atrioventricular and ventricular interdependence and possibly heart failure preventive strategies.

Inducible cardiomyocyte injury within the atrioventricular conduction system uncovers latent regenerative capacity in mice.

The cardiac conduction system (CCS) ensures regular contractile function, and injury to any of its components can cause cardiac dysrhythmia. Although all cardiomyocytes (CMs) originate from common progenitors, the CCS is composed of biologically distinct cell types with unique functional and developmental characteristics. In contrast to ventricular cardiomyocytes, which continue to proliferate after birth, most CCS cells terminally exit the cell cycle during fetal development. Although the CCS should thus provide a poor substrate for postnatal injury repair, its regenerative capacity remains untested. Here, we describe a genetic system for ablating CMs that reside within the atrioventricular conduction system (AVCS). Adult mouse AVCS ablation resulted in regenerative failure characterized by persistent atrioventricular conduction defects and contractile dysfunction. In contrast, AVCS injury in neonatal mice led to recovery in a subset of these mice, thus providing evidence for CCS plasticity. Furthermore, CM proliferation did not appear to completely account for the observed functional recovery, suggesting that mechanisms regulating recovery from dysrhythmia are likely to be distinct from cardiac regeneration associated with ventricular injury. Taken together, we anticipate that our results will motivate further mechanistic studies of CCS plasticity and enable the exploration of rhythm restoration as an alternative therapeutic strategy.

Countering reproducibility issues in mathematical models with software engineering techniques: A case study using a one-dimensional mathematical model of the atrioventricular node.

One should assume that in silico experiments in systems biology are less susceptible to reproducibility issues than their wet-lab counterparts, because they are free from natural biological variations and their environment can be fully controlled. However, recent studies show that only half of the published mathematical models of biological systems can be reproduced without substantial effort. In this article we examine the potential causes for failed or cumbersome reproductions in a case study of a one-dimensional mathematical model of the atrioventricular node, which took us four months to reproduce. The model demonstrates that even otherwise rigorous studies can be hard to reproduce due to missing information, errors in equations and parameters, a lack in available data files, non-executable code, missing or incomplete experiment protocols, and missing rationales behind equations. Many of these issues seem similar to problems that have been solved in software engineering using techniques such as unit testing, regression tests, continuous integration, version control, archival services, and a thorough modular design with extensive documentation. Applying these techniques, we reimplement the examined model using the modeling language Modelica. The resulting workflow is independent of the model and can be translated to SBML, CellML, and other languages. It guarantees methods reproducibility by executing automated tests in a virtual machine on a server that is physically separated from the development environment. Additionally, it facilitates results reproducibility, because the model is more understandable and because the complete model code, experiment protocols, and simulation data are published and can be accessed in the exact version that was used in this article. We found the additional design and documentation effort well justified, even just considering the immediate benefits during development such as easier and faster debugging, increased understandability of equations, and a reduced requirement for looking up details from the literature.

Distinct circadian mechanisms govern cardiac rhythms and susceptibility to arrhythmia.

Electrical activity in the heart exhibits 24-hour rhythmicity, and potentially fatal arrhythmias are more likely to occur at specific times of day. Here, we demonstrate that circadian clocks within the brain and heart set daily rhythms in sinoatrial (SA) and atrioventricular (AV) node activity, and impose a time-of-day dependent susceptibility to ventricular arrhythmia. Critically, the balance of circadian inputs from the autonomic nervous system and cardiomyocyte clock to the SA and AV nodes differ, and this renders the cardiac conduction system sensitive to decoupling during abrupt shifts in behavioural routine and sleep-wake timing. Our findings reveal a functional segregation of circadian control across the heart's conduction system and inherent susceptibility to arrhythmia.

Innervation and Neuronal Control of the Mammalian Sinoatrial Node a Comprehensive Atlas.

[Figure: see text].

T-box transcription factor 3 governs a transcriptional program for the function of the mouse atrioventricular conduction system.

Genome-wide association studies have identified noncoding variants near TBX3 that are associated with PR interval and QRS duration, suggesting that subtle changes in TBX3 expression affect atrioventricular conduction system function. To explore whether and to what extent the atrioventricular conduction system is affected by Tbx3 dose reduction, we first characterized electrophysiological properties and morphology of heterozygous Tbx3 mutant (Tbx3+/-) mouse hearts. We found PR interval shortening and prolonged QRS duration, as well as atrioventricular bundle hypoplasia after birth in heterozygous mice. The atrioventricular node size was unaffected. Transcriptomic analysis of atrioventricular nodes isolated by laser capture microdissection revealed hundreds of deregulated genes in Tbx3+/- mutants. Notably, Tbx3+/- atrioventricular nodes showed increased expression of working myocardial gene programs (mitochondrial and metabolic processes, muscle contractility) and reduced expression of pacemaker gene programs (neuronal, Wnt signaling, calcium/ion channel activity). By integrating chromatin accessibility profiles (ATAC sequencing) of atrioventricular tissue and other epigenetic data, we identified Tbx3-dependent atrioventricular regulatory DNA elements (REs) on a genome-wide scale. We used transgenic reporter assays to determine the functionality of candidate REs near Ryr2, an up-regulated chamber-enriched gene, and in Cacna1g, a down-regulated conduction system-specific gene. Using genome editing to delete candidate REs, we showed that a strong intronic bipartite RE selectively governs Cacna1g expression in the conduction system in vivo. Our data provide insights into the multifactorial Tbx3-dependent transcriptional network that regulates the structure and function of the cardiac conduction system, which may underlie the differences in PR duration and QRS interval between individuals carrying variants in the TBX3 locus.

Cardiomyocyte Expression of ZO-1 Is Essential for Normal Atrioventricular Conduction but Does Not Alter Ventricular Function.

RATIONALE: ZO-1 (Zonula occludens-1), a plasma membrane-associated scaffolding protein regulates signal transduction, transcription, and cellular communication. Global deletion of ZO-1 in the mouse is lethal by embryonic day 11.5. The function of ZO-1 in cardiac myocytes (CM) is largely unknown. OBJECTIVE: To determine the function of CM ZO-1 in the intact heart, given its binding to other CM proteins that have been shown instrumental in normal cardiac conduction and function. METHODS AND RESULTS: We generated ZO-1 CM-specific knockout (KO) mice using α-Myosin Heavy Chain-nuclear Cre (ZO-1cKO) and investigated physiological and electrophysiological function by echocardiography, surface ECG and conscious telemetry, intracardiac electrograms and pacing, and optical mapping studies. ZO-1cKO mice were viable, had normal Mendelian ratios, and had a normal lifespan. Ventricular morphometry and function were not significantly different between the ZO-1cKO versus control (CTL) mice, basally in young or aged mice, or even when hearts were subjected to hemodynamic loading. Atrial mass was increased in ZO-1cKO. Electrophysiological and optical mapping studies indicated high-grade atrioventricular (A-V) block in ZO-1cKO comparing to CTL hearts. While ZO-1-associated proteins such as vinculin, connexin 43, N-cadherin, and α-catenin showed no significant change with the loss of ZO-1, Connexin-45 and Coxsackie-adenovirus (CAR) proteins were reduced in atria of ZO-1cKO. Further, with loss of ZO-1, ZO-2 protein was increased significantly in ventricular CM in a presumed compensatory manner but was still not detected in the AV nodal myocytes. Importantly, the expression of the sodium channel protein NaV1.5 was altered in AV nodal cells of the ZO-1cKO versus CTL. CONCLUSIONS: ZO-1 protein has a unique physiological role in cardiac nodal tissue. This is in alignment with its known interaction with CAR and Cx45, and a new function in regulating the expression of NaV1.5 in AV node. Uniquely, ZO-1 is dispensable for function of the working myocardium.

ZO-1 Regulates Intercalated Disc Composition and Atrioventricular Node Conduction.

RATIONALE: ZO-1 (Zona occludens 1), encoded by the tight junction protein 1 (TJP1) gene, is a regulator of paracellular permeability in epithelia and endothelia. ZO-1 interacts with the actin cytoskeleton, gap, and adherens junction proteins and localizes to intercalated discs in cardiomyocytes. However, the contribution of ZO-1 to cardiac physiology remains poorly defined. OBJECTIVE: We aim to determine the role of ZO-1 in cardiac function. METHODS AND RESULTS: Inducible cardiomyocyte-specific Tjp1 deletion mice (Tjp1fl/fl; Myh6Cre/Esr1*) were generated by crossing the Tjp1 floxed mice and Myh6Cre/Esr1* transgenic mice. Tamoxifen-induced loss of ZO-1 led to atrioventricular (AV) block without changes in heart rate, as measured by ECG and ex vivo optical mapping. Mice with tamoxifen-induced conduction system-specific deletion of Tjp1 (Tjp1fl/fl; Hcn4CreERt2) developed AV block while tamoxifen-induced conduction system deletion of Tjp1 distal to the AV node (Tjp1fl/fl; Kcne1CreERt2) did not demonstrate conduction defects. Western blot and immunostaining analyses of AV nodes showed that ZO-1 loss decreased Cx (connexin) 40 expression and intercalated disc localization. Consistent with the mouse model study, immunohistochemical staining showed that ZO-1 is abundantly expressed in the human AV node and colocalizes with Cx40. Ventricular conduction was not altered despite decreased localization of ZO-1 and Cx43 at the ventricular intercalated disc and modestly decreased left ventricular ejection fraction, suggesting ZO-1 is differentially required for AV node and ventricular conduction. CONCLUSIONS: ZO-1 is a key protein responsible for maintaining appropriate AV node conduction through maintaining gap junction protein localization.

Permanent His bundle pacing and atrioventricular node ablation in a case with drug-refractory atrial fibrillation.

Permanent His bundle pacing (HBP) activates the ventricles through the normal conduction system and has become a useful technique for patients with a high ventricular pacing rate. Presently described is a case of drug-refractory atrial fibrillation (AF) with a high ventricular rate that was treated with atrioventricular (AV) node ablation and permanent HBP. A 62-year-old woman with persistent AF and a drug-refractory high ventricular response was referred for exercise intolerance and palpitation. She had a history of failed catheter ablation attempts and amiodarone toxicity. Permanent HBP and AV node ablation was planned to achieve rate control with a stepwise approach. Initially, implantation of a permanent pacemaker was performed. The His lead and right ventricular back-up leads were implanted successfully, in the manner described previously. The His lead was connected to the atrial channel of the pacemaker battery and programmed to AAI pacing mode. The AV node was ablated successfully 3 weeks later without any threshold changes in the His lead. No His lead threshold changes were observed during or after AV node ablation and the patient was subsequently asymptomatic with twice daily apixaban 5 mg. Permanent HBP after AV node ablation can be a beneficial treatment option to prevent pacing-induced ventricular dyssynchrony and heart failure in patients who are not eligible for cardiac resynchronization therapy.

[The Discovery of the Cardiac Atrioventricular Node by Sunao Tawara and Ludwig Aschoff]. / Die Entdeckung des AV-Knotens des Herzens durch Sunao Tawara und Ludwig Aschoff.

Already in 1664, the Danish anatomist and naturalist Niels Stensen proved that the heart is a muscle. But for a long time it remained unclear what triggered the heart contractions.The Dutch physiologist Willem Einthoven registered the electrical processes in the contraction of the heart muscle and thus provided the first electrophysiological basis of cardiac muscle activity. Since 1903, Sunao Tawara was assistant to Ludwig Aschoff in Marburg. Both left Marburg in 1906: Tawara went back to Japan and Aschoff to Freiburg. In 1905, Tawara discovered the connections of the His' bundle to the AV node and the Purkinje fibers. At that time, there was no thought of a functional interpretation. Tawara discovered a kind of "knot" that linked to the adjacent myocardial cells, as well as the "Tawara thighs", which frayed and went into structures known as Purkinje fibers. Tawara detected the tree-like structure he had discovered as a muscle-fiber system that controlled the arousal of the heart's musculature. Thus the old dispute between myogenic and neurogenic arousal of the heart was decided in favor of the myogenic excitation conduction. The atrioventricular node described by Tawara was given the eponym "Aschoff-Tawara node". Tawara's groundbreaking work on the conduction system was the basis for the discovery of the sinus node and the interpretation of the heart's electrophysiology.

Histological and morphometric study of the components of the sinus and atrioventricular nodes in horses and dogs.

The cardiac nodes are the source of the electrical impulse that is transmitted to the heart, the aim of this work is study the histological and morphometric characteristics of the different components of the sinus and atrioventricular nodes in horses and dogs that help to know the physiopathology of these nodes. A group of ten horse hearts and five dog hearts were used. The region of the sinus and atrioventricular nodes was sectioned serially, and the block of tissue removed for study. The samples were assessed using a morphometric analysis with the Image-Pro Plus 7.1 software and the acquisition of the images using a Leica DMD108 optic microscope. The shape of the horse's sinus node is oblong and its P cells are large. The shape of the dog's sinus is rounded or oblong. The P cells are large and pale. The area of P cells in horses was 976 (SD 223.7) µm2 and in dogs the area for P cells was 106 (SD 30.4) µm2, which indicates that the value for P cells in horses are significantly higher than in dogs (p = .001). The horse atrioventricular node presented an oblong shape and in dogs, presents a spindle shape. The lower cell density in any of the cardiac nodes, especially in P cells of sinus node, can decrease electrical conduction within the nodes and in the internodal tracts, which would reflect the presence of cardiac arrhythmias derived from poor conduction, even in morphologically normal hearts.

An integrated overview of AV node physiology.

The atrioventricular (AV) node generates half of the AV delay needed for blood pumping and filters atrial impulses that could otherwise induce life-threatening ventricular arrhythmias. It is also a pacemaker and a key target in the treatment of cardiac arrhythmias. The special roles of the AV node primarily arise from its slow conduction, long refractory period, and cellular automaticity. However, efforts to establish the dynamics of these properties and their interaction led to many controversies. In fact, the AV node's behavior is so complex that it seems to escape broadly applicable rules. This review summarizes progresses made in resolving these issues and in integrating the multiple roles of the AV node within a common functional model. Presented evidence shows that the rate-dependent conduction and refractory properties of the AV node can be reliably characterized and reconciled from nodal responses to S1 S2 S3 protocols. It also supports the concept that dual pathways constitute a feature of the normal AV node and account for its overall conduction and refractory properties. In this model, the posterior extension and compact node provide the core of the slow and fast pathway, respectively. The transitional tissues and lower nodal bundle provide a common proximal and distal pathway, respectively. These pathways would also support bidirectional conduction. The dual pathway involvement can also be extended to widely variable AV nodal responses, such as Wenckebach cycles, hysteresis, and ventricular response to atrial fibrillation. In brief, the intricate AV nodal behavior may obey a limited set of accessible and definable rules.

Using Gjd3-CreEGFP mice to examine atrioventricular node morphology and composition.

The atrioventricular node (AVN) coordinates the timing of atrial and ventricular contraction to optimize cardiac performance. To study this critical function using mouse genetics, however, new reagents are needed that allow AVN-specific manipulation. Here we describe a novel Gjd3-CreEGFP mouse line that successfully recombines floxed alleles within the AVN beginning at E12.5. These mice have been engineered to express CreEGFP under the control of endogenous Gjd3 regulatory elements without perturbing native protein expression. Detailed histological analysis of Gjd3-CreEGFP mice reveals specific labeling of AVN cardiomyocytes and a subset of cardiac endothelial cells. Importantly, we show that Gjd3-CreEGFP mice have preserved cardiac mechanical and electrical function. In one application of our newly described mouse line, we provide a three-dimensional (3D) view of the AVN using tissue clearing combined with confocal microscopy. With this 3D model as a reference, we identify specific AVN sub-structures based on marker staining characteristics. In addition, we use our Gjd3-CreEGFP mice to guide microdissection of the AVN and construction of a single-cell atlas. Thus, our results establish a new transgenic tool for AVN-specific recombination, provide an updated model of AVN morphology, and describe a roadmap for exploring AVN cellular heterogeneity.

PR interval genome-wide association meta-analysis identifies 50 loci associated with atrial and atrioventricular electrical activity.

Electrocardiographic PR interval measures atrio-ventricular depolarization and conduction, and abnormal PR interval is a risk factor for atrial fibrillation and heart block. Our genome-wide association study of over 92,000 European-descent individuals identifies 44 PR interval loci (34 novel). Examination of these loci reveals known and previously not-yet-reported biological processes involved in cardiac atrial electrical activity. Genes in these loci are over-represented in cardiac disease processes including heart block and atrial fibrillation. Variants in over half of the 44 loci were associated with atrial or blood transcript expression levels, or were in high linkage disequilibrium with missense variants. Six additional loci were identified either by meta-analysis of ~105,000 African and European-descent individuals and/or by pleiotropic analyses combining PR interval with heart rate, QRS interval, and atrial fibrillation. These findings implicate developmental pathways, and identify transcription factors, ion-channel genes, and cell-junction/cell-signaling proteins in atrio-ventricular conduction, identifying potential targets for drug development.