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

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

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.

Strengthening cardiac therapy pipelines using human pluripotent stem cell-derived cardiomyocytes

Advances in hiPSC isolation and reprogramming and hPSC-CM differentiation have prompted their therapeutic application and utilization for evaluating potential cardiovascular safety liabilities. In this perspective, we showcase key efforts toward the large-scale production of hiPSC-CMs, implementation of hiPSC-CMs in industry settings, and recent clinical applications of this technology. The key observations are a need for traceable gender and ethnically diverse hiPSC lines, approaches to reduce cost of scale-up, accessible clinical trial datasets, and transparent guidelines surrounding the safety and efficacy of hiPSC-based therapies.

Creating cell-specific computational models of stem cell-derived cardiomyocytes using optical experiments

Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) have gained traction as a powerful model in cardiac disease and therapeutics research, since iPSCs are self-renewing and can be derived from healthy and diseased patients without invasive surgery. However, current iPSC-CM differentiation methods produce cardiomyocytes with immature, fetal-like electrophysiological phenotypes, and the variety of maturation protocols in the literature results in phenotypic differences between labs. Heterogeneity of iPSC donor genetic backgrounds contributes to additional phenotypic variability. Several mathematical models of iPSC-CM electrophysiology have been developed to help understand the ionic underpinnings of, and to simulate, various cell responses, but these models individually do not capture the phenotypic variability observed in iPSC-CMs. Here, we tackle these limitations by developing a computational pipeline to calibrate cell preparation-specific iPSC-CM electrophysiological parameters. We used the genetic algorithm (GA), a heuristic parameter calibration method, to tune ion channel parameters in a mathematical model of iPSC-CM physiology. To systematically optimize an experimental protocol that generates sufficient data for parameter calibration, we created simulated datasets by applying various protocols to a population of in silico cells with known conductance variations, and we fitted to those datasets. We found that calibrating models to voltage and calcium transient data under 3 varied experimental conditions, including electrical pacing combined with ion channel blockade and changing buffer ion concentrations, improved model parameter estimates and model predictions of unseen channel block responses. This observation held regardless of whether the fitted data were normalized, suggesting that normalized fluorescence recordings, which are more accessible and higher throughput than patch clamp recordings, could sufficiently inform conductance parameters. Therefore, this computational pipeline can be applied to different iPSC-CM preparations to determine cell line-specific ion channel properties and understand the mechanisms behind variability in perturbation responses.

Injection of IK1 through dynamic clamp can make all the difference in patch-clamp studies on hiPSC-derived cardiomyocytes

Human-induced stem cell-derived cardiomyocytes (hiPSC-CMs) are a valuable tool for studying development, pharmacology, and (inherited) arrhythmias. Unfortunately, hiPSC-CMs are depolarized and spontaneously active, even the working cardiomyocyte subtypes such as atrial- and ventricular-like hiPSC-CMs, in contrast to the situation in the atria and ventricles of adult human hearts. Great efforts have been made, using many different strategies, to generate more mature, quiescent hiPSC-CMs with more close-to-physiological resting membrane potentials, but despite promising results, it is still difficult to obtain hiPSC-CMs with such properties. The dynamic clamp technique allows to inject a current with characteristics of the inward rectifier potassium current (IK1), computed in real time according to the actual membrane potential, into patch-clamped hiPSC-CMs during action potential measurements. This results in quiescent hiPSC-CMs with a close-to-physiological resting membrane potential. As a result, action potential measurements can be performed with normal ion channel availability, which is particularly important for the physiological functioning of the cardiac SCN5A-encoded fast sodium current (INa). We performed in vitro and in silico experiments to assess the beneficial effects of the dynamic clamp technique in dissecting the functional consequences of the SCN5A-1795insD+/- mutation. In two separate sets of patch-clamp experiments on control hiPSC-CMs and on hiPSC-CMs with mutations in ACADVL and GNB5, we assessed the value of dynamic clamp in detecting delayed afterdepolarizations and in investigating factors that modulate the resting membrane potential. We conclude that the dynamic clamp technique has highly beneficial effects in all of the aforementioned settings and should be widely used in patch-clamp studies on hiPSC-CMs while waiting for the ultimate fully mature hiPSC-CMs.

Model-driven optimal experimental design for calibrating cardiac electrophysiology models

BACKGROUND AND OBJECTIVE

Models of the cardiomyocyte action potential have contributed immensely to the understanding of heart function, pathophysiology, and the origin of heart rhythm disturbances. However, action potential models are highly nonlinear, making them difficult to parameterise and limiting to describing 'average cell' dynamics, when cell-specific models would be ideal to uncover inter-cell variability but are too experimentally challenging to be achieved. Here, we focus on automatically designing experimental protocols that allow us to better identify cell-specific maximum conductance values for each major current type.

METHODS AND RESULTS

We developed an approach that applies optimal experimental designs to patch-clamp experiments, including both voltage-clamp and current-clamp experiments. We assessed the models calibrated to these new optimal designs by comparing them to the models calibrated to some of the commonly used designs in the literature. We showed that optimal designs are not only overall shorter in duration but also able to perform better than many of the existing experiment designs in terms of identifying model parameters and hence model predictive power.

CONCLUSIONS

For cardiac cellular electrophysiology, this approach will allow researchers to define their hypothesis of the dynamics of the system and automatically design experimental protocols that will result in theoretically optimal designs.

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.

Electrophysiological and calcium-handling development during long-term culture of human-induced pluripotent stem cell-derived cardiomyocytes

