Predicting how pharmaceuticals may affect heart rhythm is a crucial step in drug development and requires a deep understanding of a compound’s action on ion channels. In vitro hERG channel current recordings are an important step in evaluating the proarrhythmic potential of small molecules and are now routinely performed using automated high-throughput patch-clamp platforms. These machines can execute traditional voltage-clamp protocols aimed at specific gating processes, but the array of protocols needed to fully characterize a current is typically too long to be applied in a single cell. Shorter high-information protocols have recently been introduced that have this capability, but they are not typically compatible with high-throughput platforms. We present a new 15 second protocol to characterize hERG (Kv11.1) kinetics, suitable for both manual and high-throughput systems. We demonstrate its use on the Nanion SyncroPatch 384PE, a 384-well automated patch-clamp platform, by applying it to Chinese hamster ovary cells stably expressing hERG1a. From these recordings, we construct 124 cell-specific variants/parameterizations of a hERG model at 25°C. A further eight independent protocols are run in each cell and are used to validate the model predictions. We then combine the experimental recordings using a hierarchical Bayesian model, which we use to quantify the uncertainty in the model parameters, and their variability from cell-to-cell; we use this model to suggest reasons for the variability. This study demonstrates a robust method to measure and quantify uncertainty and shows that it is possible and practical to use high-throughput systems to capture full hERG channel kinetics quantitatively and rapidly.

Cardiac electrophysiology models are among the most mature and well‐studied mathematical models of biological systems. This maturity is bringing new challenges as models are being used increasingly to make quantitative rather than qualitative predictions. As such, calibrating the parameters within ion current and action potential (AP) models to experimental data sets is a crucial step in constructing a predictive model. This review highlights some of the fundamental concepts in cardiac model calibration and is intended to be readily understood by computational and mathematical modelers working in other fields of biology. We discuss the classic and latest approaches to calibration in the electrophysiology field, at both the ion channel and cellular AP scales. We end with a discussion of the many challenges that work to date has raised and the need for reproducible descriptions of the calibration process to enable models to be recalibrated to new data sets and built upon for new studies. This article is categorized under: Analytical and Computational Methods > Computational Methods Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Cellular Models

Time series models are ubiquitous in science, arising in any situation where researchers seek to understand how a system's behaviour changes over time. A key problem in time series modelling is inference; determining properties of the underlying system based on observed time series. For both statistical and mechanistic models, inference involves finding parameter values, or distributions of parameters values, which produce outputs consistent with observations. A wide variety of inference techniques are available and different approaches are suitable for different classes of problems. This variety presents a challenge for researchers, who may not have the resources or expertise to implement and experiment with these methods. PINTS (Probabilistic Inference on Noisy Time Series-https://github.com/pints-team/pints) is an open-source (BSD 3-clause license) Python library that provides researchers with a broad suite of non-linear optimisation and sampling methods. It allows users to wrap a model and data in a transparent and straightforward interface, which can then be used with custom or pre-defined error measures for optimisation, or with likelihood functions for Bayesian inference or maximum-likelihood estimation. Derivative-free optimisation algorithms-which work without harder-to-obtain gradient information-are included, as well as inference algorithms such as adaptive Markov chain Monte Carlo and nested sampling, which estimate distributions over parameter values. By making these statistical techniques available in an open and easy-to-use framework, PINTS brings the power of these modern methods to a wider scientific audience.

Ion channel behavior can depend strongly on temperature, with faster kinetics at physiological temperatures leading to considerable changes in currents relative to room temperature. These temperature-dependent changes in voltage-dependent ion channel kinetics (rates of opening, closing, inactivating, and recovery) are commonly represented with Q 10 coefficients or an Eyring relationship. In this article, we assess the validity of these representations by characterizing channel kinetics at multiple temperatures. We focus on the human Ether-à-go-go-Related Gene (hERG) channel, which is important in drug safety assessment and commonly screened at room temperature so that results require extrapolation to physiological temperature. In Part I of this study, we established a reliable method for high-throughput characterization of hERG1a (Kv11.1) kinetics, using a 15-second information-rich optimized protocol. In this Part II, we use this protocol to study the temperature dependence of hERG kinetics using Chinese hamster ovary cells overexpressing hERG1a on the Nanion SyncroPatch 384PE, a 384-well automated patch-clamp platform, with temperature control. We characterize the temperature dependence of hERG gating by fitting the parameters of a mathematical model of hERG kinetics to data obtained at five distinct temperatures between 25 and 37°C and validate the models using different protocols. Our models reveal that activation is far more temperature sensitive than inactivation, and we observe that the temperature dependency of the kinetic parameters is not represented well by Q 10 coefficients; it broadly follows a generalized, but not the standardly-used, Eyring relationship. We also demonstrate that experimental estimations of Q 10 coefficients are protocol dependent. Our results show that a direct fit using our 15-s protocol best represents hERG kinetics at any given temperature and suggests that using the Generalized Eyring theory is preferable if no experimental data are available to derive model parameters at a given temperature.

