The zebrafish grime mutant uncovers an evolutionarily conserved role for Tmem161b in the control of cardiac rhythm.

The establishment of cardiac function in the developing embryo is essential to ensure blood flow and, therefore, growth and survival of the animal. The molecular mechanisms controlling normal cardiac rhythm remain to be fully elucidated. From a forward genetic screen, we identified a unique mutant, grime, that displayed a specific cardiac arrhythmia phenotype. We show that loss-of-function mutations in tmem161b are responsible for the phenotype, identifying Tmem161b as a regulator of cardiac rhythm in zebrafish. To examine the evolutionary conservation of this function, we generated knockout mice for Tmem161b. Tmem161b knockout mice are neonatal lethal and cardiomyocytes exhibit arrhythmic calcium oscillations. Mechanistically, we find that Tmem161b is expressed at the cell membrane of excitable cells and live imaging shows it is required for action potential repolarization in the developing heart. Electrophysiology on isolated cardiomyocytes demonstrates that Tmem161b is essential to inhibit Ca2+ and K+ currents in cardiomyocytes. Importantly, Tmem161b haploinsufficiency leads to cardiac rhythm phenotypes, implicating it as a candidate gene in heritable cardiac arrhythmia. Overall, these data describe Tmem161b as a highly conserved regulator of cardiac rhythm that functions to modulate ion channel activity in zebrafish and mice.

Leak current, even with gigaohm seals, can cause misinterpretation of stem cell-derived cardiomyocyte action potential recordings.

AIMS: Human-induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) have become an essential tool to study arrhythmia mechanisms. Much of the foundational work on these cells, as well as the computational models built from the resultant data, has overlooked the contribution of seal-leak current on the immature and heterogeneous phenotype that has come to define these cells. The aim of this study is to understand the effect of seal-leak current on recordings of action potential (AP) morphology. METHODS AND RESULTS: Action potentials were recorded in human iPSC-CMs using patch clamp and simulated using previously published mathematical models. Our in silico and in vitro studies demonstrate how seal-leak current depolarizes APs, substantially affecting their morphology, even with seal resistances (Rseal) above 1 GΩ. We show that compensation of this leak current is difficult due to challenges with obtaining accurate measures of Rseal during an experiment. Using simulation, we show that Rseal measures (i) change during an experiment, invalidating the use of pre-rupture values, and (ii) are polluted by the presence of transmembrane currents at every voltage. Finally, we posit that the background sodium current in baseline iPSC-CM models imitates the effects of seal-leak current and is increased to a level that masks the effects of seal-leak current on iPSC-CMs. CONCLUSION: Based on these findings, we make recommendations to improve iPSC-CM AP data acquisition, interpretation, and model-building. Taking these recommendations into account will improve our understanding of iPSC-CM physiology and the descriptive ability of models built from such data.

Is zebrafish heart regeneration "complete"? Lineage-restricted cardiomyocytes proliferate to pre-injury numbers but some fail to differentiate in fibrotic hearts.

Adult zebrafish are frequently described to be able to "completely" regenerate the heart. Yet, the extent to which cardiomyocytes lost to injury are replaced is unknown, since existing evidence for cardiomyocyte proliferation is indirect or non-quantitative. We established stereological methods to quantify the number of cardiomyocytes at several time-points post cryoinjury. Intriguingly, after cryoinjuries that killed about 1/3 of the ventricular cardiomyocytes, pre-injury cardiomyocyte numbers were restored already within 30 days. Yet, many hearts retained small residual scars, and a subset of cardiomyocytes bordering these fibrotic areas remained smaller, lacked differentiated sarcomeric structures, and displayed defective calcium signaling. Thus, a subset of regenerated cardiomyocytes failed to fully mature. While lineage-tracing experiments have shown that regenerating cardiomyocytes are derived from differentiated cardiomyocytes, technical limitations have previously made it impossible to test whether cardiomyocyte trans-differentiation contributes to regeneration of non-myocyte cell lineages. Using Cre responder lines that are expressed in all major cell types of the heart, we found no evidence for cardiomyocyte transdifferentiation into endothelial, epicardial, fibroblast or immune cell lineages. Overall, our results imply a refined answer to the question whether zebrafish can completely regenerate the heart: in response to cryoinjury, preinjury cardiomyocyte numbers are indeed completely regenerated by proliferation of lineage-restricted cardiomyocytes, while restoration of cardiomyocyte differentiation and function, as well as resorption of scar tissue, is less robustly achieved.

