Control of NFAT Isoform Activation and NFAT-Dependent Gene Expression through Two Coincident and Spatially Segregated Intracellular Ca2+ Signals.

Excitation-transcription coupling, linking stimulation at the cell surface to changes in nuclear gene expression, is conserved throughout eukaryotes. How closely related coexpressed transcription factors are differentially activated remains unclear. Here, we show that two Ca2+-dependent transcription factor isoforms, NFAT1 and NFAT4, require distinct sub-cellular InsP3 and Ca2+ signals for physiologically sustained activation. NFAT1 is stimulated by sub-plasmalemmal Ca2+ microdomains, whereas NFAT4 additionally requires Ca2+ mobilization from the inner nuclear envelope by nuclear InsP3 receptors. NFAT1 is rephosphorylated (deactivated) more slowly than NFAT4 in both cytoplasm and nucleus, enabling a more prolonged activation phase. Oscillations in cytoplasmic Ca2+, long considered the physiological form of Ca2+ signaling, play no role in activating either NFAT protein. Instead, effective sustained physiological activation of NFAT4 is tightly linked to oscillations in nuclear Ca2+. Our results show how gene expression can be controlled by coincident yet geographically distinct Ca2+ signals, generated by a freely diffusible InsP3 message.

Uncertainty and error in SARS-CoV-2 epidemiological parameters inferred from population-level epidemic models.

During the SARS-CoV-2 pandemic, epidemic models have been central to policy-making. Public health responses have been shaped by model-based projections and inferences, especially related to the impact of various non-pharmaceutical interventions. Accompanying this has been increased scrutiny over model performance, model assumptions, and the way that uncertainty is incorporated and presented. Here we consider a population-level model, focusing on how distributions representing host infectiousness and the infection-to-death times are modelled, and particularly on the impact of inferred epidemic characteristics if these distributions are mis-specified. We introduce an SIR-type model with the infected population structured by 'infected age', i.e. the number of days since first being infected, a formulation that enables distributions to be incorporated that are consistent with clinical data. We show that inference based on simpler models without infected age, which implicitly mis-specify these distributions, leads to substantial errors in inferred quantities relevant to policy-making, such as the reproduction number and the impact of interventions. We consider uncertainty quantification via a Bayesian approach, implementing this for both synthetic and real data focusing on UK data in the period 15 Feb-14 Jul 2020, and emphasising circumstances where it is misleading to neglect uncertainty. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

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.

Empirical Quantification of Predictive Uncertainty Due to Model Discrepancy by Training with an Ensemble of Experimental Designs: An Application to Ion Channel Kinetics.

When using mathematical models to make quantitative predictions for clinical or industrial use, it is important that predictions come with a reliable estimate of their accuracy (uncertainty quantification). Because models of complex biological systems are always large simplifications, model discrepancy arises-models fail to perfectly recapitulate the true data generating process. This presents a particular challenge for making accurate predictions, and especially for accurately quantifying uncertainty in these predictions. Experimentalists and modellers must choose which experimental procedures (protocols) are used to produce data used to train models. We propose to characterise uncertainty owing to model discrepancy with an ensemble of parameter sets, each of which results from training to data from a different protocol. The variability in predictions from this ensemble provides an empirical estimate of predictive uncertainty owing to model discrepancy, even for unseen protocols. We use the example of electrophysiology experiments that investigate the properties of hERG potassium channels. Here, 'information-rich' protocols allow mathematical models to be trained using numerous short experiments performed on the same cell. In this case, we simulate data with one model and fit it with a different (discrepant) one. For any individual experimental protocol, parameter estimates vary little under repeated samples from the assumed additive independent Gaussian noise model. Yet parameter sets arising from the same model applied to different experiments conflict-highlighting model discrepancy. Our methods will help select more suitable ion channel models for future studies, and will be widely applicable to a range of biological modelling problems.

A Quantitative Systems Pharmacology Perspective on the Importance of Parameter Identifiability.

There is an inherent tension in Quantitative Systems Pharmacology (QSP) between the need to incorporate mathematical descriptions of complex physiology and drug targets with the necessity of developing robust, predictive and well-constrained models. In addition to this, there is no "gold standard" for model development and assessment in QSP. Moreover, there can be confusion over terminology such as model and parameter identifiability; complex and simple models; virtual populations; and other concepts, which leads to potential miscommunication and misapplication of methodologies within modeling communities, both the QSP community and related disciplines. This perspective article highlights the pros and cons of using simple (often identifiable) vs. complex (more physiologically detailed but often non-identifiable) models, as well as aspects of parameter identifiability, sensitivity and inference methodologies for model development and analysis. The paper distills the central themes of the issue of identifiability and optimal model size and discusses open challenges.

