Mathematical modelling of the action potential of human embryonic stem cell derived cardiomyocytes

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.

Advances in induced pluripotent stem cell-derived cardiac myocytes: technological breakthroughs, key discoveries and new applications

A transformation is underway in precision and patient-specific medicine. Rapid progress has been enabled by multiple new technologies including induced pluripotent stem cell-derived cardiac myocytes (iPSC-CMs). Here, we delve into these advancements and their future promise, focusing on the efficiency of reprogramming techniques, the fidelity of differentiation into the cardiac lineage, the functional characterization of the resulting cardiac myocytes, and the many applications of in silico models to understand general and patient-specific mechanisms controlling excitation-contraction coupling in health and disease. Furthermore, we explore the current and potential applications of iPSC-CMs in both research and clinical settings, underscoring the far-reaching implications of this rapidly evolving field.

Arrhythmia Mechanisms in Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes

ABSTRACT

Despite major efforts by clinicians and researchers, cardiac arrhythmia remains a leading cause of morbidity and mortality in the world. Experimental work has relied on combining high-throughput strategies with standard molecular and electrophysiological studies, which are, to a great extent, based on the use of animal models. Because this poses major challenges for translation, the progress in the development of novel antiarrhythmic agents and clinical care has been mostly disappointing. Recently, the advent of human induced pluripotent stem cell-derived cardiomyocytes has opened new avenues for both basic cardiac research and drug discovery; now, there is an unlimited source of cardiomyocytes of human origin, both from healthy individuals and patients with cardiac diseases. Understanding arrhythmic mechanisms is one of the main use cases of human induced pluripotent stem cell-derived cardiomyocytes, in addition to pharmacological cardiotoxicity and efficacy testing, in vitro disease modeling, developing patient-specific models and personalized drugs, and regenerative medicine. Here, we review the advances that the human induced pluripotent stem cell-derived-based modeling systems have brought so far regarding the understanding of both arrhythmogenic triggers and substrates, while also briefly speculating about the possibilities in the future.

A computational model of induced pluripotent stem-cell derived cardiomyocytes incorporating experimental variability from multiple data sources

Key points

Induced pluripotent stem cell‐derived cardiomyocytes (iPSC‐CMs) capture patient‐specific genotype–phenotype relationships, as well as cell‐to‐cell variability of cardiac electrical activity

Computational modelling and simulation provide a high throughput approach to reconcile multiple datasets describing physiological variability, and also identify vulnerable parameter regimes

We have developed a whole‐cell model of iPSC‐CMs, composed of single exponential voltage‐dependent gating variable rate constants, parameterized to fit experimental iPSC‐CM outputs

We have utilized experimental data across multiple laboratories to model experimental variability and investigate subcellular phenotypic mechanisms in iPSC‐CMs

This framework links molecular mechanisms to cellular‐level outputs by revealing unique subsets of model parameters linked to known iPSC‐CM phenotypes

Abstract

There is a profound need to develop a strategy for predicting patient‐to‐patient vulnerability in the emergence of cardiac arrhythmia. A promising in vitro method to address patient‐specific proclivity to cardiac disease utilizes induced pluripotent stem cell‐derived cardiomyocytes (iPSC‐CMs). A major strength of this approach is that iPSC‐CMs contain donor genetic information and therefore capture patient‐specific genotype–phenotype relationships. A cited detriment of iPSC‐CMs is the cell‐to‐cell variability observed in electrical activity. We postulated, however, that cell‐to‐cell variability may constitute a strength when appropriately utilized in a computational framework to build cell populations that can be employed to identify phenotypic mechanisms and pinpoint key sensitive parameters. Thus, we have exploited variation in experimental data across multiple laboratories to develop a computational framework for investigating subcellular phenotypic mechanisms. We have developed a whole‐cell model of iPSC‐CMs composed of simple model components comprising ion channel models with single exponential voltage‐dependent gating variable rate constants, parameterized to fit experimental iPSC‐CM data for all major ionic currents. By optimizing ionic current model parameters to multiple experimental datasets, we incorporate experimentally‐observed variability in the ionic currents. The resulting population of cellular models predicts robust inter‐subject variability in iPSC‐CMs. This approach links molecular mechanisms to known cellular‐level iPSC‐CM phenotypes, as shown by comparing immature and mature subpopulations of models to analyse the contributing factors underlying each phenotype. In the future, the presented models can be readily expanded to include genetic mutations and pharmacological interventions for studying the mechanisms of rare events, such as arrhythmia triggers.

