A computational model of induced pluripotent stem-cell derived cardiomyocytes for high throughput risk stratification of KCNQ1 genetic variants

Abstract

In the last decade, there has been tremendous progress in identifying genetic anomalies linked to clinical disease. New experimental platforms have connected genetic variants to mechanisms underlying disruption of cellular and organ behavior and the emergence of proarrhythmic cardiac phenotypes. The development of induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) signifies an important advance in the study of genetic disease in a patient-specific context. However, considerable limitations of iPSC-CM technologies have not been addressed: 1) phenotypic variability in apparently identical genotype perturbations, 2) low-throughput electrophysiological measurements, and 3) an immature phenotype which may impact translation to adult cardiac response. We have developed a computational approach intended to address these problems. We applied our recent iPSC-CM computational model to predict the proarrhythmic risk of 40 KCNQ1 genetic variants. An I _Ks computational model was fit to experimental data for each mutation, and the impact of each mutation was simulated in a population of iPSC-CM models. Using a test set of 15 KCNQ1 mutations with known clinical long QT phenotypes, we developed a method to stratify the effects of KCNQ1 mutations based on proarrhythmic markers. We utilized this method to predict the severity of the remaining 25 KCNQ1 mutations with unknown clinical significance. Tremendous phenotypic variability was observed in the iPSC-CM model population following mutant perturbations. A key novelty is our reporting of the impact of individual KCNQ1 mutant models on adult ventricular cardiomyocyte electrophysiology, allowing for prediction of mutant impact across the continuum of aging. This serves as a first step toward translating predicted response in the iPSC-CM model to predicted response of the adult ventricular myocyte given the same genetic mutation. As a whole, this study presents a new computational framework that serves as a high throughput method to evaluate risk of genetic mutations based-on proarrhythmic behavior in phenotypically variable populations.

Author summary

In the last decade, there has been tremendous progress in identifying genetic mutations linked to clinical diseases, such as cardiac arrhythmia. Many experimental platforms have been developed to study this link, including induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). IPSC-CMs are patient-derived cardiac cells which allow for the study of genetic variants within a patient-specific context. However, experimentally iPSC-CMs have certain limitations, including: (1) they exhibit variability in behavior within cells that are apparently genetically identical, and (2) they are immature compared to adult cardiac cells. In our study, we have developed a computational approach to model 40 genetic variants in the KCNQ1 gene and predict the proarrhythmic risk of each variant. To do this, we modeled the ionic current determined by KCNQ1, I _Ks, to fit experimental data for each mutation. We then simulated the impact of each mutation in a population of iPSC-CMs, incorporating variability across the population. We also simulated each variant in an adult cardiac cell model, providing a link between iPSC-CM response to mutants and adult cardiac cell response to the same mutants. Overall, this study provides a new computational framework to evaluate risk of genetic mutations based-on proarrhythmic behavior diverse populations of iPSC-CM models.