Presenilin-1 ΔE9 mutation associated sarcoplasmic reticulum leak alters [Ca2+]i distribution in human iPSC-derived cardiomyocytes.

Mutations in ubiquitously expressed presenilin genes (PSENs) lead to early-onset familial Alzheimer's disease (FAD), but patients carrying the mutation also suffer from heart diseases. To elucidate the cardiac myocyte specific effects of PSEN ΔE9, we studied cardiomyocytes derived from induced pluripotent stem cells (iPSC-CMs) from patients carrying AD-causing PSEN1 exon 9 deletion (PSEN1 ΔE9). When compared with their isogenic controls, PSEN1 ΔE9 cardiomyocytes showed increased sarcoplasmic reticulum (SR) Ca2+ leak that was resistant to blockage of ryanodine receptors (RyRs) by tetracaine or inositol-3-reseceptors (IP3Rs) by 2-ABP. The SR Ca2+ leak did not affect electrophysiological properties of the hiPSC-CMs, but according to experiments and in silico simulations the leak induces a diastolic buildup of [Ca2+] near the perinuclear SR and reduces the releasable Ca2+ during systole. This demonstrates that PSEN1 ΔE9 induced SR Ca2+ leak has specific effects in iPSC-CMs, reflecting their unique structural and calcium signaling features. The results shed light on the physiological and pathological mechanisms of PSEN1 in cardiac myocytes and explain the intricacies of comorbidity associated with AD-causing mutations in PSEN1.

Translational investigation of electrophysiology in hypertrophic cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) patients are at increased risk of ventricular arrhythmias and sudden cardiac death, which can occur even in the absence of structural changes of the heart. HCM mouse models suggest mutations in myofilament components to affect Ca2+ homeostasis and thereby favor arrhythmia development. Additionally, some of them show indications of pro-arrhythmic changes in cardiac electrophysiology. In this study, we explored arrhythmia mechanisms in mice carrying a HCM mutation in Mybpc3 (Mybpc3-KI) and tested the translatability of our findings in human engineered heart tissues (EHTs) derived from CRISPR/Cas9-generated homozygous MYBPC3 mutant (MYBPC3hom) in induced pluripotent stem cells (iPSC) and to left ventricular septum samples obtained from HCM patients. We observed higher arrhythmia susceptibility in contractility measurements of field-stimulated intact cardiomyocytes and ventricular muscle strips as well as in electromyogram recordings of Langendorff-perfused hearts from adult Mybpc3-KI mice than in wild-type (WT) controls. The latter only occurred in homozygous (Hom-KI) but not in heterozygous (Het-KI) mouse hearts. Both Het- and Hom-KI are known to display pro-arrhythmic increased Ca2+ myofilament sensitivity as a direct consequence of the mutation. In the electrophysiological characterization of the model, we observed smaller repolarizing K+ currents in single cell patch clamp, longer ventricular action potentials in sharp microelectrode recordings and longer ventricular refractory periods in Langendorff-perfused hearts in Hom-KI, but not Het-KI. Interestingly, reduced K+ channel subunit transcript levels and prolonged action potentials were already detectable in newborn, pre-hypertrophic Hom-KI mice. Human iPSC-derived MYBPC3hom EHTs, which genetically mimicked the Hom-KI mice, did exhibit lower mutant mRNA and protein levels, lower force, beating frequency and relaxation time, but no significant alteration of the force-Ca2+ relation in skinned EHTs. Furthermore, MYBPC3hom EHTs did show higher spontaneous arrhythmic behavior, whereas action potentials measured by sharp microelectrode did not differ to isogenic controls. Action potentials measured in septal myectomy samples did not differ between patients with HCM and patients with aortic stenosis, except for the only sample with a MYBPC3 mutation. The data demonstrate that increased myofilament Ca2+ sensitivity is not sufficient to induce arrhythmias in the Mybpc3-KI mouse model and suggest that reduced K+ currents can be a pro-arrhythmic trigger in Hom-KI mice, probably already in early disease stages. However, neither data from EHTs nor from left ventricular samples indicate relevant reduction of K+ currents in human HCM. Therefore, our study highlights the species difference between mouse and human and emphasizes the importance of research in human samples and human-like models.

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.