Induced pluripotent stem cell‐derived cardiomyocytes (iPSC‐CMs) capture patient‐specific genotype–phenotype relationships, as well as cell‐to‐cell variability of cardiac electrical activity

Computational modelling and simulation provide a high throughput approach to reconcile multiple datasets describing physiological variability, and also identify vulnerable parameter regimes

We have developed a whole‐cell model of iPSC‐CMs, composed of single exponential voltage‐dependent gating variable rate constants, parameterized to fit experimental iPSC‐CM outputs

We have utilized experimental data across multiple laboratories to model experimental variability and investigate subcellular phenotypic mechanisms in iPSC‐CMs

This framework links molecular mechanisms to cellular‐level outputs by revealing unique subsets of model parameters linked to known iPSC‐CM phenotypes

Large-Scale Simulation of the Phenotypical Variability Induced by Loss-of-Function Long QT Mutations in Human Induced Pluripotent Stem Cell Cardiomyocytes

Loss-of-function long QT (LQT) mutations inducing LQT1 and LQT2 syndromes have been successfully translated to human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) used as disease-specific models. However, their in vitro investigation mainly relies on experiments using small numbers of cells. This is especially critical when working with cells as heterogeneous as hiPSC-CMs. We aim (i) to investigate in silico the ionic mechanisms underlying LQT1 and LQT2 hiPSC-CM phenotypic variability, and (ii) to enable massive in silico drug tests on mutant hiPSC-CMs. We combined (i) data of control and mutant slow and rapid delayed rectifying K⁺ currents, I_Kr and I_Ks respectively, (ii) a recent in silico hiPSC-CM model, and (iii) the population of models paradigm to generate control and mutant populations for LQT1 and LQT2 cardiomyocytes. Our four populations contain from 1008 to 3584 models. In line with the experimental in vitro data, mutant in silico hiPSC-CMs showed prolonged action potential (AP) duration (LQT1: +14%, LQT2: +39%) and large electrophysiological variability. Finally, the mutant populations were split into normal-like hiPSC-CMs (with action potential duration similar to control) and at risk hiPSC-CMs (with clearly prolonged action potential duration). At risk mutant hiPSC-CMs carried higher expression of L-type Ca²⁺, lower expression of I_Kr and increased sensitivity to quinidine as compared to mutant normal-like hiPSC-CMs, resulting in AP abnormalities. In conclusion, we were able to reproduce the two most common LQT syndromes with large-scale simulations, which enable investigating biophysical mechanisms difficult to assess in vitro, e.g., how variations of ion current expressions in a physiological range can impact on AP properties of mutant hiPSC-CMs.

Automatic Optimization of an in Silico Model of Human iPSC Derived Cardiomyocytes Recapitulating Calcium Handling Abnormalities

The growing importance of human induced pluripotent stem cell-derived cardiomyoyctes (hiPSC-CMs), as patient-specific and disease-specific models for studying cellular cardiac electrophysiology or for preliminary cardiotoxicity tests, generated better understanding of hiPSC-CM biophysical mechanisms and great amount of action potential and calcium transient data. In this paper, we propose a new hiPSC-CM in silico model, with particular attention to Ca²⁺ handling. We used (i) the hiPSC-CM Paci2013 model as starting point, (ii) a new dataset of Ca²⁺ transient measurements to tune the parameters of the inward and outward Ca²⁺ fluxes of sarcoplasmic reticulum, and (iii) an automatic parameter optimization to fit action potentials and Ca²⁺ transients. The Paci2018 model simulates, together with the typical hiPSC-CM spontaneous action potentials, more refined Ca²⁺ transients and delayed afterdepolarizations-like abnormalities, which the old Paci2013 was not able to predict due to its mathematical formulation. The Paci2018 model was validated against (i) the same current blocking experiments used to validate the Paci2013 model, and (ii) recently published data about effects of different extracellular ionic concentrations. In conclusion, we present a new and more versatile in silico model, which will provide a platform for modeling the effects of drugs or mutations that affect Ca²⁺ handling in hiPSC-CMs.