Mathematical models of ion channels, which constitute indispensable components of action potential models, are commonly constructed by fitting to whole-cell patch-clamp data. In a previous study, we fitted cell-specific models to hERG1a (Kv11.1) recordings simultaneously measured using an automated high-throughput system, and studied cell-cell variability by inspecting the resulting model parameters. However, the origin of the observed variability was not identified. Here, we study the source of variability by constructing a model that describes not just ion current dynamics, but the entire voltage-clamp experiment. The experimental artefact components of the model include: series resistance, membrane and pipette capacitance, voltage offsets, imperfect compensations made by the amplifier for these phenomena, and leak current. In this model, variability in the observations can be explained by either cell properties, measurement artefacts, or both. Remarkably, by assuming that variability arises exclusively from measurement artefacts, it is possible to explain a larger amount of the observed variability than when assuming cell-specific ion current kinetics. This assumption also leads to a smaller number of model parameters. This result suggests that most of the observed variability in patch-clamp data measured under the same conditions is caused by experimental artefacts, and hence can be compensated for in post-processing by using our model for the patch-clamp experiment. This study has implications for the question of the extent to which cell-cell variability in ion channel kinetics exists, and opens up routes for better correction of artefacts in patch-clamp data. This article is part of the theme issue ‘Uncertainty quantification in cardiac and cardiovascular modelling and simulation’.

1 Affiliations removed for initial submission as per guidelines. 6 ABSTRACT Predicting how pharmaceuticals may affect heart rhythm is a crucial step in drug-development, and requires a 7 deep understanding of a compound's action on ion channels. In vitro hERG-channel current recordings are an important step in 8 evaluating the pro-arrhythmic potential of small molecules, and are now routinely performed using automated high-throughput 9 patch clamp platforms. These machines can execute traditional voltage clamp protocols aimed at specific gating processes, 10 but the array of protocols needed to fully characterise a current is typically too long to be applied in a single cell. Shorter 11 high-information protocols have recently been introduced which have this capability, but they are not typically compatible with 12 high-throughput platforms. We present a new high-information 15 s protocol to characterise hERG (Kv11.1) kinetics, suitable for 13 both manual and high-throughput systems. We demonstrate its use on the Nanion SyncroPatch 384PE, a 384 well automated 14 patch clamp platform, by applying it to CHO cells stably expressing hERG1a. From these recordings we construct 124 cell-specific 15 variants/parameterisations of a hERG model at 25 • C. A further 8 independent protocols are run in each cell, and are used to 16 validate the model predictions. We then combine the experimental recordings using a hierarchical Bayesian model, which we 17 use to quantify the uncertainty in the model parameters, and their variability from cell to cell, which we use to suggest reasons 18 for the variability. This study demonstrates a robust method to measure and quantify uncertainty, and shows that it is possible 19 and practical to use high-throughput systems to capture full hERG channel kinetics quantitatively and rapidly. 20 21 We present a method for high-throughput characterisation of hERG potassium channel kinetics, via fitting a mathematical 22 model to results of over one hundred single cell patch clamp measurements collected simultaneously on an automated voltage 23 clamp platform. The automated patch clamp data are used to parameterise a mathematical ion channel model fully, opening a 24 new era of automated and rapid development of mathematical models from quick and cheap experiments. The method also 25 allows ample data for independent validation of the models and enables us to study experimental variability and propose its 26 origins. In future the method can be applied to characterise changes to hERG currents in different conditions, for instance 27 at different temperatures (see Part II of the study) or under mutations or the action of pharmaceuticals; and should be easily 28 adapted to study many other currents. Statement of Significance 30The human Ether-à-go-go-Related Gene (hERG) is of great interest in the field of cardiac safety pharmacology. hERG encodes 31 the pore-forming alpha subunit of the ion channel Kv11.1 which conducts the rapid delayed rectifier potassium current, I Kr 32(1)...

Uncertainty quantification (UQ) is a vital step in using mathematical models and simulations to take decisions. The field of cardiac simulation has begun to explore and adopt UQ methods to characterize uncertainty in model inputs and how that propagates through to outputs or predictions; examples of this can be seen in the papers of this issue. In this review and perspective piece, we draw attention to an important and under-addressed source of uncertainty in our predictions—that of uncertainty in the model structure or the equations themselves. The difference between imperfect models and reality is termed model discrepancy , and we are often uncertain as to the size and consequences of this discrepancy. Here, we provide two examples of the consequences of discrepancy when calibrating models at the ion channel and action potential scales. Furthermore, we attempt to account for this discrepancy when calibrating and validating an ion channel model using different methods, based on modelling the discrepancy using Gaussian processes and autoregressive-moving-average models, then highlight the advantages and shortcomings of each approach. Finally, suggestions and lines of enquiry for future work are provided. This article is part of the theme issue ‘Uncertainty quantification in cardiac and cardiovascular modelling and simulation’.

Human induced pluripotent stem cell derived cardiomyocytes (iPSC-CMs) have applications in disease modeling, cell therapy, drug screening and personalized medicine. Computational models can be used to interpret experimental findings in iPSC-CMs, provide mechanistic insights, and translate these findings to adult cardiomyocyte (CM) electrophysiology. However, different cell lines display different expression of ion channels, pumps and receptors, and show differences in electrophysiology. In this exploratory study, we use a mathematical model based on iPSC-CMs from Cellular Dynamic International (CDI, iCell), and compare its predictions to novel experimental recordings made with the Axiogenesis Cor.4U line. We show that tailoring this model to the specific cell line, even using limited data and a relatively simple approach, leads to improved predictions of baseline behavior and response to drugs. This demonstrates the need and the feasibility to tailor models to individual cell lines, although a more refined approach will be needed to characterize individual currents, address differences in ion current kinetics, and further improve these results.