A tale of two dogs: analyzing two models of canine ventricular electrophysiology

Uncovering the dynamics of cardiac systems using stochastic pacing and frequency domain analyses

Alternans of cardiac action potential duration (APD) is a well-known arrhythmogenic mechanism which results from dynamical instabilities. The propensity to alternans is classically investigated by examining APD restitution and by deriving APD restitution slopes as predictive markers. However, experiments have shown that such markers are not always accurate for the prediction of alternans. Using a mathematical ventricular cell model known to exhibit unstable dynamics of both membrane potential and Ca²⁺ cycling, we demonstrate that an accurate marker can be obtained by pacing at cycle lengths (CLs) varying randomly around a basic CL (BCL) and by evaluating the transfer function between the time series of CLs and APDs using an autoregressive-moving-average (ARMA) model. The first pole of this transfer function corresponds to the eigenvalue (λ(alt)) of the dominant eigenmode of the cardiac system, which predicts that alternans occurs when λ(alt) ≤ -1. For different BCLs, control values of λ(alt) were obtained using eigenmode analysis and compared to the first pole of the transfer function estimated using ARMA model fitting in simulations of random pacing protocols. In all versions of the cell model, this pole provided an accurate estimation of λ(alt). Furthermore, during slow ramp decreases of BCL or simulated drug application, this approach predicted the onset of alternans by extrapolating the time course of the estimated λ(alt). In conclusion, stochastic pacing and ARMA model identification represents a novel approach to predict alternans without making any assumptions about its ionic mechanisms. It should therefore be applicable experimentally for any type of myocardial cell.

Parameter fitting using multiple datasets in cardiac action potential modeling

A multiple dataset model fitting approach for improving parameter reliability in action potential modeling is presented. A robust generic cardiac ionic model employing membrane currents based on two-gate Hodgkin-Huxley kinetics is described. Its generic nature allows it to accurately reproduce action potential waveforms in heterogeneous cardiac tissue by optimizing parameters governing ion channel kinetics and magnitudes. The model allows a user-defined number of voltage and time-dependent ion currents to be incorporated, in order to reproduce and predict multiple action potential waveforms recorded in intact cardiac myocyte. In total 12 N(c)+2 parameters were optimized using a curvilinear gradient method, where N(c) is the user-specified number of time-dependent currents. Given appropriate experimental datasets, many of the known physiological membrane currents could be effectively reconstructed. Also, the optimized models were able to predict additional experimental action potential recordings that were not used in the optimization process.

Verification of cardiac tissue electrophysiology simulators using an N-version benchmark

Ongoing developments in cardiac modelling have resulted, in particular, in the development of advanced and increasingly complex computational frameworks for simulating cardiac tissue electrophysiology. The goal of these simulations is often to represent the detailed physiology and pathologies of the heart using codes that exploit the computational potential of high-performance computing architectures. These developments have rapidly progressed the simulation capacity of cardiac virtual physiological human style models; however, they have also made it increasingly challenging to verify that a given code provides a faithful representation of the purported governing equations and corresponding solution techniques. This study provides the first cardiac tissue electrophysiology simulation benchmark to allow these codes to be verified. The benchmark was successfully evaluated on 11 simulation platforms to generate a consensus gold-standard converged solution. The benchmark definition in combination with the gold-standard solution can now be used to verify new simulation codes and numerical methods in the future.

Multi-scale computational models of familial hypertrophic cardiomyopathy: genotype to phenotype