Human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are increasingly used for personalised medicine and preclinical cardiotoxicity testing. Reports on hiPSC-CM commonly describe heterogenous functional readouts and underdeveloped or immature phenotypical properties. Cost-effective, fully defined monolayer culture is approaching mainstream adoption; however, the optimal age at which to utilise hiPSC-CM is unknown. In this study, we identify, track and model the dynamic developmental behaviour of key ionic currents and Ca2+-handling properties in hiPSC-CM over long-term culture (30-80 days). hiPSC-CMs > 50 days post differentiation show significantly larger ICa,L density along with an increased ICa,L-triggered Ca2+-transient. INa and IK1 densities significantly increase in late-stage cells, contributing to increased upstroke velocity and reduced action potential duration, respectively. Importantly, our in silico model of hiPSC-CM electrophysiological age dependence confirmed IK1 as the key ionic determinant of action potential shortening in older cells. We have made this model available through an open source software interface that easily allows users to simulate hiPSC-CM electrophysiology and Ca2+-handling and select the appropriate age range for their parameter of interest. This tool, together with the insights from our comprehensive experimental characterisation, could be useful in future optimisation of the culture-to-characterisation pipeline in the field of hiPSC-CM research.

Assessing the arrhythmogenic risk of engineered heart tissue patches through in silico application on infarcted ventricle models

BACKGROUND

Post myocardial infarction (MI) ventricles contain fibrotic tissue and may have disrupted electrical properties, both of which predispose to an increased risk of life-threatening arrhythmias. Application of epicardial patches obtained from human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are a potential long-term therapy to treat heart failure resulting from post MI remodelling. However, whether the introduction of these patches is anti- or pro-arrhythmic has not been studied.

METHODS

We studied arrhythmic risk using in silico engineered heart tissue (EHT) patch engraftment on human post-MI ventricular models. Two patient models were studied, including one with a large dense scar and one with an apparent channel of preserved viability bordered on both sides by scar. In each heart model a virtual EHT patch was introduced as a layer of viable tissue overlying the scarred area, with hiPSC-CMs electrophysiological properties. The incidence of re-entrant and sustained activation in simulations with and without EHT patches was assessed and the arrhythmia inducibility compared in the context of different EHT patch properties (conduction velocity (CV) and action potential duration (APD)). The impact of the EHT patch on the likelihood of focal ectopic impulse propagation was estimated by assessing the minimum stimulus strength and duration required to generate a propagating impulse in the scar border zone (BZ) with and without patch.

RESULTS

We uncovered two main mechanisms by which ventricular tachycardia (VT) risk could be either augmented or attenuated by the interaction of the patch with the tissue. In the case of isthmus-related VT, our simulations predict that EHT patches can prevent the induction of VT when the, generally longer, hiPSC-CMs APD is reduced towards more physiological values. In the case of large dense scar, we found that, an EHT patch with CV similar to the host myocardium does not promote VT, while EHT patches with lower CV increase the risk of VT, by promoting both non-sustained and sustained re-entry. Finally, our simulations indicate that electrically coupled EHT patches reduce the likelihood of propagation of focal ectopic impulses.

CONCLUSIONS

The introduction of EHT patches as a treatment for heart failure has the potential to augment or attenuate the risk of ventricular arrhythmias, and variations in the anatomic configuration of the substrate, the functional properties of the BZ and the electrophysiologic properties of the patch itself will determine the overall impact. Planning for delivery of this therapy will need to consider the possible impact on arrhythmia.

Polymer Kernels as Compact Carriers for Suspended Cardiomyocytes

Induced pluripotent stem cells (iPSCs) constitute a potential source of patient-specific human cardiomyocytes for a cardiac cell replacement therapy via intramyocardial injections, providing a major benefit over other cell sources in terms of immune rejection. However, intramyocardial injection of the cardiomyocytes has substantial challenges related to cell survival and electrophysiological coupling with recipient tissue. Current methods of manipulating cell suspensions do not allow one to control the processes of adhesion of injected cells to the tissue and electrophysiological coupling with surrounding cells. In this article, we documented the possibility of influencing these processes using polymer kernels: biocompatible fiber fragments of subcellular size that can be adsorbed to a cell, thereby creating the minimum necessary adhesion foci to shape the cell and provide support for the organization of the cytoskeleton and the contractile apparatus prior to adhesion to the recipient tissue. Using optical excitation markers, the restoration of the excitability of cardiomyocytes in suspension upon adsorption of polymer kernels was shown. It increased the likelihood of the formation of a stable electrophysiological coupling in vitro. The obtained results may be considered as a proof of concept that the stochastic engraftment process of injected suspension cells can be controlled by smart biomaterials.

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 Ca²⁺ 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-Ca²⁺ 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.

An in silico-in vitro pipeline for drug cardiotoxicity screening identifies ionic pro-arrhythmia mechanisms

BACKGROUND AND PURPOSE

Before advancing to clinical trials, new drugs are screened for their pro-arrhythmic potential using a method that is overly conservative and provides limited mechanistic insight. The shortcomings of this approach can lead to the mis-classification of beneficial drugs as pro-arrhythmic.

EXPERIMENTAL APPROACH

An in silico-in vitro pipeline was developed to circumvent these shortcomings. A computational human induced pluripotent stem cell-derived cardiomyocyte (iPSC-CM) model was used as part of a genetic algorithm to design experiments, specifically electrophysiological voltage clamp (VC) protocols, to identify which of several cardiac ion channels were blocked during in vitro drug studies. Such VC data, along with dynamically clamped action potentials (AP), were acquired from iPSC-CMs before and after treatment with a control solution or a low- (verapamil), intermediate- (cisapride or quinine) or high-risk (quinidine) drug.

KEY RESULTS

Significant AP prolongation (a pro-arrhythmia marker) was seen in response to quinidine and quinine. The VC protocol identified block of IKr (a source of arrhythmias) by all strong IKr blockers, including cisapride, quinidine and quinine. The protocol also detected block of ICaL by verapamil and Ito by quinidine. Further demonstrating the power of the approach, the VC data uncovered a previously unidentified If block by quinine, which was confirmed with experiments using a HEK-293 expression system and automated patch-clamp.

CONCLUSION AND IMPLICATIONS

We developed an in silico-in vitro pipeline that simultaneously identifies pro-arrhythmia risk and mechanism of ion channel-blocking drugs. The approach offers a new tool for evaluating cardiotoxicity during preclinical drug screening.

Modelling the interaction between stem cells derived cardiomyocytes patches and host myocardium to aid non-arrhythmic engineered heart tissue design

