Required GK1 to Suppress Automaticity of iPSC-CMs Depends Strongly on IK1 Model Structure.

Human-induced pluripotent stem cells derived cardiomyocytes (hiPSC-CMs) are a virtually endless source of human cardiomyocytes that may become a great tool for safety pharmacology; however, their electrical phenotype is immature: they show spontaneous action potentials (APs) and an unstable and depolarized resting membrane potential (RMP) because of lack of IK1. Such immaturity hampers their application in assessing drug safety. The electronic overexpression of IK1 (e.g., through the dynamic clamp (DC) technique) is an option to overcome this deficit. In this computational study, we aim to estimate how much IK1 is needed to bring hiPSC-CMs to a stable and hyperpolarized RMP and which mathematical description of IK1 is most suitable for DC experiments. We compared five mature IK1 formulations (Bett, Dhamoon, Ishihara, O'Hara-Rudy, and ten Tusscher) with the native one (Paci), evaluating the main properties (outward peak, degree of rectification), and we quantified their effects on AP features (RMP, V˙max, APD50, APD90 (AP duration at 50 and 90% of repolarization), and APD50/APD90) after including them in the hiPSC-CM mathematical model by Paci. Then, we automatically identified the critical conductance for IK1 ( GK1, critical), the minimally required amount of IK1 suppressing spontaneous activity. Preconditioning the cell model with depolarizing/hyperpolarizing prepulses allowed us to highlight time dependency of the IK1 formulations. Simulations showed that inclusion of mature IK1 formulations resulted in hyperpolarized RMP and higher V˙max, and observed GK1, critical and the effect on AP duration strongly depended on IK1 formulation. Finally, the Ishihara IK1 led to shorter (-16.3%) and prolonged (+6.5%) APD90 in response to hyperpolarizing and depolarizing prepulses, respectively, whereas other models showed negligible effects. Fine-tuning of GK1 is an important step in DC experiments. Our computational work proposes a procedure to automatically identify how much IK1 current is required to inject to stop the spontaneous activity and suggests the use of the Ishihara IK1 model to perform DC experiments in hiPSC-CMs.

Hypocalcemia-Induced Slowing of Human Sinus Node Pacemaking.

Each heartbeat is initiated by cyclic spontaneous depolarization of cardiomyocytes in the sinus node forming the primary natural pacemaker. In patients with end-stage renal disease undergoing hemodialysis, it was recently shown that the heart rate drops to very low values before they suffer from sudden cardiac death with an unexplained high incidence. We hypothesize that the electrolyte changes commonly occurring in these patients affect sinus node beating rate and could be responsible for severe bradycardia. To test this hypothesis, we extended the Fabbri et al. computational model of human sinus node cells to account for the dynamic intracellular balance of ion concentrations. Using this model, we systematically tested the effect of altered extracellular potassium, calcium, and sodium concentrations. Although sodium changes had negligible (0.15 bpm/mM) and potassium changes mild effects (8 bpm/mM), calcium changes markedly affected the beating rate (46 bpm/mM ionized calcium without autonomic control). This pronounced bradycardic effect of hypocalcemia was mediated primarily by ICaL attenuation due to reduced driving force, particularly during late depolarization. This, in turn, caused secondary reduction of calcium concentration in the intracellular compartments and subsequent attenuation of inward INaCa and reduction of intracellular sodium. Our in silico findings are complemented and substantiated by an empirical database study comprising 22,501 pairs of blood samples and in vivo heart rate measurements in hemodialysis patients and healthy individuals. A reduction of extracellular calcium was correlated with a decrease of heartrate by 9.9 bpm/mM total serum calcium (p < 0.001) with intact autonomic control in the cross-sectional population. In conclusion, we present mechanistic in silico and empirical in vivo data supporting the so far neglected but experimentally testable and potentially important mechanism of hypocalcemia-induced bradycardia and asystole, potentially responsible for the highly increased and so far unexplained risk of sudden cardiac death in the hemodialysis patient population.

Computational analysis of the human sinus node action potential: model development and effects of mutations.

KEY POINTS: We constructed a comprehensive mathematical model of the spontaneous electrical activity of a human sinoatrial node (SAN) pacemaker cell, starting from the recent Severi-DiFrancesco model of rabbit SAN cells. Our model is based on electrophysiological data from isolated human SAN pacemaker cells and closely matches the action potentials and calcium transient that were recorded experimentally. Simulated ion channelopathies explain the clinically observed changes in heart rate in corresponding mutation carriers, providing an independent qualitative validation of the model. The model shows that the modulatory role of the 'funny current' (If ) in the pacing rate of human SAN pacemaker cells is highly similar to that of rabbit SAN cells, despite its considerably lower amplitude. The model may prove useful in the design of experiments and the development of heart-rate modulating drugs. ABSTRACT: The sinoatrial node (SAN) is the normal pacemaker of the mammalian heart.  Over several decades, a large amount of data on the ionic mechanisms underlying the spontaneous electrical activity of SAN pacemaker cells has been obtained, mostly in experiments on single cells isolated from rabbit SAN. This wealth of data has allowed the development of mathematical models of the electrical activity of rabbit SAN pacemaker cells. The present study aimed to construct a comprehensive model of the electrical activity of a human SAN pacemaker cell using recently obtained electrophysiological data from human SAN pacemaker cells.  We based our model on the recent Severi-DiFrancesco model of a rabbit SAN pacemaker cell. The action potential and calcium transient of the resulting model are close to the experimentally recorded values. The model has a much smaller 'funny current' (If ) than do rabbit cells, although its modulatory role is highly similar. Changes in pacing rate upon the implementation of mutations associated with sinus node dysfunction agree with the clinical observations. This agreement holds for both loss-of-function and gain-of-function mutations in the HCN4, SCN5A and KCNQ1 genes, underlying ion channelopathies in If , fast sodium current and slow delayed rectifier potassium current, respectively. We conclude that our human SAN cell model can be a useful tool in the design of experiments and the development of drugs that aim to modulate heart rate.

Automated Dynamic Clamp for Simulation of IK1 in Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes in Real Time Using Patchliner Dynamite8.

Current in vitro assays typically poorly predict cardiac liability as they focus on single ion channels overexpressed in cell lines. Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), on the other hand, provide a unique opportunity for drug testing on human cardiomyocytes using high-throughput systems. However, these cells can differ from adult cardiomyocytes in their ion channel expression and, therefore, electrophysiologic properties. One of the main challenges of hiPSC-CMs is the physiologic expression of ion channels such as the inward rectifiers (e.g., Kir2.1-2.3), which conduct the cardiac inward rectifier potassium current (IK1 ). IK1 is one of the primary contributors in maintaining a stable resting membrane potential in cardiac cells, which is essential for excitability. This is only expressed in low levels, or sometimes not at all, in hiPSC-CMs as shown by patch clamp studies. Dynamic clamp is a method of electronically introducing ion currents (e.g., IK1 ) into cells to compensate for the lack of endogenous expression, thus offering the potential to record more stable action potentials in hiPSC-CMs. In this article, we describe the method of using hiPSC-CMs on an automated patch clamp device (Patchliner) coupled with the automated dynamic clamp add-on (Dynamite8 ). We describe protocols for optimized cell handling and harvesting for use on the Patchliner and the steps required for automated execution of experiments and data analysis in dynamic clamp mode. © 2019 by John Wiley & Sons, Inc. Basic Protocol: Recording action potential pharmacology from human induced pluripotent stem cell-derived cardiomyocytes in automated patch clamp combined with dynamic clamp to introduce simulated IK1 and compensate seal resistance Support Protocol 1: Cardiomyocyte plating and culture Support Protocol 2: Cell harvesting and dissociation Alternate Protocol: Recording action potential pharmacology at physiologic temperatures.

