ParamAP: Standardized Parameterization of Sinoatrial Node Myocyte Action Potentials

A modified method for isolating sinoatrial node myocytes from adult mice

Sinoatrial node (SAN) is the pacemaker of the heart in charge of initiating spontaneous electronical activity and controlling heart rate. Myocytes from SAN can generate spontaneous rhythmic action potentials, which propagate through the myocardium, thereby triggering cardiac myocyte contraction. Acutely, the method for isolating sinoatrial node myocytes (SAMs) is critical in studying the protein expression and function of myocytes in SAN. Currently, the SAMs were isolated by transferring SAN tissue directly into the digestion solution, but it is difficult to judge the degree of digestion, and the system was unstable. Here, we present a modified protocol for the isolation of SAMs in mice, based on the collagenase II and protease perfusion of the heart using a Langendorff apparatus and subsequent dissociation of SAMs. The appearance and droplet flow rate of the heart could be significantly changed during enzymatic digestion via perfusion, which allowed us to easily judge the degree of digestion and avoid incomplete or excessive digestion. The SAMs with stable yield and viability achieved from our optimized approach would facilitate the follow-up experiments.

SanPy: Software for the analysis and visualization of whole-cell current-clamp recordings

The analysis of action potentials and other membrane voltage fluctuations provides a powerful approach for interrogating the function of excitable cells. However, a major bottleneck in the interpretation of this critical data is the lack of intuitive, agreed-upon software tools for its analysis. Here, we present SanPy, an open-source and freely available software package for the analysis and exploration of whole-cell current-clamp recordings written in Python. SanPy provides a robust computational engine with an application programming interface. Using this, we have developed a cross-platform desktop application with a graphical user interface that does not require programming. SanPy is designed to extract common parameters from action potentials, including threshold time and voltage, peak, half-width, and interval statistics. In addition, several cardiac parameters are measured, including the early diastolic duration and rate. SanPy is built to be fully extensible by providing a plugin architecture for the addition of new file loaders, analysis, and visualizations. A key feature of SanPy is its focus on quality control and data exploration. In the desktop interface, all plots of the data and analysis are linked, allowing simultaneous data visualization from different dimensions with the goal of obtaining ground-truth analysis. We provide documentation for all aspects of SanPy, including several use cases and examples. To test SanPy, we performed analysis on current-clamp recordings from heart and brain cells. Taken together, SanPy is a powerful tool for whole-cell current-clamp analysis and lays the foundation for future extension by the scientific community.

SanPy: A whole-cell electrophysiology analysis pipeline

The analysis of action potentials and other membrane voltage fluctuations provide a powerful approach for interrogating the function of excitable cells. Yet, a major bottleneck in the interpretation of this critical data is the lack of intuitive, agreed upon software tools for its analysis. Here, we present SanPy, a Python-based open-source and freely available software pipeline for the analysis and exploration of whole-cell current-clamp recordings. SanPy provides a robust computational engine with an application programming interface. Using this, we have developed a cross-platform graphical user interface that does not require programming. SanPy is designed to extract common parameters from action potentials including threshold time and voltage, peak, half-width, and interval statistics. In addition, several cardiac parameters are measured including the early diastolic duration and rate. SanPy is built to be fully extensible by providing a plugin architecture for the addition of new file loaders, analysis, and visualizations. A key feature of SanPy is its focus on quality control and data exploration. In the desktop interface, all plots of the data and analysis are linked allowing simultaneous data visualization from different dimensions with the goal of obtaining ground truth analysis. We provide documentation for all aspects of SanPy including several use cases and examples. To test SanPy, we have performed analysis on current-clamp recordings from heart and brain cells. Taken together, SanPy is a powerful tool for whole-cell current-clamp analysis and lays the foundation for future extension by the scientific community.

ActionPytential: An open source tool for analyzing and visualizing cardiac action potential data