Refractoriness in human atria: Time and voltage dependence of sodium channel availability.

BACKGROUND: Refractoriness of cardiac cells limits maximum frequency of electrical activity and protects the heart from tonic contractions. Short refractory periods support major arrhythmogenic substrates and augmentation of refractoriness is therefore seen as a main mechanism of antiarrhythmic drugs. Cardiomyocyte excitability depends on availability of sodium channels, which involves both time- and voltage-dependent recovery from inactivation. This study therefore aims to characterise how sodium channel inactivation affects refractoriness in human atria. METHODS AND RESULTS: Steady-state activation and inactivation parameters of sodium channels measured in vitro in isolated human atrial cardiomyocytes were used to parameterise a mathematical human atrial cell model. Action potential data were acquired from human atrial trabeculae of patients in either sinus rhythm or chronic atrial fibrillation. The ex vivo measurements of action potential duration, effective refractory period and resting membrane potential were well-replicated in simulations using this new in silico model. Notably, the voltage threshold potential at which refractoriness was observed was not different between sinus rhythm and chronic atrial fibrillation tissues and was neither affected by changes in frequency (1 vs. 3Hz). CONCLUSIONS: Our results suggest a preferentially voltage-dependent, rather than time-dependent, effect with respect to refractoriness at physiologically relevant rates in human atria. However, as the resting membrane potential is hyperpolarized in chronic atrial fibrillation, the voltage-dependence of excitability dominates, profoundly increasing the risk for arrhythmia re-initiation and maintenance in fibrillating atria. Our results thereby highlight resting membrane potential as a potential target in pharmacological management of chronic atrial fibrillation.

Calcium signalling in developing cardiomyocytes: implications for model systems and disease.

Adult cardiomyocytes exhibit complex Ca(2+) homeostasis, enabling tight control of contraction and relaxation. This intricate regulatory system develops gradually, with progressive maturation of specialized structures and increasing capacity of Ca(2+) sources and sinks. In this review, we outline current understanding of these developmental processes, and draw parallels to pathophysiological conditions where cardiomyocytes exhibit a striking regression to an immature state of Ca(2+) homeostasis. We further highlight the importance of understanding developmental physiology when employing immature cardiomyocyte models such as cultured neonatal cells and stem cells.

Preservation of cardiac function by prolonged action potentials in mice deficient of KChIP2.

Inherited ion channelopathies and electrical remodeling in heart disease alter the cardiac action potential with important consequences for excitation-contraction coupling. Potassium channel-interacting protein 2 (KChIP2) is reduced in heart failure and interacts under physiological conditions with both Kv4 to conduct the fast-recovering transient outward K(+) current (Ito,f) and with CaV1.2 to mediate the inward L-type Ca(2+) current (ICa,L). Anesthetized KChIP2(-/-) mice have normal cardiac contraction despite the lower ICa,L, and we hypothesized that the delayed repolarization could contribute to the preservation of contractile function. Detailed analysis of current kinetics shows that only ICa,L density is reduced, and immunoblots demonstrate unaltered CaV1.2 and CaVß2 protein levels. Computer modeling suggests that delayed repolarization would prolong the period of Ca(2+) entry into the cell, thereby augmenting Ca(2+)-induced Ca(2+) release. Ca(2+) transients in disaggregated KChIP2(-/-) cardiomyocytes are indeed comparable to wild-type transients, corroborating the preserved contractile function and suggesting that the compensatory mechanism lies in the Ca(2+)-induced Ca(2+) release event. We next functionally probed dyad structure, ryanodine receptor Ca(2+) sensitivity, and sarcoplasmic reticulum Ca(2+) load and found that increased temporal synchronicity of the Ca(2+) release in KChIP2(-/-) cardiomyocytes may reflect improved dyad structure aiding the compensatory mechanisms in preserving cardiac contractile force. Thus the bimodal effect of KChIP2 on Ito,f and ICa,L constitutes an important regulatory effect of KChIP2 on cardiac contractility, and we conclude that delayed repolarization and improved dyad structure function together to preserve cardiac contraction in KChIP2(-/-) mice.