Application of epicardial patches constructed from human-induced pluripotent stem cell- derived cardiomyocytes (hiPSC-CMs) has been proposed as a long-term therapy to treat scarred hearts post myocardial infarction (MI). Understanding electrical interaction between engineered heart tissue patches (EHT) and host myocardium represents a key step toward a successful patch engraftment. EHT retain different electrical properties with respect to the host heart tissue due to the hiPSC-CMs immature phenotype, which may lead to increased arrhythmia risk. We developed a modelling framework to examine the influence of patch design on electrical activation at the engraftment site. We performed an in silico investigation of different patch design approaches to restore pre-MI activation properties and evaluated the associated arrhythmic risk. We developed an in silico cardiac electrophysiology model of a transmural cross section of host myocardium. The model featured an infarct region, an epicardial patch spanning the infarct region and a bath region. The patch is modelled as a layer of hiPSC-CM, combined with a layer of conductive polymer (CP). Tissue and patch geometrical dimensions and conductivities were incorporated through 10 modifiable model parameters. We validated our model against 4 independent experimental studies and showed that it can qualitatively reproduce their findings. We performed a global sensitivity analysis (GSA) to isolate the most important parameters, showing that the stimulus propagation is mainly governed by the scar depth, radius and conductivity when the scar is not transmural, and by the EHT patch conductivity when the scar is transmural. We assessed the relevance of small animal studies to humans by comparing simulations of rat, rabbit and human myocardium. We found that stimulus propagation paths and GSA sensitivity indices are consistent across species. We explored which EHT design variables have the potential to restore physiological propagation. Simulations predict that increasing EHT conductivity from 0.28 to 1-1.1 S/m recovered physiological activation in rat, rabbit and human. Finally, we assessed arrhythmia risk related to increasing EHT conductivity and tested increasing the EHT Na+ channel density as an alternative strategy to match healthy activation. Our results revealed a greater arrhythmia risk linked to increased EHT conductivity compared to increased Na+ channel density. We demonstrated that our modeling framework could capture the interaction between host and EHT patches observed in in vitro experiments. We showed that large (patch and tissue dimensions) and small (cardiac myocyte electrophysiology) scale differences between small animals and humans do not alter EHT patch effect on infarcted tissue. Our model revealed that only when the scar is transmural do EHT properties impact activation times and isolated the EHT conductivity as the main parameter influencing propagation. We predicted that restoring physiological activation by tuning EHT conductivity is possible but may promote arrhythmic behavior. Finally, our model suggests that acting on hiPSC-CMs low action potential upstroke velocity and lack of IK1 may restore pre-MI activation while not promoting arrhythmia.

Computational modeling of aberrant electrical activity following remuscularization with intramyocardially injected pluripotent stem cell-derived cardiomyocytes

Acute engraftment arrhythmias (EAs) remain a serious complication of remuscularization therapy. Preliminary evidence suggests that a focal source underlies these EAs stemming from the automaticity of immature pluripotent stem cell-derived cardiomyocytes (PSC-CMs) in nascent myocardial grafts. How these EAs arise though during early engraftment remains unclear. In a series of in silico experiments, we probed the origin of EAs-exploring aspects of altered impulse formation and altered impulse propagation within nascent PSC-CM grafts and at the host-graft interface. To account for poor gap junctional coupling during early PSC-CM engraftment, the voltage dependence of gap junctions and the possibility of ephaptic coupling were incorporated. Inspired by cardiac development, we also studied the contributions of another feature of immature PSC-CMs, circumferential sodium channel (NaCh) distribution in PSC-CMs. Ectopic propagations emerged from nascent grafts of immature PSC-CMs at a rate of <96 bpm. Source-sink effects dictated this rate and contributed to intermittent capture between host and graft. Moreover, ectopic beats emerged from dynamically changing sites along the host-graft interface. The latter arose in part because circumferential NaCh distribution in PSC-CMs contributed to preferential conduction slowing and block of electrical impulses from host to graft myocardium. We conclude that additional mechanisms, in addition to focal ones, contribute to EAs and recognize that their relative contributions are dynamic across the engraftment process.

Quantitative cross-species translators of cardiac myocyte electrophysiology: Model training, experimental validation, and applications

Animal experimentation is key in the evaluation of cardiac efficacy and safety of novel therapeutic compounds. However, interspecies differences in the mechanisms regulating excitation-contraction coupling can limit the translation of experimental findings from animal models to human physiology and undermine the assessment of drugs’ efficacy and safety. Here, we built a suite of translators for quantitatively mapping electrophysiological responses in ventricular myocytes across species. We trained these statistical operators using a broad dataset obtained by simulating populations of our biophysically detailed computational models of action potential and Ca2+ transient in mouse, rabbit, and human. We then tested our translators against experimental data describing the response to stimuli, such as ion channel block, change in beating rate, and β-adrenergic challenge. We demonstrate that this approach is well suited to predicting the effects of perturbations across different species or experimental conditions and suggest its integration into mechanistic studies and drug development pipelines.

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.

Analysis of the response of human iPSC-derived cardiomyocyte tissue to ICaL block. A combined in vitro and in silico approach

The high incidence of cardiac arrythmias underlines the need for the assessment of pharmacological therapies. In this field of drug efficacy, as in the field of drug safety highlighted by the Comprehensive in Vitro Proarrhythmia Assay initiative, new pillars for research have become crucial: firstly, the integration of in-silico experiments, and secondly the evaluation of fully integrated biological systems, such as human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). In this study, we therefore aimed to combine in-vitro experiments and in-silico simulations to evaluate the antiarrhythmic effect of L-type calcium current (I_CaL) block in hiPSC-CMs. For this, hiPSC-CM preparations were cultured and an equivalent virtual tissue was modeled. Re-entry patterns of electrical activation were induced and several biomarkers were obtained before and after I_CaL block. The virtual hiPSC-CM simulations were also reproduced using a tissue composed of adult ventricular cardiomyocytes (hAdultV-CMs). The analysis of phases, currents and safety factor for propagation showed an increased size of the re-entry core when I_CaL was blocked as a result of depressed cellular excitability. The bigger wavefront curvature yielded reductions of 12.2%, 6.9%, and 4.2% in the frequency of the re-entry for hiPSC-CM cultures, virtual hiPSC-CM, and hAdultV-CM tissues, respectively. Furthermore, I_CaL block led to a 47.8% shortening of the vulnerable window for re-entry in the virtual hiPSC-CM tissue and to re-entry vanishment in hAdultV-CM tissue. The consistent behavior between in-vitro and in-silico hiPSC-CMs and between in-silico hiPSC-CMs and hAdultV-CMs evidences that virtual hiPSC-CM tissues are suitable for assessing cardiac efficacy, as done in the present study through the analysis of I_CaL block.