A systematic strategy for estimating hERG block potency and its implications in a new cardiac safety paradigm.

INTRODUCTION: hERG block potency is widely used to calculate a drug's safety margin against its torsadogenic potential. Previous studies are confounded by use of different patch clamp electrophysiology protocols and a lack of statistical quantification of experimental variability. Since the new cardiac safety paradigm being discussed by the International Council for Harmonisation promotes a tighter integration of nonclinical and clinical data for torsadogenic risk assessment, a more systematic approach to estimate the hERG block potency and safety margin is needed. METHODS: A cross-industry study was performed to collect hERG data on 28 drugs with known torsadogenic risk using a standardized experimental protocol. A Bayesian hierarchical modeling (BHM) approach was used to assess the hERG block potency of these drugs by quantifying both the inter-site and intra-site variability. A modeling and simulation study was also done to evaluate protocol-dependent changes in hERG potency estimates. RESULTS: A systematic approach to estimate hERG block potency is established. The impact of choosing a safety margin threshold on torsadogenic risk evaluation is explored based on the posterior distributions of hERG potency estimated by this method. The modeling and simulation results suggest any potency estimate is specific to the protocol used. DISCUSSION: This methodology can estimate hERG block potency specific to a given voltage protocol. The relationship between safety margin thresholds and torsadogenic risk predictivity suggests the threshold should be tailored to each specific context of use, and safety margin evaluation may need to be integrated with other information to form a more comprehensive risk assessment.

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'.

Considering discrepancy when calibrating a mechanistic electrophysiology model.

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'.

Four Ways to Fit an Ion Channel Model.

Mathematical models of ionic currents are used to study the electrophysiology of the heart, brain, gut, and several other organs. Increasingly, these models are being used predictively in the clinic, for example, to predict the risks and results of genetic mutations, pharmacological treatments, or surgical procedures. These safety-critical applications depend on accurate characterization of the underlying ionic currents. Four different methods can be found in the literature to fit voltage-sensitive ion channel models to whole-cell current measurements: method 1, fitting model equations directly to time-constant, steady-state, and I-V summary curves; method 2, fitting by comparing simulated versions of these summary curves to their experimental counterparts; method 3, fitting to the current traces themselves from a range of protocols; and method 4, fitting to a single current trace from a short and rapidly fluctuating voltage-clamp protocol. We compare these methods using a set of experiments in which hERG1a current was measured in nine Chinese hamster ovary cells. In each cell, the same sequence of fitting protocols was applied, as well as an independent validation protocol. We show that methods 3 and 4 provide the best predictions on the independent validation set and that short, rapidly fluctuating protocols like that used in method 4 can replace much longer conventional protocols without loss of predictive ability. Although data for method 2 are most readily available from the literature, we find it performs poorly compared to methods 3 and 4 both in accuracy of predictions and computational efficiency. Our results demonstrate how novel experimental and computational approaches can improve the quality of model predictions in safety-critical applications.

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 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.

Rapid Characterization of hERG Channel Kinetics II: Temperature Dependence.

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 Q10 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 Q10 coefficients; it broadly follows a generalized, but not the standardly-used, Eyring relationship. We also demonstrate that experimental estimations of Q10 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.

Representation of Multiple Cellular Phenotypes Within Tissue-Level Simulations of Cardiac Electrophysiology.

Distinct electrophysiological phenotypes are exhibited by biological cells that have differentiated into particular cell types. The usual approach when simulating the cardiac electrophysiology of tissue that includes different cell types is to model the different cell types as occupying spatially distinct yet coupled regions. Instead, we model the electrophysiology of well-mixed cells by using homogenisation to derive an extension to the commonly used monodomain or bidomain equations. These new equations permit spatial variations in the distribution of the different subtypes of cells and will reduce the computational demands of solving the governing equations. We validate the homogenisation computationally, and then use the new model to explain some experimental observations from stem cell-derived cardiomyocyte monolayers.

Sinusoidal voltage protocols for rapid characterisation of ion channel kinetics.