Accounting for variability in ion current recordings using a mathematical model of artefacts in voltage-clamp experiments.

Mathematical models of ion channels, which constitute indispensable components of action potential models, are commonly constructed by fitting to whole-cell patch-clamp data. In a previous study, we fitted cell-specific models to hERG1a (Kv11.1) recordings simultaneously measured using an automated high-throughput system, and studied cell-cell variability by inspecting the resulting model parameters. However, the origin of the observed variability was not identified. Here, we study the source of variability by constructing a model that describes not just ion current dynamics, but the entire voltage-clamp experiment. The experimental artefact components of the model include: series resistance, membrane and pipette capacitance, voltage offsets, imperfect compensations made by the amplifier for these phenomena, and leak current. In this model, variability in the observations can be explained by either cell properties, measurement artefacts, or both. Remarkably, by assuming that variability arises exclusively from measurement artefacts, it is possible to explain a larger amount of the observed variability than when assuming cell-specific ion current kinetics. This assumption also leads to a smaller number of model parameters. This result suggests that most of the observed variability in patch-clamp data measured under the same conditions is caused by experimental artefacts, and hence can be compensated for in post-processing by using our model for the patch-clamp experiment. This study has implications for the question of the extent to which cell-cell variability in ion channel kinetics exists, and opens up routes for better correction of artefacts in patch-clamp data. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

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.

Flotillins in the intercalated disc are potential modulators of cardiac excitability.

BACKGROUND: The intercalated disc (ID) is important for cardiac remodeling and has become a subject of intensive research efforts. However, as yet the composition of the ID has still not been conclusively resolved and the role of many proteins identified in the ID, like Flotillin-2, is often unknown. The Flotillin proteins are known to be involved in the stabilization of cadherins and desmosomes in the epidermis and upon cancer development. However, their role in the heart has so far not been investigated. Therefore, in this study, we aimed at identifying the role of Flotillin-1 and Flotillin-2 in the cardiac ID. METHODS: Location of Flotillins in human and murine cardiac tissue was evaluated by fluorescent immunolabeling and co-immunoprecipitation. In addition, the effect of Flotillin knockout (KO) on proteins of the ID and in electrical excitation and conduction was investigated in cardiac samples of wildtype (WT), Flotillin-1 KO, Flotilin-2 KO and Flotilin-1/2 double KO mice. Consequences of Flotillin knockdown (KD) on cardiac function were studied (patch clamp and Multi Electrode Array (MEA)) in neonatal rat cardiomyocytes (NRCMs) transfected with siRNAs against Flotillin-1 and/or Flotillin-2. RESULTS: First, we confirmed presence in the ID and mutual binding of Flotillin-1 and Flotillin-2 in murine and human cardiac tissue. Flotillin KO mice did not show cardiac fibrosis, nor hypertrophy or changes in expression of the desmosomal ID proteins. However, protein expression of the cardiac sodium channel NaV1.5 was significantly decreased in Flotillin-1 and Flotillin-1/2 KO mice compared to WT mice. In addition, sodium current density showed a significant decrease upon Flotillin-1/2 KD in NRCMs as compared to scrambled siRNA-transfected NRCMs. MEA recordings of Flotillin-2 KD NRCM cultures showed a significantly decreased spike amplitude and a tendency of a reduced spike slope when compared to control and scrambled siRNA-transfected cultures. CONCLUSIONS: In this study, we demonstrate the presence of Flotillin-1, in addition to Flotillin-2 in the cardiac ID. Our findings indicate a modulatory role of Flotillins on NaV1.5 expression at the ID, with potential consequences for cardiac excitation.

Sinusoidal voltage protocols for rapid characterisation of ion channel kinetics.

