Computational prediction of drug response in short QT syndrome type 1 based on measurements of compound effect in stem cell-derived cardiomyocytes

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.

A novel ionic model for matured and paced atrial-like human iPSC-CMs integrating IKur and IKCa currents

This work introduces the first atrial-specific in-silico human induced pluripotent stem cells-derived cardiomyocytes (hiPSC-CMs) model, based on a set of phenotype-specific IKur,IKCa and IK1 membrane currents. This model is built on novel in-vitro experimental data recently published by some of the co-authors to simulate the paced action potential of matured atrial-like hiPSC-CMs. The model consists of a system of stiff ordinary differential equations depending on several parameters, which have been tuned by automatic optimization techniques to closely match selected experimental biomarkers. The new model effectively simulates the electronic in-vitro hiPSC-CMs maturation process, transitioning from an unstable depolarized membrane diastolic potential to a stable hyperpolarized resting potential, and exhibits spontaneous firing activity in unpaced conditions. Moreover, our model accurately reflects the experimental rate dependence data at different cycle length and demonstrates the expected response to a specific current blocker. This atrial-specific in-silico model provides a novel computational tool for electrophysiological studies of cardiac stem cells and their applications to drug evaluation and atrial fibrillation treatment.

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.

Impacts of gene variants on drug effects-the foundation of genotype-guided pharmacologic therapy for long QT syndrome and short QT syndrome

The clinical significance of optimal pharmacotherapy for inherited arrhythmias such as short QT syndrome (SQTS) and long QT syndrome (LQTS) has been increasingly recognised. The advancement of gene technology has opened up new possibilities for identifying genetic variations and investigating the pathophysiological roles and mechanisms of genetic arrhythmias. Numerous variants in various genes have been proven to be causative in genetic arrhythmias. Studies have demonstrated that the effectiveness of certain drugs is specific to the patient or genotype, indicating the important role of gene-variants in drug response. This review aims to summarize the reported data on the impact of different gene-variants on drug response in SQTS and LQTS, as well as discuss the potential mechanisms by which gene-variants alter drug response. These findings may provide valuable information for future studies on the influence of gene variants on drug efficacy and the development of genotype-guided or precision treatment for these diseases.

The simplified Kirchhoff network model (SKNM): a cell-based reaction-diffusion model of excitable tissue

Cell-based models of excitable tissues offer the advantage of cell-level precision, which cannot be achieved using traditional homogenized electrophysiological models. However, this enhanced accuracy comes at the cost of increased computational demands, necessitating the development of efficient cell-based models. The widely-accepted bidomain model serves as the standard in computational cardiac electrophysiology, and under certain anisotropy ratio conditions, it is well known that it can be reduced to the simpler monodomain model. Recently, the Kirchhoff Network Model (KNM) was developed as a cell-based counterpart to the bidomain model. In this paper, we aim to demonstrate that KNM can be simplified using the same steps employed to derive the monodomain model from the bidomain model. We present the cell-based Simplified Kirchhoff Network Model (SKNM), which produces results closely aligned with those of KNM while requiring significantly less computational resources.

Leak current, even with gigaohm seals, can cause misinterpretation of stem cell-derived cardiomyocyte action potential recordings

AIMS

Human-induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) have become an essential tool to study arrhythmia mechanisms. Much of the foundational work on these cells, as well as the computational models built from the resultant data, has overlooked the contribution of seal-leak current on the immature and heterogeneous phenotype that has come to define these cells. The aim of this study is to understand the effect of seal-leak current on recordings of action potential (AP) morphology.

METHODS AND RESULTS

Action potentials were recorded in human iPSC-CMs using patch clamp and simulated using previously published mathematical models. Our in silico and in vitro studies demonstrate how seal-leak current depolarizes APs, substantially affecting their morphology, even with seal resistances (Rseal) above 1 GΩ. We show that compensation of this leak current is difficult due to challenges with obtaining accurate measures of Rseal during an experiment. Using simulation, we show that Rseal measures (i) change during an experiment, invalidating the use of pre-rupture values, and (ii) are polluted by the presence of transmembrane currents at every voltage. Finally, we posit that the background sodium current in baseline iPSC-CM models imitates the effects of seal-leak current and is increased to a level that masks the effects of seal-leak current on iPSC-CMs.

CONCLUSION

Based on these findings, we make recommendations to improve iPSC-CM AP data acquisition, interpretation, and model-building. Taking these recommendations into account will improve our understanding of iPSC-CM physiology and the descriptive ability of models built from such data.

Computational analysis of arrhythmogenesis in KCNH2 T618I mutation-associated short QT syndrome and the pharmacological effects of quinidine and sotalol