KEY POINTS: Ion current kinetics are commonly represented by current-voltage relationships, time constant-voltage relationships and subsequently mathematical models fitted to these. These experiments take substantial time, which means they are rarely performed in the same cell. Rather than traditional square-wave voltage clamps, we fitted a model to the current evoked by a novel sum-of-sinusoids voltage clamp that was only 8 s long. Short protocols that can be performed multiple times within a single cell will offer many new opportunities to measure how ion current kinetics are affected by changing conditions. The new model predicts the current under traditional square-wave protocols well, with better predictions of underlying currents than literature models. The current under a novel physiologically relevant series of action potential clamps is predicted extremely well. The short sinusoidal protocols allow a model to be fully fitted to individual cells, allowing us to examine cell-cell variability in current kinetics for the first time. ABSTRACT: Understanding the roles of ion currents is crucial to predict the action of pharmaceuticals and mutations in different scenarios, and thereby to guide clinical interventions in the heart, brain and other electrophysiological systems. Our ability to predict how ion currents contribute to cellular electrophysiology is in turn critically dependent on our characterisation of ion channel kinetics - the voltage-dependent rates of transition between open, closed and inactivated channel states. We present a new method for rapidly exploring and characterising ion channel kinetics, applying it to the hERG potassium channel as an example, with the aim of generating a quantitatively predictive representation of the ion current. We fitted a mathematical model to currents evoked by a novel 8 second sinusoidal voltage clamp in CHO cells overexpressing hERG1a. The model was then used to predict over 5 minutes of recordings in the same cell in response to further protocols: a series of traditional square step voltage clamps, and also a novel voltage clamp comprising a collection of physiologically relevant action potentials. We demonstrate that we can make predictive cell-specific models that outperform the use of averaged data from a number of different cells, and thereby examine which changes in gating are responsible for cell-cell variability in current kinetics. Our technique allows rapid collection of consistent and high quality data, from single cells, and produces more predictive mathematical ion channel models than traditional approaches.

Systems Toxicology: Real World Applications and Opportunities.

Systems Toxicology aims to change the basis of how adverse biological effects of xenobiotics are characterized from empirical end points to describing modes of action as adverse outcome pathways and perturbed networks. Toward this aim, Systems Toxicology entails the integration of in vitro and in vivo toxicity data with computational modeling. This evolving approach depends critically on data reliability and relevance, which in turn depends on the quality of experimental models and bioanalysis techniques used to generate toxicological data. Systems Toxicology involves the use of large-scale data streams ("big data"), such as those derived from omics measurements that require computational means for obtaining informative results. Thus, integrative analysis of multiple molecular measurements, particularly acquired by omics strategies, is a key approach in Systems Toxicology. In recent years, there have been significant advances centered on in vitro test systems and bioanalytical strategies, yet a frontier challenge concerns linking observed network perturbations to phenotypes, which will require understanding pathways and networks that give rise to adverse responses. This summary perspective from a 2016 Systems Toxicology meeting, an international conference held in the Alps of Switzerland, describes the limitations and opportunities of selected emerging applications in this rapidly advancing field. Systems Toxicology aims to change the basis of how adverse biological effects of xenobiotics are characterized, from empirical end points to pathways of toxicity. This requires the integration of in vitro and in vivo data with computational modeling. Test systems and bioanalytical technologies have made significant advances, but ensuring data reliability and relevance is an ongoing concern. The major challenge facing the new pathway approach is determining how to link observed network perturbations to phenotypic toxicity.

The Cardiac Electrophysiology Web Lab.