KEY POINTS: Ion current kinetics are commonly represented by current-voltage relationships, time constant-voltage relationships and subsequently mathematical models fitted to these. These experiments take substantial time, which means they are rarely performed in the same cell. Rather than traditional square-wave voltage clamps, we fitted a model to the current evoked by a novel sum-of-sinusoids voltage clamp that was only 8 s long. Short protocols that can be performed multiple times within a single cell will offer many new opportunities to measure how ion current kinetics are affected by changing conditions. The new model predicts the current under traditional square-wave protocols well, with better predictions of underlying currents than literature models. The current under a novel physiologically relevant series of action potential clamps is predicted extremely well. The short sinusoidal protocols allow a model to be fully fitted to individual cells, allowing us to examine cell-cell variability in current kinetics for the first time. ABSTRACT: Understanding the roles of ion currents is crucial to predict the action of pharmaceuticals and mutations in different scenarios, and thereby to guide clinical interventions in the heart, brain and other electrophysiological systems. Our ability to predict how ion currents contribute to cellular electrophysiology is in turn critically dependent on our characterisation of ion channel kinetics - the voltage-dependent rates of transition between open, closed and inactivated channel states. We present a new method for rapidly exploring and characterising ion channel kinetics, applying it to the hERG potassium channel as an example, with the aim of generating a quantitatively predictive representation of the ion current. We fitted a mathematical model to currents evoked by a novel 8 second sinusoidal voltage clamp in CHO cells overexpressing hERG1a. The model was then used to predict over 5 minutes of recordings in the same cell in response to further protocols: a series of traditional square step voltage clamps, and also a novel voltage clamp comprising a collection of physiologically relevant action potentials. We demonstrate that we can make predictive cell-specific models that outperform the use of averaged data from a number of different cells, and thereby examine which changes in gating are responsible for cell-cell variability in current kinetics. Our technique allows rapid collection of consistent and high quality data, from single cells, and produces more predictive mathematical ion channel models than traditional approaches.

Cardiac optogenetics: using light to monitor cardiac physiology.

Our current understanding of cardiac excitation and its coupling to contraction is largely based on ex vivo studies utilising fluorescent organic dyes to assess cardiac action potentials and signal transduction. Recent advances in optogenetic sensors open exciting new possibilities for cardiac research and allow us to answer research questions that cannot be addressed using the classic organic dyes. Especially thrilling is the possibility to use optogenetic sensors to record parameters of cardiac excitation and contraction in vivo. In addition, optogenetics provide a high spatial resolution, as sensors can be coupled to motifs and targeted to specific cell types and subcellular domains of the heart. In this review, we will give a comprehensive overview of relevant optogenetic sensors, how they can be utilised in cardiac research and how they have been applied in cardiac research up to now.

Sensing Cardiac Electrical Activity With a Cardiac Myocyte--Targeted Optogenetic Voltage Indicator.

RATIONALE: Monitoring and controlling cardiac myocyte activity with optogenetic tools offer exciting possibilities for fundamental and translational cardiovascular research. Genetically encoded voltage indicators may be particularly attractive for minimal invasive and repeated assessments of cardiac excitation from the cellular to the whole heart level. OBJECTIVE: To test the hypothesis that cardiac myocyte-targeted voltage-sensitive fluorescence protein 2.3 (VSFP2.3) can be exploited as optogenetic tool for the monitoring of electric activity in isolated cardiac myocytes and the whole heart as well as function and maturity in induced pluripotent stem cell-derived cardiac myocytes. METHODS AND RESULTS: We first generated mice with cardiac myocyte-restricted expression of VSFP2.3 and demonstrated distinct localization of VSFP2.3 at the t-tubulus/junctional sarcoplasmic reticulum microdomain without any signs for associated pathologies (assessed by echocardiography, RNA-sequencing, and patch clamping). Optically recorded VSFP2.3 signals correlated well with membrane voltage measured simultaneously by patch clamping. The use of VSFP2.3 for human action potential recordings was confirmed by simulation of immature and mature action potentials in murine VSFP2.3 cardiac myocytes. Optical cardiograms could be monitored in whole hearts ex vivo and minimally invasively in vivo via fiber optics at physiological heart rate (10 Hz) and under pacing-induced arrhythmia. Finally, we reprogrammed tail-tip fibroblasts from transgenic mice and used the VSFP2.3 sensor for benchmarking functional and structural maturation in induced pluripotent stem cell-derived cardiac myocytes. CONCLUSIONS: We introduce a novel transgenic voltage-sensor model as a new method in cardiovascular research and provide proof of concept for its use in optogenetic sensing of physiological and pathological excitation in mature and immature cardiac myocytes in vitro and in vivo.

