Cell-to-cell heterogeneity in ion channel conductance impacts substrate vulnerability to arrhythmia.

Determining whether an ectopic depolarization will lead to a self-perpetuating arrhythmia is of critical importance in determining arrhythmia risk, so it is necessary to understand what factors impact substrate vulnerability. This study sought to explore the impact of cell-to-cell heterogeneity in ion channel conductance on substrate vulnerability to arrhythmia by measuring the duration of the vulnerable window in computational models of one-dimensional cables of ventricular cardiomyocytes. We began by using a population of uniform cable models to determine the mechanisms underlying the vulnerable window phenomenon. We found that in addition to the known importance of GNa, the conductances GCa,L and GKr also play a minor role in determining the vulnerable window duration. We also found that a steeper slope of the repolarizing action potential during the vulnerable window correlated with a shorter vulnerable window duration in uniform cables. We applied our understanding from these initial simulations to an investigation of the vulnerable window in heterogeneous cable models. The heterogeneous cables displayed a great deal of intra-cable variation in vulnerable window duration, highly sensitive to the cardiomyocytes in the local environment of the ectopic stimulus. Coupling strength modulated not only the magnitude of the vulnerable window duration but also the extent of intra-tissue variability in vulnerable window duration.NEW & NOTEWORTHY We investigate the impact of cell-to-cell heterogeneity in ion channel conductance on substrate vulnerability to arrhythmia by measuring the vulnerable window duration in computational cardiomyocyte cable models. We demonstrate a wide range of intra-cable variability in vulnerable window duration (VWD) and show how this is changed by ion channel block and coupling strength perturbations.

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.

Investigational Treatments for COVID-19 may Increase Ventricular Arrhythmia Risk Through Drug Interactions.

Many drugs that have been proposed for treatment of coronavirus disease 2019 (COVID-19) are reported to cause cardiac adverse events, including ventricular arrhythmias. In order to properly weigh risks against potential benefits, particularly when decisions must be made quickly, mathematical modeling of both drug disposition and drug action can be useful for predicting patient response and making informed decisions. Here, we explored the potential effects on cardiac electrophysiology of four drugs proposed to treat COVID-19: lopinavir, ritonavir, chloroquine, and azithromycin, as well as combination therapy involving these drugs. Our study combined simulations of pharmacokinetics (PKs) with quantitative systems pharmacology (QSP) modeling of ventricular myocytes to predict potential cardiac adverse events caused by these treatments. Simulation results predicted that drug combinations can lead to greater cellular action potential prolongation, analogous to QT prolongation, compared with drugs given in isolation. The combination effect can result from both PK and pharmacodynamic drug interactions. Importantly, simulations of different patient groups predicted that women with pre-existing heart disease are especially susceptible to drug-induced arrhythmias, compared with diseased men or healthy individuals of either sex. Statistical analysis of population simulations revealed the molecular factors that make certain women with heart failure especially susceptible to arrhythmias. Overall, the results illustrate how PK and QSP modeling may be combined to more precisely predict cardiac risks of COVID-19 therapies.

Investigational treatments for COVID-19 may increase ventricular arrhythmia risk through drug interactions.

Many drugs that have been proposed for treatment of COVID-19 are reported to cause cardiac adverse events, including ventricular arrhythmias. In order to properly weigh risks against potential benefits, particularly when decisions must be made quickly, mathematical modeling of both drug disposition and drug action can be useful for predicting patient response and making informed decisions. Here we explored the potential effects on cardiac electrophysiology of 4 drugs proposed to treat COVID-19: lopinavir, ritonavir, chloroquine, and azithromycin, as well as combination therapy involving these drugs. Our study combined simulations of pharmacokinetics (PK) with quantitative systems pharmacology (QSP) modeling of ventricular myocytes to predict potential cardiac adverse events caused by these treatments. Simulation results predicted that drug combinations can lead to greater cellular action potential prolongation, analogous to QT prolongation, compared with drugs given in isolation. The combination effect can result from both pharmacokinetic and pharmacodynamic drug interactions. Importantly, simulations of different patient groups predicted that females with pre-existing heart disease are especially susceptible to drug-induced arrhythmias, compared males with disease or healthy individuals of either sex. Overall, the results illustrate how PK and QSP modeling may be combined to more precisely predict cardiac risks of COVID-19 therapies.