Electrophysiology of hiPSC-Cardiomyocytes Co-Cultured with HEK Cells Expressing the Inward Rectifier Channel

The immature electrophysiology of human-induced pluripotent stem cell-derived cardiomyocytes (hiCMs) complicates their use for therapeutic and pharmacological purposes. An insufficient inward rectifying current (I_K1) and the presence of a funny current (if) cause spontaneous electrical activity. This study tests the hypothesis that the co-culturing of hiCMs with a human embryonic kidney (HEK) cell-line expressing the Kir2.1 channel (HEK-I_K1) can generate an electrical syncytium with an adult-like cardiac electrophysiology. The mechanical activity of co-cultures using different HEK-I_K1:hiCM ratios was compared with co-cultures using wildtype (HEK–WT:hiCM) or hiCM alone on days 3–8 after plating. Only ratios of 1:3 and 1:1 showed a significant reduction in spontaneous rate at days 4 and 6, suggesting that I_K1 was influencing the electrophysiology. Detailed analysis at day 4 revealed an increased incidence of quiescent wells or sub-areas. Electrical activity showed a decreased action potential duration (APD) at 20% and 50%, but not at 90%, alongside a reduced amplitude of the aggregate AP signal. A computational model of the 1:1 co-culture replicates the electrophysiological effects of HEK–WT. The addition of the I_K1 conductance reduced the spontaneous rate and APD20, 50 and 90, and minor variation in the intercellular conductance caused quiescence. In conclusion, a 1:1 co-culture HEK-I_K1:hiCM caused changes in electrophysiology and spontaneous activity consistent with the integration of I_K1 into the electrical syncytium. However, the additional electrical effects of the HEK cell at 1:1 increased the possibility of electrical quiescence before sufficient I_K1 was integrated into the syncytium.

An in silico hiPSC-Derived Cardiomyocyte Model Built With Genetic Algorithm

The formulation of in silico biophysical models generally requires optimization strategies for reproducing experimentally observed phenomena. In electrophysiological modeling, robust nonlinear regressive methods are often crucial for guaranteeing high fidelity models. Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), though nascent, have proven to be useful in cardiac safety pharmacology, regenerative medicine, and in the implementation of patient-specific test benches for investigating inherited cardiac disorders. This study demonstrates the potency of heuristic techniques at formulating biophysical models, with emphasis on a hiPSC-CM model using a novel genetic algorithm (GA) recipe we proposed. The proposed GA protocol was used to develop a hiPSC-CM biophysical computer model by fitting mathematical formulations to experimental data for five ionic currents recorded in hiPSC-CMs. The maximum conductances of the remaining ionic channels were scaled based on recommendations from literature to accurately reproduce the experimentally observed hiPSC-CM action potential (AP) metrics. Near-optimal parameter fitting was achieved for the GA-fitted ionic currents. The resulting model recapitulated experimental AP parameters such as AP durations (APD₅₀, APD₇₅, and APD₉₀), maximum diastolic potential, and frequency of automaticity. The outcome of this work has implications for validating the biophysics of hiPSC-CMs in their use as viable substitutes for human cardiomyocytes, particularly in cardiac safety pharmacology and in the study of inherited cardiac disorders. This study presents a novel GA protocol useful for formulating robust numerical biophysical models. The proposed protocol is used to develop a hiPSC-CM model with implications for cardiac safety pharmacology.

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 Ca²⁺ transient duration (CTD₉₀) and occurrence of arrhythmic events. Then, the sensitivity of each model to different degrees of blockage of Na⁺ (I_Na), L-type Ca²⁺ (I_CaL), and rapid delayed rectifying K⁺ (I_Kr) currents was quantified. Finally, we compared the drug-induced effects on CTD₉₀ against the corresponding in vitro experiments.

Results: The observed CTD₉₀ 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% CTD₉₀ 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 I_Kr. 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 I_CaL and I_Kr in Koivumäki2018 and Paci2020, while only by strong I_CaL blockers, e.g. nisoldipine, in Kernik2019. These results were confirmed by the sensitivity analysis.

Conclusion: To conclude, The CTD₉₀ 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.

Clinical Trial in a Dish: Using Patient-Derived Induced Pluripotent Stem Cells to Identify Risks of Drug-Induced Cardiotoxicity

Drug-induced cardiotoxicity is a significant clinical issue, with many drugs in the market being labeled with warnings on cardiovascular adverse effects. Treatments are often prematurely halted when cardiotoxicity is observed, which limits their therapeutic potential. Moreover, cardiotoxicity is a major reason for abandonment during drug development, reducing available treatment options for diseases and creating a significant financial burden and disincentive for drug developers. Thus, it is important to minimize the cardiotoxic effects of medications that are in use or in development. To this end, identifying patients at a higher risk of developing cardiovascular adverse effects for the drug of interest may be an effective strategy. The discovery of human induced pluripotent stem cells has enabled researchers to generate relevant cell types that retain a patient's own genome and examine patient-specific disease mechanisms, paving the way for precision medicine. Combined with the rapid development of pharmacogenomic analysis, the ability of induced pluripotent stem cell-derivatives to recapitulate patient-specific drug responses provides a powerful platform to identify subsets of patients who are particularly vulnerable to drug-induced cardiotoxicity. In this review, we will discuss the current use of patient-specific induced pluripotent stem cells in identifying populations who are at risk to drug-induced cardiotoxicity and their potential applications in future precision medicine practice. Graphic Abstract: A graphic abstract is available for this article.

