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.

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' (I ) 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 I , 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.

Key points• Computational models of the electrical activity of sinoatrial cells (SANCs) have been proposed to gain a deeper understanding of the cellular basis of cardiac pacemaking.• However, they fail to reproduce a number of experimental data, among which are effects measured after modifications of the 'funny' (I f ) current.• We developed a novel SANC mathematical model by updating the description of membrane currents and intracellular mechanisms on the basis of experimental acquisitions, in an attempt to reproduce pacemaker activity and its physiological and pharmacological modulation.• Our model describes satisfactorily experimental data on pacemaking regulation due to neural modulation, I f block and inhibition of the intracellular Ca 2+ handling.• Computer simulation results suggest that a detailed description of the intracellular Ca 2+ fluxes is fully compatible with the observation that I f is a major component of pacemaking and heart rate modulation.Abstract The cellular basis of cardiac pacemaking is still debated. Reliable computational models of the sinoatrial node (SAN) action potential (AP) may help gain a deeper understanding of the phenomenon. Recently, novel models incorporating detailed Ca 2+ -handling dynamics have been proposed, but they fail to reproduce a number of experimental data, and more specifically effects of 'funny' (I f ) current modifications. We therefore developed a SAN AP model, based on available experimental data, in an attempt to reproduce physiological and pharmacological heart rate modulation. Cell compartmentalization and intracellular Ca 2+ -handling mechanisms were formulated as in the Maltsev-Lakatta model, focusing on Ca 2+ -cycling processes. Membrane current equations were revised on the basis of published experimental data. Modifications of the formulation of currents/pumps/exchangers to simulate I f blockers, autonomic modulators and Ca 2+ -dependent mechanisms (ivabradine, caesium, acetylcholine, isoprenaline, BAPTA) were derived from experimental data. The model generates AP waveforms typical of rabbit SAN cells, whose parameters fall within the experimental ranges: 352 ms cycle length, 80 mV AP amplitude, −58 mV maximum diastolic potential (MDP), 108 ms APD 50 , and 7.1 V s −1 maximum upstroke velocity. Rate modulation by I f -blocking drugs agrees with experimental findings: 20% and 22% caesium-induced (5 mM) and ivabradine-induced (3 μM) rate reductions, respectively, due to changes in diastolic depolarization (DD) slope, with no changes in either MDP or take-off potential (TOP). The model consistently reproduces the effects of autonomic modulation: 20% rate decrease with 10 nM acetylcholine and 28% increase with 1 μM isoprenaline, again entirely due to increase in the DD slope, with no changes in either MDP or TOP. Model testing of BAPTA effects showed slowing of rate, −26%, without cessation of beating. Our up-to-date model describes satisfactorily experimental data concerning autonomic stimulation, funny-channel blockade and inhibition of the Ca 2+ -rela...

Human-based modelling and simulations are becoming ubiquitous in biomedical science due to their ability to augment experimental and clinical investigations. Cardiac electrophysiology is one of the most advanced areas, with cardiac modelling and simulation being considered for virtual testing of pharmacological therapies and medical devices. Current models present inconsistencies with experimental data, which limit further progress. In this study, we present the design, development, calibration and independent validation of a human-based ventricular model (ToR-ORd) for simulations of electrophysiology and excitation-contraction coupling, from ionic to whole-organ dynamics, including the electrocardiogram. Validation based on substantial multiscale simulations supports the credibility of the ToR-ORd model under healthy and key disease conditions, as well as drug blockade. In addition, the process uncovers new theoretical insights into the biophysical properties of the L-type calcium current, which are critical for sodium and calcium dynamics. These insights enable the reformulation of L-type calcium current, as well as replacement of the hERG current model.

