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.Author summary: Many drug treatments or environmental factors can trigger cardiac arrhythmias, which are dangerous and often unpredictable. Human cardiomyocytes derived from donor stem cells have proven to be a promising model for studying these events, but variability in donor genetic background and cell maturation methods, as well as overall immaturity of stem cell-derived cardiomyocytes relative to the adult heart, have hindered reproducibility and reliability of these studies. Mathematical models of these cells can aid in understanding the underlying electrophysiological contributors to this variability, but determining these models’ parameters for multiple cell preparations is challenging. In this study, we tackle these limitations by developing a computational method to simultaneously estimate multiple model parameters using data from imaging-based experiments, which can be easily scaled to rapidly characterize multiple cell lines. This method can generate many personalized models of individual cell preparations, improving drug response predictions and revealing specific differences in electrophysiological properties that contribute to variability in cardiac maturity and arrhythmia susceptibility.