A computational model of induced pluripotent stem‐cell derived cardiomyocytes incorporating experimental variability from multiple data sources

KEY POINTS: Induced pluripotent stem cell‐derived cardiomyocytes (iPSC‐CMs) capture patient‐specific genotype–phenotype relationships, as well as cell‐to‐cell variability of cardiac electrical activity. Computational modelling and simulation provide a high throughput approach to reconcile multiple datasets describing physiological variability, and also identify vulnerable parameter regimes. We have developed a whole‐cell model of iPSC‐CMs, composed of single exponential voltage‐dependent gating variable rate constants, parameterized to fit experimental iPSC‐CM outputs. We have utilized experimental data across multiple laboratories to model experimental variability and investigate subcellular phenotypic mechanisms in iPSC‐CMs. This framework links molecular mechanisms to cellular‐level outputs by revealing unique subsets of model parameters linked to known iPSC‐CM phenotypes. ABSTRACT: There is a profound need to develop a strategy for predicting patient‐to‐patient vulnerability in the emergence of cardiac arrhythmia. A promising in vitro method to address patient‐specific proclivity to cardiac disease utilizes induced pluripotent stem cell‐derived cardiomyocytes (iPSC‐CMs). A major strength of this approach is that iPSC‐CMs contain donor genetic information and therefore capture patient‐specific genotype–phenotype relationships. A cited detriment of iPSC‐CMs is the cell‐to‐cell variability observed in electrical activity. We postulated, however, that cell‐to‐cell variability may constitute a strength when appropriately utilized in a computational framework to build cell populations that can be employed to identify phenotypic mechanisms and pinpoint key sensitive parameters. Thus, we have exploited variation in experimental data across multiple laboratories to develop a computational framework for investigating subcellular phenotypic mechanisms. We have developed a whole‐cell model of iPSC‐CMs composed of simple model components comprising ion channel models with single exponential voltage‐dependent gating variable rate constants, parameterized to fit experimental iPSC‐CM data for all major ionic currents. By optimizing ionic current model parameters to multiple experimental datasets, we incorporate experimentally‐observed variability in the ionic currents. The resulting population of cellular models predicts robust inter‐subject variability in iPSC‐CMs. This approach links molecular mechanisms to known cellular‐level iPSC‐CM phenotypes, as shown by comparing immature and mature subpopulations of models to analyse the contributing factors underlying each phenotype. In the future, the presented models can be readily expanded to include genetic mutations and pharmacological interventions for studying the mechanisms of rare events, such as arrhythmia triggers.