Cardiac cell modelling: observations from the heart of the cardiac physiome project

Reconstructing ventricular cardiomyocyte dynamics and parameter estimation using data assimilation

Cardiac ventricular myocyte action potential dynamics are regulated by intricate and nonlinear interactions between the cell transmembrane potential and ionic currents and concentrations. Present technology limits the ability to measure transmembrane potential and multiple ionic currents simultaneously, which narrows the scope of experiments to provide a complete snapshot of the cardiac myocyte state. This limitation presents an obstacle for understanding how perturbations can trigger arrhythmias and more broadly how the myocyte responds to different conditions, such as changes in pacing rate or responses to drug treatment. In this study, we demonstrate that a data-assimilation approach can successfully reconstruct and predict the dynamics of a heterogeneous virtual cardiac ventricular myocyte population in the presence of parameter uncertainty. A population of heterogeneous cardiac ventricular myocytes is generated by varying ionic current conductance parameters, and additional observational uncertainty is mimicked by the addition of Gaussian noise to the transmembrane potential. We demonstrate that the data-assimilation approach accurately reconstructs transmembrane potential, with error less than the magnitude of the noise. Further, the data-assimilation approach successfully estimates the conductances of ionic currents generally with high accuracy and requiring low computational time. As a proof of concept, we apply the data-assimilation approach to reconstruct action potential dynamics from optical mapping experiments in an ex vivo isolated guinea pig heart. Critically, we demonstrate that the ionic conductance parameters estimated from a recording at one pacing frequency can accurately predict action potential dynamics at different rates.

Fast interactive simulations of cardiac electrical activity in anatomically accurate heart structures by compressing sparse uniform cartesian grids

BACKGROUND AND OBJECTIVE

Numerical simulations are valuable tools for studying cardiac arrhythmias. Not only do they complement experimental studies, but there is also an increasing expectation for their use in clinical applications to guide patient-specific procedures. However, numerical studies that solve the reaction-diffusion equations describing cardiac electrical activity remain challenging to set up, are time-consuming, and in many cases, are prohibitively computationally expensive for long studies. The computational cost of cardiac simulations of complex models on anatomically accurate structures necessitates parallel computing. Graphics processing units (GPUs), which have thousands of cores, have been introduced as a viable technology for carrying out fast cardiac simulations, sometimes including real-time interactivity. Our main objective is to increase the performance and accuracy of such GPU implementations while conserving computational resources.

METHODS

In this work, we present a compression algorithm that can be used to conserve GPU memory and improve efficiency by managing the sparsity that is inherent in using Cartesian grids to represent cardiac structures directly obtained from high-resolution MRI and mCT scans. Furthermore, we present a discretization scheme that includes the cross-diagonal terms in the computational cell to increase numerical accuracy, which is especially important for simulating thin tissue sections without the need for costly mesh refinement.

RESULTS

Interactive WebGL simulations of atrial/ventricular structures (on PCs, laptops, tablets, and phones) demonstrate the algorithm's ability to reduce memory demand by an order of magnitude and achieve calculations up to 20x faster. We further showcase its superiority in slender tissues and validate results against experiments performed in live explanted human hearts.

CONCLUSIONS

In this work, we present a compression algorithm that accelerates electrical activity simulations on realistic anatomies by an order of magnitude (up to 20x), thereby allowing the use of finer grid resolutions while conserving GPU memory. Additionally, improved accuracy is achieved through cross-diagonal terms, which are essential for thin tissues, often found in heart structures such as pectinate muscles and trabeculae, as well as Purkinje fibers. Our method enables interactive simulations with even interactive domain boundary manipulation (unlike finite element/volume methods). Finally, agreement with experiments and ease of mesh import into WebGL paves the way for virtual cohorts and digital twins, aiding arrhythmia analysis and personalized therapies.

Ionic Mechanisms of Propagated Repolarization in a One-Dimensional Strand of Human Ventricular Myocyte Model

Although repolarization has been suggested to propagate in cardiac tissue both theoretically and experimentally, it has been challenging to estimate how and to what extent the propagation of repolarization contributes to relaxation because repolarization only occurs in the course of membrane excitation in normal hearts. We established a mathematical model of a 1D strand of 600 myocytes stabilized at an equilibrium potential near the plateau potential level by introducing a sustained component of the late sodium current (INaL). By applying a hyperpolarizing stimulus to a small part of the strand, we succeeded in inducing repolarization which propagated along the strand at a velocity of 1~2 cm/s. The ionic mechanisms responsible for repolarization at the myocyte level, i.e., the deactivation of both the INaL and the L-type calcium current (ICaL), and the activation of the rapid component of delayed rectifier potassium current (IKr) and the inward rectifier potassium channel (IK1), were found to be important for the propagation of repolarization in the myocyte strand. Using an analogy with progressive activation of the sodium current (INa) in the propagation of excitation, regenerative activation of the predominant magnitude of IK1 makes the myocytes at the wave front start repolarization in succession through the electrical coupling via gap junction channels.

Diversity of cells and signals in the cardiovascular system

This white paper is the outcome of the seventh UC Davis Cardiovascular Research Symposium on Systems Approach to Understanding Cardiovascular Disease and Arrhythmia. This biannual meeting aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2022 Symposium was 'Cell Diversity in the Cardiovascular System, cell-autonomous and cell-cell signalling'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies, and challenges in examining cell and signal diversity, co-ordination and interrelationships involved in cardiovascular function. This paper originates from the topics of formal presentations and informal discussions from the Symposium, which aimed to develop a holistic view of how the multiple cell types in the cardiovascular system integrate to influence cardiovascular function, disease progression and therapeutic strategies. The first section describes the major cell types (e.g. cardiomyocytes, vascular smooth muscle and endothelial cells, fibroblasts, neurons, immune cells, etc.) and the signals involved in cardiovascular function. The second section emphasizes the complexity at the subcellular, cellular and system levels in the context of cardiovascular development, ageing and disease. Finally, the third section surveys the technological innovations that allow the interrogation of this diversity and advancing our understanding of the integrated cardiovascular function and dysfunction.

