The Ca2+ transient as a feedback sensor controlling cardiomyocyte ionic conductances in mouse populations.

Conductances of ion channels and transporters controlling cardiac excitation may vary in a population of subjects with different cardiac gene expression patterns. However, the amount of variability and its origin are not quantitatively known. We propose a new conceptual approach to predict this variability that consists of finding combinations of conductances generating a normal intracellular Ca2+ transient without any constraint on the action potential. Furthermore, we validate experimentally its predictions using the Hybrid Mouse Diversity Panel, a model system of genetically diverse mouse strains that allows us to quantify inter-subject versus intra-subject variability. The method predicts that conductances of inward Ca2+ and outward K+ currents compensate each other to generate a normal Ca2+ transient in good quantitative agreement with current measurements in ventricular myocytes from hearts of different isogenic strains. Our results suggest that a feedback mechanism sensing the aggregate Ca2+ transient of the heart suffices to regulate ionic conductances

Early afterdepolarisations and ventricular arrhythmias in cardiac tissue: a computational study

Afterdepolarisations are associated with arrhythmias in the heart, but are difficult to study experimentally. In this study we used a simplified computational model of 1D and 2D cardiac ventricular tissue, where we could control the size of the region generating afterdepolarisations, as well as the properties of the afterdepolarisation waveform. Provided the size of the afterdepolarisation region was greater than around 1 mm, propagating extrasystoles were produced in both 1D and 2D. The number of extrasystoles produced depended on the amplitude, period, and duration of the oscillatory EAD waveform. In 2D, re-entry was also initiated for specific combinations of EAD amplitude, period, and duration, with the afterdepolarisation region acting as a common pathway. The main finding from this modelling study is therefore that afterdepolarisations can act as potent sources of propagating extrasystoles, as well as a source of re-entrant activation

A comparative study of early afterdepolarization-mediated fibrillation in two mathematical models for human ventricular cells

Early afterdepolarizations (EADs), which are abnormal oscillations of the membrane potential at the plateau phase of an action potential, are implicated in the development of cardiac arrhythmias like Torsade de Pointes. We carry out extensive numerical simulations of the TP06 and ORd mathematical models for human ventricular cells with EADs. We investigate the different regimes in both these models, namely, the parameter regimes where they exhibit (1) a normal action potential (AP) with no EADs, (2) an AP with EADs, and (3) an AP with EADs that does not go back to the resting potential. We also study the dependence of EADs on the rate of at which we pace a cell, with the specific goal of elucidating EADs that are induced by slow or fast rate pacing. In our simulations in two-and three-dimensional domains, in the presence of EADs, we find the following wave types: (A) waves driven by the fast sodium current and the L-type calcium current (Na-Ca-mediated waves); (B) waves driven only by the L-type calcium current (Ca-mediated waves); (C) phase waves, which are pseudo-travelling waves. Furthermore, we compare the wave patterns of the various wave-types (Na-Ca-mediated, Ca-mediated, and phase waves) in both these models. We find that the two models produce qualitatively similar results in terms of exhibiting Na-Ca-mediated wave patterns that are more chaotic than those for the Ca-mediated and phase waves. However, there are quantitative differences in the wave patterns of each wave type. The Na-Ca-mediated waves in the ORd model show short-lived spirals but the TP06 model does not. The TP06 model supports more Ca-mediated spirals than those in the ORd model, and the TP06 model exhibits more phase-wave patterns than does the ORd model

Bridging the gap between computation and clinical biology: validation of cable theory in humans

