Corrigendum: The Electrogenic Na+/K+ Pump Is a Key Determinant of Repolarization Abnormality Susceptibility in Human Ventricular Cardiomyocytes: A Population-Based Simulation Study

[This corrects the article on p. 278 in vol. 8, PMID: 28529489.].

Human In Silico Drug Trials Demonstrate Higher Accuracy than Animal Models in Predicting Clinical Pro-Arrhythmic Cardiotoxicity

Early prediction of cardiotoxicity is critical for drug development. Current animal models raise ethical and translational questions, and have limited accuracy in clinical risk prediction. Human-based computer models constitute a fast, cheap and potentially effective alternative to experimental assays, also facilitating translation to human. Key challenges include consideration of inter-cellular variability in drug responses and integration of computational and experimental methods in safety pharmacology. Our aim is to evaluate the ability of in silico drug trials in populations of human action potential (AP) models to predict clinical risk of drug-induced arrhythmias based on ion channel information, and to compare simulation results against experimental assays commonly used for drug testing. A control population of 1,213 human ventricular AP models in agreement with experimental recordings was constructed. In silico drug trials were performed for 62 reference compounds at multiple concentrations, using pore-block drug models (IC₅₀/Hill coefficient). Drug-induced changes in AP biomarkers were quantified, together with occurrence of repolarization/depolarization abnormalities. Simulation results were used to predict clinical risk based on reports of Torsade de Pointes arrhythmias, and further evaluated in a subset of compounds through comparison with electrocardiograms from rabbit wedge preparations and Ca²⁺-transient recordings in human induced pluripotent stem cell-derived cardiomyocytes (hiPS-CMs). Drug-induced changes in silico vary in magnitude depending on the specific ionic profile of each model in the population, thus allowing to identify cell sub-populations at higher risk of developing abnormal AP phenotypes. Models with low repolarization reserve (increased Ca²⁺/late Na⁺ currents and Na⁺/Ca²⁺-exchanger, reduced Na⁺/K⁺-pump) are highly vulnerable to drug-induced repolarization abnormalities, while those with reduced inward current density (fast/late Na⁺ and Ca²⁺ currents) exhibit high susceptibility to depolarization abnormalities. Repolarization abnormalities in silico predict clinical risk for all compounds with 89% accuracy. Drug-induced changes in biomarkers are in overall agreement across different assays: in silico AP duration changes reflect the ones observed in rabbit QT interval and hiPS-CMs Ca²⁺-transient, and simulated upstroke velocity captures variations in rabbit QRS complex. Our results demonstrate that human in silico drug trials constitute a powerful methodology for prediction of clinical pro-arrhythmic cardiotoxicity, ready for integration in the existing drug safety assessment pipelines.

Quantitative Comparison of Effects of Dofetilide, Sotalol, Quinidine, and Verapamil between Human Ex vivo Trabeculae and In silico Ventricular Models Incorporating Inter-Individual Action Potential Variability

Background:In silico modeling could soon become a mainstream method of pro-arrhythmic risk assessment in drug development. However, a lack of human-specific data and appropriate modeling techniques has previously prevented quantitative comparison of drug effects between in silico models and recordings from human cardiac preparations. Here, we directly compare changes in repolarization biomarkers caused by dofetilide, dl-sotalol, quinidine, and verapamil, between in silico populations of human ventricular cell models and ex vivo human ventricular trabeculae. Methods and Results:Ex vivo recordings from human ventricular trabeculae in control conditions were used to develop populations of in silico human ventricular cell models that integrated intra- and inter-individual variability in action potential (AP) biomarker values. Models were based on the O'Hara-Rudy ventricular cardiomyocyte model, but integrated experimental AP variability through variation in underlying ionic conductances. Changes to AP duration, triangulation and early after-depolarization occurrence from application of the four drugs at multiple concentrations and pacing frequencies were compared between simulations and experiments. To assess the impact of variability in IC50 measurements, and the effects of including state-dependent drug binding dynamics, each drug simulation was repeated with two different IC50 datasets, and with both the original O'Hara-Rudy hERG model and a recently published state-dependent model of hERG and hERG block. For the selective hERG blockers dofetilide and sotalol, simulation predictions of AP prolongation and repolarization abnormality occurrence showed overall good agreement with experiments. However, for multichannel blockers quinidine and verapamil, simulations were not in agreement with experiments across all IC50 datasets and IKr block models tested. Quinidine simulations resulted in overprolonged APs and high incidence of repolarization abnormalities, which were not observed in experiments. Verapamil simulations showed substantial AP prolongation while experiments showed mild AP shortening. Conclusions: Results for dofetilide and sotalol show good agreement between experiments and simulations for selective compounds, however lack of agreement from simulations of quinidine and verapamil suggest further work is needed to understand the more complex electrophysiological effects of these multichannel blocking drugs.