A nonlinear and time-dependent leak current in the presence of calcium fluoride patch-clamp seal enhancer.

Automated patch-clamp platforms are widely used and vital tools in both academia and industry to enable high-throughput studies such as drug screening. A leak current to ground occurs whenever the seal between a pipette and cell (or internal solution and cell in high-throughput machines) is not perfectly insulated from the bath (extracellular) solution. Over 1 GΩ seal resistance between pipette and bath solutions is commonly used as a quality standard for manual patch work. With automated platforms it can be difficult to obtain such a high seal resistance between the intra- and extra-cellular solutions. One suggested method to alleviate this problem is using an F - containing internal solution together with a Ca 2+ containing external solution - so that a CaF 2 crystal forms when the two solutions meet which 'plugs the holes' to enhance the seal resistance. However, we observed an unexpected nonlinear-in-voltage and time-dependent current using these solutions on an automated patch-clamp platform. We performed manual patch-clamp experiments with the automated patch-clamp solutions, but no biological cell, and observed the same nonlinear time-dependent leak current. The current could be completely removed by washing out F - ions to leave a conventional leak current that was linear and not time-dependent. We therefore conclude fluoride ions interacting with the CaF 2 crystal are the origin of the nonlinear time-dependent leak current. The consequences of such a nonlinear and time-dependent leak current polluting measurements should be considered carefully if it cannot be isolated and subtracted.

Inter-Species Differences in the Response of Sinus Node Cellular Pacemaking to Changes of Extracellular Calcium.

Changes of serum and extracellular ion concentrations occur regularly in patients with chronic kidney disease (CKD). Recently, hypocalcaemia, i.e. a decrease of the extracellular calcium concentration [Ca2+]o, has been suggested as potential pathomechanism contributing to the unexplained high rate of sudden cardiac death (SCD) in CKD patients. In particular, there is a hypothesis that hypocalcaemia could slow down natural pacemaking in the human sinus node to fatal degrees. Here, we address the question whether there are inter-species differences in the response of cellular sinus node pacemaking to changes of [Ca2+]o. Towards this end, we employ computational models of mouse, rabbit and human sinus node cells. The Fabbri et al. human model was updated to consider changes of intracellular ion concentrations. We identified crucial inter-species differences in the response of cellular pacemaking in the sinus node to changes of [Ca2+]o with little changes of cycle length in mouse and rabbit models (<; 83 ms) in contrast to a pronounced bradycardic effect in the human model (up to >1000 ms). Our results suggest that experiments with human sinus node cells are required to investigate the potential mechanism of hypocalcaemia-induced bradycardic SCD in CKD patients and small animal models are not well suited.