Familial hypertrophic cardiomyopathy (FHC) is an inherited disorder affecting roughly one in 500 people. Its hallmark is abnormal thickening of the ventricular wall, leading to serious complications that include heart failure and sudden cardiac death. Treatment is complicated by variation in the severity, symptoms and risks for sudden death within the patient population. Nearly all of the genetic lesions associated with FHC occur in genes encoding sarcomeric proteins, indicating that defects in cardiac muscle contraction underlie the condition. Detailed biophysical data are increasingly available for computational analyses that could be used to predict heart phenotypes based on genotype. These models must integrate the dynamic processes occurring in cardiac cells with properties of myocardial tissue, heart geometry and haemodynamic load in order to predict strain and stress in the ventricular walls and overall pump function. Recent advances have increased the biophysical detail in these models at the myofilament level, which will allow properties of FHC-linked mutant proteins to be accurately represented in simulations of whole heart function. The short-term impact of these models will be detailed descriptions of contractile dysfunction and altered myocardial strain patterns at the earliest stages of the disease-predictions that could be validated in genetically modified animals. Long term, these multi-scale models have the potential to improve clinical management of FHC through genotype-based risk stratification and personalized therapy.

3D current-voltage-time surfaces unveil critical repolarization differences underlying similar cardiac action potentials: A model study

The number of mathematical models of cardiac cellular excitability is rapidly growing, and compact graphical representations of their properties can make new acquisitions available for a broader range of scientists in cardiac field. Particularly, the intrinsic over-determination of the model equations systems when fitted only to action potential (AP) waveform and the fact that they are frequently tuned on data covering only a relatively narrow range of dynamic conditions, often lead modellers to compare very similar AP profiles, which underlie though quite different excitable properties. In this study I discuss a novel compact 3D representation of the cardiac cellular AP, where the third dimension represents the instantaneous current-voltage profile of the membrane, measured as repolarization proceeds. Measurements of this type have been used previously for in vivo experiments, and are adopted here iteratively at a very high time, voltage, current-resolution on (i) the same human ventricular model, endowed with two different parameters sets which generate the same AP waveform, and on (ii) three different models of the same human ventricular cell type. In these 3D representations, the AP waveforms lie at the intersection between instantaneous time-voltage-current surfaces and the zero-current plane. Different surfaces can share the same intersection and therefore the same AP; in these cases, the morphology of the current surface provides a compact view of important differences within corresponding repolarization dynamics. Refractory period, supernormal excitability window, and extent of repolarization reserve can be visualized at once. Two pivotal dynamical properties can be precisely assessed, i.e. all-or-nothing repolarization window and membrane resistance during recovery. I discuss differences in these properties among the membranes under study, and show relevant implications for cardiac cellular repolarization.

Systematic characterization of the ionic basis of rabbit cellular electrophysiology using two ventricular models

Several mathematical models of rabbit ventricular action potential (AP) have been proposed to investigate mechanisms of arrhythmias and excitation-contraction coupling. Our study aims at systematically characterizing how ionic current properties modulate the main cellular biomarkers of arrhythmic risk using two widely-used rabbit ventricular models, and comparing simulation results using the two models with experimental data available for rabbit. A sensitivity analysis of AP properties, Ca²⁺ and Na⁺ dynamics, and their rate dependence to variations (±15% and ±30%) in the main transmembrane current conductances and kinetics was performed using the Shannon et al. (2004) and the Mahajan et al. (2008a,b) AP rabbit models. The effects of severe transmembrane current blocks (up to 100%) on steady-state AP and calcium transients, and AP duration (APD) restitution curves were also simulated using both models. Our simulations show that, in both virtual rabbit cardiomyocytes, APD is significantly modified by most repolarization currents, AP triangulation is regulated mostly by the inward rectifier K⁺ current (I(K1)) whereas APD rate adaptation as well as [Na⁺](i) rate dependence is influenced by the Na⁺/K⁺ pump current (I(NaK)). In addition, steady-state [Ca²⁺](i) levels, APD restitution properties and [Ca²⁺](i) rate dependence are strongly dependent on I(NaK), the L-Type Ca²⁺ current (I(CaL)) and the Na⁺/Ca²⁺ exchanger current (I(NaCa)), although the relative role of these currents is markedly model dependent. Furthermore, our results show that simulations using both models agree with many experimentally-reported electrophysiological characteristics. However, our study shows that the Shannon et al. model mimics rabbit electrophysiology more accurately at normal pacing rates, whereas Mahajan et al. model behaves more appropriately at faster rates. Our results reinforce the usefulness of sensitivity analysis for further understanding of cellular electrophysiology and validation of cardiac AP models.