Introduction: Computerized simulations of cardiac activity have significantly contributed to our understanding of cardiac electrophysiology, but techniques of simulations based on patient-acquired data remain in their infancy. We sought to integrate data acquired from human electrophysiological studies into patient-specific models, and validated this approach by testing whether electrophysiological responses to sequential premature stimuli could be predicted in a quantitatively accurate manner. Methods: Eleven patients with structurally normal hearts underwent electrophysiological studies. Semi-automated analysis was used to reconstruct activation and repolarization dynamics for each electrode. This S(2) extrastimuli data was used to inform individualized models of cardiac conduction, including a novel derivation of conduction velocity restitution. Activation dynamics of multiple premature extrastimuli were then predicted from this model and compared against measured patient data as well as data derived from the ten-Tusscher cell-ionic model. Results: Activation dynamics following a premature S(3) were significantly different from those after an S(2). Patient specific models demonstrated accurate prediction of the S(3) activation wave, (Pearson's R(2) = 0.90, median error 4%). Examination of the modeled conduction dynamics allowed inferences into the spatial dispersion of activation delay. Further validation was performed against data from the ten-Tusscher cell-ionic model, with our model accurately recapitulating predictions of repolarization times (R(2) = 0.99). Conclusions: Simulations based on clinically acquired data can be used to successfully predict complex activation patterns following sequential extrastimuli. Such modeling techniques may be useful as a method of incorporation of clinical data into predictive models

Endogenous driving and synchronization in cardiac and uterine virtual tissues: bifurcations and local coupling

Cardiac and uterine muscle cells and tissue can be either autorhythmic or excitable. These behaviours exchange stability at bifurcations produced by changes in parameters, which if spatially localized can produce an ectopic pacemaking focus. The effects of these parameters on cell dynamics have been identified and quantified using continuation algorithms and by numerical solutions of virtual cells. The ability of a compact pacemaker to drive the surrounding excitable tissues depends on both the size of the pacemaker and the strength of electrotonic coupling between cells within, between, and outside the pacemaking region. We investigate an ectopic pacemaker surrounded by normal excitable tissue. Cell–cell coupling is simulated by the diffusion coefficient for voltage. For uniformly coupled tissues, the behaviour of the hybrid tissue can take one of the three forms: (i) the surrounding tissue electrotonically suppresses the pacemaker; (ii) depressed rate oscillatory activity in the pacemaker but no propagation; and (iii) pacemaker driving propagations into the excitable region. However, real tissues are heterogeneous with spatial changes in cell–cell coupling. In the gravid uterus during early pregnancy, cells are weakly coupled, with the cell–cell coupling increasing during late pregnancy, allowing synchronous contractions during labour. These effects are investigated for a caricature uterine tissue by allowing both excitability and diffusion coefficient to vary stochastically with space, and for cardiac tissues by spatial gradients in the diffusion coefficient

High speed sCMOS-based oblique plane microscopy applied to the study of calcium dynamics in cardiac myocytes

blique plane microscopy (OPM) is a form of light sheet microscopy that uses a single high numerical aperture microscope objective for both fluorescence excitation and collection. In this paper, measurements of the relative collection efficiency of OPM are presented. An OPM system incorporating two sCMOS cameras is then introduced that enables single isolated cardiac myocytes to be studied continuously for 22 seconds in two dimensions at 667 frames per second with 960 × 200 pixels and for 30 seconds with 960 × 200 × 20 voxels at 25 volumes per second. In both cases OPM is able to record in two spectral channels, enabling intracellular calcium to be studied via the probe Fluo-4 AM simultaneously with the sarcolemma and transverse tubule network via the membrane dye Cellmask Orange. The OPM system was then applied to determine the spatial origin of spontaneous calcium waves for the first time and to measure the cell transverse tubule structure at their point of origin. Further results are presented to demonstrate that the OPM system can also be used to study calcium spark parameters depending on their relationship to the transverse tubule structure

Effect of myocyte-fibroblast coupling on the onset of pathological dynamics in a model of ventricular tissue

Managing lethal cardiac arrhythmias is one of the biggest challenges in modern cardiology, and hence it is very important to understand the factors underlying such arrhythmias. While early afterdepolarizations (EAD) of cardiac cells is known to be one such arrhythmogenic factor, the mechanisms underlying the emergence of tissue level arrhythmias from cellular level EADs is not fully understood. Another known arrhythmogenic condition is fibrosis of cardiac tissue that occurs both due to aging and in many types of heart diseases. In this paper we describe the results of a systematic insilico study, using the TNNP model of human cardiac cells and MacCannell model for (myo) fibroblasts, on the possible effects of diffuse fibrosis on arrhythmias occurring via EADs. We find that depending on the resting potential of fibroblasts (VFR), M-F coupling can either increase or decrease the region of parameters showing EADs. Fibrosis increases the probability of occurrence of arrhythmias after a single focal stimulation and this effect increases with the strength of the M-F coupling. While in our simulations, arrhythmias occur due to fibrosis induced ectopic activity, we do not observe any specific fibrotic pattern that promotes the occurrence of these ectopic sources

Examination of the Effects of Heterogeneous Organization of RyR Clusters, Myofibrils and Mitochondria on Ca2+ Release Patterns in Cardiomyocytes

Spatio-temporal dynamics of intracellular calcium, [Ca2+]i, regulate the contractile function of cardiac muscle cells. Measuring [Ca2+]i flux is central to the study of mechanisms that underlie both normal cardiac function and calcium-dependent etiologies in heart disease. However, current imaging techniques are limited in the spatial resolution to which changes in [Ca2+]i can be detected. Using spatial point process statistics techniques we developed a novel method to simulate the spatial distribution of RyR clusters, which act as the major mediators of contractile Ca2+ release, upon a physiologically-realistic cellular landscape composed of tightly-packed mitochondria and myofibrils.We applied this method to computationally combine confocal-scale (~ 200 nm) data of RyR clusters with 3D electron microscopy data (~ 30 nm) of myofibrils and mitochondria, both collected from adult rat left ventricular myocytes. Using this hybrid-scale spatial model, we simulated reaction-diffusion of [Ca2+]i during the rising phase of the transient (first 30 ms after initiation). At 30 ms, the average peak of the simulated [Ca2+]i transient and of the simulated fluorescence intensity signal, F/F0, reached values similar to that found in the literature ([Ca2+]i 1 μM; F/F0 5.5). However, our model predicted the variation in [Ca2+]i to be between 0.3 and 12.7 μM (~3 to 100 fold from resting value of 0.1 μM) and the corresponding F/F0 signal ranging from 3 to 9.5. We demonstrate in this study that: (i) heterogeneities in the [Ca2+]i transient are due not only to heterogeneous distribution and clustering of mitochondria; (ii) but also to heterogeneous local densities of RyR clusters. Further, we show that: (iii) these structureinduced heterogeneities in [Ca2+]i can appear in line scan data. Finally, using our unique method for generating RyR cluster distributions, we demonstrate the robustness in the [Ca2+]i transient to differences in RyR cluster distributions measured between rat and human cardiomyocytes