Dynamic Model for Characterizing Contractile Behaviors and Mechanical Properties of a Cardiomyocyte

Studies on the contractile dynamics of heart cells have attracted broad attention for the development of both heart disease therapies and cardiomyocyte-actuated micro-robotics. In this study, a linear dynamic model of a single cardiomyocyte cell was proposed at the subcellular scale to characterize the contractile behaviors of heart cells, with system parameters representing the mechanical properties of the subcellular components of living cardiomyocytes. The system parameters of the dynamic model were identified with the cellular beating pattern measured by a scanning ion conductance microscope. The experiments were implemented with cardiomyocytes in one control group and two experimental groups with the drugs cytochalasin-D or nocodazole, to identify the system parameters of the model based on scanning ion conductance microscope measurements, measurement of the cellular Young's modulus with atomic force microscopy indentation, measurement of cellular contraction forces using the micro-pillar technique, and immunofluorescence staining and imaging of the cytoskeleton. The proposed mathematical model was both indirectly and qualitatively verified by the variation in cytoskeleton, beating amplitude, and contractility of cardiomyocytes among the control and the experimental groups, as well as directly and quantitatively validated by the simulation and the significant consistency of 90.5% in the comparison between the ratios of the Young's modulus and the equivalent comprehensive cellular elasticities of cells in the experimental groups to those in the control group. Apart from mechanical properties (mass, elasticity, and viscosity) of subcellular structures, other properties of cardiomyocytes have also been studied, such as the properties of the relative action potential pattern and cellular beating frequency. This work has potential implications for research on cytobiology, drug screening, mechanisms of the heart, and cardiomyocyte-based bio-syncretic robotics.

Tailoring Mathematical Models to Stem-Cell Derived Cardiomyocyte Lines Can Improve Predictions of Drug-Induced Changes to Their Electrophysiology

Human induced pluripotent stem cell derived cardiomyocytes (iPSC-CMs) have applications in disease modeling, cell therapy, drug screening and personalized medicine. Computational models can be used to interpret experimental findings in iPSC-CMs, provide mechanistic insights, and translate these findings to adult cardiomyocyte (CM) electrophysiology. However, different cell lines display different expression of ion channels, pumps and receptors, and show differences in electrophysiology. In this exploratory study, we use a mathematical model based on iPSC-CMs from Cellular Dynamic International (CDI, iCell), and compare its predictions to novel experimental recordings made with the Axiogenesis Cor.4U line. We show that tailoring this model to the specific cell line, even using limited data and a relatively simple approach, leads to improved predictions of baseline behavior and response to drugs. This demonstrates the need and the feasibility to tailor models to individual cell lines, although a more refined approach will be needed to characterize individual currents, address differences in ion current kinetics, and further improve these results.

The Hyperpolarization-Activated Cyclic Nucleotide-Gated Channels: from Biophysics to Pharmacology of a Unique Family of Ion Channels

Hyperpolarization-activated, cyclic nucleotide-gated (HCN) channels are important members of the voltage-gated pore loop channels family. They show unique features: they open at hyperpolarizing potential, carry a mixed Na/K current, and are regulated by cyclic nucleotides. Four different isoforms have been cloned (HCN1-4) that can assemble to form homo- or heterotetramers, characterized by different biophysical properties. These proteins are widely distributed throughout the body and involved in different physiologic processes, the most important being the generation of spontaneous electrical activity in the heart and the regulation of synaptic transmission in the brain. Their role in heart rate, neuronal pacemaking, dendritic integration, learning and memory, and visual and pain perceptions has been extensively studied; these channels have been found also in some peripheral tissues, where their functions still need to be fully elucidated. Genetic defects and altered expression of HCN channels are linked to several pathologies, which makes these proteins attractive targets for translational research; at the moment only one drug (ivabradine), which specifically blocks the hyperpolarization-activated current, is clinically available. This review discusses current knowledge about HCN channels, starting from their biophysical properties, origin, and developmental features, to (patho)physiologic role in different tissues and pharmacological modulation, ending with their present and future relevance as drug targets.