High-throughput functional curation of cellular electrophysiology models

Effective reuse of a quantitative mathematical model requires not just access to curated versions of the model equations, but also an understanding of the functional capabilities of the model, and the advisable scope of its application. To enable this "functional curation" we have developed a simulation environment that provides high-throughput evaluation of a mathematical model's functional response to an arbitrary user-defined protocol, and optionally compares the results against experimental data. In this study we demonstrate the efficacy of this simulation environment on 31 cardiac electrophysiology cell models using two test cases. The S1-S2 response is evaluated to characterise the models' restitution curves, and their L-type calcium channel current-voltage curves are evaluated. The significant variation in the response of these models, even when the models represent the same species and temperature, demonstrates the importance of knowing the functional characteristics of a model prior to its reuse. We also discuss the wider implications for this approach, in improving the selection of models for reuse, enabling the identification of models that exhibit particular experimentally observed phenomena, and making the incremental development of models more robust.

Quantification of repolarization reserve to understand interpatient variability in the response to proarrhythmic drugs: a computational analysis

BACKGROUND

"Repolarization reserve" is frequently invoked to explain why potentially proarrhythmic drugs cause, across a population, a range of changes to cardiac action potentials (APs). However, the mechanisms underlying this interindividual variability are not understood quantitatively.

OBJECTIVE

The purpose of this study was to perform a novel analysis of mathematical models of ventricular myocytes to quantify repolarization reserve and gain insight into the factors responsible for variability in the response to proarrhythmic drugs.

METHODS/RESULTS

In several models of human or canine ventricular myocytes, variability was simulated by randomizing model parameters and running repeated simulations. With each randomly generated set of parameters, APs before and after simulated 75% block of the rapid delayed rectifier current (I(Kr)) were calculated. Multivariable regression was performed to determine how much each model parameter attenuated or exacerbated the AP prolongation caused by the I(Kr)-blocking drug. Simulations with a human ventricular myocyte model suggest that drug response is influenced most strongly by (1) the density of I(Kr), (2) the density of slow delayed rectifier current I(Ks), (3) the voltage dependence of I(Kr) inactivation, (4) the density of L-type Ca2+ current, and (5) the kinetics of I(Ks) activation. The analysis also identified mechanisms underlying nonintuitive behavior, such as ionic currents that prolong baseline APs but decrease drug-induced AP prolongation. Finally, the simulations provided quantitative insight into conditions that aggravate the drug response, such as silent ion channel mutations and heart failure.

CONCLUSION

These modeling results provide the first thorough quantification of repolarization reserve and improve our understanding of interindividual variability in adverse drug reactions.

Reentrant excitation in an analog-digital hybrid circuit model of cardiac tissue

We propose an analog-digital hybrid circuit model of one-dimensional cardiac tissue with hardware implementation that allows us to perform real-time simulations of spatially conducting cardiac action potentials. Each active nodal compartment of the tissue model is designed using analog circuits and a dsPIC microcontroller, by which the time-dependent and time-independent nonlinear current-voltage relationships of six types of ion channel currents employed in the Luo-Rudy phase I (LR-I) model for a single mammalian cardiac ventricular cell can be reproduced quantitatively. Here, we perform real-time simulations of reentrant excitation conduction in a ring-shaped tissue model that includes eighty nodal compartments. In particular, we show that the hybrid tissue model can exhibit real-time dynamics for initiation of reentries induced by uni-directional block, as well as those for phase resetting that leads to annihilation of the reentry in response to impulsive current stimulations at appropriate nodes and timings. The dynamics of the hybrid model are comparable to those of a spatially distributed tissue model with LR-I compartments. Thus, it is conceivable that the hybrid model might be a useful tool for large scale simulations of cardiac tissue dynamics, as an alternative to numerical simulations, leading toward further understanding of the reentrant mechanisms.