The lost art of mathematical modelling

We provide a critique of mathematical biology in light of rapid developments in modern machine learning. We argue that out of the three modelling activities - (1) formulating models; (2) analysing models; and (3) fitting or comparing models to data - inherent to mathematical biology, researchers currently focus too much on activity (2) at the cost of (1). This trend, we propose, can be reversed by realising that any given biological phenomenon can be modelled in an infinite number of different ways, through the adoption of an pluralistic approach, where we view a system from multiple, different points of view. We explain this pluralistic approach using fish locomotion as a case study and illustrate some of the pitfalls - universalism, creating models of models, etc. - that hinder mathematical biology. We then ask how we might rediscover a lost art: that of creative mathematical modelling.

Assay for evaluation of proarrhythmic effects of herbal products: Case study with 12 Evodia preparations

Guidelines for preclinical drug development reduce the occurrence of arrhythmia-related side effects. Besides ample evidence for the presence of arrhythmogenic substances in plants, there is no consensus on a research strategy for the evaluation of proarrhythmic effects of herbal products. Here, we propose a cardiac safety assay for the detection of proarrhythmic effects of plant extracts based on the experimental approaches described in the Comprehensive In vitro Proarrhythmia Assay (CiPA). Microelectrode array studies (MEAs) and voltage sensing optical technique on human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) were combined with ionic current measurements in mammalian cell lines, In-silico simulations of cardiac action potentials (APs) and statistic regression analysis. Proarrhythmic effects of 12 Evodia preparations, containing different amounts of the hERG inhibitors dehydroevodiamine (DHE) and hortiamine were analysed. Extracts produced different prolongation of the AP, occurrence of early after depolarisations and triangulation of the AP in hiPSC-CMs depending on the contents of the hERG inhibitors. DHE and hortiamine dose-dependently prolonged the field potential duration in hiPSC-CMs studied with MEAs. In-silico simulations of ventricular AP support a scenario where proarrhythmic effects of Evodia extracts are predominantly caused by the content of the selective hERG inhibitors. Statistic regression analysis revealed a high torsadogenic risk for both compounds that was comparable to drugs assigned to the high-risk category in a CiPA study.

Computational cardiac physiology for new modelers: Origins, foundations, and future

Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.

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.

Bayesian Calibration of Electrophysiology Models Using Restitution Curve Emulators

Calibration of cardiac electrophysiology models is a fundamental aspect of model personalization for predicting the outcomes of cardiac therapies, simulation testing of device performance for a range of phenotypes, and for fundamental research into cardiac function. Restitution curves provide information on tissue function and can be measured using clinically feasible measurement protocols. We introduce novel "restitution curve emulators" as probabilistic models for performing model exploration, sensitivity analysis, and Bayesian calibration to noisy data. These emulators are built by decomposing restitution curves using principal component analysis and modeling the resulting coordinates with respect to model parameters using Gaussian processes. Restitution curve emulators can be used to study parameter identifiability via sensitivity analysis of restitution curve components and rapid inference of the posterior distribution of model parameters given noisy measurements. Posterior uncertainty about parameters is critical for making predictions from calibrated models, since many parameter settings can be consistent with measured data and yet produce very different model behaviors under conditions not effectively probed by the measurement protocols. Restitution curve emulators are therefore promising probabilistic tools for calibrating electrophysiology models.

Accelerating the convergence to a limit cycle in 3D cardiac electromechanical simulations through a data-driven 0D emulator

The results of numerical simulations of cardiac electromechanics are typically characterized by a long transient before reaching a periodic solution known as limit cycle. This yields a serious computational overhead, as the only clinically relevant output is associated with such limit cycle. To accelerate the convergence to the limit cycle, we propose a strategy based on a surrogate model, wherein the computationally demanding 3D components are replaced by a 0D emulator, built through an automated data-driven algorithm on the basis of pressure-volume transients of as few as three heartbeats simulated with the 3D model. The 0D emulator, consisting of a time-dependent pressure-volume relationship, can provide the 3D model with an initial guess, such that in just two heartbeats a solution is reached that is as close to the limit cycle as the one obtained after more than 20 heartbeats with the 3D model. The 0D emulator is also recommended in many-query settings (e.g. when performing sensitivity analysis, parameter estimation and uncertainty quantification), that call for the repeated solution of the model for different values of the parameters. Indeed, the construction of the emulator does not have to be repeated when the parameters of the circulation model it is coupled with vary. Finally, should the parameters of the 3D electromechanical model vary as well, we propose a parametric emulator, obtained by interpolation of emulators constructed for given values of the parameters. This paper is accompanied by a Python library implementing the proposed algorithm, open to integration with existing cardiac solvers.

Data integration for the numerical simulation of cardiac electrophysiology

The increasing availability of extensive and accurate clinical data is rapidly shaping cardiovascular care by improving the understanding of physiological and pathological mechanisms of the cardiovascular system and opening new frontiers in designing therapies and interventions. In this direction, mathematical and numerical models provide a complementary relevant tool, able not only to reproduce patient-specific clinical indicators but also to predict and explore unseen scenarios. With this goal, clinical data are processed and provided as inputs to the mathematical model, which quantitatively describes the physical processes that occur in the cardiac tissue. In this paper, the process of integration of clinical data and mathematical models is discussed. Some challenges and contributions in the field of cardiac electrophysiology are reported.

Novel Approach for EKG Signals Analysis Based on Markovian and Non-Markovian Fractalization Type in Scale Relativity Theory

Two distinct operational procedures are proposed for diagnosis and tracking of heart disease evolution (in particular atrial fibrillations). The first procedure, based on the application of non-linear dynamic methods (strange attractors, skewness, kurtosis, histograms, Lyapunov exponent, etc.) analyzes the electrical activity of the heart (electrocardiogram signals). The second procedure, based on multifractalization through Markovian and non-Markovian-type stochasticizations in the framework of the scale relativity theory, reconstructs any type of EKG signal by means of harmonic mappings from the usual space to the hyperbolic one. These mappings mime various scale transitions by differential geometries, in Riemann spaces with symmetries of SL(2R)-type. Then, the two operational procedures are not mutually exclusive, but rather become complementary, through their finality, which is gaining valuable information concerning fibrillation crises. As such, the author’s proposed method could be used for developing new models for medical diagnosis and evolution tracking of heart diseases (patterns dynamics, signal reconstruction, etc.).