Integrating cardiomyocytes from human pluripotent stem cells in safety pharmacology: has the time come?

Cardiotoxicity is a severe side effect of drugs that induce structural or electrophysiological changes in heart muscle cells. As a result, the heart undergoes failure and potentially lethal arrhythmias. It is still a major reason for drug failure in preclinical and clinical phases of drug discovery. Current methods for predicting cardiotoxicity are based on guidelines that combine electrophysiological analysis of cell lines expressing ion channels ectopically in vitro with animal models and clinical trials. Although no new cases of drugs linked to lethal arrhythmias have been reported since the introduction of these guidelines in 2005, their limited predictive power likely means that potentially valuable drugs may not reach clinical practice. Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) are now emerging as potentially more predictive alternatives, particularly for the early phases of preclinical research. However, these cells are phenotypically immature and culture and assay methods not standardized, which could be a hurdle to the development of predictive computational models and their implementation into the drug discovery pipeline, in contrast to the ambitions of the comprehensive pro-arrhythmia in vitro assay (CiPA) initiative. Here, we review present and future preclinical cardiotoxicity screening and suggest possible hPSC-CM-based strategies that may help to move the field forward. Coordinated efforts by basic scientists, companies and hPSC banks to standardize experimental conditions for generating reliable and reproducible safety indices will be helpful not only for cardiotoxicity prediction but also for precision medicine.

LINKED ARTICLES

This article is part of a themed section on New Insights into Cardiotoxicity Caused by Chemotherapeutic Agents. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v174.21/issuetoc.

Correlation between human ether-a-go-go-related gene channel inhibition and action potential prolongation

BACKGROUND AND PURPOSE

Human ether-a-go-go-related gene (hERG; Kv 11.1) channel inhibition is a widely accepted predictor of cardiac arrhythmia. hERG channel inhibition alone is often insufficient to predict pro-arrhythmic drug effects. This study used a library of dofetilide derivatives to investigate the relationship between standard measures of hERG current block in an expression system and changes in action potential duration (APD) in human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). The interference from accompanying block of Cav 1.2 and Nav 1.5 channels was investigated along with an in silico AP model.

EXPERIMENTAL APPROACH

Drug-induced changes in APD were assessed in hiPSC-CMs using voltage-sensitive dyes. The IC50 values for dofetilide and 13 derivatives on hERG current were estimated in an HEK293 expression system. The relative potency of each drug on APD was estimated by calculating the dose (D150 ) required to prolong the APD at 90% (APD90 ) repolarization by 50%.

KEY RESULTS

The D150 in hiPSC-CMs was linearly correlated with IC50 of hERG current. In silico simulations supported this finding. Three derivatives inhibited hERG without prolonging APD, and these compounds also inhibited Cav 1.2 and/or Nav 1.5 in a channel state-dependent manner. Adding Cav 1.2 and Nav 1.2 block to the in silico model recapitulated the direction but not the extent of the APD change.

CONCLUSIONS AND IMPLICATIONS

Potency of hERG current inhibition correlates linearly with an index of APD in hiPSC-CMs. The compounds that do not correlate have additional effects including concomitant block of Cav 1.2 and/or Nav 1.5 channels. In silico simulations of hiPSC-CMs APs confirm the principle of the multiple ion channel effects.

Developmental changes in the balance of glycolytic ATP production and oxidative phosphorylation in ventricular cells: A simulation study