How different two almost identical action potentials can be: a model study on cardiac repolarization

Spatial heterogeneity in the properties of ion channels generates spatial dispersion of ventricular repolarization, which is modulated by gap junctional coupling. However, it is possible to simulate conditions in which local differences in excitation properties are electrophysiologically silent and only play a role in pathological states. We use a numerical procedure on the Luo-Rudy phase 1 model of the ventricular action potential (AP1) in order to find a modified set of model parameters which generates an action potential profile (AP2) almost identical to AP1. We show that, although the two waveforms elicited from resting conditions as a single AP are very similar and belong to membranes sharing similar passive electrical properties, the modified membrane generating AP2 is a weaker current source than the one generating AP1, has different sensitivity to up/down-regulation of ion channels and to extracellular potassium, and a different electrical restitution profile. We study electrotonic interaction of AP1- and AP2- type membranes in cell pairs and in cable conduction, and find differences in source-sink properties which are masked in physiological conditions and become manifest during intercellular uncoupling or partial block of ion channels, leading to unidirectional block and spatial repolarization gradients. We provide contour plot representations that summarize differences and similarities. The present report characterizes an inverse problem in cardiac cells, and strengthen the recently emergent notion that a comprehensive characterization and validation of cell models and their components are necessary in order to correctly understand simulation results at higher levels of complexity.

Models of cardiac tissue electrophysiology: progress, challenges and open questions

Models of cardiac tissue electrophysiology are an important component of the Cardiac Physiome Project, which is an international effort to build biophysically based multi-scale mathematical models of the heart. Models of tissue electrophysiology can provide a bridge between electrophysiological cell models at smaller scales, and tissue mechanics, metabolism and blood flow at larger scales. This paper is a critical review of cardiac tissue electrophysiology models, focussing on the micro-structure of cardiac tissue, generic behaviours of action potential propagation, different models of cardiac tissue electrophysiology, the choice of parameter values and tissue geometry, emergent properties in tissue models, numerical techniques and computational issues. We propose a tentative list of information that could be included in published descriptions of tissue electrophysiology models, and used to support interpretation and evaluation of simulation results. We conclude with a discussion of challenges and open questions.

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

In this manuscript we review the state of cardiac cell modelling in the context of international initiatives such as the IUPS Physiome and Virtual Physiological Human Projects, which aim to integrate computational models across scales and physics. In particular we focus on the relationship between experimental data and model parameterisation across a range of model types and cellular physiological systems. Finally, in the context of parameter identification and model reuse within the Cardiac Physiome, we suggest some future priority areas for this field.

Multistability property in cardiac ionic models of mammalian and human ventricular cells

The underlying mechanisms of irregular cardiac rhythms are still poorly understood. Many experimental and modeling studies are aimed at identifying factors which cause cardiac arrhythmias. However, a lack of understanding of heart rhythm dynamical properties makes it difficult to uncover precise mechanisms of electrical instabilities, and hence to predict the onset of heart rhythm disorders. We review and compare the existing methods of studying cardiac dynamics, including restitution protocol (S1-S2), dynamic restitution protocol and multistability test protocol (S1-CI-S2). We focus on cardiac cell dynamics to elucidate regularities of heart rhythm. We demonstrate the advantages of our newly proposed systematic approach of analysis of cardiac cell dynamics using mammalian Luo Rudy 1991 and human ventricular Ten Tusscher 2006 single cell models under healthy and diseased conditions such as altered K(+) or Ca(2+) related currents. We investigate the role of ionic properties and the shape of an action potential on the nonlinear dynamics of electrical processes in periodically stimulated cardiac cells. We show the existence of multistability property for human ventricular cells. Moreover, the multistability is proposed to be an intrinsic property of cardiac cells, and is also suggested to be one of the mechanisms which could underlie the sudden triggering of life-threatening ventricular arrhythmias in the human heart.