Modeling and Analysis of Cardiac Hybrid Cellular Automata via GPU-Accelerated Monte Carlo Simulation

The heart consists of a complex network of billions of cells. Under physiological conditions, cardiac cells propagate electrical signals in space, generating the heartbeat in a synchronous and coordinated manner. When such a synchronization fails, life-threatening events can arise. The inherent complexity of the underlying nonlinear dynamics and the large number of biological components involved make the modeling and the analysis of electrophysiological properties in cardiac tissue still an open challenge. We consider here a Hybrid Cellular Automata (HCA) approach modeling the cardiac cell-cell membrane resistance with a free variable. We show that the modeling approach can reproduce important and complex spatiotemporal properties paving the ground for promising future applications. We show how GPU-based technology can considerably accelerate the simulation and the analysis. Furthermore, we study the cardiac behavior within a unidimensional domain considering inhomogeneous resistance and we perform a Monte Carlo analysis to evaluate our approach.

Biophysically detailed mathematical models of multiscale cardiac active mechanics

We propose four novel mathematical models, describing the microscopic mechanisms of force generation in the cardiac muscle tissue, which are suitable for multiscale numerical simulations of cardiac electromechanics. Such models are based on a biophysically accurate representation of the regulatory and contractile proteins in the sarcomeres. Our models, unlike most of the sarcomere dynamics models that are available in the literature and that feature a comparable richness of detail, do not require the time-consuming Monte Carlo method for their numerical approximation. Conversely, the models that we propose only require the solution of a system of PDEs and/or ODEs (the most reduced of the four only involving 20 ODEs), thus entailing a significant computational efficiency. By focusing on the two models that feature the best trade-off between detail of description and identifiability of parameters, we propose a pipeline to calibrate such parameters starting from experimental measurements available in literature. Thanks to this pipeline, we calibrate these models for room-temperature rat and for body-temperature human cells. We show, by means of numerical simulations, that the proposed models correctly predict the main features of force generation, including the steady-state force-calcium and force-length relationships, the length-dependent prolongation of twitches and increase of peak force, the force-velocity relationship. Moreover, they correctly reproduce the Frank-Starling effect, when employed in multiscale 3D numerical simulation of cardiac electromechanics.

Active Force Generation in Cardiac Muscle Cells: Mathematical Modeling and Numerical Simulation of the Actin-Myosin Interaction

Cardiac in silico numerical simulations are based on mathematical models describing the physical processes involved in the heart function. In this review paper, we critically survey biophysically-detailed mathematical models describing the subcellular mechanisms behind the generation of active force, that is the process by which the chemical energy of ATP (adenosine triphosphate) is transformed into mechanical work, thus making the muscle tissue contract. While presenting these models, that feature different levels of biophysical detail, we analyze the trade-off between the accuracy in the description of the subcellular mechanisms and the number of parameters that need to be estimated from experiments. Then, we focus on a generalized version of the classic Huxley model, built on the basis of models available in the literature, that is able to reproduce the main experimental characterizations associated to the time scales typical of a heartbeat-such as the force-velocity relationship and the tissue stiffness in response to small steps-featuring only four independent parameters. Finally, we show how those parameters can be calibrated starting from macroscopic measurements available from experiments.

Calibration of ionic and cellular cardiac electrophysiology models

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.

Reducing complexity and unidentifiability when modelling human atrial cells

Mathematical models of a cellular action potential (AP) in cardiac modelling have become increasingly complex, particularly in gating kinetics, which control the opening and closing of individual ion channel currents. As cardiac models advance towards use in personalized medicine to inform clinical decision-making, it is critical to understand the uncertainty hidden in parameter estimates from their calibration to experimental data. This study applies approximate Bayesian computation to re-calibrate the gating kinetics of four ion channels in two existing human atrial cell models to their original datasets, providing a measure of uncertainty and indication of potential issues with selecting a single unique value given the available experimental data. Two approaches are investigated to reduce the uncertainty present: re-calibrating the models to a more complete dataset and using a less complex formulation with fewer parameters to constrain. The re-calibrated models are inserted back into the full cell model to study the overall effect on the AP. The use of more complete datasets does not eliminate uncertainty present in parameter estimates. The less complex model, particularly for the fast sodium current, gave a better fit to experimental data alongside lower parameter uncertainty and improved computational speed. 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'.

An audit of uncertainty in multi-scale cardiac electrophysiology models

Models of electrical activation and recovery in cardiac cells and tissue have become valuable research tools, and are beginning to be used in safety-critical applications including guidance for clinical procedures and for drug safety assessment. As a consequence, there is an urgent need for a more detailed and quantitative understanding of the ways that uncertainty and variability influence model predictions. In this paper, we review the sources of uncertainty in these models at different spatial scales, discuss how uncertainties are communicated across scales, and begin to assess their relative importance. We conclude by highlighting important challenges that continue to face the cardiac modelling community, identifying open questions, and making recommendations for future studies. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Reducing complexity and unidentifiability when modelling human atrial cells

Mathematical models of a cellular action potential (AP) in cardiac modelling have become increasingly complex, particularly in gating kinetics, which control the opening and closing of individual ion channel currents. As cardiac models advance towards use in personalized medicine to inform clinical decision-making, it is critical to understand the uncertainty hidden in parameter estimates from their calibration to experimental data. This study applies approximate Bayesian computation to re-calibrate the gating kinetics of four ion channels in two existing human atrial cell models to their original datasets, providing a measure of uncertainty and indication of potential issues with selecting a single unique value given the available experimental data. Two approaches are investigated to reduce the uncertainty present: re-calibrating the models to a more complete dataset and using a less complex formulation with fewer parameters to constrain. The re-calibrated models are inserted back into the full cell model to study the overall effect on the AP. The use of more complete datasets does not eliminate uncertainty present in parameter estimates. The less complex model, particularly for the fast sodium current, gave a better fit to experimental data alongside lower parameter uncertainty and improved computational speed. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Sensitivity and Uncertainty Analysis of Two Human Atrial Cardiac Cell Models Using Gaussian Process Emulators