The developmental program of the heart requires accurate regulation to ensure continuous circulation and simultaneous cardiac morphogenesis, because any functional abnormalities may progress to congenital heart malformation. Notably, energy metabolism in fetal ventricular cells is regulated in a manner that differs from adult ventricular cells: fetal cardiomyocytes generally have immature mitochondria and fetal ventricular cells show greater dependence on glycolytic ATP production. However, although various characteristics of energy metabolism in fetal ventricular cells have been reported, to our knowledge, a quantitative description of the contributions of these factors to fetal ventricular cell functions has not yet been established. Here, we constructed a mathematical model to integrate various characteristics of fetal ventricular cells and predicted the contribution of each characteristic to the maintenance of intracellular ATP concentration and sarcomere contraction under anoxic conditions. Our simulation results demonstrated that higher glycogen content, higher hexokinase activity, and lower creatine concentration helped prolong the time for which ventricular cell contraction was maintained under anoxic conditions. The integrated model also enabled us to quantitatively assess the contributions of factors related to energy metabolism in ventricular cells. Because fetal cardiomyocytes exhibit similar energy metabolic profiles to stem cell-derived cardiomyocytes and those in the failing heart, an improved understanding of these fetal ventricular cells will contribute to a better comprehension of the processes in stem cell-derived cardiomyocytes or under pathological conditions.

Human induced pluripotent stem cell-derived versus adult cardiomyocytes: an in silico electrophysiological study on effects of ionic current block

Background and Purpose

Two new technologies are likely to revolutionize cardiac safety and drug development: in vitro experiments on human‐induced pluripotent stem cell‐derived cardiomyocytes (hiPSC‐CMs) and in silico human adult ventricular cardiomyocyte (hAdultV‐CM) models. Their combination was recently proposed as a potential replacement for the present hERG‐based QT study for pharmacological safety assessments. Here, we systematically compared in silico the effects of selective ionic current block on hiPSC‐CM and hAdultV‐CM action potentials (APs), to identify similarities/differences and to illustrate the potential of computational models as supportive tools for evaluating new in vitro technologies.

Experimental Approach

In silico AP models of ventricular‐like and atrial‐like hiPSC‐CMs and hAdultV‐CM were used to simulate the main effects of four degrees of block of the main cardiac transmembrane currents.

Key Results

Qualitatively, hiPSC‐CM and hAdultV‐CM APs showed similar responses to current block, consistent with results from experiments. However, quantitatively, hiPSC‐CMs were more sensitive to block of (i) L‐type Ca²⁺ currents due to the overexpression of the Na⁺/Ca²⁺ exchanger (leading to shorter APs) and (ii) the inward rectifier K⁺ current due to reduced repolarization reserve (inducing diastolic potential depolarization and repolarization failure).

Conclusions and Implications

In silico hiPSC‐CMs and hAdultV‐CMs exhibit a similar response to selective current blocks. However, overall hiPSC‐CMs show greater sensitivity to block, which may facilitate in vitro identification of drug‐induced effects. Extrapolation of drug effects from hiPSC‐CM to hAdultV‐CM and pro‐arrhythmic risk assessment can be facilitated by in silico predictions using biophysically‐based computational models.

Instantaneous current-voltage relationships during the course of the human cardiac ventricular action potential: new computational insights into repolarization dynamics

AIMS

To adopt a novel three-dimensional (3D) representation of cardiac action potential (AP) to compactly visualize dynamical properties of human cellular ventricular repolarization.

METHODS AND RESULTS

We have recently established a novel 3D representation of cardiac AP, which is based on the iterative measurement of instantaneous ion current-voltage profiles during the course of an AP. Such an approach has been originally developed on real patch-clamped ventricular cells, and subsequently improved in silico on several cardiac ventricular AP models of different mammals, and on models of different AP types of the human heart. We apply it here on two different models of human ventricular AP, and show that it compactly provides further insights into repolarization dynamics. The 3D representation of the AP includes equilibrium points during repolarization, and can be screened in terms of what we have shown to be a region, during late repolarization, when membrane conductance becomes negative and repolarization therefore auto-regenerative. We have called this time window auto-regenerative-repolarization-phase (ARRP).

CONCLUSION