Parameter sensitivity analysis in electrophysiological models using multivariable regression

Computational models of electrical activity and calcium signaling in cardiac myocytes are important tools for understanding physiology. The sensitivity of these models to changes in parameters is often not well-understood, however, because parameter evaluation can be a time-consuming, tedious process. I demonstrate here what I believe is a novel method for rapidly determining how changes in parameters affect outputs. In three models of the ventricular action potential, parameters were randomized, repeated simulations were run, important outputs were calculated, and multivariable regression was performed on the collected results. Random parameters included both maximal rates of ion transport and gating variable characteristics. The procedure generated simplified, empirical models that predicted outputs resulting from new sets of input parameters. The linear regression models were quite accurate, despite nonlinearities in the mechanistic models. Moreover, the regression coefficients, which represent parameter sensitivities, were robust, even when parameters were varied over a wide range. Most importantly, a side-by-side comparison of two similar models identified fundamental differences in model behavior, and revealed model predictions that were both consistent with, and inconsistent with, experimental data. This new method therefore shows promise as a tool for the characterization and assessment of computational models. The general strategy may also suggest methods for integrating traditional quantitative models with large-scale data sets obtained using high-throughput technologies.

Properties and ionic mechanisms of action potential adaptation, restitution, and accommodation in canine epicardium

Computational models of cardiac myocytes are important tools for understanding ionic mechanisms of arrhythmia. This work presents a new model of the canine epicardial myocyte that reproduces a wide range of experimentally observed rate-dependent behaviors in cardiac cell and tissue, including action potential (AP) duration (APD) adaptation, restitution, and accommodation. Model behavior depends on updated formulations for the 4-aminopyridine-sensitive transient outward current (I(to1)), the slow component of the delayed rectifier K(+) current (I(Ks)), the L-type Ca(2+) channel current (I(Ca,L)), and the Na(+)-K(+) pump current (I(NaK)) fit to data from canine ventricular myocytes. We found that I(to1) plays a limited role in potentiating peak I(Ca,L) and sarcoplasmic reticulum Ca(2+) release for propagated APs but modulates the time course of APD restitution. I(Ks) plays an important role in APD shortening at short diastolic intervals, despite a limited role in AP repolarization at longer cycle lengths. In addition, we found that I(Ca,L) plays a critical role in APD accommodation and rate dependence of APD restitution. Ca(2+) entry via I(Ca,L) at fast rate drives increased Na(+)-Ca(2+) exchanger Ca(2+) extrusion and Na(+) entry, which in turn increases Na(+) extrusion via outward I(NaK). APD accommodation results from this increased outward I(NaK). Our simulation results provide valuable insight into the mechanistic basis of rate-dependent phenomena important for determining the heart's response to rapid and irregular pacing rates (e.g., arrhythmia). Accurate simulation of rate-dependent phenomena and increased understanding of their mechanistic basis will lead to more realistic multicellular simulations of arrhythmia and identification of molecular therapeutic targets.

A meta-analysis of cardiac electrophysiology computational models

Computational models of cardiac electrophysiology are exemplar demonstrations of the integration of multiple data sets into a consistent biophysical framework. These models encapsulate physiological understanding to provide quantitative predictions of function. The combination or extension of existing models within a common framework allows integrative phenomena in larger systems to be investigated. This methodology is now routinely applied, as demonstrated by the increasing number of studies which use or extend previously developed models. In this study, we present a meta-analysis of this model re-use for two leading models of cardiac electrophysiology in the form of parameter inheritance trees, a sensitivity analysis and a comparison of the functional significance of the sodium potassium pump for defining restitution curves. These results indicate that even though the models aim to represent the same physiological system, both the sources of parameter values and the function of equivalent components are significantly different.

Termination of equine atrial fibrillation by quinidine: an optical mapping study