In silico screening of the key cellular remodeling targets in chronic atrial fibrillation.

Chronic atrial fibrillation (AF) is a complex disease with underlying changes in electrophysiology, calcium signaling and the structure of atrial myocytes. How these individual remodeling targets and their emergent interactions contribute to cell physiology in chronic AF is not well understood. To approach this problem, we performed in silico experiments in a computational model of the human atrial myocyte. The remodeled function of cellular components was based on a broad literature review of in vitro findings in chronic AF, and these were integrated into the model to define a cohort of virtual cells. Simulation results indicate that while the altered function of calcium and potassium ion channels alone causes a pronounced decrease in action potential duration, remodeling of intracellular calcium handling also has a substantial impact on the chronic AF phenotype. We additionally found that the reduction in amplitude of the calcium transient in chronic AF as compared to normal sinus rhythm is primarily due to the remodeling of calcium channel function, calcium handling and cellular geometry. Finally, we found that decreased electrical resistance of the membrane together with remodeled calcium handling synergistically decreased cellular excitability and the subsequent inducibility of repolarization abnormalities in the human atrial myocyte in chronic AF. We conclude that the presented results highlight the complexity of both intrinsic cellular interactions and emergent properties of human atrial myocytes in chronic AF. Therefore, reversing remodeling for a single remodeled component does little to restore the normal sinus rhythm phenotype. These findings may have important implications for developing novel therapeutic approaches for chronic AF.

Variable t-tubule organization and Ca2+ homeostasis across the atria.

Although t-tubules have traditionally been thought to be absent in atrial cardiomyocytes, recent studies have suggested that t-tubules exist in the atria of large mammals. However, it is unclear whether regional differences in t-tubule organization exist that define cardiomyocyte function across the atria. We sought to investigate regional t-tubule density in pig and rat atria and the consequences for cardiomyocyte Ca(2+) homeostasis. We observed t-tubules in approximately one-third of rat atrial cardiomyocytes, in both tissue cryosections and isolated cardiomyocytes. In a minority (≈10%) of atrial cardiomyocytes, the t-tubular network was well organized, with a transverse structure resembling that of ventricular cardiomyocytes. In both rat and pig atrial tissue, we observed higher t-tubule density in the epicardium than in the endocardium. Consistent with high variability in the distribution of t-tubules and Ca(2+) channels among cells, L-type Ca(2+) current amplitude was also highly variable and steeply dependent on capacitance and t-tubule density. Accordingly, Ca(2+) transients showed great variability in Ca(2+) release synchrony. Simultaneous imaging of the cell membrane and Ca(2+) transients confirmed t-tubule functionality. Results from mathematical modeling indicated that a transmural gradient in t-tubule organization and Ca(2+) release kinetics supports synchronization of contraction across the atrial wall and may underlie transmural differences in the refractory period. In conclusion, our results indicate that t-tubule density is highly variable across the atria. We propose that higher t-tubule density in cells localized in the epicardium may promote synchronization of contraction across the atrial wall.

Investigations of the Navß1b sodium channel subunit in human ventricle; functional characterization of the H162P Brugada syndrome mutant.