In addition to previous findings obtained through the same procedure, we show here that 3D current-voltage-time representations of human ventricular AP allow compact visualization of dynamical properties, which are relevant for the physiology and pathology of ventricular repolarization. In particular, we suggest that the volume under the current surface corresponding to the ARRP might be used as a predictor of safety of repolarization, in single cells and during AP conduction in cell pairs.

Stem cells and stem cell-derived tissues and their use in safety assessment

Toxicology has long relied on animal models in a tedious approach to understanding risk of exposure to an uncharacterized molecule. Stem cell-derived tissues can be made in high purity, quality, and quantity to enable a new approach to this problem. Currently, stem cell-derived tissues are primarily "generic" genetic backgrounds; the future will see the integration of various genetic backgrounds and complex three-dimensional models to create truly unique in vitro organoids. This minireview focuses on the state of the art of a number of stem cell-derived tissues and details their application in toxicology.

Computational models of ventricular- and atrial-like human induced pluripotent stem cell derived cardiomyocytes

The clear importance of human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) as an in-vitro model highlights the relevance of studying these cells and their function also in-silico. Moreover, the phenotypical differences between the hiPSC-CM and adult myocyte action potentials (APs) call for understanding of how hiPSC-CMs are maturing towards adult myocytes. Using recently published experimental data, we developed two computational models of the hiPSC-CM AP, distinguishing between the ventricular-like and atrial-like phenotypes, emerging during the differentiation process of hiPSC-CMs. Also, we used the computational approach to quantitatively assess the role of ionic mechanisms which are likely responsible for the not completely mature phenotype of hiPSC-CMs. Our models reproduce the typical hiPSC-CM ventricular-like and atrial-like spontaneous APs and the response to prototypical current blockers, namely tetrodotoxine, nifedipine, E4041 and 3R4S-Chromanol 293B. Moreover, simulations using our ventricular-like model suggest that the interplay of immature I Na, I f and I K1 currents has a fundamental role in the hiPSC-CM spontaneous beating whereas a negative shift in I CaL activation causes the observed long lasting AP. In conclusion, this work provides two novel tools useful in investigating the electrophysiological features of hiPSC-CMs, whose importance is growing fast as in-vitro models for pharmacological studies.

Molecular and functional evidence of HCN4 and caveolin-3 interaction during cardiomyocyte differentiation from human embryonic stem cells

Maturation of human embryonic stem cell-derived cardiomyocytes (hESC-CM) is accompanied by changes in ion channel expression, with relevant electrophysiological consequences. In rodent CM, the properties of hyperpolarization-activated cyclic nucleotide-gated channel (HCN)4, a major f-channel isoform, depends on the association with caveolin-3 (Cav3). To date, no information exists on changes in Cav3 expression and its associative relationship with HCN4 upon hESC-CM maturation. We hypothesize that Cav3 expression and its compartmentalization with HCN4 channels during hESC-CM maturation accounts for the progression of f-current properties toward adult phenotypes. To address this, hESC were differentiated into spontaneously beating CM and examined at ∼30, ∼60, and ∼110 days of differentiation. Human adult and fetal CM served as references. HCN4 and Cav3 expression and localization were analyzed by real time PCR and immunocyto/histochemistry. F-current was measured in patch-clamped single cells. HCN4 and Cav3 colocalize in adult human atrial and ventricular CM, but not in fetal CM. Proteins and mRNA for Cav3 were not detected in undifferentiated hESC, but expression increased during hESC-CM maturation. At 110 days, HCN4 appeared to be colocalized with Cav3. Voltage-dependent activation of the f-current was significantly more positive in fetal CM and 60-day hESC-CM (midpoint activation, V1/2, ∼ -82 mV) than in 110-day hESC-CM or adult CM (V1/2∼-100 mV). In the latter cells, caveolae disruption reversed voltage dependence toward a more positive or an immature phenotype, with V1/2 at -75 mV, while in fetal CM voltage dependence was not affected. Our data show, for the first time, a developmental change in HCN4-Cav3 association in hESC-CM. Cav3 expression and its association with ionic channels likely represent a crucial step of cardiac maturation.