OBJECTIVE

To perform the first optical mapping studies of equine atrium to assess the spatiotemporal dynamics of atrial fibrillation (AF) and of its termination by quinidine.

ANIMALS

Intact, perfused atrial preparations obtained from four horses with normal cardiovascular examinations.

MATERIALS AND METHODS

AF was induced by a rapid pacing protocol with or without acetylcholine perfusion, and optical mapping was used to determine spatial dominant frequency distributions, electrical activity maps, and single-pixel optical signals. Following induction of AF, quinidine gluconate was perfused into the preparation and these parameters were monitored during quinidine-induced termination of AF.

RESULTS

Equine AF develops in the context of spatial gradients in action potential duration (APD) and diastolic interval (DI) that produce alternans, conduction block, and Wenckebach conduction in different regions at fast pacing rates. Quinidine terminates AF and prevents subsequent reinduction by reducing the maximal frequency and increasing frequency homogeneity.

CONCLUSIONS

Heterogeneity of APD and DI promote alternans and conduction block at fast pacing rates in the equine atrium, predisposing to the development of AF. Quinidine terminates AF by reducing maximum frequency and increasing frequency homogeneity. Our results are consistent with the hypothesis that quinidine increases effective refractory period, thereby decreasing frequency.

Properties of two human atrial cell models in tissue: restitution, memory, propagation, and reentry

To date, two detailed ionic models of human atrial cell electrophysiology have been developed, the Nygren et al. model (NM) and the Courtemanche et al. model (CM). Although both models draw from similar experimental data, they have vastly different properties. This paper provides the first systematic analysis and comparison of the dynamics of these models in spatially extended systems including one-dimensional cables and rings, two-dimensional sheets, and a realistic three-dimensional human atrial geometry. We observe that, as in single cells, the CM adapts to rate changes primarily by changes in action potential duration (APD) and morphology, while for the NM rate changes affect resting membrane potential (RMP) more than APD. The models also exhibit different memory properties as assessed through S1-S2 APD and conduction velocity (CV) restitution curves with different S1 cycle lengths. Reentrant wave dynamics also differ, with the NM exhibiting stable, non-breaking spirals and the CM exhibiting frequent transient wave breaks. The realistic atrial geometry modifies dynamics in some cases through drift, transient pinning, and breakup. Previously proposed modifications to represent atrial fibrillation-remodeled electrophysiology produce altered dynamics, including reduced rate adaptation and memory for both models and conversion to stable reentry for the CM. Furthermore, proposed variations to the NM to reproduce action potentials more closely resembling those of the CM do not substantially alter the underlying dynamics of the model, so that tissue simulations using these modifications still behave more like the unmodified NM. Finally, interchanging the transmembrane current formulations of the two models suggests that currents contribute more strongly to RMP and CV, intracellular calcium dynamics primarily determine reentrant wave dynamics, and both are important in APD restitution and memory in these models. This finding implies that the formulation of intracellular calcium processes is as important to producing realistic models as transmembrane currents.

Indeterminacy of spatiotemporal cardiac alternans

Cardiac alternans, a beat-to-beat alternation in action potential duration (at the cellular level) or in electrocardiogram morphology (at the whole heart level), is a marker of ventricular fibrillation, a fatal heart rhythm that kills hundreds of thousands of people in the United States each year. Investigating cardiac alternans may lead to a better understanding of the mechanisms of cardiac arrhythmias and eventually better algorithms for the prediction and prevention of such dreadful diseases. In paced cardiac tissue, alternans develops under increasingly shorter pacing period. Existing experimental and theoretical studies adopt the assumption that alternans in homogeneous cardiac tissue is exclusively determined by the pacing period. In contrast, we find that, when calcium-driven alternans develops in cardiac fibers, it may take different spatiotemporal patterns depending on the pacing history. Because there coexist multiple alternans solutions for a given pacing period, the alternans pattern on a fiber becomes unpredictable. Using numerical simulation and theoretical analysis, we show that the coexistence of multiple alternans patterns is induced by the interaction between electrotonic coupling and an instability in calcium cycling.