In vitro and In silico Models to Study SARS-CoV-2 Infection: Integrating Experimental and Computational Tools to Mimic "COVID-19 Cardiomyocyte"

The rapid dissemination of SARS-CoV-2 has made COVID-19 a tremendous social, economic, and health burden. Despite the efforts to understand the virus and treat the disease, many questions remain unanswered about COVID-19 mechanisms of infection and progression. Severe Acute Respiratory Syndrome (SARS) infection can affect several organs in the body including the heart, which can result in thromboembolism, myocardial injury, acute coronary syndromes, and arrhythmias. Numerous cardiac adverse events, from cardiomyocyte death to secondary effects caused by exaggerated immunological response against the virus, have been clinically reported. In addition to the disease itself, repurposing of treatments by using "off label" drugs can also contribute to cardiotoxicity. Over the past several decades, animal models and more recently, stem cell-derived cardiomyocytes have been proposed for studying diseases and testing treatments in vitro. In addition, mechanistic in silico models have been widely used for disease and drug studies. In these models, several characteristics such as gender, electrolyte imbalance, and comorbidities can be implemented to study pathophysiology of cardiac diseases and to predict cardiotoxicity of drug treatments. In this Mini Review, we (1) present the state of the art of in vitro and in silico cardiomyocyte modeling currently in use to study COVID-19, (2) review in vitro and in silico models that can be adopted to mimic the effects of SARS-CoV-2 infection on cardiac function, and (3) provide a perspective on how to combine some of these models to mimic "COVID-19 cardiomyocytes environment.".

Evolution of mathematical models of cardiomyocyte electrophysiology

Advanced computational techniques and mathematical modeling have become more and more important to the study of cardiac electrophysiology. In this review, we provide a brief history of the evolution of cardiomyocyte electrophysiology models and highlight some of the most important ones that had a major impact on our understanding of the electrical activity of the myocardium and associated transmembrane ion fluxes in normal and pathological states. We also present the use of these models in the study of various arrhythmogenesis mechanisms, particularly the integration of experimental pharmacology data into advanced humanized models for in silico proarrhythmogenic risk prediction as an essential component of the Comprehensive in vitro Proarrhythmia Assay (CiPA) drug safety paradigm.

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.

Cardiac transmembrane ion channels and action potentials: cellular physiology and arrhythmogenic behavior

Cardiac arrhythmias are among the leading causes of mortality. They often arise from alterations in the electrophysiological properties of cardiac cells and their underlying ionic mechanisms. It is therefore critical to further unravel the pathophysiology of the ionic basis of human cardiac electrophysiology in health and disease. In the first part of this review, current knowledge on the differences in ion channel expression and properties of the ionic processes that determine the morphology and properties of cardiac action potentials and calcium dynamics from cardiomyocytes in different regions of the heart are described. Then the cellular mechanisms promoting arrhythmias in congenital or acquired conditions of ion channel function (electrical remodeling) are discussed. The focus is on human-relevant findings obtained with clinical, experimental, and computational studies, given that interspecies differences make the extrapolation from animal experiments to human clinical settings difficult. Deepening the understanding of the diverse pathophysiology of human cellular electrophysiology will help in developing novel and effective antiarrhythmic strategies for specific subpopulations and disease conditions.

Using cardiac ionic cell models to interpret clinical data

For over 100 years cardiac electrophysiology has been measured in the clinic. The electrical signals that can be measured span from noninvasive ECG and body surface potentials measurements through to detailed invasive measurements of local tissue electrophysiology. These electrophysiological measurements form a crucial component of patient diagnosis and monitoring; however, it remains challenging to quantitatively link changes in clinical electrophysiology measurements to biophysical cellular function. Multi-scale biophysical computational models represent one solution to this problem. These models provide a formal framework for linking cellular function through to emergent whole organ function and routine clinical diagnostic signals. In this review, we describe recent work on the use of computational models to interpret clinical electrophysiology signals. We review the simulation of human cardiac myocyte electrophysiology in the atria and the ventricles and how these models are being used to link organ scale function to patient disease mechanisms and therapy response in patients receiving implanted defibrillators, \cardiac resynchronisation therapy or suffering from atrial fibrillation and ventricular tachycardia. There is a growing use of multi-scale biophysical models to interpret clinical data. This allows cardiologists to link clinical observations with cellular mechanisms to better understand cardiopathophysiology and identify novel treatment strategies. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.

A greedy classifier optimization strategy to assess ion channel blocking activity and pro-arrhythmia in hiPSC-cardiomyocytes

Novel studies conducting cardiac safety assessment using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are promising but might be limited by their specificity and predictivity. It is often challenging to correctly classify ion channel blockers or to sufficiently predict the risk for Torsade de Pointes (TdP). In this study, we developed a method combining in vitro and in silico experiments to improve machine learning approaches in delivering fast and reliable prediction of drug-induced ion-channel blockade and proarrhythmic behaviour. The algorithm is based on the construction of a dictionary and a greedy optimization, leading to the definition of optimal classifiers. Finally, we present a numerical tool that can accurately predict compound-induced pro-arrhythmic risk and involvement of sodium, calcium and potassium channels, based on hiPSC-CM field potential data.

Combining tissue engineering and optical imaging approaches to explore interactions along the neuro-cardiac axis