IntroductionHypertrophic cardiomyopathy (HCM) is a cause of sudden arrhythmic death, but the understanding of its pro-arrhythmic mechanisms and an effective pharmacological treatment are lacking. HCM electrophysiological remodelling includes both increased inward and reduced outward currents, but their role in promoting repolarisation abnormalities remains unknown. The goal of this study is to identify key ionic mechanisms driving repolarisation abnormalities in human HCM, and to evaluate anti-arrhythmic effects of single and multichannel inward current blocks.MethodsExperimental ionic current, action potential (AP) and Ca2 +-transient (CaT) recordings were used to construct populations of human non-diseased and HCM AP models (n = 9118), accounting for inter-subject variability. Simulations were conducted for several degrees of selective and combined inward current block.ResultsSimulated HCM cardiomyocytes exhibited prolonged AP and CaT, diastolic Ca2 + overload and decreased CaT amplitude, in agreement with experiments. Repolarisation abnormalities in HCM models were consistently driven by L-type Ca2 + current (ICaL) re-activation, and ICaL block was the most effective intervention to normalise repolarisation and diastolic Ca2 +, but compromised CaT amplitude. Late Na+ current (INaL) block partially abolished repolarisation abnormalities, with small impact on CaT. Na+/Ca2 + exchanger (INCX) block effectively restored repolarisation and CaT amplitude, but increased Ca2 + overload. Multichannel block increased efficacy in normalising repolarisation, AP biomarkers and CaT amplitude compared to selective block.ConclusionsExperimentally-calibrated populations of human AP models identify ICaL re-activation as the key mechanism for repolarisation abnormalities in HCM, and combined INCX, INaL and ICaL block as effective anti-arrhythmic therapies also able to partially reverse the HCM electrophysiological phenotype.

Ca-calmodulin-dependent protein kinase II (CaMKII) was recently shown to alter Na(+) channel gating and recapitulate a human Na(+) channel genetic mutation that causes an unusual combined arrhythmogenic phenotype in patients: simultaneous long QT syndrome and Brugada syndrome. CaMKII is upregulated in heart failure where arrhythmias are common, and CaMKII inhibition can reduce arrhythmias. Thus, CaMKII-dependent channel modulation may contribute to acquired arrhythmic disease. We developed a Markovian Na(+) channel model including CaMKII-dependent changes, and incorporated it into a comprehensive myocyte action potential (AP) model with Na(+) and Ca(2+) transport. CaMKII shifts Na(+) current (I(Na)) availability to more negative voltage, enhances intermediate inactivation, and slows recovery from inactivation (all loss-of-function effects), but also enhances late noninactivating I(Na) (gain of function). At slow heart rates, with long diastolic time for I(Na) recovery, late I(Na) is the predominant effect, leading to AP prolongation (long QT syndrome). At fast heart rates, where recovery time is limited and APs are shorter, there is little effect on AP duration, but reduced availability decreases I(Na), AP upstroke velocity, and conduction (Brugada syndrome). CaMKII also increases cardiac Ca(2+) and K(+) currents (I(Ca) and I(to)), complicating CaMKII-dependent AP changes. Incorporating I(Ca) and I(to) effects individually prolongs and shortens AP duration. Combining I(Na), I(Ca), and I(to) effects results in shortening of AP duration with CaMKII. With transmural heterogeneity of I(to) and I(to) downregulation in heart failure, CaMKII may accentuate dispersion of repolarization. This provides a useful initial framework to consider pathways by which CaMKII may contribute to arrhythmogenesis.

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 Ca2+ handling. We used (i) the hiPSC-CM Paci2013 model as starting point, (ii) a new dataset of Ca2+ transient measurements to tune the parameters of the inward and outward Ca2+ fluxes of sarcoplasmic reticulum, and (iii) an automatic parameter optimization to fit action potentials and Ca2+ transients. The Paci2018 model simulates, together with the typical hiPSC-CM spontaneous action potentials, more refined Ca2+ 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 Ca2+ handling in hiPSC-CMs.

The main functional derangement in CALM1-F142L was prolonged repolarization with altered rate-dependency and sensitivity to β-adrenergic stimulation. Impaired CDI of ICaL underlined the electrical abnormality, which was sensitive to ICaL blockade. High mutation penetrance was confirmed in the presence of the native genotype, implying strong dominance of effects.