An in silico-in vitro pipeline for drug cardiotoxicity screening identifies ionic pro-arrhythmia mechanisms

Resolving artefacts in voltage-clamp experiments with computational modelling: an application to fast sodium current recordings

Cellular electrophysiology is the foundation of many fields, from basic science in neurology, cardiology, oncology to safety critical applications for drug safety testing, clinical phenotyping, etc. Patch-clamp voltage clamp is the gold standard technique for studying cellular electrophysiology. Yet, the quality of these experiments is not always transparent, which may lead to erroneous conclusions for studies and applications. Here, we have developed a new computational approach that allows us to explain and predict the experimental artefacts in voltage-clamp experiments. The computational model captures the experimental procedure and its inadequacies, including: voltage offset, series resistance, membrane capacitance and (imperfect) amplifier compensations, such as series resistance compensation and supercharging. The computational model was validated through a series of electrical model cell experiments. Using this computational approach, the artefacts in voltage-clamp experiments of cardiac fast sodium current, one of the most challenging currents to voltage clamp, were able to be resolved and explained through coupling the observed current and the simulated membrane voltage, including some typically observed shifts and delays in the recorded currents. We further demonstrated that the typical way of averaging data for current-voltage relationships would lead to biases in the peak current and shifts in the peak voltage, and such biases can be in the same order of magnitude as those differences reported for disease-causing mutations. Therefore, the presented new computational pipeline will provide a new standard of assessing the voltage-clamp experiments and interpreting the experimental data, which may be able to rectify and provide a better understanding of ion channel mutations and other related applications.

Single-cell ionic current phenotyping elucidates non-canonical features and predictive potential of cardiomyocytes during automated drug experiments

All new drugs must go through preclinical screening tests to determine their proarrhythmic potential. While these assays effectively filter out dangerous drugs, they are too conservative, often misclassifying safe compounds as proarrhythmic. In this study, we attempt to address this shortcoming with a novel, medium-throughput drug-screening approach: we use an automated patch-clamp system to acquire optimized voltage clamp (VC) and action potential (AP) data from human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) at several drug concentrations (baseline, 3×, 10× and 20× the effective free plasma concentrations). With our novel method, we show correlations between INa block and upstroke slowing after treatment with flecainide or quinine. Additionally, after quinine treatment, we identify significant reductions in current during voltage steps designed to isolate If and IKs. However, we do not detect any IKr block by either drug, and upon further investigation, do not see any IKr present in the iPSC-CMs when prepared for automated patch experiments (i.e. in suspension) - this is in contrast to similar experiments we have conducted with these cells using the manual patch setup. In this study, we: (1) present a proof-of-concept demonstration of a single-cell medium-throughput drug study, and (2) characterize the non-canonical electrophysiology of iPSC-CMs when prepared for experiments in a medium-throughput setting. KEY POINTS: Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) offer potential as an in vitro model to study the proarrhythmic potential of drugs, but insights from these cells are often limited by the low throughput of manual patch-clamp. In this study, we use a medium-throughput automated patch-clamp system to acquire action potential (AP) and complex voltage clamp (VC) data from single iPSC-CMs at multiple drug concentrations. A correlation between AP upstroke and INa transients was identified and drug-induced changes in ionic currents found. We also characterize the substantially altered physiology of iPSC-CMs when patched in an automated system, suggesting the need to investigate differences between manual and automated patch experiments.

Single-cell ionic current phenotyping explains stem cell-derived cardiomyocyte action potential morphology

Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) are a promising tool to study arrhythmia-related factors, but the variability of action potential (AP) recordings from these cells limits their use as an in vitro model. In this study, we use recently published brief (10 s), dynamic voltage-clamp (VC) data to provide mechanistic insights into the ionic currents contributing to AP heterogeneity; we call this approach rapid ionic current phenotyping (RICP). Features of this VC data were correlated to AP recordings from the same cells, and we used computational models to generate mechanistic insights into cellular heterogeneity. This analysis uncovered several interesting links between AP morphology and ionic current density: both L-type calcium and sodium currents contribute to upstroke velocity, rapid delayed rectifier K+ current is the main determinant of the maximal diastolic potential, and an outward current in the activation range of slow delayed rectifier K+ is the main determinant of AP duration. Our analysis also identified an outward current in several cells at 6 mV that is not reproduced by iPSC-CM mathematical models but contributes to determining AP duration. RICP can be used to explain how cell-to-cell variability in ionic currents gives rise to AP heterogeneity. Because of its brief duration (10 s) and ease of data interpretation, we recommend the use of RICP for single-cell patch-clamp experiments that include the acquisition of APs.NEW & NOTEWORTHY We present rapid ionic current phenotyping (RICP), a current quantification approach based on an optimized voltage-clamp protocol. The method captures a rich snapshot of the ionic current dynamics, providing quantitative information about multiple currents (e.g., ICa,L, IKr) in the same cell. The protocol helped to identify key ionic determinants of cellular action potential heterogeneity in iPSC-CMs. This included unexpected results, such as the critical role of IKr in establishing the maximum diastolic potential.