Short QT syndrome (SQTS) is a rare but dangerous genetic disease. In this research, we conducted a comprehensive in silico investigation into the arrhythmogenesis in KCNH2 T618I-associated SQTS using a multi-scale human ventricle model. A Markov chain model of IKr was developed firstly to reproduce the experimental observations. It was then incorporated into cell, tissue, and organ models to explore how the mutation provided substrates for ventricular arrhythmias. Using this T618I Markov model, we explicitly revealed the subcellular level functional alterations by T618I mutation, particularly the changes of ion channel states that are difficult to demonstrate in wet experiments. The following tissue and organ models also successfully reproduced the changed dynamics of reentrant spiral waves and impaired rate adaptions in hearts of T618I mutation. In terms of pharmacotherapy, we replicated the different effects of a drug under various conditions using identical mathematical descriptions for drugs. This study not only simulated the actions of an effective drug (quinidine) at various physiological levels, but also elucidated why the IKr inhibitor sotalol failed in SQT1 patients through profoundly analyzing its mutation-dependent actions.

Metabolically driven maturation of human-induced-pluripotent-stem-cell-derived cardiac microtissues on microfluidic chips

The immature physiology of cardiomyocytes derived from human induced pluripotent stem cells (hiPSCs) limits their utility for drug screening and disease modelling. Here we show that suitable combinations of mechanical stimuli and metabolic cues can enhance the maturation of hiPSC-derived cardiomyocytes, and that the maturation-inducing cues have phenotype-dependent effects on the cells' action-potential morphology and calcium handling. By using microfluidic chips that enhanced the alignment and extracellular-matrix production of cardiac microtissues derived from genetically distinct sources of hiPSC-derived cardiomyocytes, we identified fatty-acid-enriched maturation media that improved the cells' mitochondrial structure and calcium handling, and observed divergent cell-source-dependent effects on action-potential duration (APD). Specifically, in the presence of maturation media, tissues with abnormally prolonged APDs exhibited shorter APDs, and tissues with aberrantly short APDs displayed prolonged APDs. Regardless of cell source, tissue maturation reduced variabilities in spontaneous beat rate and in APD, and led to converging cell phenotypes (with APDs within the 300-450 ms range characteristic of human left ventricular cardiomyocytes) that improved the modelling of the effects of pro-arrhythmic drugs on cardiac tissue.

Arrhythmogenic influence of mutations in a myocyte-based computational model of the pulmonary vein sleeve

In the heart, electrophysiological dysregulation arises from defects at many biological levels (from point mutations in ion channel proteins to gross structural abnormalities). These defects disrupt the normal pattern of electrical activation, producing ectopic activity and reentrant arrhythmia. To interrogate mechanisms that link these primary biological defects to macroscopic electrophysiologic dysregulation most prior computational studies have utilized either (i) detailed models of myocyte ion channel dynamics at limited spatial scales, or (ii) homogenized models of action potential conduction that reproduce arrhythmic activity at tissue and organ levels. Here we apply our recent model (EMI), which integrates electrical activation and propagation across these scales, to study human atrial arrhythmias originating in the pulmonary vein (PV) sleeves. These small structures initiate most supraventricular arrhythmias and include pronounced myocyte-to-myocyte heterogeneities in ion channel expression and intercellular coupling. To test EMI's cell-based architecture in this physiological context we asked whether ion channel mutations known to underlie atrial fibrillation are capable of initiating arrhythmogenic behavior via increased excitability or reentry in a schematic PV sleeve geometry. Our results illustrate that EMI's improved spatial resolution can directly interrogate how electrophysiological changes at the individual myocyte level manifest in tissue and as arrhythmia in the PV sleeve.

Integrative Computational Modeling of Cardiomyocyte Calcium Handling and Cardiac Arrhythmias: Current Status and Future Challenges

Cardiomyocyte calcium-handling is the key mediator of cardiac excitation-contraction coupling. In the healthy heart, calcium controls both electrical impulse propagation and myofilament cross-bridge cycling, providing synchronous and adequate contraction of cardiac muscles. However, calcium-handling abnormalities are increasingly implicated as a cause of cardiac arrhythmias. Due to the complex, dynamic and localized interactions between calcium and other molecules within a cardiomyocyte, it remains experimentally challenging to study the exact contributions of calcium-handling abnormalities to arrhythmogenesis. Therefore, multiscale computational modeling is increasingly being used together with laboratory experiments to unravel the exact mechanisms of calcium-mediated arrhythmogenesis. This article describes various examples of how integrative computational modeling makes it possible to unravel the arrhythmogenic consequences of alterations to cardiac calcium handling at subcellular, cellular and tissue levels, and discusses the future challenges on the integration and interpretation of such computational data.