Biophysically detailed cardiac cell models reconstruct the action potential and calcium dynamics of cardiac myocytes. They aim to capture the biophysics of current flow through ion channels, pumps, and exchangers in the cell membrane, and are highly detailed. However, the relationship between model parameters and model outputs is difficult to establish because the models are both complex and non-linear. The consequences of uncertainty and variability in model parameters are therefore difficult to determine without undertaking large numbers of model evaluations. The aim of the present study was to demonstrate how sensitivity and uncertainty analysis using Gaussian process emulators can be used for a systematic and quantitive analysis of biophysically detailed cardiac cell models. We selected the Courtemanche and Maleckar models of the human atrial action potential for analysis because these models describe a similar set of currents, with different formulations. In our approach Gaussian processes emulate the main features of the action potential and calcium transient. The emulators were trained with a set of design data comprising samples from parameter space and corresponding model outputs, initially obtained from 300 model evaluations. Variance based sensitivity indices were calculated using the emulators, and first order and total effect indices were calculated for each combination of parameter and output. The differences between the first order and total effect indices indicated that the effect of interactions between parameters was small. A second set of emulators were then trained using a new set of design data with a subset of the model parameters with a sensitivity index of more than 0.1 (10%). This second stage analysis enabled comparison of mechanisms in the two models. The second stage sensitivity indices enabled the relationship between the L-type Ca 2+ current and the action potential plateau to be quantified in each model. Our quantitative analysis predicted that changes in maximum conductance of the ultra-rapid K + channel I Kur would have opposite effects on action potential duration in the two models, and this prediction was confirmed by additional simulations. This study has demonstrated that Gaussian process emulators are an effective tool for sensitivity and uncertainty analysis of biophysically detailed cardiac cell models.

Predicting critical drug concentrations and torsadogenic risk using a multiscale exposure-response simulator

Torsades de pointes is a serious side effect of many drugs that can trigger sudden cardiac death, even in patients with structurally normal hearts. Torsadogenic risk has traditionally been correlated with the blockage of a specific potassium channel and a prolonged recovery period in the electrocardiogram. However, the precise mechanisms by which single channel block translates into heart rhythm disorders remain incompletely understood. Here we establish a multiscale exposure-response simulator that converts block-concentration characteristics from single cell recordings into three-dimensional excitation profiles and electrocardiograms to rapidly assess torsadogenic risk. For the drug dofetilide, we characterize the QT interval and heart rate at different drug concentrations and identify the critical concentration at the onset of torsades de pointes: For dofetilide concentrations of 2x, 3x, and 4x, as multiples of the free plasma concentration Cmax = 2.1 nM, the QT interval increased by +62.0%, +71.2%, and +82.3% compared to baseline, and the heart rate changed by -21.7%, -23.3%, and +88.3%. The last number indicates that, at the critical concentration of 4x, the heart spontaneously developed an episode of a torsades-like arrhythmia. Strikingly, this critical drug concentration is higher than the concentration estimated from early afterdepolarizations in single cells and lower than in one-dimensional cable models. Our results highlight the importance of whole heart modeling and explain, at least in part, why current regulatory paradigms often fail to accurately quantify the pro-arrhythmic potential of a drug. Our exposure-response simulator could provide a more mechanistic assessment of pro-arrhythmic risk and help establish science-based guidelines to reduce rhythm disorders, design safer drugs, and accelerate drug development.

Understanding the Beat-to-Beat Variations of P-Waves Morphologies in AF Patients During Sinus Rhythm: A Scoping Review of the Atrial Simulation Studies

The remarkable advances in high-performance computing and the resulting increase of the computational power have the potential to leverage computational cardiology toward improving our understanding of the pathophysiological mechanisms of arrhythmias, such as Atrial Fibrillation (AF). In AF, a complex interaction between various triggers and the atrial substrate is considered to be the leading cause of AF initiation and perpetuation. In electrocardiography (ECG), P-wave is supposed to reflect atrial depolarization. It has been found that even during sinus rhythm (SR), multiple P-wave morphologies are present in AF patients with a history of AF, suggesting a higher dispersion of the conduction route in this population. In this scoping review, we focused on the mechanisms which modify the electrical substrate of the atria in AF patients, while investigating the existence of computational models that simulate the propagation of the electrical signal through different routes. The adopted review methodology is based on a structured analytical framework which includes the extraction of the keywords based on an initial limited bibliographic search, the extensive literature search and finally the identification of relevant articles based on the reference list of the studies. The leading mechanisms identified were classified according to their scale, spanning from mechanisms in the cell, tissue or organ level, and the produced outputs. The computational modeling approaches for each of the factors that influence the initiation and the perpetuation of AF are presented here to provide a clear overview of the existing literature. Several levels of categorization were adopted while the studies which aim to translate their findings to ECG phenotyping are highlighted. The results denote the availability of multiple models, which are appropriate under specific conditions. However, the consideration of complex scenarios taking into account multiple spatiotemporal scales, personalization of electrophysiological and anatomical models and the reproducibility in terms of ECG phenotyping has only partially been tackled so far.

Chaotic tip trajectories of a single spiral wave in the presence of heterogeneities

Spiral waves have been observed in a variety of physical, chemical, and biological systems. They play a major role in cardiac arrhythmias, including fibrillation, where the observed irregular activation patterns are generally thought to arise from the continuous breakup of multiple unstable spiral waves. Using spatially extended simulations of different electrophysiological models of cardiac tissue, we show that a single spiral wave in the presence of heterogeneities can display chaotic tip trajectories, consistent with fibrillation. We also show that the simulated spiral tip dynamics, including chaotic trajectories, can be captured by a simple particle model which only describes the dynamics of the spiral tip. This shows that spiral wave breakup, or interactions with other waves, are not necessary to initiate chaos in spiral waves.