The action potential forms the basis of cardiac pacemaking, conduction, and contraction. Action potentials can be recorded from numerous preparation types, including ventricular or atrial trabecules, Purkinje fibers, isolated cardiac myocytes. Numerous techniques are also available as well, such as the conventional microelectrode and the single-cell current clamp techniques, optical mapping, or in silico modeling. With such a vast array of electrophysiological methods comes an array of available hardware and software solutions. In this work, we present a software with an intuitive graphical user interface, ActionPytential, that enables the analysis of any type of cardiac action potential, regardless of acquisition method or tissue type. In most available software tools, the analysis of continuous (gap-free) recordings often requires manual user interaction to segment the individual action potentials. We provide an automated solution for this, both for slow-response and for externally paced action potentials. As of now, ActionPytential calculates 34 parameters from each action potential. The most often utilized ones, including amplitude, maximal rate of depolarization, and action potential duration values, were validated on 1200 action potentials from human, dog, rabbit, guinea pig, and rat cardiac preparations. We also provide new parameters that were previously only measurable manually, including the position and the depth of the notch in potentials showing a spike-and-dome morphology. Further notable features include a Butterworth-type low-pass filter, the averaging of multiple potentials, automated corrections for baseline drifting, aided manual analysis, high-quality plots, and batch processing for any number of potentials. ActionPytential is available for all major platforms (Windows, MacOS, GNU + Linux, BSD).

The virtual sinoatrial node: What did computational models tell us about cardiac pacemaking?

Since its discovery, the sinoatrial node (SAN) has represented a fascinating and complex matter of research. Despite over a century of discoveries, a full comprehension of pacemaking has still to be achieved. Experiments often produced conflicting evidence that was used either in support or against alternative theories, originating intense debates. In this context, mathematical descriptions of the phenomena underlying the heartbeat have grown in importance in the last decades since they helped in gaining insights where experimental evaluation could not reach. This review presents the most updated SAN computational models and discusses their contribution to our understanding of cardiac pacemaking. Electrophysiological, structural and pathological aspects - as well as the autonomic control over the SAN - are taken into consideration to reach a holistic view of SAN activity.

MicroRNA-dependent suppression of biological pacemaker activity induced by TBX18

Chemically modified mRNA (CMmRNA) with selectively altered nucleotides are used to deliver transgenes, but translation efficiency is variable. We have transfected CMmRNA encoding human T-box transcription factor 18 (CMmTBX18) into heart cells or the left ventricle of rats with atrioventricular block. TBX18 protein expression from CMmTBX18 is weak and transient, but Acriflavine, an Argonaute 2 inhibitor, boosts TBX18 levels. Small RNA sequencing identified two upregulated microRNAs (miRs) in CMmTBX18-transfected cells. Co-administration of miR-1-3p and miR-1b antagomiRs with CMmTBX18 prolongs TBX18 expression in vitro and in vivo and is sufficient to generate electrical stimuli capable of pacing the heart. Different suppressive miRs likewise limit the expression of VEGF-A CMmRNA. Cells therefore resist translation of CMmRNA therapeutic transgenes by upregulating suppressive miRs. Blockade of suppressive miRs enhances CMmRNA expression of genes driving biological pacing or angiogenesis. Such counterstrategies constitute an approach to boost the efficacy and efficiency of CMmRNA therapies.

The funny current If is essential for the fight-or-flight response in cardiac pacemaker cells

The sympathetic nervous system fight-or-flight response is characterized by a rapid increase in heart rate, which is mediated by an increase in the spontaneous action potential (AP) firing rate of pacemaker cells in the sinoatrial node. Sympathetic neurons stimulate sinoatrial myocytes (SAMs) by activating β adrenergic receptors (βARs) and increasing cAMP. The funny current (If) is among the cAMP-sensitive currents in SAMs. If is critical for pacemaker activity, however, its role in the fight-or-flight response remains controversial. In this study, we used AP waveform analysis, machine learning, and dynamic clamp experiments in acutely isolated SAMs from mice to quantitatively define the AP waveform changes and role of If in the fight-or-flight increase in AP firing rate. We found that while βAR stimulation significantly altered nearly all AP waveform parameters, the increase in firing rate was only correlated with changes in a subset of parameters (diastolic duration, late AP duration, and diastolic depolarization rate). Dynamic clamp injection of the βAR-sensitive component of If showed that it accounts for ∼41% of the fight-or-flight increase in AP firing rate and 60% of the decrease in the interval between APs. Thus, If is an essential contributor to the fight-or-flight increase in heart rate.