The concept of triple wavefront fusion during biventricular pacing: Using the EGM to produce the best acute hemodynamic improvement in CRT.

BACKGROUND: Previous reports suggest that biventricular pacing (BiVp) fused with intrinsic conduction (BiVp-fusion, triple wavefront fusion) is associated with improved resynchronization compared to pure-BiVp in cardiac resynchronization therapy (CRT). This study aimed to assess the association between acute hemodynamic benefit of CRT and signs of BiVp-fusion by using a novel electrogram (EGM)-based method. METHODS: In 17 patients undergoing CRT implantation, 28 combinations of atrioventricular (AV) and interventricular (VV) delays were applied while invasively measuring acute hemodynamic response based on maximum rate of left ventricular (LV) pressure rise (LV dP/dtmax ) to assess optimal BiVp settings. BiVp-fusion was noted if farfield signal (caused by first intrinsic ventricular depolarization) was seen prior to right ventricular (RV) pacing (RVp) artifact on integrated bipolar RV EGM, or QRS morphology changed compared to pure-BiVp (short AV-delay) as seen on electrocardiogram (ECG). RESULTS: Mean optimal RVp timing was at 98 ± 17% of intrinsic right atrial (RA)-RVfarfield (interval from right atrial pace or sense to RV farfield signal) interval, while preactivating the LV at 50 ± 11% of RA-RVsense (interval from right atrial pace or sense to RV sense interval) interval. BiVp-fusion was noted in 16 of 17 (94%) patients on ECG during optimal BiVp. Eight of these patients showed intrinsic farfield signal prior to RVp artifact on RV EGM. In the remaining eight, the RVp was paced just within the RA-RVfarfield interval with a mean of 25 ± 14 ms prior to the onset; therefore, the intrinsic farfield was masked. CONCLUSION: Optimal hemodynamic BiVp facilitates triple wavefront fusion, by pacing the RV around the onset of intrinsic farfield signal on RV EGM, while preactivating the LV. Aiming at BiVp-fusion could be a target for noninvasive EGM-based CRT device setting optimization.

Uncertainty and variability in models of the cardiac action potential: Can we build trustworthy models?

Cardiac electrophysiology models have been developed for over 50years, and now include detailed descriptions of individual ion currents and sub-cellular calcium handling. It is commonly accepted that there are many uncertainties in these systems, with quantities such as ion channel kinetics or expression levels being difficult to measure or variable between samples. Until recently, the original approach of describing model parameters using single values has been retained, and consequently the majority of mathematical models in use today provide point predictions, with no associated uncertainty. In recent years, statistical techniques have been developed and applied in many scientific areas to capture uncertainties in the quantities that determine model behaviour, and to provide a distribution of predictions which accounts for this uncertainty. In this paper we discuss this concept, which is termed uncertainty quantification, and consider how it might be applied to cardiac electrophysiology models. We present two case studies in which probability distributions, instead of individual numbers, are inferred from data to describe quantities such as maximal current densities. Then we show how these probabilistic representations of model parameters enable probabilities to be placed on predicted behaviours. We demonstrate how changes in these probability distributions across data sets offer insight into which currents cause beat-to-beat variability in canine APs. We conclude with a discussion of the challenges that this approach entails, and how it provides opportunities to improve our understanding of electrophysiology.

Perpetuation of torsade de pointes in heterogeneous hearts: competing foci or re-entry?