Resonant model-A new paradigm for modeling an action potential of biological cells

Organ level simulation of bioelectric behavior in the body benefits from flexible and efficient models of cellular membrane potential. These computational organ and cell models can be used to study the impact of pharmaceutical drugs, test hypotheses, assess risk and for closed-loop validation of medical devices. To move closer to the real-time requirements of this modeling a new flexible Fourier based general membrane potential model, called as a Resonant model, is developed that is computationally inexpensive. The new model accurately reproduces non-linear potential morphologies for a variety of cell types. Specifically, the method is used to model human and rabbit sinoatrial node, human ventricular myocyte and squid giant axon electrophysiology. The Resonant models are validated with experimental data and with other published models. Dynamic changes in biological conditions are modeled with changing model coefficients and this approach enables ionic channel alterations to be captured. The Resonant model is used to simulate entrainment between competing sinoatrial node cells. These models can be easily implemented in low-cost digital hardware and an alternative, resource-efficient implementations of sine and cosine functions are presented and it is shown that a Fourier term is produced with two additions and a binary shift.

Real-time interactive simulations of large-scale systems on personal computers and cell phones: Toward patient-specific heart modeling and other applications

Cardiac dynamics modeling has been useful for studying and treating arrhythmias. However, it is a multiscale problem requiring the solution of billions of differential equations describing the complex electrophysiology of interconnected cells. Therefore, large-scale cardiac modeling has been limited to groups with access to supercomputers and clusters. Many areas of computational science face similar problems where computational costs are too high for personal computers so that supercomputers or clusters currently are necessary. Here, we introduce a new approach that makes high-performance simulation of cardiac dynamics and other large-scale systems like fluid flow and crystal growth accessible to virtually anyone with a modest computer. For cardiac dynamics, this approach will allow not only scientists and students but also physicians to use physiologically accurate modeling and simulation tools that are interactive in real time, thereby making diagnostics, research, and education available to a broader audience and pushing the boundaries of cardiac science.

Computational models in cardiology

The treatment of individual patients in cardiology practice increasingly relies on advanced imaging, genetic screening and devices. As the amount of imaging and other diagnostic data increases, paralleled by the greater capacity to personalize treatment, the difficulty of using the full array of measurements of a patient to determine an optimal treatment seems also to be paradoxically increasing. Computational models are progressively addressing this issue by providing a common framework for integrating multiple data sets from individual patients. These models, which are based on physiology and physics rather than on population statistics, enable computational simulations to reveal diagnostic information that would have otherwise remained concealed and to predict treatment outcomes for individual patients. The inherent need for patient-specific models in cardiology is clear and is driving the rapid development of tools and techniques for creating personalized methods to guide pharmaceutical therapy, deployment of devices and surgical interventions.

A short history of the development of mathematical models of cardiac mechanics

Cardiac mechanics plays a crucial role in atrial and ventricular function, in the regulation of growth and remodelling, in the progression of disease, and the response to treatment. The spatial scale of the critical mechanisms ranges from nm (molecules) to cm (hearts) with the fastest events occurring in milliseconds (molecular events) and the slowest requiring months (growth and remodelling). Due to its complexity and importance, cardiac mechanics has been studied extensively both experimentally and through mathematical models and simulation. Models of cardiac mechanics evolved from seminal studies in skeletal muscle, and developed into cardiac specific, species specific, human specific and finally patient specific calculations. These models provide a formal framework to link multiple experimental assays recorded over nearly 100 years into a single unified representation of cardiac function. This review first provides a summary of the proteins, physiology and anatomy involved in the generation of cardiac pump function. We then describe the evolution of models of cardiac mechanics starting with the early theoretical frameworks describing the link between sarcomeres and muscle contraction, transitioning through myosin-level models to calcium-driven systems, and ending with whole heart patient-specific models.

Early Afterdepolarisations Induced by an Enhancement in the Calcium Current

Excitable biological cells, such as cardiac muscle cells, can exhibit complex patterns of oscillations such as spiking and bursting. Moreover, it is well known that an enhancement in calcium currents may yield certain kind of cardiac arrhythmia, so-called early afterdepolarisations (EADs). The presence of EADs strongly correlates with the onset of dangerous cardiac arrhythmia. In this paper we study mathematically and numerically the dynamics of a cardiac muscle cell with respect to the calcium current by investigating a simplistic system of differential equations. For the study of this phenomena, we use bifurcation theory, numerical bifurcation analysis, geometric singular perturbation theory and computational methods to investigate a nonlinear multiple time scales system. It will turn out that EADs related to an enhanced calcium current are canard–induced and that we have to combine these theories to derive a better understanding of the dynamics behind EADs. Moreover, a suitable time scale separation argument determines the important and sensitive system parameters which are related to the occurrence of EADs.

Extinction dynamics of spiral defect chaos

Spatially extended excitable systems can exhibit spiral defect chaos (SDC) during which spiral waves continuously form and disappear. To address how this dynamical state terminates using simulations can be computationally challenging, especially for large systems. To circumvent this limitation, we treat the number of spiral waves as a stochastic population with a corresponding birth-death equation and use techniques from statistical physics to determine the mean episode duration of SDC. Motivated by cardiac fibrillation, during which the heart's electrical activity becomes disorganized and shows fragmenting spiral waves, we use generic models of cardiac electrophysiology. We show that the duration can be computed in minimal computational time and that it depends exponentially on domain size. Therefore, the approach can result in efficient and accurate predictions of mean episode duration which may be extended to more complex geometries and models.

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.

Fully coupled fluid-electro-mechanical model of the human heart for supercomputers

