Neural Network Differential Equations For Ion Channel Modelling

Mathematical models of cardiac ion channels have been widely used to study and predict the behaviour of ion currents. Typically models are built using biophysically-based mechanistic principles such as Hodgkin-Huxley or Markov state transitions. These models provide an abstract description of the underlying conformational changes of the ion channels. However, due to the abstracted conformation states and assumptions for the rates of transition between them, there are differences between the models and reality—termed model discrepancy or misspecification. In this paper, we demonstrate the feasibility of using a mechanistically-inspired neural network differential equation model, a hybrid non-parametric model, to model ion channel kinetics. We apply it to the hERG potassium ion channel as an example, with the aim of providing an alternative modelling approach that could alleviate certain limitations of the traditional approach. We compare and discuss multiple ways of using a neural network to approximate extra hidden states or alternative transition rates. In particular we assess their ability to learn the missing dynamics, and ask whether we can use these models to handle model discrepancy. Finally, we discuss the practicality and limitations of using neural networks and their potential applications

Importance of modelling hERG binding in predicting drug-induced action potential prolongations for drug safety assessment

Reduction of the rapid delayed rectifier potassium current (IKr) via drug binding to the human Ether-à-go-go-Related Gene (hERG) channel is a well recognised mechanism that can contribute to an increased risk of Torsades de Pointes. Mathematical models have been created to replicate the effects of channel blockers, such as reducing the ionic conductance of the channel. Here, we study the impact of including state-dependent drug binding in a mathematical model of hERG when translating hERG inhibition to action potential changes. We show that the difference in action potential predictions when modelling drug binding of hERG using a state-dependent model versus a conductance scaling model depends not only on the properties of the drug and whether the experiment achieves steady state, but also on the experimental protocols. Furthermore, through exploring the model parameter space, we demonstrate that the state-dependent model and the conductance scaling model generally predict different action potential prolongations and are not interchangeable, while at high binding and unbinding rates, the conductance scaling model tends to predict shorter action potential prolongations. Finally, we observe that the difference in simulated action potentials between the models is determined by the binding and unbinding rate, rather than the trapping mechanism. This study demonstrates the importance of modelling drug binding and highlights the need for improved understanding of drug trapping which can have implications for the uses in drug safety assessment

A Bayesian nonparametric method for detecting rapid changes in disease transmission

Whether an outbreak of infectious disease is likely to grow or dissipate is determined through the time-varying reproduction number, Rt. Real-time or retrospective identification of changes in Rt following the imposition or relaxation of interventions can thus contribute important evidence about disease transmission dynamics which can inform policymaking. Here, we present a method for estimating shifts in Rt within a renewal model framework. Our method, which we call EpiCluster, is a Bayesian nonparametric model based on the Pitman-Yor process. We assume that Rt is piecewise-constant, and the incidence data and priors determine when or whether Rt should change and how many times it should do so throughout the series. We also introduce a prior which induces sparsity over the number of changepoints. Being Bayesian, our approach yields a measure of uncertainty in Rt and its changepoints. EpiCluster is fast, straightforward to use, and we demonstrate that it provides automated detection of rapid changes in transmission, either in real-time or retrospectively, for synthetic data series where the Rt profile is known. We illustrate the practical utility of our method by fitting it to case data of outbreaks of COVID-19 in Australia and Hong Kong, where it finds changepoints coinciding with the imposition of non-pharmaceutical interventions. Bayesian nonparametric methods, such as ours, allow the volume and complexity of the data to dictate the number of parameters required to approximate the process and should find wide application in epidemiology

Rapid Characterization of hERG Channel Kinetics I: Using an Automated High-Throughput System

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 pro-arrhythmic 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 characterise a current is typically too long to be applied in a single cell. Shorter high-information protocols have recently been introduced which have this capability, but they are not typically compatible with high-throughput platforms. We present a new 15 second protocol to characterise 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 CHO cells stably expressing hERG1a. From these recordings we construct 124 cell-specific variants/parameterisations of a hERG model at 25C. A further 8 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, which we use 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

Model-driven optimal experimental design for calibrating cardiac electrophysiology models

Background and Objective: Models of the cardiomyocyte action potential have contributed immensely to the understanding of heart function, pathophysiology, and the origin of heart rhythm disturbances. However, action potential models are highly nonlinear, making them difficult to parameterise and limiting to describing ‘average cell’ dynamics, when cell-specific models would be ideal to uncover inter-cell variability but are too experimentally challenging to be achieved. Here, we focus on automatically designing experimental protocols that allow us to better identify cell-specific maximum conductance values for each major current type.Methods and Results: We developed an approach that applies optimal experimental designs to patch-clamp experiments, including both voltage-clamp and current-clamp experiments. We assessed the models calibrated to these new optimal designs by comparing them to the models calibrated to some of the commonly used designs in the literature. We showed that optimal designs are not only overall shorter in duration but also able to perform better than many of the existing experiment designs in terms of identifying model parameters and hence model predictive power.Conclusions: For cardiac cellular electrophysiology, this approach will allow researchers to define their hypothesis of the dynamics of the system and automatically design experimental protocols that will result in theoretically optimal designs

Ten simple rules for teaching sustainable software engineering

Computational methods and associated software implementations are central to every field of scientific investigation. Modern biological research, particularly within systems biology, has relied heavily on the development of software tools to process and organize increasingly large datasets, simulate complex mechanistic models, provide tools for the analysis and management of data, and visualize and organize outputs. However, developing high-quality research software requires scientists to develop a host of software development skills, and teaching these skills to students is challenging. There has been a growing importance placed on ensuring reproducibility and good development practices in computational research. However, less attention has been devoted to informing the specific teaching strategies which are effective at nurturing in researchers the complex skillset required to produce high-quality software that, increasingly, is required to underpin both academic and industrial biomedical research. Recent articles in the Ten Simple Rules collection have discussed the teaching of foundational computer science and coding techniques to biology students. We advance this discussion by describing the specific steps for effectively teaching the necessary skills scientists need to develop sustainable software packages which are fit for (re-)use in academic research or more widely. Although our advice is likely to be applicable to all students and researchers hoping to improve their software development skills, our guidelines are directed towards an audience of students that have some programming literacy but little formal training in software development or engineering, typical of early doctoral students. These practices are also applicable outside of doctoral training environments, and we believe they should form a key part of postgraduate training schemes more generally in the life sciences.Comment: Prepared for submission to PLOS Computational Biology's 10 Simple Rules collectio

Autocorrelated measurement processes and inference for ordinary differential equation models of biological systems

Ordinary differential equation models are used to describe dynamic processes across biology. To perform likelihood-based parameter inference on these models, it is necessary to specify a statistical process representing the contribution of factors not explicitly included in the mathematical model. For this, independent Gaussian noise is commonly chosen, with its use so widespread that researchers typically provide no explicit justification for this choice. This noise model assumes `random' latent factors affect the system in ephemeral fashion resulting in unsystematic deviation of observables from their modelled counterparts. However, like the deterministically modelled parts of a system, these latent factors can have persistent effects on observables. Here, we use experimental data from dynamical systems drawn from cardiac physiology and electrochemistry to demonstrate that highly persistent differences between observations and modelled quantities can occur. Considering the case when persistent noise arises due only to measurement imperfections, we use the Fisher information matrix to quantify how uncertainty in parameter estimates is artificially reduced when erroneously assuming independent noise. We present a workflow to diagnose persistent noise from model fits and describe how to remodel accounting for correlated errors

Accounting for variability in ion current recordings using a mathematical model of artefacts in voltage-clamp experiments

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'