Brugada syndrome (BrS) is a rare inherited disease that can give rise to ventricular arrhythmia and ultimately sudden cardiac death. Numerous loss-of-function mutations in the cardiac sodium channel Nav1.5 have been associated with BrS. However, few mutations in the auxiliary Navß1-4 subunits have been linked to this disease. Here we investigated differences in expression and function between Navß1 and Navß1b and whether the H162P/Navß1b mutation found in a BrS patient is likely to be the underlying cause of disease. The impact of Navß subunits was investigated by patch-clamp electrophysiology, and the obtained in vitro values were used for subsequent in silico modeling. We found that Navß1b transcripts were expressed at higher levels than Navß1 transcripts in the human heart. Navß1 and Navß1b coexpressed with Nav1.5 induced a negative shift on steady state of activation and inactivation compared with Nav1.5 alone. Furthermore, Navß1b was found to increase the current level when coexpressed with Nav1.5, Navß1b/H162P mutated subunit peak current density was reduced by 48% (-645 ± 151 vs. -334 ± 71 pA/pF), V1/2 steady-state inactivation shifted by -6.7 mV (-70.3 ± 1.5 vs. -77.0 ± 2.8 mV), and time-dependent recovery from inactivation slowed by >50% compared with coexpression with Navß1b wild type. Computer simulations revealed that these electrophysiological changes resulted in a reduction in both action potential amplitude and maximum upstroke velocity. The experimental data thereby indicate that Navß1b/H162P results in reduced sodium channel activity functionally affecting the ventricular action potential. This result is an important replication to support the notion that BrS can be linked to the function of Navß1b and is associated with loss-of-function of the cardiac sodium channel.

NS5806 partially restores action potential duration but fails to ameliorate calcium transient dysfunction in a computational model of canine heart failure.

AIMS: The study investigates how increased Ito, as mediated by the activator NS5806, affects excitation-contraction coupling in chronic heart failure (HF). We hypothesized that restoring spike-and-dome morphology of the action potential (AP) to a healthy phenotype would be insufficient to restore the intracellular Ca(2) (+) transient (CaT), due to HF-induced remodelling of Ca(2+) handling. METHODS AND RESULTS: An existing mathematical model of the canine ventricular myocyte was modified to incorporate recent experimental data from healthy and failing myocytes, resulting in models of both healthy and HF epicardial, midmyocardial, and endocardial cell variants. Affects of NS5806 were also included in HF models through its direct interaction with Kv4.3 and Kv1.4. Single-cell simulations performed in all models (control, HF, and HF + drug) and variants (epi, mid, and endo) assessed AP morphology and underlying ionic processes with a focus on calcium transients (CaT), how these were altered in HF across the ventricular wall, and the subsequent effects of varying compound concentration in HF. Heart failure model variants recapitulated a characteristic increase in AP duration (APD) in the disease. The qualitative effects of application of half-maximal effective concentration (EC50) of NS5806 on APs and CaT are heterogeneous and non-linear. Deepening in the AP notch with drug is a direct effect of the activation of Ito; both Ito and consequent alteration of IK1 kinetics cause decrease in AP plateau potential. Decreased APD50 and APD90 are both due to altered IK1. Analysis revealed that drug effects depend on transmurality. Ca(2+) transient morphology changes-increased amplitude and shorter time to peak-are due to direct increase in ICa,L and indirect larger SR Ca(2+) release subsequent to Ito activation. CONCLUSIONS: Downstream effects of a compound acting exclusively on sarcolemmal ion channels are difficult to predict. Remediation of APD to pre-failing state does not ameliorate dysfunction in CaT; however, restoration of notch depth appears to impart modest benefit and a likelihood of therapeutic value in modulating early repolarization.

In silico study of the mechanisms of hypoxia and contractile dysfunction during ischemia and reperfusion of hiPSC cardiomyocytes.

Interconnected mechanisms of ischemia and reperfusion (IR) has increased the interest in IR in vitro experiments using human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). We developed a whole-cell computational model of hiPSC-CMs including the electromechanics, a metabolite-sensitive sarcoplasmic reticulum Ca2+-ATPase (SERCA) and an oxygen dynamics formulation to investigate IR mechanisms. Moreover, we simulated the effect and action mechanism of levosimendan, which recently showed promising anti-arrhythmic effects in hiPSC-CMs in hypoxia. The model was validated using hiPSC-CM and in vitro animal data. The role of SERCA in causing relaxation dysfunction in IR was anticipated to be comparable to its function in sepsis-induced heart failure. Drug simulations showed that levosimendan counteracts the relaxation dysfunction by utilizing a particular Ca2+-sensitizing mechanism involving Ca2+-bound troponin C and Ca2+ flux to the myofilament, rather than inhibiting SERCA phosphorylation. The model demonstrates extensive characterization and promise for drug development, making it suitable for evaluating IR therapy strategies based on the changing levels of cardiac metabolites, oxygen and molecular pathways.