In this work, we present a fully coupled fluid-electro-mechanical model of a 50th percentile human heart. The model is implemented on Alya, the BSC multi-physics parallel code, capable of running efficiently in supercomputers. Blood in the cardiac cavities is modeled by the incompressible Navier-Stokes equations and an arbitrary Lagrangian-Eulerian (ALE) scheme. Electrophysiology is modeled with a monodomain scheme and the O'Hara-Rudy cell model. Solid mechanics is modeled with a total Lagrangian formulation for discrete strains using the Holzapfel-Ogden cardiac tissue material model. The three problems are simultaneously and bidirectionally coupled through an electromechanical feedback and a fluid-structure interaction scheme. In this paper, we present the scheme in detail and propose it as a computational cardiac workbench.

Modeling effects of voltage dependent properties of the cardiac muscarinic receptor on human sinus node function

The cardiac muscarinic receptor (M2R) regulates heart rate, in part, by modulating the acetylcholine (ACh) activated K⁺ current I_K,ACh through dissociation of G-proteins, that in turn activate K_ACh channels. Recently, M2Rs were noted to exhibit intrinsic voltage sensitivity, i.e. their affinity for ligands varies in a voltage dependent manner. The voltage sensitivity of M2R implies that the affinity for ACh (and thus the ACh effect) varies throughout the time course of a cardiac electrical cycle. The aim of this study was to investigate the contribution of M2R voltage sensitivity to the rate and shape of the human sinus node action potentials in physiological and pathophysiological conditions. We developed a Markovian model of the I_K,ACh modulation by voltage and integrated it into a computational model of human sinus node. We performed simulations with the integrated model varying ACh concentration and voltage sensitivity. Low ACh exerted a larger effect on I_K,ACh at hyperpolarized versus depolarized membrane voltages. This led to a slowing of the pacemaker rate due to an attenuated slope of phase 4 depolarization with only marginal effect on action potential duration and amplitude. We also simulated the theoretical effects of genetic variants that alter the voltage sensitivity of M2R. Modest negative shifts in voltage sensitivity, predicted to increase the affinity of the receptor for ACh, slowed the rate of phase 4 depolarization and slowed heart rate, while modest positive shifts increased heart rate. These simulations support our hypothesis that altered M2R voltage sensitivity contributes to disease and provide a novel mechanistic foundation to study clinical disorders such as atrial fibrillation and inappropriate sinus tachycardia.

Heart rate regulation is dependent upon a delicate interplay between parasympathetic and sympathetic nerve activity at the level of the sinus node. Acetylcholine slows the heart rate by activating the M2 muscarinic receptor (M2R) that, in turn, opens the acetylcholine-activated potassium channel (I_K,ACh) to slow the firing of the sinus node. Surprisingly, the M2R is sensitive to membrane potential and undergoes conformational changes throughout the cardiac action potential that alter the affinity for acetylcholine, with secondary consequences for I_K,ACh activity. Here, we investigated the contribution of M2R voltage sensitivity to the rate and shape of the human sinus node action potential in physiological and pathophysiological conditions, using a Markovian model of the I_K,ACh channel integrated into a computational model of human sinus node. The computational model allowed us to assess the effects of potential genetic variants that alter specific properties of voltage sensitivity. Our results indicate that alterations in the M2R voltage sensitivity play a significant role in the physiology and pathophysiology of the human sinus node and atria. Our computational model is relevant for future studies aimed at the design and development of anti-arrhythmic agents that specifically target the unique voltage-sensitive properties of M2R.

A thermodynamic framework for modelling membrane transporters

Membrane transporters contribute to the regulation of the internal environment of cells by translocating substrates across cell membranes. Like all physical systems, the behaviour of membrane transporters is constrained by the laws of thermodynamics. However, many mathematical models of transporters, especially those incorporated into whole-cell models, are not thermodynamically consistent, leading to unrealistic behaviour. In this paper we use a physics-based modelling framework, in which the transfer of energy is explicitly accounted for, to develop thermodynamically consistent models of transporters. We then apply this methodology to model two specific transporters: the cardiac sarcoplasmic/endoplasmic Ca²⁺ ATPase (SERCA) and the cardiac Na⁺/K⁺ ATPase.

Fitting two human atrial cell models to experimental data using Bayesian history matching

Cardiac cell models are potentially valuable tools for applications such as quantitative safety pharmacology, but have many parameters. Action potentials in real cardiac cells also vary from beat to beat, and from one cell to another. Calibrating cardiac cell models to experimental observations is difficult, because the parameter space is large and high-dimensional. In this study we have demonstrated the use of history matching to calibrate the maximum conductance of ion channels and exchangers in two detailed models of the human atrial action potential against measurements of action potential biomarkers. History matching is an approach developed in other modelling communities, based on constructing fast-running Gaussian process emulators of the model. Emulators were constructed from a small number of model runs (around 10²), and then run many times (>10⁶) at low computational cost, each time with a different set of model parameters. Emulator outputs were compared with experimental biomarkers using an implausibility measure, which took into account experimental variance as well as emulator variance. By repeating this process, the region of non-implausible parameter space was iteratively reduced. Both cardiac cell models were successfully calibrated to experimental datasets, resulting in sets of parameters that could be sampled to produce variable action potentials. However, model parameters did not occupy a small range of values. Instead, the history matching process exposed inputs that can co-vary across a wide range and still be consistent with a particular biomarker. We also found correlations between some biomarkers, indicating a need for better descriptors of action potential shape.

A Heart for Diversity: Simulating Variability in Cardiac Arrhythmia Research

In cardiac electrophysiology, there exist many sources of inter- and intra-personal variability. These include variability in conditions and environment, and genotypic and molecular diversity, including differences in expression and behavior of ion channels and transporters, which lead to phenotypic diversity (e.g., variable integrated responses at the cell, tissue, and organ levels). These variabilities play an important role in progression of heart disease and arrhythmia syndromes and outcomes of therapeutic interventions. Yet, the traditional in silico framework for investigating cardiac arrhythmias is built upon a parameter/property-averaging approach that typically overlooks the physiological diversity. Inspired by work done in genetics and neuroscience, new modeling frameworks of cardiac electrophysiology have been recently developed that take advantage of modern computational capabilities and approaches, and account for the variance in the biological data they are intended to illuminate. In this review, we outline the recent advances in statistical and computational techniques that take into account physiological variability, and move beyond the traditional cardiac model-building scheme that involves averaging over samples from many individuals in the construction of a highly tuned composite model. We discuss how these advanced methods have harnessed the power of big (simulated) data to study the mechanisms of cardiac arrhythmias, with a special emphasis on atrial fibrillation, and improve the assessment of proarrhythmic risk and drug response. The challenges of using in silico approaches with variability are also addressed and future directions are proposed.