Strengthening cardiac therapy pipelines using human pluripotent stem cell-derived cardiomyocytes

Advances in hiPSC isolation and reprogramming and hPSC-CM differentiation have prompted their therapeutic application and utilization for evaluating potential cardiovascular safety liabilities. In this perspective, we showcase key efforts toward the large-scale production of hiPSC-CMs, implementation of hiPSC-CMs in industry settings, and recent clinical applications of this technology. The key observations are a need for traceable gender and ethnically diverse hiPSC lines, approaches to reduce cost of scale-up, accessible clinical trial datasets, and transparent guidelines surrounding the safety and efficacy of hiPSC-based therapies.

Optimization of a cardiomyocyte model illuminates role of increased INa,L in repolarization reserve

Cardiac ion currents may compensate for each other when one is compromised by a congenital or drug-induced defect. Such redundancy contributes to a robust repolarization reserve that can prevent the development of lethal arrhythmias. Most efforts made to describe this phenomenon have quantified contributions by individual ion currents. However, it is important to understand the interplay between all major ion-channel conductances, as repolarization reserve is dependent on the balance between all ion currents in a cardiomyocyte. Here, a genetic algorithm was designed to derive profiles of nine ion-channel conductances that optimize repolarization reserve in a mathematical cardiomyocyte model. Repolarization reserve was quantified using a previously defined metric, repolarization reserve current, i.e., the minimum constant current to prevent normal action potential repolarization in a cell. The optimization improved repolarization reserve current up to 84% compared to baseline in a human adult ventricular myocyte model and increased resistance to arrhythmogenic insult. The optimized conductance profiles were not only characterized by increased repolarizing current conductances but also uncovered a previously unreported behavior by the late sodium current. Simulations demonstrated that upregulated late sodium increased action potential duration, without compromising repolarization reserve current. The finding was generalized to multiple models. Ultimately, this computational approach, in which multiple currents were studied simultaneously, illuminated mechanistic insights into how the metric's magnitude could be increased and allowed for the unexpected role of late sodium to be elucidated.NEW & NOTEWORTHY Genetic algorithms are typically used to fit models or extract desired parameters from data. Here, we use the tool to produce a ventricular cardiomyocyte model with increased repolarization reserve. Since arrhythmia mitigation is dependent on multiple cardiac ion-channel conductances, study using a comprehensive, unbiased, and systems-level approach is important. The use of this optimization strategy allowed us to find robust profiles that illuminated unexpected mechanistic determinants of key ion-channel conductances in repolarization reserve.

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.

The Dynamic Clamp Technique: A Robust Toolkit for Investigating Potassium Channel Function

The dynamic clamp technique has emerged as a powerful tool in the field of cardiac electrophysiology, enabling researchers to investigate the intricate dynamics of ion currents in cardiac cells. Potassium channels play a critical role in the functioning of cardiac cells and the overall electrical stability of the heart. This chapter provides a comprehensive overview of the methods and applications of dynamic clamp in the study of key potassium currents in cardiac cells. A step-by-step guide is presented, detailing the experimental setup and protocols required for implementing the dynamic clamp technique in cardiac cell studies. Special attention is given to the design and construction of a dynamic clamp setup with Real Time eXperimental Interface, configurations, and the incorporation of mathematical models to mimic ion channel behavior. The chapter's core focuses on applying dynamic clamp to elucidate the properties of various potassium channels in cardiac cells. It discusses how dynamic clamp can be used to investigate channel kinetics, voltage-dependent properties, and the impact of different potassium channel subtypes on cardiac electrophysiology. The chapter will also include examples of specific dynamic clamp experiments that studied potassium currents or their applications in cardiac cells.

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.

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.

Diversity of cells and signals in the cardiovascular system