Building blocks of microphysiological system to model physiology and pathophysiology of human heart.

Microphysiological systems (MPS) are drawing increasing interest from academia and from biomedical industry due to their improved capability to capture human physiology. MPS offer an advanced in vitro platform that can be used to study human organ and tissue level functions in health and in diseased states more accurately than traditional single cell cultures or even animal models. Key features in MPS include microenvironmental control and monitoring as well as high biological complexity of the target tissue. To reach these qualities, cross-disciplinary collaboration from multiple fields of science is required to build MPS. Here, we review different areas of expertise and describe essential building blocks of heart MPS including relevant cardiac cell types, supporting matrix, mechanical stimulation, functional measurements, and computational modelling. The review presents current methods in cardiac MPS and provides insights for future MPS development with improved recapitulation of human physiology.

Human induced pluripotent stem cell-derived cardiomyocytes as an electrophysiological model: Opportunities and challenges-The Hamburg perspective.

Models based on human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) are proposed in almost any field of physiology and pharmacology. The development of human induced pluripotent stem cell-derived cardiomyocytes is expected to become a step forward to increase the translational power of cardiovascular research. Importantly they should allow to study genetic effects on an electrophysiological background close to the human situation. However, biological and methodological issues revealed when human induced pluripotent stem cell-derived cardiomyocytes were used in experimental electrophysiology. We will discuss some of the challenges that should be considered when human induced pluripotent stem cell-derived cardiomyocytes will be used as a physiological model.

Impact of sarcoplasmic reticulum calcium release on calcium dynamics and action potential morphology in human atrial myocytes: a computational study.

Electrophysiological studies of the human heart face the fundamental challenge that experimental data can be acquired only from patients with underlying heart disease. Regarding human atria, there exist sizable gaps in the understanding of the functional role of cellular Ca²+ dynamics, which differ crucially from that of ventricular cells, in the modulation of excitation-contraction coupling. Accordingly, the objective of this study was to develop a mathematical model of the human atrial myocyte that, in addition to the sarcolemmal (SL) ion currents, accounts for the heterogeneity of intracellular Ca²+ dynamics emerging from a structurally detailed sarcoplasmic reticulum (SR). Based on the simulation results, our model convincingly reproduces the principal characteristics of Ca²+ dynamics: 1) the biphasic increment during the upstroke of the Ca²+ transient resulting from the delay between the peripheral and central SR Ca²+ release, and 2) the relative contribution of SL Ca²+ current and SR Ca²+ release to the Ca²+ transient. In line with experimental findings, the model also replicates the strong impact of intracellular Ca²+ dynamics on the shape of the action potential. The simulation results suggest that the peripheral SR Ca²+ release sites define the interface between Ca²+ and AP, whereas the central release sites are important for the fire-diffuse-fire propagation of Ca²+ diffusion. Furthermore, our analysis predicts that the modulation of the action potential duration due to increasing heart rate is largely mediated by changes in the intracellular Na+ concentration. Finally, the results indicate that the SR Ca²+ release is a strong modulator of AP duration and, consequently, myocyte refractoriness/excitability. We conclude that the developed model is robust and reproduces many fundamental aspects of the tight coupling between SL ion currents and intracellular Ca²+ signaling. Thus, the model provides a useful framework for future studies of excitation-contraction coupling in human atrial myocytes.

Altered contractility in mutation-specific hypertrophic cardiomyopathy: A mechano-energetic in silico study with pharmacological insights.