The Electrogenic Na+/K+ Pump Is a Key Determinant of Repolarization Abnormality Susceptibility in Human Ventricular Cardiomyocytes: A Population-Based Simulation Study

Background: Cellular repolarization abnormalities occur unpredictably due to disease and drug effects, and can occur even in cardiomyocytes that exhibit normal action potentials (AP) under control conditions. Variability in ion channel densities may explain differences in this susceptibility to repolarization abnormalities. Here, we quantify the importance of key ionic mechanisms determining repolarization abnormalities following ionic block in human cardiomyocytes yielding normal APs under control conditions.

Methods and Results: Sixty two AP recordings from non-diseased human heart preparations were used to construct a population of human ventricular models with normal APs and a wide range of ion channel densities. Multichannel ionic block was applied to investigate susceptibility to repolarization abnormalities. I_Kr block was necessary for the development of repolarization abnormalities. Models that developed repolarization abnormalities over the widest range of blocks possessed low Na⁺/K⁺ pump conductance below 50% of baseline, and I_CaL conductance above 70% of baseline. Furthermore, I_NaK made the second largest contribution to repolarizing current in control simulations and the largest contribution under 75% I_Kr block. Reversing intracellular Na⁺ overload caused by reduced I_NaK was not sufficient to prevent abnormalities in models with low Na⁺/K⁺ pump conductance, while returning Na⁺/K⁺ pump conductance to normal substantially reduced abnormality occurrence, indicating I_NaK is an important repolarization current.

Conclusions: I_NaK is an important determinant of repolarization abnormality susceptibility in human ventricular cardiomyocytes, through its contribution to repolarization current rather than homeostasis. While we found I_Kr block to be necessary for repolarization abnormalities to occur, I_NaK decrease, as in disease, may amplify the pro-arrhythmic risk of drug-induced I_Kr block in humans.

Variability in cardiac electrophysiology: Using experimentally-calibrated populations of models to move beyond the single virtual physiological human paradigm

Physiological variability manifests itself via differences in physiological function between individuals of the same species, and has crucial implications in disease progression and treatment. Despite its importance, physiological variability has traditionally been ignored in experimental and computational investigations due to averaging over samples from multiple individuals. Recently, modelling frameworks have been devised for studying mechanisms underlying physiological variability in cardiac electrophysiology and pro-arrhythmic risk under a variety of conditions and for several animal species as well as human. One such methodology exploits populations of cardiac cell models constrained with experimental data, or experimentally-calibrated populations of models. In this review, we outline the considerations behind constructing an experimentally-calibrated population of models and review the studies that have employed this approach to investigate variability in cardiac electrophysiology in physiological and pathological conditions, as well as under drug action. We also describe the methodology and compare it with alternative approaches for studying variability in cardiac electrophysiology, including cell-specific modelling approaches, sensitivity-analysis based methods, and populations-of-models frameworks that do not consider the experimental calibration step. We conclude with an outlook for the future, predicting the potential of new methodologies for patient-specific modelling extending beyond the single virtual physiological human paradigm.

Defining the limits of single-molecule FRET resolution in TIRF microscopy

Single-molecule FRET (smFRET) has long been used as a molecular ruler for the study of biology on the nanoscale (∼2-10 nm); smFRET in total-internal reflection fluorescence (TIRF) Förster resonance energy transfer (TIRF-FRET) microscopy allows multiple biomolecules to be simultaneously studied with high temporal and spatial resolution. To operate at the limits of resolution of the technique, it is essential to investigate and rigorously quantify the major sources of noise and error; we used theoretical predictions, simulations, advanced image analysis, and detailed characterization of DNA standards to quantify the limits of TIRF-FRET resolution. We present a theoretical description of the major sources of noise, which was in excellent agreement with results for short-timescale smFRET measurements (<200 ms) on individual molecules (as opposed to measurements on an ensemble of single molecules). For longer timescales (>200 ms) on individual molecules, and for FRET distributions obtained from an ensemble of single molecules, we observed significant broadening beyond theoretical predictions; we investigated the causes of this broadening. For measurements on individual molecules, analysis of the experimental noise allows us to predict a maximum resolution of a FRET change of 0.08 with 20-ms temporal resolution, sufficient to directly resolve distance differences equivalent to one DNA basepair separation (0.34 nm). For measurements on ensembles of single molecules, we demonstrate resolution of distance differences of one basepair with 1000-ms temporal resolution, and differences of two basepairs with 80-ms temporal resolution. Our work paves the way for ultra-high-resolution TIRF-FRET studies on many biomolecules, including DNA processing machinery (DNA and RNA polymerases, helicases, etc.), the mechanisms of which are often characterized by distance changes on the scale of one DNA basepair.