Heart rate variability (beat-to-beat changes in the RR interval) has attracted considerable attention over the last 30+ years (PubMed currently lists >17,000 publications). Clinically, a decrease in heart rate variability is correlated to higher morbidity and mortality in diverse conditions, from heart disease to foetal distress. It is usually attributed to fluctuation in cardiac autonomic nerve activity. We calculated heart rate variability parameters from a variety of cardiac preparations (including humans, living animals, Langendorff-perfused heart and single sinoatrial nodal cell) in diverse species, combining this with data from previously published papers. We show that regardless of conditions, there is a universal exponential decay-like relationship between heart rate variability and heart rate. Using two biophysical models, we develop a theory for this, and confirm that heart rate variability is primarily dependent on heart rate and cannot be used in any simple way to assess autonomic nerve activity to the heart. We suggest that the correlation between a change in heart rate variability and altered morbidity and mortality is substantially attributable to the concurrent change in heart rate. This calls for re-evaluation of the findings from many papers that have not adjusted properly or at all for heart rate differences when comparing heart rate variability in multiple circumstances.

Chronic atrial fibrillation (AF) is associated with structural and electrical remodelling in the atria, which are associated with a high recurrence of AF. Through biophysically detailed computer modelling, this study investigated mechanisms by which AF-induced electrical remodelling promotes and perpetuates AF. A family of Courtemanche–Ramirez–Nattel variant models of human atrial cell action potentials (APs), taking into account of intrinsic atrial electrophysiological properties, was modified to incorporate various experimental data sets on AF-induced changes of major ionic channel currents (ICaL, IKur, Ito, IK1, IKs, INaCa) and on intracellular Ca2+ handling. The single cell models for control and AF-remodelled conditions were incorporated into multicellular three-dimensional (3D) atrial tissue models. Effects of the AF-induced electrical remodelling were quantified as the changes of AP profile, AP duration (APD) and its dispersion across the atria, and the vulnerability of atrial tissue to the initiation of re-entry. The dynamic behaviour of re-entrant excitation waves in the 3D models was characterised. In our simulations, AF-induced electrical remodelling abbreviated atrial APD non-uniformly across the atria; this resulted in relatively short APDs co-existing with marked regional differences in the APD at junctions of the crista terminalis/pectinate muscle, pulmonary veins/left atrium. As a result, the measured tissue vulnerability to re-entry initiation at these tissue junctions was increased. The AF-induced electrical remodelling also stabilized and accelerated re-entrant excitation waves, leading to rapid and sustained re-entry. Under the AF-remodelled condition, re-entrant scroll waves in the 3D model degenerated into persistent and erratic wavelets, leading to fibrillation. In conclusion, realistic 3D atrial tissue models indicate that AF-induced electrical remodelling produces regionally heterogeneous and shortened APD; these respectively facilitate initiation and maintenance of re-entrant excitation waves.

Investigating the mechanisms underlying the genesis and conduction of electrical excitation in the atria at physiological and pathological states is of great importance. To provide knowledge concerning the mechanisms of excitation, we constructed a biophysical detailed and anatomically accurate computer model of human atria that incorporates both structural and electrophysiological heterogeneities. The three-dimensional geometry was extracted from the visible female dataset. The sinoatrial node (SAN) and atrium, including crista terminalis (CT), pectinate muscles (PM), appendages (APG) and Bachmann's bundle (BB) were segmented in this work. Fibre orientation in CT, PM and BB was set to local longitudinal direction. Descriptions for all used cell types were based on modifications of the Courtemanche et al. model of a human atrial cell. Maximum conductances of Ito, IKr and ICa,L were modified for PM, CT, APG and atrioventricular ring to reproduce measured action potentials (AP). Pacemaker activity in the human SAN was reproduced by removing IK1, but including If, ICa,T, and gradients of channel conductances as described in previous studies for heterogeneous rabbit SAN. Anisotropic conduction was computed with a monodomain model using the finite element method. The transversal to longitudinal ratio of conductivity for PM, CT and BB was 1:9. Atrial working myocardium (AWM) was set to be isotropic. Simulation of atrial electrophysiology showed initiation of APs in the SAN centre. The excitation spread afterwards to the periphery near to the region of the CT and preferentially towards the atrioventricular region. The excitation extends over the right atrium along PM. Both CT and PM activated the right AWM. Earliest activation of the left atrium was through BB and excitation spread over to the APG. The conduction velocities were 0.6ms-1 for AWM, 1.2ms-1 for CT, 1.6ms-1 for PM and 1.1ms-1 for BB at a rate of 63bpm. The simulations revealed that bundles form dominant pathways for atrial conduction. The preferential conduction towards CT and along PM is comparable with clinical mapping. Repolarization is more homogeneous than excitation due to the heterogeneous distribution of electrophysiological properties and hence the action potential duration.