Introduction: Mavacamten (MAVA), Blebbistatin (BLEB), and Omecamtiv mecarbil (OM) are promising drugs directly targeting sarcomere dynamics, with demonstrated efficacy against hypertrophic cardiomyopathy (HCM) in (pre)clinical trials. However, the molecular mechanism affecting cardiac contractility regulation, and the diseased cell mechano-energetics are not fully understood yet. Methods: We present a new metabolite-sensitive computational model of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) electromechanics to investigate the pathology of R403Q HCM mutation and the effect of MAVA, BLEB, and OM on the cell mechano-energetics. Results: We offer a mechano-energetic HCM calibration of the model, capturing the prolonged contractile relaxation due to R403Q mutation (∼33%), without assuming any further modifications such as an additional Ca2+ flux to the thin filaments. The HCM model variant correctly predicts the negligible alteration in ATPase activity in R403Q HCM condition compared to normal hiPSC-CMs. The simulated inotropic effects of MAVA, OM, and BLEB, along with the ATPase activities in the control and HCM model variant agree with in vitro results from different labs. The proposed model recapitulates the tension-Ca2+ relationship and action potential duration change due to 1 µM OM and 5 µM BLEB, consistently with in vitro data. Finally, our model replicates the experimental dose-dependent effect of OM and BLEB on the normalized isometric tension. Conclusion: This work is a step toward deep-phenotyping the mutation-specific HCM pathophysiology, manifesting as altered interfilament kinetics. Accordingly, the modeling efforts lend original insights into the MAVA, BLEB, and OM contributions to a new interfilament balance resulting in a cardioprotective effect.

Computational cardiac physiology for new modelers: Origins, foundations, and future.

Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.

Regulation of APD and Force by the Na+/Ca2+ Exchanger in Human-Induced Pluripotent Stem Cell-Derived Engineered Heart Tissue.

The physiological importance of NCX in human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) is not well characterized but may depend on the relative strength of the current, compared to adult cardiomyocytes, and on the exact spatial arrangement of proteins involved in Ca2+ extrusion. Here, we determined NCX currents and its contribution to action potential and force in hiPSC-CMs cultured in engineered heart tissue (EHT). The results were compared with data from rat and human left ventricular tissue. The NCX currents in hiPSC-CMs were larger than in ventricular cardiomyocytes isolated from human left ventricles (1.3 ± 0.2 pA/pF and 3.2 ± 0.2 pA/pF for human ventricle and EHT, respectively, p < 0.05). SEA0400 (10 µM) markedly shortened the APD90 in EHT (by 26.6 ± 5%, p < 0.05) and, to a lesser extent, in rat ventricular tissue (by 10.7 ± 1.6%, p < 0.05). Shortening in human left ventricular preparations was small and not different from time-matched controls (TMCs; p > 0.05). Force was increased by the NCX block in rat ventricle (by 31 ± 5.4%, p < 0.05) and EHT (by 20.8 ± 3.9%, p < 0.05), but not in human left ventricular preparations. In conclusion, hiPSC-CMs possess NCX currents not smaller than human left ventricular tissue. Robust NCX block-induced APD shortening and inotropy makes EHT an attractive pharmacological model.

Ca2+-calmodulin-dependent protein kinase II represses cardiac transcription of the L-type calcium channel alpha(1C)-subunit gene (Cacna1c) by DREAM translocation.

Recent studies have demonstrated that changes in the activity of calcium-calmodulin-dependent protein kinase II (CaMKII) induce a unique cardiomyocyte phenotype through the regulation of specific genes involved in excitation-contraction (E-C)-coupling. To explain the transcriptional effects of CaMKII we identified a novel CaMKII-dependent pathway for controlling the expression of the pore-forming α-subunit (Cav1.2) of the L-type calcium channel (LTCC) in cardiac myocytes. We show that overexpression of either cytosolic (δC) or nuclear (δB) CaMKII isoforms selectively downregulate the expression of the Cav1.2. Pharmacological inhibition of CaMKII activity induced measurable changes in LTCC current density and subsequent changes in cardiomyocyte calcium signalling in less than 24 h. The effect of CaMKII on the α1C-subunit gene (Cacna1c) promoter was abolished by deletion of the downstream regulatory element (DRE), which binds transcriptional repressor DREAM/calsenilin/KChIP3. Imaging DREAM-GFP (green fluorescent protein)-expressing cardiomyocytes showed that CaMKII potentiates the calcium-induced nuclear translocation of DREAM. Thereby CaMKII increases DREAM binding to the DRE consensus sequence of the endogenous Cacna1c gene. By mathematical modelling we demonstrate that the LTCC downregulation through the Ca2+-CaMKII-DREAM cascade constitutes a physiological feedback mechanism enabling cardiomyocytes to adjust the calcium intrusion through LTCCs to the amount of intracellular calcium detected by CaMKII.