KEY POINTS: The underlying mechanism of torsade de pointes (TdP) remains of debate: perpetuation may be due to (1) focal activity or (2) re-entrant activity. The onset of TdP correlates with action potential heterogeneities in different regions of the heart. We studied the mechanism of perpetuation of TdP in silico using a 2D model of human cardiac tissue and an anatomically accurate model of the ventricles of the human heart. We found that the mechanism of perpetuation TdP depends on the degree of heterogeneity. If the degree of heterogeneity is large, focal activity alone can sustain a TdP, otherwise re-entrant activity emerges. This result can help to understand the relationship between the mechanisms of TdP and tissue properties and may help in developing new drugs against it. ABSTRACT: Torsade de pointes (TdP) can be the consequence of cardiac remodelling, drug effects or a combination of both. The mechanism underlying TdP is unclear, and may involve triggered focal activity or re-entry. Recent work by our group has indicated that both cases may exist, i.e. TdPs induced in the chronic atrioventricular block (CAVB) dog model may have a focal origin or are due to re-entry. Also it was found that heterogeneities might play an important role. In the current study we have used computational modelling to further investigate the mechanisms involved in TdP initiation and perpetuation, especially in the CAVB dog model, by the addition of heterogeneities with reduced repolarization reserve in comparison with the surrounding tissue. For this, the TNNP computer model was used for computations. We demonstrated in 2D and 3D simulations that ECGs with the typical TdP morphology can be caused by both multiple competing foci and re-entry circuits as a result of introduction of heterogeneities, depending on whether the heterogeneities have a large or a smaller reduced repolarization reserve in comparison with the surrounding tissue. Large heterogeneities can produce ectopic TdP, while smaller heterogeneities will produce re-entry-type TdP.

Recording ten-fold larger IKr conductances with automated patch clamping using equimolar Cs+ solutions.

Background: The rapid delayed rectifier potassium current (IKr) is important for cardiac repolarization and is most often involved in drug-induced arrhythmias. However, accurately measuring this current can be challenging in human-induced pluripotent stem cell (hiPSC)-derived cardiomyocytes because of its small current density. Interestingly, the ion channel conducting IKr, hERG channel, is not only permeable to K+ ions but also to Cs+ ions when present in equimolar concentrations inside and outside of the cell. Methods: In this study, IhERG was measured from Chinese hamster ovary (CHO)-hERG cells and hiPSC-CM using either Cs+ or K+ as the charge carrier. Equimolar Cs+ has been used in the literature in manual patch-clamp experiments, and here, we apply this approach using automated patch-clamp systems. Four different (pre)clinical drugs were tested to compare their effects on Cs+- and K+-based currents. Results: Using equimolar Cs+ solutions gave rise to approximately ten-fold larger hERG conductances. Comparison of Cs+- and K+-mediated currents upon application of dofetilide, desipramine, moxifloxacin, or LUF7244 revealed many similarities in inhibition or activation properties of the drugs studied. Using equimolar Cs+ solutions gave rise to approximately ten-fold larger hERG conductances. In hiPSC-CM, the Cs+-based conductance is larger compared to the known K+-based conductance, and the Cs+ hERG conductance can be inhibited similarly to the K+-based conductance. Conclusion: Using equimolar Cs+ instead of K+ for IhERG measurements in an automated patch-clamp system gives rise to a new method by which, for example, quick scans can be performed on effects of drugs on hERG currents. This application is specifically relevant when such experiments are performed using cells which express small IKr current densities in combination with small membrane capacitances.

Inflammation-induced mitochondrial and metabolic disturbances in sensory neurons control the switch from acute to chronic pain.

Pain often persists in patients with an inflammatory disease, even when inflammation has subsided. The molecular mechanisms leading to this failure in pain resolution and the transition to chronic pain are poorly understood. Mitochondrial dysfunction in sensory neurons links to chronic pain, but its role in resolution of inflammatory pain is unclear. Transient inflammation causes neuronal plasticity, called hyperalgesic priming, which impairs resolution of pain induced by a subsequent inflammatory stimulus. We identify that hyperalgesic priming in mice increases the expression of a mitochondrial protein (ATPSc-KMT) and causes mitochondrial and metabolic disturbances in sensory neurons. Inhibition of mitochondrial respiration, knockdown of ATPSCKMT expression, or supplementation of the affected metabolite is sufficient to restore resolution of inflammatory pain and prevents chronic pain development. Thus, inflammation-induced mitochondrial-dependent disturbances in sensory neurons predispose to a failure in resolution of inflammatory pain and development of chronic pain.