This white paper is the outcome of the seventh UC Davis Cardiovascular Research Symposium on Systems Approach to Understanding Cardiovascular Disease and Arrhythmia. This biannual meeting aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2022 Symposium was 'Cell Diversity in the Cardiovascular System, cell-autonomous and cell-cell signalling'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies, and challenges in examining cell and signal diversity, co-ordination and interrelationships involved in cardiovascular function. This paper originates from the topics of formal presentations and informal discussions from the Symposium, which aimed to develop a holistic view of how the multiple cell types in the cardiovascular system integrate to influence cardiovascular function, disease progression and therapeutic strategies. The first section describes the major cell types (e.g. cardiomyocytes, vascular smooth muscle and endothelial cells, fibroblasts, neurons, immune cells, etc.) and the signals involved in cardiovascular function. The second section emphasizes the complexity at the subcellular, cellular and system levels in the context of cardiovascular development, ageing and disease. Finally, the third section surveys the technological innovations that allow the interrogation of this diversity and advancing our understanding of the integrated cardiovascular function and dysfunction.

Emerging methods to model cardiac ion channel and myocyte electrophysiology

In the field of cardiac electrophysiology, modeling has played a central role for many decades. However, even though the effort is well-established, it has recently seen a rapid and sustained evolution in the complexity and predictive power of the models being created. In particular, new approaches to modeling have allowed the tracking of parallel and interconnected processes that span from the nanometers and femtoseconds that determine ion channel gating to the centimeters and minutes needed to describe an arrhythmia. The connection between scales has brought unprecedented insight into cardiac arrhythmia mechanisms and drug therapies. This review focuses on the generation of these models from first principles, generation of detailed models to describe ion channel kinetics, algorithms to create and numerically solve kinetic models, and new approaches toward data gathering that parameterize these models. While we focus on application of these models for cardiac arrhythmia, these concepts are widely applicable to model the physiology and pathophysiology of any excitable cell.

Human-Induced Pluripotent Stem Cell-Based Differentiation of Cardiomyocyte Subtypes for Drug Discovery and Cell Therapy

Drug attrition rates have increased over the past few years, accompanied with growing costs for the pharmaceutical industry and consumers. Lack of in vitro models connecting the results of toxicity screening assays with clinical outcomes accounts for this high attrition rate. The emergence of cardiomyocytes derived from human pluripotent stem cells provides an amenable source of cells for disease modeling, drug discovery, and cardiotoxicity screening. Functionally similar to to embryonic stem cells, but with fewer ethical concerns, induced pluripotent stem cells (iPSCs) can recapitulate patient-specific genetic backgrounds, which would be a huge revolution for personalized medicine. The generated iPSC-derived cardiomyocytes (iPSC-CMs) represent different subtypes including ventricular-, atrial-, and nodal-like cardiomyocytes. Purifying these subtypes for chamber-specific drug screening presents opportunities and challenges. In this chapter, we discuss the strategies for the purification of iPSC-CMs, the use of iPSC-CMs for drug discovery and cardiotoxicity test, and the current limitations of iPSC-CMs that should be overcome for wider and more precise cardiovascular applications.

Inhibitory Effects of Nobiletin on Voltage-Gated Na+ Channel in Rat Ventricular Myocytes Based on Electrophysiological Analysis and Molecular Docking Method

Nobiletin (NOB) has attracted much attention owing to its outstanding bioactivities. This study aimed to investigate its anti-arrhythmic effect through electrophysiological and molecular docking studies. We assessed the anti-arrhythmic effects of NOB using aconitine-induced ventricular arrhythmia in a rat model and the electrophysiological effects of NOB on rat cardiomyocytes utilizing whole-cell patch-clamp techniques. Moreover, we investigated the binding characters of NOB with rNav1.5, rNav1.5/QQQ, and hNa_V1.5 via docking analysis, comparing them with amiodarone and aconitine. NOB pretreatment delayed susceptibility to ventricular premature and ventricular tachycardia and decreased the incidence of fatal ventricular fibrillation. Whole-cell patch-clamp assays demonstrated that the peak current density of the voltage-gated Na⁺ channel current was reversibly reduced by NOB in a concentration-dependent manner. The steady-state activation and recovery curves were shifted in the positive direction along the voltage axis, and the steady-state inactivation curve was shifted in the negative direction along the voltage axis, as shown by gating kinetics. The molecular docking study showed NOB formed a π-π stacking interaction with rNav1.5 and rNav1.5/QQQ upon Phe-1762, which is the homolog to Phe-1760 in hNa_V1.5 and plays an important role in antiarrhythmic action This study reveals that NOB may act as a class I sodium channel anti-arrhythmia agent.