A computational method for identifying an optimal combination of existing drugs to repair the action potentials of SQT1 ventricular myocytes

Mutations are known to cause perturbations in essential functional features of integral membrane proteins, including ion channels. Even restricted or point mutations can result in substantially changed properties of ion currents. The additive effect of these alterations for a specific ion channel can result in significantly changed properties of the action potential (AP). Both AP shortening and AP prolongation can result from known mutations, and the consequences can be life-threatening. Here, we present a computational method for identifying new drugs utilizing combinations of existing drugs. Based on the knowledge of theoretical effects of existing drugs on individual ion currents, our aim is to compute optimal combinations that can ‘repair’ the mutant AP waveforms so that the baseline AP-properties are restored. More specifically, we compute optimal, combined, drug concentrations such that the waveforms of the transmembrane potential and the cytosolic calcium concentration of the mutant cardiomyocytes (CMs) becomes as similar as possible to their wild type counterparts after the drug has been applied. In order to demonstrate the utility of this method, we address the question of computing an optimal drug for the short QT syndrome type 1 (SQT1). For the SQT1 mutation N588K, there are available data sets that describe the effect of various drugs on the mutated K⁺ channel. These published findings are the basis for our computational analysis which can identify optimal compounds in the sense that the AP of the mutant CMs resembles essential biomarkers of the wild type CMs. Using recently developed insights regarding electrophysiological properties among myocytes from different species, we compute optimal drug combinations for hiPSC-CMs, rabbit ventricular CMs and adult human ventricular CMs with the SQT1 mutation. Since the ‘composition’ of ion channels that form the AP is different for the three types of myocytes under consideration, so is the composition of the optimal drug.

Poly-pharmacology (using multiple drugs to treat disease) has been proposed for improving cardiac anti-arrhythmic therapy for at least two decades. However, the specific arrhythmia contexts in which polytherapy is likely to be both safe and effective have remained elusive. Type 1 short QT syndrome (SQT1) is a rare form of cardiac arrhythmia that results from mutations to the human Ether-á-go-go Related Gene (hERG) potassium channel. Functionally, these mutations are remarkably consistent in that they permit the channel to open earlier during each heart beat. While hundreds of compounds are known to inhibit hERG channels, the specific effect of SQT1 mutations that allows for early channel opening also limits the ability of most of those compounds to correct SQT1 dysfunction. Here, we have applied a suite of ventricular cardiomyocyte computational models to ask whether polytherapy may offer a more effective therapeutic strategy in SQT1, and if so, what the likely characteristics of that strategy are. Our analyses suggest that simultaneous induction of late sodium current and partial hERG blockade offers a promising strategy. While no activators of late sodium current have been clinically approved, several experimental compounds are available and may provide a basis for interrogating this strategy. The method presented here can be used to compute optimal drug combinations provided that the effect of each drug on every relevant ion channel is known.

Effects of Antiarrhythmic Drugs on hERG Gating in Human-Induced Pluripotent Stem Cell-Derived Cardiomyocytes From a Patient With Short QT Syndrome Type 1

Aims: The short QT syndrome type 1 (SQT1) is linked to hERG channel mutations (e.g., N588K). Drug effects on hERG channel gating kinetics in SQT1-cells have not been investigated. Methods: This study used hiPSC-CMs of a healthy donor and a SQT1-patient carrying the N588K mutation and patch clamp to examine the drug effects on hERG channel gating kinetics. Results: Ajmaline, amiodarone, ivabradine, flecainide, quinidine, mexiletine and ranolazine inhibited the hERG channel current (I_Kr) less strongly in hiPSC-CMs from the SQTS1-patient (SQT1-hiPSC-CMs) comparing with cells from the healthy donor (donor-hiPSC-CMs). Quinidine and mexiletine reduced, but ajmaline, amiodarone, ivabradine and ranolazine increased the time to peak of I_Kr similarly in SQT1-hiPSC-CMs and donor-hiPSC-CMs. Although regarding the shift of activation and inactivation curves, tested drugs showed differential effects in donor- and SQT1-hiPSC-CMs, quinidine, ajmaline, ivabradine and mexiletine but not amiodarone, flecainide and ranolazine reduced the window current in SQT1-hiPSC-CMs. Quinidine, ajmaline, ivabradine and mexiletine differentially changed the time constant of recovery from inactivation, but all of them increased the time constant of deactivation in SQT1-hiPSC-CMs. Conclusion: The window current-reducing and deactivation-slowing effects may be important for the antiarrhythmic effect of ajmaline, ivabradine, quinidine and mexiletine in SQT1-cells. This information may be helpful for selecting drugs for treating SQT1-patients with hERG channel mutation.