Frequency-Dependent Properties of the Hyperpolarization-Activated Cation Current, If, in Adult Mouse Heart Primary Pacemaker Myocytes

A number of distinct electrophysiological mechanisms that modulate the myogenic spontaneous pacemaker activity in the sinoatrial node (SAN) of the mammalian heart have been investigated extensively. There is agreement that several (3 or 4) different transmembrane ionic current changes (referred to as the voltage clock) are involved; and that the resulting net current interacts with direct and indirect effects of changes in intracellular Ca²⁺ (the calcium clock). However, significant uncertainties, and important knowledge gaps, remain concerning the functional roles in SAN spontaneous pacing of many of the individual ion channel- or exchanger-mediated transmembrane current changes. We report results from patch clamp studies and mathematical modeling of the hyperpolarization-activated current, I_f, in the generation/modulation of the diastolic depolarization, or pacemaker potential, produced by individual myocytes that were enzymatically isolated from the adult mouse sinoatrial node (SAN). Amphotericin-mediated patch microelectrode recordings at 35 °C were made under control conditions and in the presence of 5 or 10 nM isoproterenol (ISO). These sets of results were complemented and integrated with mathematical modeling of the current changes that take place in the range of membrane potentials (−70 to −50 mV), which corresponds to the ‘pacemaker depolarization’ in the adult mouse SAN. Our results reveal a very small, but functionally important, approximately steady-state or time-independent current generated by residual activation of I_f channels that are expressed in these pacemaker myocytes. Recordings of the pacemaker depolarization and action potential, combined with measurements of changes in I_f, and the well-known increases in the L-type Ca²⁺ current, I_CaL, demonstrated that I_CaL activation, is essential for myogenic pacing. Moreover, after being enhanced (approximately 3-fold) by 5 or 10 nM ISO, I_CaL contributes significantly to the positive chronotropic effect. Our mathematical model has been developed in an attempt to better understand the underlying mechanisms for the pacemaker depolarization and action potential in adult mouse SAN myocytes. After being updated with our new experimental data describing I_f, our simulations reveal a novel functional component of I_f in adult mouse SAN. Computational work carried out with this model also confirms that in the presence of ISO the residual activation of I_f and opening of I_CaL channels combine to generate a net current change during the slow diastolic depolarization phase that is essential for the observed accelerated pacemaking rate of these SAN myocytes.

Hypoglycaemia combined with mild hypokalaemia reduces the heart rate and causes abnormal pacemaker activity in a computational model of a human sinoatrial cell

Low blood glucose, hypoglycaemia, has been implicated as a possible contributing factor to sudden cardiac death (SCD) in people with diabetes but it is challenging to investigate in clinical studies. We hypothesized the effects of hypoglycaemia on the sinoatrial node (SAN) in the heart to be a candidate mechanism and adapted a computational model of the human SAN action potential developed by Fabbri et al., to investigate the effects of hypoglycaemia on the pacemaker rate. Using Latin hypercube sampling, we combined the effects of low glucose (LG) on the human ether-a-go-go-related gene channel with reduced blood potassium, hypokalaemia, and added sympathetic and parasympathetic stimulus. We showed that hypoglycaemia on its own causes a small decrease in heart rate but there was also a marked decrease in heart rate when combined with hypokalaemia. The effect of the sympathetic stimulus was diminished, causing a smaller increase in heart rate, with LG and hypokalaemia compared to normoglycaemia. By contrast, the effect of the parasympathetic stimulus was enhanced, causing a greater decrease in heart rate. We therefore demonstrate a potential mechanistic explanation for hypoglycaemia-induced bradycardia and show that sinus arrest is a plausible mechanism for SCD in people with diabetes.

An OpenCV-Based Approach for Automated Cardiac Rhythm Measurement in Zebrafish from Video Datasets