Dynamics of human atrial cell models: restitution, memory, and intracellular calcium dynamics in single cells

Mathematical models of cardiac cells have become important tools for investigating the electrophysiological properties and behavior of the heart. As the number of published models increases, it becomes more difficult to choose a model appropriate for the conditions to be studied, especially when multiple models describing the species and region of the heart of interest are available. In this paper, we will review and compare two detailed ionic models of human atrial myocytes, the Nygren et al. model (NM) and the Courtemanche et al. model (CM). Although both models include the same transmembrane currents and are largely based on the same experimental data from human atrial cells, the two models exhibit vastly different properties, especially in their dynamical behavior, including restitution and memory effects. The CM produces pronounced rate adaptation of action potential duration (APD) with limited memory effects, while the NM exhibits strong rate dependence of resting membrane potential (RMP), limited APD restitution, and stronger memory, as well as delayed afterdepolarizations and auto-oscillatory behavior upon cessation of rapid pacing. Channel conductance modifications based on experimentally measured changes during atrial fibrillation modify rate adaptation and memory in both models, but do not change the primary rate-dependent properties of APD and RMP for the CM and NM, respectively. Two sets of proposed changes to the NM that yield a spike-and-dome action potential morphology qualitatively similar to the CM at slow pacing rates similarly do not change the underlying dynamics of the model. Moreover, interchanging the formulations of all transmembrane currents between the two models while leaving calcium handling and ionic concentrations intact indicates that the currents strongly influence memory and the rate adaptation of RMP, while intracellular calcium dynamics primarily determine APD rate adaptation. Our results suggest that differences in intracellular calcium handling between the two human atrial myocyte models are responsible for marked dynamical differences and may prevent reconciliation between the models by straightforward channel conductance modifications.

Using Physiome standards to couple cellular functions for rat cardiac excitation-contraction

Scientific endeavour is reliant upon the extension and reuse of previous knowledge. The formalization of this process for computational modelling is facilitated by the use of accepted standards with which to describe and simulate models, ensuring consistency between the models and thus reducing the development and propagation of errors. CellML 1.1, an XML-based programming language, has been designed as a modelling standard which, by virtue of its import and grouping functions, facilitates model combination and reuse. Using CellML 1.1, we demonstrate the process of formalized model reuse by combining three separate models of rat cardiomyocyte function (an electrophysiology model, a model of cellular calcium dynamics and a mechanics model) which together make up the Pandit-Hinch-Niederer et al. cell model. Not only is this integrative model of rat electromechanics a useful tool for cardiac modelling but it is also an ideal framework with which to demonstrate both the power of model reuse and the challenges associated with this process. We highlight and classify a number of these issues associated with combining models and provide some suggested solutions.

Effect of bundle branch block on cardiac output: a whole heart simulation study

The heart is an electrically controlled fluid pump which operates by mechanical contraction. Whole heart modelling is a computationally daunting task which must incorporate several subsystems: mechanical, electrical, and fluidic. Numerous feedback mechanisms on many levels, and operating at different scales, exist to finely control behaviour. Understanding these interactions is necessary to understand heart operation, as well as pathologies and therapies. A review of the components in such a model is given. The authors then present a framework for their electro-mechano-fluidic whole heart model based on cable methods. The model incorporates atria and ventricles, and has functioning valves with papillary muscles. The effect of altered propagation due to left and right bundle branch block on cardiac output is examined using the cable-based model. Results are compared to clinically observed phenomena. Good agreement was obtained, but tighter coupling of mechanical and electrical events is needed to fully account for behaviour. Cable-based models offer an alternative to continuum models.

Asymptotic analysis and analytical solutions of a model of cardiac excitation

