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

Enhanced drug classification using machine learning with multiplexed cardiac contractility assays

Cardiac screening of newly discovered drugs remains a longstanding challenge for the pharmaceutical industry. While therapeutic efficacy and cardiotoxicity are evaluated through preclinical biochemical and animal testing, 90 % of lead compounds fail to meet safety and efficacy benchmarks during human clinical trials. A preclinical model more representative of the human cardiac response is needed; heart tissue engineered from human pluripotent stem cell derived cardiomyocytes offers such a platform. In this study, three functionally distinct and independently validated engineered cardiac tissue assays are exposed to increasing concentrations of known compounds representing 5 classes of mechanistic action, creating a robust electrophysiology and contractility dataset. Combining results from six individual models, the resulting ensemble algorithm can classify the mechanistic action of unknown compounds with 86.2 % predictive accuracy. This outperforms single-assay models and offers a strategy to enhance future clinical trial success aligned with the recent FDA Modernization Act 2.0.

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 to predict 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 in silico datasets by simulating various protocols applied to a population of models with known conductance variations, and then fitted parameters to those datasets. We found that calibrating 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 also held when 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.

The Current State of Realistic Heart Models for Disease Modelling and Cardiotoxicity

One of the many unresolved obstacles in the field of cardiovascular research is an uncompromising in vitro cardiac model. While primary cell sources from animal models offer both advantages and disadvantages, efforts over the past half-century have aimed to reduce their use. Additionally, obtaining a sufficient quantity of human primary cardiomyocytes faces ethical and legal challenges. As the practically unlimited source of human cardiomyocytes from induced pluripotent stem cells (hiPSC-CM) is now mostly resolved, there are great efforts to improve their quality and applicability by overcoming their intrinsic limitations. The greatest bottleneck in the field is the in vitro ageing of hiPSC-CMs to reach a maturity status that closely resembles that of the adult heart, thereby allowing for more appropriate drug developmental procedures as there is a clear correlation between ageing and developing cardiovascular diseases. Here, we review the current state-of-the-art techniques in the most realistic heart models used in disease modelling and toxicity evaluations from hiPSC-CM maturation through heart-on-a-chip platforms and in silico models to the in vitro models of certain cardiovascular diseases.

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.

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.

Identification and characterization of two novel KCNH2 mutations contributing to long QT syndrome

We identified two different inherited mutations in KCNH2 gene, or human ether-a-go-go related gene (hERG), which are linked to Long QT Syndrome. The first mutation was in a 1-day-old infant, whereas the second was in a 14-year-old girl. The two KCNH2 mutations were transiently transfected into either human embryonic kidney (HEK) cells or human induced pluripotent stem-cell derived cardiomyocytes. We performed associated multiscale computer simulations to elucidate the arrhythmogenic potentials of the KCNH2 mutations. Genetic screening of the first and second index patients revealed a heterozygous missense mutation in KCNH2, resulting in an amino acid change (P632L) in the outer loop of the channel and substitution at position 428 from serine to proline (S428P), respectively. Heterologous expression of P632L and S428P into HEK cells produced no hERG current compared to the wild type (WT). Moreover, the co-transfection of WT and P632L yielded no hERG current; however, the co-transfection of WT and S428P yielded partial hERG current. Action potentials were prolonged in a complete or partial blockade of hERG current from computer simulations which was more severe in Purkinje than ventricular myocytes. Three dimensional simulations revealed a higher susceptibility to reentry in the presence of hERG current blockade. Our experimental findings suggest that both P632L and S428P mutations may impair the KCNH2 gene. The Purkinje cells exhibit a more severe phenotype than ventricular myocytes, and the hERG current blockade renders the ventricles an arrhythmogenic substrate from computer modeling.

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.

Rapid ionic current phenotyping (RICP) identifies mechanistic underpinnings of iPSC-CM AP heterogeneity

As a renewable, easily accessible, human-derived in vitro model, human induced pluripotent stem cell derived cardiomyocytes (iPSC-CMs) are a promising tool for studying arrhythmia-related factors, including cardiotoxicity and congenital proarrhythmia risks. An oft-mentioned limitation of iPSC-CMs is the abundant cell-to-cell variability in recordings of their electrical activity. Here, we develop a new method, rapid ionic current phenotyping (RICP), that utilizes a short (10 s) voltage clamp protocol to quantify cell-to-cell heterogeneity in key ionic currents. We correlate these ionic current dynamics to action potential recordings from the same cells and produce mechanistic insights into cellular heterogeneity. We present evidence that the L-type calcium current is the main determinant of upstroke velocity, rapid delayed rectifier K+ current is the main determinant of the maximal diastolic potential, and an outward current in the excitable range of slow delayed rectifier K+ is the main determinant of action potential duration. We measure an unidentified outward current in several cells at 6 mV that is not recapitulated by iPSC-CM mathematical models but contributes to determining action potential duration. In this way, our study both quantifies cell-to-cell variability in membrane potential and ionic currents, and demonstrates how the ionic current variability gives rise to action potential heterogeneity. Based on these results, we argue that iPSC-CM heterogeneity should not be viewed simply as a problem to be solved but as a model system to understand the mechanistic underpinnings of cellular variability.

Computational design of custom therapeutic cells to correct failing human cardiomyocytes

Background

Myocardial delivery of non-excitable cells-namely human mesenchymal stem cells (hMSCs) and c-kit⁺ cardiac interstitial cells (hCICs)-remains a promising approach for treating the failing heart. Recent empirical studies attempt to improve such therapies by genetically engineering cells to express specific ion channels, or by creating hybrid cells with combined channel expression. This study uses a computational modeling approach to test the hypothesis that custom hypothetical cells can be rationally designed to restore a healthy phenotype when coupled to human heart failure (HF) cardiomyocytes.

Methods

Candidate custom cells were simulated with a combination of ion channels from non-excitable cells and healthy human cardiomyocytes (hCMs). Using a genetic algorithm-based optimization approach, candidate cells were accepted if a root mean square error (RMSE) of less than 50% relative to healthy hCM was achieved for both action potential and calcium transient waveforms for the cell-treated HF cardiomyocyte, normalized to the untreated HF cardiomyocyte.

Results

Custom cells expressing only non-excitable ion channels were inadequate to restore a healthy cardiac phenotype when coupled to either fibrotic or non-fibrotic HF cardiomyocytes. In contrast, custom cells also expressing cardiac ion channels led to acceptable restoration of a healthy cardiomyocyte phenotype when coupled to fibrotic, but not non-fibrotic, HF cardiomyocytes. Incorporating the cardiomyocyte inward rectifier K⁺ channel was critical to accomplishing this phenotypic rescue while also improving single-cell action potential metrics associated with arrhythmias, namely resting membrane potential and action potential duration. The computational approach also provided insight into the rescue mechanisms, whereby heterocellular coupling enhanced cardiomyocyte L-type calcium current and promoted calcium-induced calcium release. Finally, as a therapeutically translatable strategy, we simulated delivery of hMSCs and hCICs genetically engineered to express the cardiomyocyte inward rectifier K⁺ channel, which decreased action potential and calcium transient RMSEs by at least 24% relative to control hMSCs and hCICs, with more favorable single-cell arrhythmia metrics.

Conclusion

Computational modeling facilitates exploration of customizable engineered cell therapies. Optimized cells expressing cardiac ion channels restored healthy action potential and calcium handling phenotypes in fibrotic HF cardiomyocytes and improved single-cell arrhythmia metrics, warranting further experimental validation studies of the proposed custom therapeutic cells.