Interactions along the neuro-cardiac axis are being explored with regard to their involvement in cardiac diseases, including catecholaminergic polymorphic ventricular tachycardia, hypertension, atrial fibrillation, long QT syndrome and sudden death in epilepsy. Interrogation of the pathophysiology and pathogenesis of neuro-cardiac diseases in animal models present challenges resulting from species differences, phenotypic variation, developmental effects and limited availability of data relevant at both the tissue and cellular level. By contrast, tissue-engineered models containing cardiomyocytes and peripheral sympathetic and parasympathetic neurons afford characterization of cellular- and tissue-level behaviours while maintaining precise control over developmental conditions, cellular genotype and phenotype. Such approaches are uniquely suited to long-term, high-throughput characterization using optical recording techniques with the potential for increased translational benefit compared to more established techniques. Furthermore, tissue-engineered constructs provide an intermediary between whole animal/tissue experiments and in silico models. This paper reviews the advantages of tissue engineering methods of multiple cell types and optical imaging techniques for the characterization of neuro-cardiac diseases.

All-Optical Electrophysiology Refines Populations of In Silico Human iPSC-CMs for Drug Evaluation

High-throughput in vitro drug assays have been impacted by recent advances in human induced pluripotent stem cell-derived cardiomyocyte (hiPSC-CM) technology and by contact-free all-optical systems simultaneously measuring action potentials (APs) and Ca2+ transients (CaTrs). Parallel computational advances have shown that in silico simulations can predict drug effects with high accuracy. We combine these in vitro and in silico technologies and demonstrate the utility of high-throughput experimental data to refine in silico hiPSC-CM populations and to predict and explain drug action mechanisms. Optically obtained hiPSC-CM APs and CaTrs were used from spontaneous activity and under optical pacing in control and drug conditions at multiple doses. An updated version of the Paci2018 model was developed to refine the description of hiPSC-CM spontaneous electrical activity; a population of in silico hiPSC-CMs was constructed and calibrated using simultaneously recorded APs and CaTrs. We tested in silico five drugs (astemizole, dofetilide, ibutilide, bepridil, and diltiazem) and compared the outcomes to in vitro optical recordings. Our simulations showed that physiologically accurate population of models can be obtained by integrating AP and CaTr control records. Thus, constructed population of models correctly predicted the drug effects and occurrence of adverse episodes, even though the population was optimized only based on control data and in vitro drug testing data were not deployed during its calibration. Furthermore, the in silico investigation yielded mechanistic insights; e.g., through simulations, bepridil's more proarrhythmic action in adult cardiomyocytes compared to hiPSC-CMs could be traced to the different expression of ion currents in the two. Therefore, our work 1) supports the utility of all-optical electrophysiology in providing high-content data to refine experimentally calibrated populations of in silico hiPSC-CMs, 2) offers insights into certain limitations when translating results obtained in hiPSC-CMs to humans, and 3) shows the strength of combining high-throughput in vitro and population in silico approaches.

SMN-deficiency disrupts SERCA2 expression and intracellular Ca2+ signaling in cardiomyocytes from SMA mice and patient-derived iPSCs

Spinal muscular atrophy (SMA) is a neurodegenerative disease characterized by loss of alpha motor neurons and skeletal muscle atrophy. The disease is caused by mutations of the SMN1 gene that result in reduced functional expression of survival motor neuron (SMN) protein. SMN is ubiquitously expressed, and there have been reports of cardiovascular dysfunction in the most severe SMA patients and animal models of the disease. In this study, we directly assessed the function of cardiomyocytes isolated from a severe SMA model mouse and cardiomyocytes generated from patient-derived IPSCs. Consistent with impaired cardiovascular function at the very early disease stages in mice, heart failure markers such as brain natriuretic peptide were significantly elevated. Functionally, cardiomyocyte relaxation kinetics were markedly slowed and the T₅₀ for Ca²⁺ sequestration increased to 146 ± 4 ms in SMN-deficient cardiomyocytes from 126 ± 4 ms in wild type cells. Reducing SMN levels in cardiomyocytes from control patient IPSCs slowed calcium reuptake similar to SMA patent-derived cardiac cells. Importantly, restoring SMN increased calcium reuptake rate. Taken together, these results indicate that SMN deficiency impairs cardiomyocyte function at least partially through intracellular Ca²⁺ cycling dysregulation.

Simulation of the Effects of Extracellular Calcium Changes Leads to a Novel Computational Model of Human Ventricular Action Potential With a Revised Calcium Handling

The importance of electrolyte concentrations for cardiac function is well established. Electrolyte variations can lead to arrhythmias onset, due to their important role in the action potential (AP) genesis and in maintaining cell homeostasis. However, most of the human AP computer models available in literature were developed with constant electrolyte concentrations, and fail to simulate physiological changes induced by electrolyte variations. This is especially true for Ca2+, even in the O'Hara-Rudy model (ORd), one of the most widely used models in cardiac electrophysiology. Therefore, the present work develops a new human ventricular model (BPS2020), based on ORd, able to simulate the inverse dependence of AP duration (APD) on extracellular Ca2+ concentration ([Ca2+]o), and APD rate dependence at 4 mM extracellular K+. The main changes needed with respect to ORd are: (i) an increased sensitivity of L-type Ca2+ current inactivation to [Ca2+]o; (ii) a single compartment description of the sarcoplasmic reticulum; iii) the replacement of Ca2+ release. BPS2020 is able to simulate the physiological APD-[Ca2+]o relationship, while also retaining the well-reproduced properties of ORd (APD rate dependence, restitution, accommodation and current block effects). We also used BPS2020 to generate an experimentally-calibrated population of models to investigate: (i) the occurrence of repolarization abnormalities in response to hERG current block; (ii) the rate adaptation variability; (iii) the occurrence of alternans and delayed after-depolarizations at fast pacing. Our results indicate that we successfully developed an improved version of ORd, which can be used to investigate electrophysiological changes and pro-arrhythmic abnormalities induced by electrolyte variations and current block at multiple rates and at the population level.

Multi-lineage Human iPSC-Derived Platforms for Disease Modeling and Drug Discovery

Human induced pluripotent stem cells (hiPSCs) provide a powerful platform for disease modeling and have unlocked new possibilities for understanding the mechanisms governing human biology, physiology, and genetics. However, hiPSC-derivatives have traditionally been utilized in two-dimensional monocultures, in contrast to the multi-systemic interactions that influence cells in the body. We will discuss recent advances in generating more complex hiPSC-based systems using three-dimensional organoids, tissue-engineering, microfluidic organ-chips, and humanized animal systems. While hiPSC differentiation still requires optimization, these next-generation multi-lineage technologies can augment the biomedical researcher's toolkit and enable more realistic models of human tissue function.