We describe an asymptotic approach to gated ionic models of single-cell cardiac excitability. It has a form essentially different from the Tikhonov fast-slow form assumed in standard asymptotic reductions of excitable systems. This is of interest since the standard approaches have been previously found inadequate to describe phenomena such as the dissipation of cardiac wave fronts and the shape of action potential at repolarization. The proposed asymptotic description overcomes these deficiencies by allowing, among other non-Tikhonov features, that a dynamical variable may change its character from fast to slow within a single solution. The general asymptotic approach is best demonstrated on an example which should be both simple and generic. The classical model of Purkinje fibers (Noble in J. Physiol. 160:317-352, 1962) has the simplest functional form of all cardiac models but according to the current understanding it assigns a physiologically incorrect role to the Na current. This leads us to suggest an "Archetypal Model" with the simplicity of the Noble model but with a structure more typical to contemporary cardiac models. We demonstrate that the Archetypal Model admits a complete asymptotic solution in quadratures. To validate our asymptotic approach, we proceed to consider an exactly solvable "caricature" of the Archetypal Model and demonstrate that the asymptotic of its exact solution coincides with the solutions obtained by substituting the "caricature" right-hand sides into the asymptotic solution of the generic Archetypal Model. This is necessary, because, unlike in standard asymptotic descriptions, no general results exist which can guarantee the proximity of the non-Tikhonov asymptotic solutions to the solutions of the corresponding detailed ionic model.

Rhythms of high-grade block in an ionic model of a strand of regionally ischemic ventricular muscle

Electrical alternans, a beat-to-beat alternation in the electrocardiogram or electrogram, is frequently seen during the first few minutes of acute myocardial ischemia, and is often immediately followed by malignant cardiac arrhythmias such as ventricular tachycardia and ventricular fibrillation. As ischemia progresses, higher-order periodic rhythms (e.g., period-4) can replace the period-2 alternans rhythm. This is also seen in modelling work on a two-dimensional (2-D) sheet of regionally ischemic ventricular muscle. In addition, in the experimental work, ventricular arrhythmias are overwhelmingly seen only after the higher-order rhythms arise. We investigate an ionic model of a strand of ischemic ventricular muscle, constructed as a 3-cm-long 1-D cable with a centrally located 1-cm-long segment exposed to an elevated extracellular potassium concentration ([K(+)](o)). As [K(+)](o) is raised in this "ischemic segment" to represent one major effect of ongoing ischemia, the sequence of rhythms {1:1-->2:2 (alternans)-->2:1} is seen. With further increase in [K(+)](o), one sees higher-order periodic 2N:M rhythms {2:1-->4:2-->4:1-->6:2-->6:1-->8:2-->8:1}. In a 2N:M cycle, only M of the 2N action potentials generated at the proximal end of the cable successfully traverse the ischemic segment, with the remaining ones being blocked within the ischemic segment. Finally, there is a transition to complete block {8:1-->2:0-->1:0} (in an n:0 rhythm, all action potentials die out within the ischemic segment). Changing the length of the ischemic segment results in different rhythms and transitions being seen: e.g., when the ischemic segment is 2 cm long, the period-6 rhythms are not seen; when it is 0.5 cm long, there is a 3:1 rhythm interposed between the 2:1 and 1:0 rhythms. We discuss the relevance of our results to the experimental observations on the higher-order rhythms that presage reentrant ischemic ventricular arrhythmias.

A guide to modelling cardiac electrical activity in anatomically detailed ventricles

One of the most recent trends in cardiac electrophysiology is the development of integrative anatomically accurate models of the heart, which include description of cardiac activity from sub-cellular and cellular level to the level of the whole organ. In order to construct this type of model, a researcher needs to collect a wide range of information from books and journal articles on various aspects of biology, physiology, electrophysiology, numerical mathematics and computer programming. The aim of this methodological article is to survey recent developments in integrative modelling of electrical activity in the ventricles of the heart, and to provide a practical guide to the resources and tools that are available for work in this exciting and challenging area.