A mathematical model of hiPSC cardiomyocytes electromechanics.

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are becoming instrumental in cardiac research, human-based cell level cardiotoxicity tests, and developing patient-specific care. As one of the principal functional readouts is contractility, we propose a novel electromechanical hiPSC-CM computational model named the hiPSC-CM-CE. This model comprises a reparametrized version of contractile element (CE) by Rice et al., 2008, with a new passive force formulation, integrated into a hiPSC-CM electrophysiology formalism by Paci et al. in 2020. Our simulated results were validated against in vitro data reported for hiPSC-CMs at matching conditions from different labs. Specifically, key action potential (AP) and calcium transient (CaT) biomarkers simulated by the hiPSC-CM-CE model were within the experimental ranges. On the mechanical side, simulated cell shortening, contraction-relaxation kinetic indices (RT50 and RT25 ), and the amplitude of tension fell within the experimental intervals. Markedly, as an inter-scale analysis, correct classification of the inotropic effects due to non-cardiomyocytes in hiPSC-CM tissues was predicted on account of the passive force expression introduced to the CE. Finally, the physiological inotropic effects caused by Verapamil and Bay-K 8644 and the aftercontractions due to the early afterdepolarizations (EADs) were simulated and validated against experimental data. In the future, the presented model can be readily expanded to take in pharmacological trials and genetic mutations, such as those involved in hypertrophic cardiomyopathy, and study arrhythmia trigger mechanisms.

Comparison of the Simulated Response of Three in Silico Human Stem Cell-Derived Cardiomyocytes Models and in Vitro Data Under 15 Drug Actions.

Objectives: Improvements in human stem cell-derived cardiomyocyte (hSC-CM) technology have promoted their use for drug testing and disease investigations. Several in silico hSC-CM models have been proposed to augment interpretation of experimental findings through simulations. This work aims to assess the response of three hSC-CM in silico models (Koivumäki2018, Kernik2019, and Paci2020) to simulated drug action, and compare simulation results against in vitro data for 15 drugs. Methods: First, simulations were conducted considering 15 drugs, using a simple pore-block model and experimental data for seven ion channels. Similarities and differences were analyzed in the in silico responses of the three models to drugs, in terms of Ca2+ transient duration (CTD90) and occurrence of arrhythmic events. Then, the sensitivity of each model to different degrees of blockage of Na+ (INa), L-type Ca2+ (ICaL), and rapid delayed rectifying K+ (IKr) currents was quantified. Finally, we compared the drug-induced effects on CTD90 against the corresponding in vitro experiments. Results: The observed CTD90 changes were overall consistent among the in silico models, all three showing changes of smaller magnitudes compared to the ones measured in vitro. For example, sparfloxacin 10 µM induced +42% CTD90 prolongation in vitro, and +17% (Koivumäki2018), +6% (Kernik2019), and +9% (Paci2020) in silico. Different arrhythmic events were observed following drug application, mainly for drugs affecting IKr. Paci2020 and Kernik2019 showed only repolarization failure, while Koivumäki2018 also displayed early and delayed afterdepolarizations. The spontaneous activity was suppressed by Na+ blockers and by drugs with similar effects on ICaL and IKr in Koivumäki2018 and Paci2020, while only by strong ICaL blockers, e.g. nisoldipine, in Kernik2019. These results were confirmed by the sensitivity analysis. Conclusion: To conclude, The CTD90 changes observed in silico are qualitatively consistent with our in vitro data, although our simulations show differences in drug responses across the hSC-CM models, which could stem from variability in the experimental data used in their construction.