Reproducible model development in the cardiac electrophysiology Web Lab

The modelling of the electrophysiology of cardiac cells is one of the most mature areas of systems biology. This extended concentration of research effort brings with it new challenges, foremost among which is that of choosing which of these models is most suitable for addressing a particular scientific question. In a previous paper, we presented our initial work in developing an online resource for the characterisation and comparison of electrophysiological cell models in a wide range of experimental scenarios. In that work, we described how we had developed a novel protocol language that allowed us to separate the details of the mathematical model (the majority of cardiac cell models take the form of ordinary differential equations) from the experimental protocol being simulated. We developed a fully-open online repository (which we termed the Cardiac Electrophysiology Web Lab) which allows users to store and compare the results of applying the same experimental protocol to competing models. In the current paper we describe the most recent and planned extensions of this work, focused on supporting the process of model building from experimental data. We outline the necessary work to develop a machine-readable language to describe the process of inferring parameters from wet lab datasets, and illustrate our approach through a detailed example of fitting a model of the hERG channel using experimental data. We conclude by discussing the future challenges in making further progress in this domain towards our goal of facilitating a fully reproducible approach to the development of cardiac cell models.

Simulation of Spiral Waves and Point Sources in Atrial Fibrillation with Application to Rotor Localization

Cardiac simulations play an important role in studies involving understanding and investigating the mechanisms of cardiac arrhythmias. Today, studies of arrhythmogenesis and maintenance are largely being performed by creating simulations of a particular arrhythmia with high accuracy comparable to the results of clinical experiments. Atrial fibrillation (AF), the most common arrhythmia in the United States and many other parts of the world, is one of the major field where simulation and modeling is largely used. AF simulations not only assist in understanding its mechanisms but also help to develop, evaluate and improve the computer algorithms used in electrophysiology (EP) systems for ablation therapies. In this paper, we begin with a brief overeview of some common techniques used in simulations to simulate two major AF mechanisms - spiral waves (or rotors) and point (or focal) sources. We particularly focus on 2D simulations using Nygren et al.'s mathematical model of human atrial cell. Then, we elucidate an application of the developed AF simulation to an algorithm designed for localizing AF rotors for improving current AF ablation therapies. Our simulation methods and results, along with the other discussions presented in this paper is aimed to provide engineers and professionals with a working-knowledge of application-specific simulations of spirals and foci.

Confocal Microscopy-Based Estimation of Parameters for Computational Modeling of Electrical Conduction in the Normal and Infarcted Heart

Computational modeling is an important tool to advance our knowledge on cardiac diseases and their underlying mechanisms. Computational models of conduction in cardiac tissues require identification of parameters. Our knowledge on these parameters is limited, especially for diseased tissues. Here, we assessed and quantified parameters for computational modeling of conduction in cardiac tissues. We used a rabbit model of myocardial infarction (MI) and an imaging-based approach to derive the parameters. Left ventricular tissue samples were obtained from fixed control hearts (animals: 5) and infarcted hearts (animals: 6) within 200 μm (region 1), 250-750 μm (region 2) and 1,000-1,250 μm (region 3) of the MI border. We assessed extracellular space, fibroblasts, smooth muscle cells, nuclei and gap junctions by a multi-label staining protocol. With confocal microscopy we acquired three-dimensional (3D) image stacks with a voxel size of 200 × 200 × 200 nm. Image segmentation yielded 3D reconstructions of tissue microstructure, which were used to numerically derive extracellular conductivity tensors. Volume fractions of myocyte, extracellular, interlaminar cleft, vessel and fibroblast domains in control were (in %) 65.03 ± 3.60, 24.68 ± 3.05, 3.95 ± 4.84, 7.71 ± 2.15, and 2.48 ± 1.11, respectively. Volume fractions in regions 1 and 2 were different for myocyte, myofibroblast, vessel, and extracellular domains. Fibrosis, defined as increase in fibrotic tissue constituents, was (in %) 21.21 ± 1.73, 16.90 ± 9.86, and 3.58 ± 8.64 in MI regions 1, 2, and 3, respectively. For control tissues, image-based computation of longitudinal, transverse and normal extracellular conductivity yielded (in S/m) 0.36 ± 0.11, 0.17 ± 0.07, and 0.1 ± 0.06, respectively. Conductivities were markedly increased in regions 1 (+75, +171, and +100%), 2 (+53, +165, and +80%), and 3 (+42, +141, and +60%). Volume fractions of the extracellular space including interlaminar clefts strongly correlated with conductivities in control and MI hearts. Our study provides novel quantitative data for computational modeling of conduction in normal and MI hearts. Notably, our study introduces comprehensive statistical information on tissue composition and extracellular conductivities on a microscopic scale in the MI border zone. We suggest that the presented data fill a significant gap in modeling parameters and extend our foundation for computational modeling of cardiac conduction.

Cardiac ischemia-insights from computational models

BACKGROUND

Complementary to clinical and experimental studies, computational cardiac modeling serves to obtain a comprehensive understanding of the cardiovascular system in order to analyze dysfunction, evaluate existing, and develop novel treatment strategies.

OBJECTIVES

We describe the basics of multiscale computational modeling of cardiac electrophysiology from the molecular ion channel to the whole body scale. By modeling cardiac ischemia, we illustrate how in silico experiments can contribute to our understanding of how the pathophysiological mechanisms translate into changes observed in diagnostic tools such as the electrocardiogram (ECG).

MATERIALS AND METHODS

Quantitative in silico modeling spans a wide range of scales from ion channel biophysics to ECG signals. For each of the scales, a set of mathematical equations describes electrophysiology in relation to the other scales. Integration of ischemia-induced changes is performed on the ion channel, single-cell, and tissue level. This approach allows us to study how effects simulated at molecular scales translate to changes in the ECG.