Computational modeling of cardiac cellular electrophysiology has a long history, and many models are now available for different species, cell types, and experimental preparations. This success brings with it a challenge: how do we assess and compare the underlying hypotheses and emergent behaviors so that we can choose a model as a suitable basis for a new study or to characterize how a particular model behaves in different scenarios? We have created an online resource for the characterization and comparison of electrophysiological cell models in a wide range of experimental scenarios. The details of the mathematical model (quantitative assumptions and hypotheses formulated as ordinary differential equations) are separated from the experimental protocol being simulated. Each model and protocol is then encoded in computer-readable formats. A simulation tool runs virtual experiments on models encoded in CellML, and a website (https://chaste.cs.ox.ac.uk/WebLab) provides a friendly interface, allowing users to store and compare results. The system currently contains a sample of 36 models and 23 protocols, including current-voltage curve generation, action potential properties under steady pacing at different rates, restitution properties, block of particular channels, and hypo-/hyperkalemia. This resource is publicly available, open source, and free, and we invite the community to use it and become involved in future developments. Investigators interested in comparing competing hypotheses using models can make a more informed decision, and those developing new models can upload them for easy evaluation under the existing protocols, and even add their own protocols.

Uncertainty and variability in models of the cardiac action potential: Can we build trustworthy models?

Cardiac electrophysiology models have been developed for over 50years, and now include detailed descriptions of individual ion currents and sub-cellular calcium handling. It is commonly accepted that there are many uncertainties in these systems, with quantities such as ion channel kinetics or expression levels being difficult to measure or variable between samples. Until recently, the original approach of describing model parameters using single values has been retained, and consequently the majority of mathematical models in use today provide point predictions, with no associated uncertainty. In recent years, statistical techniques have been developed and applied in many scientific areas to capture uncertainties in the quantities that determine model behaviour, and to provide a distribution of predictions which accounts for this uncertainty. In this paper we discuss this concept, which is termed uncertainty quantification, and consider how it might be applied to cardiac electrophysiology models. We present two case studies in which probability distributions, instead of individual numbers, are inferred from data to describe quantities such as maximal current densities. Then we show how these probabilistic representations of model parameters enable probabilities to be placed on predicted behaviours. We demonstrate how changes in these probability distributions across data sets offer insight into which currents cause beat-to-beat variability in canine APs. We conclude with a discussion of the challenges that this approach entails, and how it provides opportunities to improve our understanding of electrophysiology.

Uncertainty and variability in computational and mathematical models of cardiac physiology.

KEY POINTS: Mathematical and computational models of cardiac physiology have been an integral component of cardiac electrophysiology since its inception, and are collectively known as the Cardiac Physiome. We identify and classify the numerous sources of variability and uncertainty in model formulation, parameters and other inputs that arise from both natural variation in experimental data and lack of knowledge. The impact of uncertainty on the outputs of Cardiac Physiome models is not well understood, and this limits their utility as clinical tools. We argue that incorporating variability and uncertainty should be a high priority for the future of the Cardiac Physiome. We suggest investigating the adoption of approaches developed in other areas of science and engineering while recognising unique challenges for the Cardiac Physiome; it is likely that novel methods will be necessary that require engagement with the mathematics and statistics community. ABSTRACT: The Cardiac Physiome effort is one of the most mature and successful applications of mathematical and computational modelling for describing and advancing the understanding of physiology. After five decades of development, physiological cardiac models are poised to realise the promise of translational research via clinical applications such as drug development and patient-specific approaches as well as ablation, cardiac resynchronisation and contractility modulation therapies. For models to be included as a vital component of the decision process in safety-critical applications, rigorous assessment of model credibility will be required. This White Paper describes one aspect of this process by identifying and classifying sources of variability and uncertainty in models as well as their implications for the application and development of cardiac models. We stress the need to understand and quantify the sources of variability and uncertainty in model inputs, and the impact of model structure and complexity and their consequences for predictive model outputs. We propose that the future of the Cardiac Physiome should include a probabilistic approach to quantify the relationship of variability and uncertainty of model inputs and outputs.

Computational cardiology and risk stratification for sudden cardiac death: one of the grand challenges for cardiology in the 21st century.

Risk stratification in the context of sudden cardiac death has been acknowledged as one of the major challenges facing cardiology for the past four decades. In recent years, the advent of high performance computing has facilitated organ-level simulation of the heart, meaning we can now examine the causes, mechanisms and impact of cardiac dysfunction in silico. As a result, computational cardiology, largely driven by the Physiome project, now stands at the threshold of clinical utility in regards to risk stratification and treatment of patients at risk of sudden cardiac death. In this white paper, we outline a roadmap of what needs to be done to make this translational step, using the relatively well-developed case of acquired or drug-induced long QT syndrome as an exemplar case.

Cardiac tissue slices: preparation, handling, and successful optical mapping.

Cardiac tissue slices are becoming increasingly popular as a model system for cardiac electrophysiology and pharmacology research and development. Here, we describe in detail the preparation, handling, and optical mapping of transmembrane potential and intracellular free calcium concentration transients (CaT) in ventricular tissue slices from guinea pigs and rabbits. Slices cut in the epicardium-tangential plane contained well-aligned in-slice myocardial cell strands ("fibers") in subepicardial and midmyocardial sections. Cut with a high-precision slow-advancing microtome at a thickness of 350 to 400 µm, tissue slices preserved essential action potential (AP) properties of the precutting Langendorff-perfused heart. We identified the need for a postcutting recovery period of 36 min (guinea pig) and 63 min (rabbit) to reach 97.5% of final steady-state values for AP duration (APD) (identified by exponential fitting). There was no significant difference between the postcutting recovery dynamics in slices obtained using 2,3-butanedione 2-monoxime or blebistatin as electromechanical uncouplers during the cutting process. A rapid increase in APD, seen after cutting, was caused by exposure to ice-cold solution during the slicing procedure, not by tissue injury, differences in uncouplers, or pH-buffers (bicarbonate; HEPES). To characterize intrinsic patterns of CaT, AP, and conduction, a combination of multipoint and field stimulation should be used to avoid misinterpretation based on source-sink effects. In summary, we describe in detail the preparation, mapping, and data analysis approaches for reproducible cardiac tissue slice-based investigations into AP and CaT dynamics.