Despite a vast amount of experimental and clinical data on the underlying ionic, cellular and tissue substrates, the mechanisms of common atrial arrhythmias (such as atrial fibrillation, AF) arising from the functional interactions at the whole atria level remain unclear. Computational modelling provides a quantitative framework for integrating such multi-scale data and understanding the arrhythmogenic behaviour that emerges from the collective spatio-temporal dynamics in all parts of the heart. In this study, we have developed a multi-scale hierarchy of biophysically detailed computational models for the human atria--the 3D virtual human atria. Primarily, diffusion tensor MRI reconstruction of the tissue geometry and fibre orientation in the human sinoatrial node (SAN) and surrounding atrial muscle was integrated into the 3D model of the whole atria dissected from the Visible Human dataset. The anatomical models were combined with the heterogeneous atrial action potential (AP) models, and used to simulate the AP conduction in the human atria under various conditions: SAN pacemaking and atrial activation in the normal rhythm, break-down of regular AP wave-fronts during rapid atrial pacing, and the genesis of multiple re-entrant wavelets characteristic of AF. Contributions of different properties of the tissue to mechanisms of the normal rhythm and arrhythmogenesis were investigated. Primarily, the simulations showed that tissue heterogeneity caused the break-down of the normal AP wave-fronts at rapid pacing rates, which initiated a pair of re-entrant spiral waves; and tissue anisotropy resulted in a further break-down of the spiral waves into multiple meandering wavelets characteristic of AF. The 3D virtual atria model itself was incorporated into the torso model to simulate the body surface ECG patterns in the normal and arrhythmic conditions. Therefore, a state-of-the-art computational platform has been developed, which can be used for studying multi-scale electrical phenomena during atrial conduction and AF arrhythmogenesis. Results of such simulations can be directly compared with electrophysiological and endocardial mapping data, as well as clinical ECG recordings. The virtual human atria can provide in-depth insights into 3D excitation propagation processes within atrial walls of a whole heart in vivo, which is beyond the current technical capabilities of experimental or clinical set-ups.

No abstract

Early development of ionic models for cardiac myocytes, from the pioneering modification of the Hodgkin-Huxley giant squid axon model by Noble to the iconic DiFrancesco-Noble model integrating voltage-gated ionic currents, ion pumps and exchangers, Ca 2C sequestration and Ca 2C -induced Ca 2C release, provided a general description for a mammalian Purkinje fibre (PF) and the framework for modern cardiac models. In the past two decades, development has focused on tissue-specific models with an emphasis on the sino-atrial (SA) node, atria and ventricles, while the PFs have largely been neglected. However, achieving the ultimate goal of creating a virtual human heart will require detailed models of all distinctive regions of the cardiac conduction system, including the PFs, which play an important role in conducting cardiac excitation and ensuring the synchronized timing and sequencing of ventricular contraction. In this paper, we present details of our newly developed model for the human PF cell including validation against experimental data. Ionic mechanisms underlying the heterogeneity between the PF and ventricular action potentials in humans and other species are analysed. The newly developed PF cell model adds a new member to the family of human cardiac cell models developed previously for the SA node, atrial and ventricular cells, which can be incorporated into an anatomical model of the human heart with details of its electrophysiological heterogeneity and anatomical complexity.

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.

We have examined sino-atrial node (SAN) function in hearts from adult mice with heterozygous targeted disruption of the Scn5a gene to clarify the role of Scn5a-encoded cardiac Na + channels in normal SAN function and the mechanism(s) by which reduced Na + channel function might cause sinus node dysfunction. Scn5a +/− mice showed depressed heart rates and occasional sino-atrial (SA) block. Their isolated peripheral SAN pacemaker cells showed a reduced Na + channel expression and slowed intrinsic pacemaker rates. Wild-type (WT) and Scn5a +/− SAN preparations exhibited similar activation patterns but with significantly slower SA conduction and frequent sino-atrial conduction block in Scn5a +/− SAN preparations. Furthermore, isolated WT and Scn5a +/− SAN cells demonstrated differing correlations between cycle length, maximum upstroke velocity and action potential amplitude, and cell size. Small myocytes showed similar, but large myocytes reduced pacemaker rates, implicating the larger peripheral SAN cells in the reduced pacemaker rate that was observed in Scn5a +/− myocytes. These findings were successfully reproduced in a model that implicated i Na directly in action potential propagation through the SAN and from SAN to atria, and in modifying heart rate through a coupling of SAN and atrial cells. Functional alterations in the SAN following heterozygous-targeted disruption of Scn5a thus closely resemble those observed in clinical sinus node dysfunction. The findings accordingly provide a basis for understanding of the role of cardiac-type Na + channels in normal SAN function and the pathophysiology of sinus node dysfunction and suggest new potential targets for its clinical management.