Improved Computational Identification of Drug Response Using Optical Measurements of Human Stem Cell Derived Cardiomyocytes in Microphysiological Systems

Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) hold great potential for drug screening applications. However, their usefulness is limited by the relative immaturity of the cells' electrophysiological properties as compared to native cardiomyocytes in the adult human heart. In this work, we extend and improve on methodology to address this limitation, building on previously introduced computational procedures which predict drug effects for adult cells based on changes in optical measurements of action potentials and Ca2+ transients made in stem cell derived cardiac microtissues. This methodology quantifies ion channel changes through the inversion of data into a mathematical model, and maps this response to an adult phenotype through the assumption of functional invariance of fundamental intracellular and membrane channels during maturation. Here, we utilize an updated action potential model to represent both hiPSC-CMs and adult cardiomyocytes, apply an IC50-based model of dose-dependent drug effects, and introduce a continuation-based optimization algorithm for analysis of dose escalation measurements using five drugs with known effects. The improved methodology can identify drug induced changes more efficiently, and quantitate important metrics such as IC50 in line with published values. Consequently, the updated methodology is a step towards employing computational procedures to elucidate drug effects in adult cardiomyocytes for new drugs using stem cell-derived experimental tissues.

General Principles for the Validation of Proarrhythmia Risk Prediction Models: An Extension of the CiPA In Silico Strategy

This white paper presents principles for validating proarrhythmia risk prediction models for regulatory use as discussed at the In Silico Breakout Session of a Cardiac Safety Research Consortium/Health and Environmental Sciences Institute/US Food and Drug Administration-sponsored Think Tank Meeting on May 22, 2018. The meeting was convened to evaluate the progress in the development of a new cardiac safety paradigm, the Comprehensive in Vitro Proarrhythmia Assay (CiPA). The opinions regarding these principles reflect the collective views of those who participated in the discussion of this topic both at and after the breakout session. Although primarily discussed in the context of in silico models, these principles describe the interface between experimental input and model-based interpretation and are intended to be general enough to be applied to other types of nonclinical models for proarrhythmia assessment. This document was developed with the intention of providing a foundation for more consistency and harmonization in developing and validating different models for proarrhythmia risk prediction using the example of the CiPA paradigm.

Mutation-specific differences in arrhythmias and drug responses in CPVT patients: simultaneous patch clamp and video imaging of iPSC derived cardiomyocytes

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is an inherited cardiac disease characterized by arrhythmias under adrenergic stress. Mutations in the cardiac ryanodine receptor (RYR2) are the leading cause for CPVT. We characterized electrophysiological properties of CPVT patient-specific induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) carrying different mutations in RYR2 and evaluated effects of carvedilol and flecainide on action potential (AP) and contractile properties of hiPSC-CMs. iPSC-CMs were generated from skin biopsies of CPVT patients carrying exon 3 deletion (E3D) and L4115F mutation in RYR2. APs and contractile movement were recorded simultaneously from the same hiPSC-CMs. Differences in AP properties of ventricular like CMs were seen in CPVT and control CMs: APD90 of both E3D (n = 20) and L4115F (n = 25) CPVT CMs was shorter than in control CMs (n = 15). E3D-CPVT CMs had shortest AP duration, lowest AP amplitude, upstroke velocity and more depolarized diastolic potential than controls. Adrenaline had positive and carvedilol and flecainide negative chronotropic effect in all hiPSC CMs. CPVT CMs had increased amount of delayed after depolarizations (DADs) and early after depolarizations (EADs) after adrenaline exposure. E3D CPVT CMs had the most DADs, EADs, and tachyarrhythmia. Discordant negatively coupled alternans was seen in L4115F CPVT CMs. Carvedilol cured almost all arrhythmias in L4115F CPVT CMs. Both drugs decreased contraction amplitude in all hiPSC CMs. E3D CPVT CMs have electrophysiological properties, which render them more prone to arrhythmias. iPSC-CMs provide a unique platform for disease modeling and drug screening for CPVT. Combining electrophysiological measurements, we can gain deeper insight into mechanisms of arrhythmias.

Stem Cells in Cardiovascular Medicine: Historical Overview and Future Prospects

Cardiovascular diseases remain the leading cause of death in the developed world, accounting for more than 30% of all deaths. In a large proportion of these patients, acute myocardial infarction is usually the first manifestation, which might further progress to heart failure. In addition, the human heart displays a low regenerative capacity, leading to a loss of cardiomyocytes and persistent tissue scaring, which entails a morbid pathologic sequela. Novel therapeutic approaches are urgently needed. Stem cells, such as induced pluripotent stem cells or embryonic stem cells, exhibit great potential for cell-replacement therapy and an excellent tool for disease modeling, as well as pharmaceutical screening of novel drugs and their cardiac side effects. This review article covers not only the origin of stem cells but tries to summarize their translational potential, as well as potential risks and clinical translation.

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 I_K1. Such immaturity hampers their application in assessing drug safety. The electronic overexpression of I_K1 (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 I_K1 is needed to bring hiPSC-CMs to a stable and hyperpolarized RMP and which mathematical description of I_K1 is most suitable for DC experiments. We compared five mature I_K1 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, APD₅₀, APD₉₀ (AP duration at 50 and 90% of repolarization), and APD₅₀/APD₉₀) after including them in the hiPSC-CM mathematical model by Paci. Then, we automatically identified the critical conductance for I_K1 ( G_{K1, critical}), the minimally required amount of I_K1 suppressing spontaneous activity. Preconditioning the cell model with depolarizing/hyperpolarizing prepulses allowed us to highlight time dependency of the I_K1 formulations. Simulations showed that inclusion of mature I_K1 formulations resulted in hyperpolarized RMP and higher V˙max, and observed G_{K1, critical} and the effect on AP duration strongly depended on I_K1 formulation. Finally, the Ishihara I_K1 led to shorter (−16.3%) and prolonged (+6.5%) APD₉₀ in response to hyperpolarizing and depolarizing prepulses, respectively, whereas other models showed negligible effects. Fine-tuning of G_K1 is an important step in DC experiments. Our computational work proposes a procedure to automatically identify how much I_K1 current is required to inject to stop the spontaneous activity and suggests the use of the Ishihara I_K1 model to perform DC experiments in hiPSC-CMs.