RESULTS

Ischemia induces action potential shortening and conduction slowing. Hence, ischemic myocardium has distinct and significant effects on propagation and repolarization of excitation, depending on the intramural extent of the ischemic region. For transmural and subendocardial ischemic regions, ST segment elevation and depression, respectively, were observed, whereas intermediate ischemic regions were found to be electrically silent (NSTEMI).

CONCLUSIONS

In silico modeling contributes quantitative and mechanistic insight into fundamental ischemia-related arrhythmogenic mechanisms. In addition, computational modeling can help to translate experimental findings at the (sub-)cellular level to the organ and body context (e. g., ECG), thereby providing a thorough understanding of this routinely used diagnostic tool that may translate into optimized applications.

Next-generation, personalised, model-based critical care medicine: a state-of-the art review of in silico virtual patient models, methods, and cohorts, and how to validation them

Critical care, like many healthcare areas, is under a dual assault from significantly increasing demographic and economic pressures. Intensive care unit (ICU) patients are highly variable in response to treatment, and increasingly aging populations mean ICUs are under increasing demand and their cohorts are increasingly ill. Equally, patient expectations are growing, while the economic ability to deliver care to all is declining. Better, more productive care is thus the big challenge. One means to that end is personalised care designed to manage the significant inter- and intra-patient variability that makes the ICU patient difficult. Thus, moving from current "one size fits all" protocolised care to adaptive, model-based "one method fits all" personalised care could deliver the required step change in the quality, and simultaneously the productivity and cost, of care. Computer models of human physiology are a unique tool to personalise care, as they can couple clinical data with mathematical methods to create subject-specific models and virtual patients to design new, personalised and more optimal protocols, as well as to guide care in real-time. They rely on identifying time varying patient-specific parameters in the model that capture inter- and intra-patient variability, the difference between patients and the evolution of patient condition. Properly validated, virtual patients represent the real patients, and can be used in silico to test different protocols or interventions, or in real-time to guide care. Hence, the underlying models and methods create the foundation for next generation care, as well as a tool for safely and rapidly developing personalised treatment protocols over large virtual cohorts using virtual trials. This review examines the models and methods used to create virtual patients. Specifically, it presents the models types and structures used and the data required. It then covers how to validate the resulting virtual patients and trials, and how these virtual trials can help design and optimise clinical trial. Links between these models and higher order, more complex physiome models are also discussed. In each section, it explores the progress reported up to date, especially on core ICU therapies in glycemic, circulatory and mechanical ventilation management, where high cost and frequency of occurrence provide a significant opportunity for model-based methods to have measurable clinical and economic impact. The outcomes are readily generalised to other areas of medical care.

Validation and Trustworthiness of Multiscale Models of Cardiac Electrophysiology

Computational models of cardiac electrophysiology have a long history in basic science applications and device design and evaluation, but have significant potential for clinical applications in all areas of cardiovascular medicine, including functional imaging and mapping, drug safety evaluation, disease diagnosis, patient selection, and therapy optimisation or personalisation. For all stakeholders to be confident in model-based clinical decisions, cardiac electrophysiological (CEP) models must be demonstrated to be trustworthy and reliable. Credibility, that is, the belief in the predictive capability, of a computational model is primarily established by performing validation, in which model predictions are compared to experimental or clinical data. However, there are numerous challenges to performing validation for highly complex multi-scale physiological models such as CEP models. As a result, credibility of CEP model predictions is usually founded upon a wide range of distinct factors, including various types of validation results, underlying theory, evidence supporting model assumptions, evidence from model calibration, all at a variety of scales from ion channel to cell to organ. Consequently, it is often unclear, or a matter for debate, the extent to which a CEP model can be trusted for a given application. The aim of this article is to clarify potential rationale for the trustworthiness of CEP models by reviewing evidence that has been (or could be) presented to support their credibility. We specifically address the complexity and multi-scale nature of CEP models which makes traditional model evaluation difficult. In addition, we make explicit some of the credibility justification that we believe is implicitly embedded in the CEP modeling literature. Overall, we provide a fresh perspective to CEP model credibility, and build a depiction and categorisation of the wide-ranging body of credibility evidence for CEP models. This paper also represents a step toward the extension of model evaluation methodologies that are currently being developed by the medical device community, to physiological models.

Mechanism of doxorubicin cardiotoxicity evaluated by integrating multiple molecular effects into a biophysical model

BACKGROUND AND PURPOSE

Doxorubicin (DOX) is an effective cancer therapeutic agent but causes therapy-limiting cardiotoxicity. The effects of DOX and its metabolite doxorubicinol (DOXL) on individual channels have been well characterized in isolation. However, it is unknown how the action and interaction of affected channels combine to generate the phenotypic cardiotoxic outcome. We sought to develop an in silico model that links drug effects on channels to action potential duration (APD) and intracellular Ca²⁺ concentration in order to address this gap in knowledge.

EXPERIMENTAL APPROACH

We first propose two methods to obtain, from published values, consensus drug effects on the currents of individual channels, transporters and pumps. Separately, we obtained equivalent values for APD and Ca²⁺ concentration (the readouts used as surrogates for cardiotoxicity). Once derived, the consensus effects on the currents were incorporated into established biophysical models of the cardiac myocyte and were refined adjusting the sarcoplasmic reticulum Ca²⁺ leak current (I_Leak ) until the consensus effects on APD and Ca²⁺ dynamics were replicated. Using factorial analysis, we then quantified the relative contribution of each channel to DOX and DOXL cardiotoxicity.

KEY RESULTS

The factorial analysis identified the rapid delayed rectifying K⁺ current, the L-type Ca²⁺ current and the sarcoplasmic reticulum I_Leak as the targets primarily responsible for the cardiotoxic effects on APD and Ca²⁺ dynamics.

CONCLUSIONS AND IMPLICATIONS

This study provides insight into the mechanisms of DOX-induced cardiotoxicity and a framework for the development of future diagnostic and therapeutic strategies.