Cardiac arrhythmia has been defined as one of the abnormal heart rhythm symptoms, which is a common problem dealt with by cardiologists. Zebrafish were established as a powerful animal model with a transparent body that enables optical observation to analyze cardiac morphology and cardiac rhythm regularity. Currently, research has observed heart-related parameters in zebrafish, which used different approaches, such as starting from the use of fluorescent transgenic zebrafish, different software, and different observation methods. In this study, we developed an innovative approach by using the OpenCV library to measure zebrafish larvae heart rate and rhythm. The program is designed in Python, with the feature of multiprocessing for simultaneous region-of-interest (ROI) detection, covering both the atrium and ventricle regions in the video, and was designed to be simple and user-friendly, having utility even for users who are unfamiliar with Python. Results were validated with our previously published method using ImageJ, which observes pixel changes. In summary, the results showed good consistency in heart rate-related parameters. In addition, the established method in this study also can be widely applied to other invertebrates (like Daphnia) for cardiac rhythm measurement.

Intracellular Na+ Modulates Pacemaking Activity in Murine Sinoatrial Node Myocytes: An In Silico Analysis

Background: The mechanisms underlying dysfunction in the sinoatrial node (SAN), the heart's primary pacemaker, are incompletely understood. Electrical and Ca2+-handling remodeling have been implicated in SAN dysfunction associated with heart failure, aging, and diabetes. Cardiomyocyte [Na+]i is also elevated in these diseases, where it contributes to arrhythmogenesis. Here, we sought to investigate the largely unexplored role of Na+ homeostasis in SAN pacemaking and test whether [Na+]i dysregulation may contribute to SAN dysfunction. Methods: We developed a dataset-specific computational model of the murine SAN myocyte and simulated alterations in the major processes of Na+ entry (Na+/Ca2+ exchanger, NCX) and removal (Na+/K+ ATPase, NKA). Results: We found that changes in intracellular Na+ homeostatic processes dynamically regulate SAN electrophysiology. Mild reductions in NKA and NCX function increase myocyte firing rate, whereas a stronger reduction causes bursting activity and loss of automaticity. These pathologic phenotypes mimic those observed experimentally in NCX- and ankyrin-B-deficient mice due to altered feedback between the Ca2+ and membrane potential clocks underlying SAN firing. Conclusions: Our study generates new testable predictions and insight linking Na+ homeostasis to Ca2+ handling and membrane potential dynamics in SAN myocytes that may advance our understanding of SAN (dys)function.

Cardiac Pacemaker Activity and Aging

A progressive decline in maximum heart rate (mHR) is a fundamental aspect of aging in humans and other mammals. This decrease in mHR is independent of gender, fitness, and lifestyle, affecting in equal measure women and men, athletes and couch potatoes, spinach eaters and fast food enthusiasts. Importantly, the decline in mHR is the major determinant of the age-dependent decline in aerobic capacity that ultimately limits functional independence for many older individuals. The gradual reduction in mHR with age reflects a slowing of the intrinsic pacemaker activity of the sinoatrial node of the heart, which results from electrical remodeling of individual pacemaker cells along with structural remodeling and a blunted β-adrenergic response. In this review, we summarize current evidence about the tissue, cellular, and molecular mechanisms that underlie the reduction in pacemaker activity with age and highlight key areas for future work.

ParamAP: Standardized Parameterization of Sinoatrial Node Myocyte Action Potentials

Sinoatrial node myocytes act as cardiac pacemaker cells by generating spontaneous action potentials (APs). Much information is encoded in sinoatrial AP waveforms, but both the analysis and the comparison of AP parameters between studies is hindered by the lack of standardized parameter definitions and the absence of automated analysis tools. Here we introduce ParamAP, a standalone cross-platform computational tool that uses a template-free detection algorithm to automatically identify and parameterize APs from text input files. ParamAP employs a graphic user interface with automatic and user-customizable input modes, and it outputs data files in text and PDF formats. ParamAP returns a total of 16 AP waveform parameters including time intervals such as the AP duration, membrane potentials such as the maximum diastolic potential, and rates of change of the membrane potential such as the diastolic depolarization rate. ParamAP provides a robust AP detection algorithm in combination with a standardized AP parameter analysis over a wide range of AP waveforms and firing rates, owing in part to the use of an iterative algorithm for the determination of the threshold potential and the diastolic depolarization rate that is independent of the maximum upstroke velocity, a parameter that can vary significantly among sinoatrial APs. Because ParamAP is implemented in Python 3, it is also highly customizable and extensible. In conclusion, ParamAP is a powerful computational tool that facilitates quantitative analysis and enables comparison of sinoatrial APs by standardizing parameter definitions and providing an automated work flow.