Comparison of in vitro and computational experiments on the relation of inter-beat interval and duration of repolarization in a specific type of human induced pluripotent stem cell-derived cardiomyocytes

We compared a published computational model of the action potential of a specific type of human induced pluripotent stem cell -derived cardiomyocytes (hiPSC-CM) with experimental field potential data with regard to their inter-beat interval and the duration of repolarization. In particular, concomitant changes in inter-beat interval and duration of repolarization were calculated after reduction and/or augmentation of specific ion channel conductances as a surrogate for pharmacological manipulation. The observed mismatches between calculations and experimental data indicate that there is information missing about the cellular test system. Based on our results we hypothesize that, among other currents, the actual I_f (“funny current”) may deviate from the prediction. We show that replacement of the I_f formulation by alternative equations causes the model predictions to change qualitatively, however, none of the available formulations is actually achieving a satisfactory match with experimental data. We suggest a strategy to clarify whether the mismatch can be completely resolved at all using single cell models and, if yes, how this goal could be reached.

The influence of hERG1a and hERG1b isoforms on drug safety screening in iPSC-CMs

The human Ether-à-go-go Related Gene (hERG) encodes the pore forming subunit of the channel that conducts the rapid delayed rectifier potassium current I_Kr. I_Kr drives repolarization in the heart and when I_Kr is dysfunctional, cardiac repolarization delays, the QT interval on the electrocardiogram (ECG) prolongs and the risk of developing lethal arrhythmias such as Torsade de Pointes (TdP) increases. TdP risk is incorporated in drug safety screening for cardiotoxicity where hERG is the main target since the I_Kr channels appear highly sensitive to blockage. hERG block is also included as an important read-out in the Comprehensive in Vitro Proarrhythmia Assay (CiPA) initiative which aims to combine in vitro and in silico experiments on induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) to screen for cardiotoxicity. However, the hERG channel has some unique features to consider for drug safety screening, which we will discuss in this study. The hERG channel consists of different isoforms, hERG1a and hERG1b, which individually influence the kinetics of the channel and the drug response in the human heart and in iPSC-CMs. hERG1b is often underappreciated in iPSC-CM studies, drug screening assays and in silico models, and the fact that its contribution might substantially differ between iPSC-CM and healthy but also diseased human heart, adds to this problem. In this study we show that the activation kinetics in iPSC-CMs resemble hERG1b kinetics using Cs⁺ as a charge carrier. Not including hERG1b in drug safety testing might underestimate the actual role of hERG1b in repolarization and drug response, and might lead to inappropriate conclusions. We stress to focus more on including hERG1b in drug safety testing concerning I_Kr.

Substrate Stiffness, Cell Anisotropy, and Cell-Cell Contact Contribute to Enhanced Structural and Calcium Handling Properties of Human Embryonic Stem Cell-Derived Cardiomyocytes

Human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs) can be utilized to understand the mechanisms underlying the development and progression of heart disease, as well as to develop better interventions and treatments for this disease. However, these cells are structurally and functionally immature, which undermines some of their adequacy in modeling adult heart tissue. Previous studies with immature cardiomyocytes have shown that altering substrate stiffness, cell anisotropy, and/or cell-cell contact can enhance the contractile and structural maturation of hPSC-CMs. In this study, the structural and calcium handling properties of human embryonic stem cell-derived cardiomyocytes (hESC-CMs) were enhanced by exposure to a downselected combination of these three maturation stimuli. First, hESC-CMs were seeded onto substrates composed of two commercial formulations of polydimethylsiloxane (PDMS), Sylgard 184 and Sylgard 527, whose stiffness ranged from 5 kPa to 101 kPa. Upon analyzing the morphological and calcium transient properties of these cells, it was concluded that a 21 kPa substrate yielded cells with the highest degree of maturation. Next, these PDMS substrates were microcontact-printed with laminin to force the cultured cells into rod-shaped geometries using line patterns that were 12, 18, or 24 μm in width. We found that cells on the 18 and 24 μm pattern widths had structural and functional properties that were superior to those on the 12 μm pattern. The hESC-CMs were then seeded onto these line-stamped surfaces at a density of 500 000 cells per 25-mm-diameter substrate, to enable the formation of cell-cell contacts at their distal ends. We discovered that this combination of culture conditions resulted in cells that were more structurally and functionally mature than those that were only exposed to one or two stimuli. Our results suggest that downselecting a combination of mechanobiological stimuli could prove to be an effective means of maturing hPSC-CMs in vitro.

Early afterdepolarizations in cardiac action potentials as mixed mode oscillations due to a folded node singularity

Early afterdepolarizations (EADs) are pathological voltage oscillations during the repolarization phase of cardiac action potentials. They are considered as potential precursors to cardiac arrhythmias and have recently gained much attention in the context of preclinical drug safety testing under the Comprehensive in vitro Proarrhythmia Assay (CiPA) paradigm. From the viewpoint of multiple time scales theory, the onset of EADs has previously been studied by means of mathematical action potential models with one slow ion channel gating variable. In this article, we for the first time associate EADs with mixed mode oscillations in dynamical systems with two slow gating variables and present a folded node singularity of the slow flow as a novel mechanism for EADs genesis. We derive regions of the pharmacology parameter space in which EADs occur using both the folded node analysis and a full system bifurcation analysis, and we suggest the normal distance to the boundary of the EADs region as a mechanism-based risk metric to computationally estimate a drug’s proarrhythmic liability.

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.