Brodie A. J. Lawson

Unlocking data sets by calibrating populations of models to data density: A study in atrial electrophysiology

We describe a statistically informed calibration of in silico populations to explore variability in complex systems. The understanding of complex physical or biological systems nearly always requires a characterization of the variability that underpins these processes. In addition, the data used to calibrate these models may also often exhibit considerable variability. A recent approach to deal with these issues has been to calibrate populations of models (POMs), multiple copies of a single mathematical model but with different parameter values, in response to experimental data. To date, this calibration has been largely limited to selecting models that produce outputs that fall within the ranges of the data set, ignoring any trends that might be present in the data. We present here a novel and general methodology for calibrating POMs to the distributions of a set of measured values in a data set. We demonstrate our technique using a data set from a cardiac electrophysiology study based on the differences in atrial action potential readings between patients exhibiting sinus rhythm (SR) or chronic atrial fibrillation (cAF) and the Courtemanche-Ramirez-Nattel model for human atrial action potentials. Not only does our approach accurately capture the variability inherent in the experimental population, but we also demonstrate how the POMs that it produces may be used to extract additional information from the data used for calibration, including improved identification of the differences underlying stratified data. We also show how our approach allows different hypotheses regarding the variability in complex systems to be quantitatively compared.

Sampling methods for exploring between-subject variability in cardiac electrophysiology experiments

Between-subject and within-subject variability is ubiquitous in biology and physiology, and understanding and dealing with this is one of the biggest challenges in medicine. At the same time, it is difficult to investigate this variability by experiments alone. A recent modelling and simulation approach, known as population of models (POM), allows this exploration to take place by building a mathematical model consisting of multiple parameter sets calibrated against experimental data. However, finding such sets within a high-dimensional parameter space of complex electrophysiological models is computationally challenging. By placing the POM approach within a statistical framework, we develop a novel and efficient algorithm based on sequential Monte Carlo (SMC). We compare the SMC approach with Latin hypercube sampling (LHS), a method commonly adopted in the literature for obtaining the POM, in terms of efficiency and output variability in the presence of a drug block through an in-depth investigation via the Beeler–Reuter cardiac electrophysiological model. We show improved efficiency for SMC that produces similar responses to LHS when making out-of-sample predictions in the presence of a simulated drug block. Finally, we show the performance of our approach on a complex atrial electrophysiological model, namely the Courtemanche–Ramirez–Nattel model.

Space-limited mitosis in the Glazier–Graner–Hogeweg Model

The Glazier–Graner–Hogeweg (GGH) model is a cellular automata framework for representing the time evolution of cellular systems, appealing because unlike many other individual-cell-based models it dynamically simulates changes in cell shape and size. Proliferation has seen some implementation into this modelling framework, but without consensus in the literature as to how this behaviour is best represented. Additionally, the majority of published GGH model implementations which feature proliferation do so in order to simulate a certain biological situation where mitosis is important, but without analysis of how these proliferation routines operate on a fundamental level. Here, a method of proliferation for the GGH model which uses separate cell phenotypes to differentiate cells which have entered or just left the mitotic phase of the cell cycle is presented and demonstrated to correctly predict logistic growth on a macroscopic scale (in accordance with experimental evidence). Comparisons between model simulations and the generalised logistic growth model provide an interpretation of the latter’s ‘shape parameter’, and the proliferation routine used here is shown to offer the modeller somewhat predictable control over the proliferation rate, important for ensuring temporal consistency between different cellular behaviours in the model. All results are found to be insensitive to the inclusion of active cell motility. The implications of these simulated proliferation assays towards problems in cell biology are also discussed.

Slow Recovery of Excitability Increases Ventricular Fibrillation Risk as Identified by Emulation

Purpose: Rotor stability and meandering are key mechanisms determining and sustaining cardiac fibrillation, with important implications for anti-arrhythmic drug development. However, little is yet known on how rotor dynamics are modulated by variability in cellular electrophysiology, particularly on kinetic properties of ion channel recovery. Methods: We propose a novel emulation approach, based on Gaussian process regression augmented with machine learning, for data enrichment, automatic detection, classification, and analysis of re-entrant biomarkers in cardiac tissue. More than 5,000 monodomain simulations of long-lasting arrhythmic episodes with Fenton-Karma ionic dynamics, further enriched by emulation to 80 million electrophysiological scenarios, were conducted to investigate the role of variability in ion channel densities and kinetics in modulating rotor-driven arrhythmic behavior. Results: Our methods predicted the class of excitation behavior with classification accuracy up to 96%, and emulation effectively predicted frequency, stability, and spatial biomarkers of functional re-entry. We demonstrate that the excitation wavelength interpretation of re-entrant behavior hides critical information about rotor persistence and devolution into fibrillation. In particular, whereas action potential duration directly modulates rotor frequency and meandering, critical windows of excitability are identified as the main determinants of breakup. Further novel electrophysiological insights of particular relevance for ventricular arrhythmias arise from our multivariate analysis, including the role of incomplete activation of slow inward currents in mediating tissue rate-dependence and dispersion of repolarization, and the emergence of slow recovery of excitability as a significant promoter of this mechanism of dispersion and increased arrhythmic risk. Conclusions: Our results mechanistically explain pro-arrhythmic effects of class Ic anti-arrhythmics in the ventricles despite their established role in the pharmacological management of atrial fibrillation. This is mediated by their slow recovery of excitability mode of action, promoting incomplete activation of slow inward currents and therefore increased dispersion of repolarization, given the larger influence of these currents in modulating the action potential in the ventricles compared to the atria. These results exemplify the potential of emulation techniques in elucidating novel mechanisms of arrhythmia and further application to cardiac electrophysiology.