PITX2 induction leads to impaired cardiomyocyte function in arrhythmogenic cardiomyopathy.

Arrhythmogenic cardiomyopathy (ACM) is an inherited progressive disease characterized by electrophysiological and structural remodeling of the ventricles. However, the disease-causing molecular pathways, as a consequence of desmosomal mutations, are poorly understood. Here, we identified a novel missense mutation within desmoplakin in a patient clinically diagnosed with ACM. Using CRISPR-Cas9, we corrected this mutation in patient-derived human induced pluripotent stem cells (hiPSCs) and generated an independent knockin hiPSC line carrying the same mutation. Mutant cardiomyocytes displayed a decline in connexin 43, NaV1.5, and desmosomal proteins, which was accompanied by a prolonged action potential duration. Interestingly, paired-like homeodomain 2 (PITX2), a transcription factor that acts a repressor of connexin 43, NaV1.5, and desmoplakin, was induced in mutant cardiomyocytes. We validated these results in control cardiomyocytes in which PITX2 was either depleted or overexpressed. Importantly, knockdown of PITX2 in patient-derived cardiomyocytes is sufficient to restore the levels of desmoplakin, connexin 43, and NaV1.5.

Interplay between calcium and sarcomeres directs cardiomyocyte maturation during regeneration.

Zebrafish hearts can regenerate by replacing damaged tissue with new cardiomyocytes. Although the steps leading up to the proliferation of surviving cardiomyocytes have been extensively studied, little is known about the mechanisms that control proliferation and redifferentiation to a mature state. We found that the cardiac dyad, a structure that regulates calcium handling and excitation-contraction coupling, played a key role in the redifferentiation process. A component of the cardiac dyad called leucine-rich repeat-containing 10 (Lrrc10) acted as a negative regulator of proliferation, prevented cardiomegaly, and induced redifferentiation. We found that its function was conserved in mammalian cardiomyocytes. This study highlights the importance of the underlying mechanisms required for heart regeneration and their application to the generation of fully functional cardiomyocytes.

Application of human stem cell-derived cardiomyocytes in safety pharmacology requires caution beyond hERG.

Human embryonic stem cell-derived cardiomyocytes (hESC-CM) have been proposed as a new model for safety pharmacology. So far, a thorough description of their basic electrophysiology and extensive testing, and mechanistic explanations, of their overall pro-arrhythmic ability is lacking. Under standardized conditions, we have evaluated the sensitivity of hESC-CM to proarrhythmic provocations by blockade of hERG and other channels. Using voltage patch clamp, some ion current densities (pA/pF) in hESC-CM were comparable to adult CM: I(Kr) (-12.5 ± 6.9), I(Ks) (0.65 ± 0.12), I(Na,peak) (-72 ± 21), I(Na,late) (-1.10 ± 0.36), and I(Ca,L) (-4.3 ± 0.6). I(f) density was larger (-10 ± 1.1) and I(K1) not existent or very small (-2.67 ± 0.3). The low I(K1) density was corroborated by low KCNJ2 mRNA levels. Effects of pro-arrhythmic compounds on action potential (AP) parameters and provocation of early afterdepolarizations (EADs) revealed that Chromanol293B (100 µmol/l) and Bay K8644 (1 µmol/l) both significantly prolonged APD(90). ATX-II (<1 µmol/l ) and BaCl(2) (10 µmol/l ) had no effect on APD. The only compound that triggered EADs was hERG blocker Cisapride. Computer simulations and AP clamp showed that the immature AP of hESC-CM prevents proper functioning of I(Na)-channels, and result in lower peak/maximal currents of several other channels, compared to the adult situation. Lack of functional I(K1) channels and shifted I(Na) channel activation cause a rather immature electrophysiological phenotype in hESC-CM, and thereby limits the potential of this model to respond accurately to pro-arrhythmic triggers other than hERG block. Maturation of the electrical phenotype is a prerequiste for future implementation of the model in arrhythmogenic safety testing.