Caloric Restriction Rejuvenates Skeletal Muscle Growth in Heart Failure With Preserved Ejection Fraction.

Heart failure with preserved ejection fraction (HFpEF) is a major clinical problem, with limited treatments. HFpEF is characterized by a distinct, but poorly understood, skeletal muscle pathology, which could offer an alternative therapeutic target. In a rat model, we identified impaired myonuclear accretion as a mechanism for low myofiber growth in HFpEF following resistance exercise. Acute caloric restriction rescued skeletal muscle pathology in HFpEF, whereas cardiac therapies had no effect. Mechanisms regulating myonuclear accretion were dysregulated in patients with HFpEF. Overall, these findings may have widespread implications in HFpEF, indicating combined dietary with exercise interventions as a beneficial approach to overcome skeletal muscle pathology.

Interactions between calcium-induced arrhythmia triggers and the electrophysiological-anatomical substrate underlying the induction of atrial fibrillation.

Atrial fibrillation (AF) is the most common cardiac arrhythmia and is sustained by spontaneous focal excitations and re-entry. Spontaneous electrical firing in the pulmonary vein (PV) sleeves is implicated in AF generation. The aim of this simulation study was to identify the mechanisms determining the localisation of AF triggers in the PVs and their contribution to the genesis of AF. A novel biophysical model of the canine atria was used that integrates stochastic, spontaneous subcellular Ca2+ release events (SCRE) with regional electrophysiological heterogeneity in ionic properties and a detailed three-dimensional model of atrial anatomy, microarchitecture and patchy fibrosis. Simulations highlighted the importance of the smaller inward rectifier potassium current (IK1 ) in PV cells compared to the surrounding atria, which enabled SCRE more readily to result in delayed-afterdepolarisations that induced triggered activity. There was a leftward shift in the dependence of the probability of triggered activity on sarcoplasmic reticulum Ca2+ load. This feature was accentuated in 3D tissue compared to single cells (Δ half-maximal [Ca2+ ]SR  = 58 µM vs. 22 µM). In 3D atria incorporating electrical heterogeneity, excitations preferentially emerged from the PV region. These triggered focal excitations resulted in transient re-entry in the left atrium. Addition of fibrotic patches promoted localised emergence of focal excitations and wavebreaks that had a more substantial impact on generating AF-like patterns than the PVs. Thus, a reduced IK1 , less negative resting membrane potential, and fibrosis-induced changes of the electrotonic load all contribute to the emergence of complex excitation patterns from spontaneous focal triggers. KEY POINTS: Focal excitations in the atria are most commonly associated with the pulmonary veins, but the mechanisms for this localisation are yet to be elucidated. We applied a multi-scale computational modelling approach to elucidate the mechanisms underlying such localisations. Myocytes in the pulmonary vein region of the atria have a less negative resting membrane potential and reduced time-independent potassium current; we demonstrate that both of these factors promote triggered activity in single cells and tissues. The less negative resting membrane potential also contributes to heterogeneous inactivation of the fast sodium current, which can enable re-entrant-like excitation patterns to emerge without traditional conduction block.

Patchy fibrosis promotes trigger-substrate interactions that both generate and maintain atrial fibrillation.

Fibrosis has been mechanistically linked to arrhythmogenesis in multiple cardiovascular conditions, including atrial fibrillation (AF). Previous studies have demonstrated that fibrosis can create functional barriers to conduction which may promote excitation wavebreak and the generation of re-entry, while also acting to pin re-entrant excitation in stable rotors during AF. However, few studies have investigated the role of fibrosis in the generation of AF triggers in detail. We apply our in-house computational framework to study the impact of fibrosis on the generation of AF triggers and trigger-substrate interactions in two- and three-dimensional atrial tissue models. Our models include a reduced and efficient description of stochastic, spontaneous cellular triggers as well as a simple model of heterogeneous inter-cellular coupling. Our results demonstrate that fibrosis promotes the emergence of focal excitations, primarily through reducing the electrotonic load on individual fibre strands. This enables excitation to robustly initiate within these single strands before spreading to neighbouring strands and inducing a full tissue focal excitation. Enhanced conduction block can allow trigger-substrate interactions that result in the emergence of complex, re-entrant excitation patterns. This study provides new insight into the mechanisms by which fibrosis promotes the triggers and substrate necessary to induce and sustain arrhythmia.

A simple approach for image-based modelling of the heart that enables robust simulation of highly heterogeneous electrical excitation.

Remodelling of cardiac tissue structure, including intercellular electrical coupling, is a major determinant of the complex and heterogeneous excitation patterns associated with cardiac arrhythmias. Evaluation of the precise mechanisms by which local tissue structure determines global arrhythmic excitation patterns is a major challenge that may be critically important for the development of effective treatment strategies. Computational modelling is a key tool in the study of cardiac arrhythmias, yet the established approaches for organ-scale modelling are unsuitable to capture the impact of local conduction heterogeneities; a novel approach is required to provide this multi-scale mechanistic insight. We present a fundamentally simple yet powerful approach to simulate electrical excitation in highly heterogeneous whole-heart models that exploits the underlying discreteness of the myocardium. Preliminary simulations demonstrate that this approach can capture lower conduction velocities and reproduce wave breakdown and the development of re-entry in a range of conditions.

Ionic current changes underlying action potential repolarization responses to physiological pacing and adrenergic stimulation in adult rat ventricular myocytes.

This study aimed to simulate ventricular responses to elevations in myocyte pacing and adrenergic stimulation using a novel electrophysiological rat model and investigate ion channel responses underlying action potential (AP) modulations. Peak ion currents and AP repolarization to 50% and 90% of full repolarization (APD50-90 ) were recorded during simulations at 1-10 Hz pacing under control and adrenergic stimulation conditions. Further simulations were performed with incremental ion current block (L-type calcium current, ICa ; transient outward current, Ito ; slow delayed rectifier potassium current, IKs ; rapid delayed rectifier potassium current, IKr ; inward rectifier potassium current, IK1 ) to identify current influence on AP response to exercise. Simulated APD50-90 closely resembled experimental findings. Rate-dependent increases in IKs (6%-101%), IKr (141%-1339%), and ICa (0%-15%) and reductions in Ito (11%-57%) and IK1 (1%-9%) were observed. Meanwhile, adrenergic stimulation triggered moderate increases in all currents (23%-67%) except IK1 . Further analyses suggest AP plateau is most sensitive to modulations in Ito and ICa while late repolarization is most sensitive to IK1 , ICa , and IKs , with alterations in IKs predominantly stimulating the greatest magnitude of influence on late repolarization (35%-846% APD90 prolongation). The modified Leeds rat model (mLR) is capable of accurately modeling APs during physiological stress. This study highlights the importance of ICa , Ito , IK1, and IKs in controlling electrophysiological responses to exercise. This work will benefit the study of cardiac dysfunction, arrythmia, and disease, though future physiologically relevant experimental studies and model development are required.

Multi-Scale Computational Modeling of Spatial Calcium Handling From Nanodomain to Whole-Heart: Overview and Perspectives.

Regulation of intracellular calcium is a critical component of cardiac electrophysiology and excitation-contraction coupling. The calcium spark, the fundamental element of the intracellular calcium transient, is initiated in specialized nanodomains which co-locate the ryanodine receptors and L-type calcium channels. However, calcium homeostasis is ultimately regulated at the cellular scale, by the interaction of spatially separated but diffusively coupled nanodomains with other sub-cellular and surface-membrane calcium transport channels with strong non-linear interactions; and cardiac electrophysiology and arrhythmia mechanisms are ultimately tissue-scale phenomena, regulated by the interaction of a heterogeneous population of coupled myocytes. Recent advances in imaging modalities and image-analysis are enabling the super-resolution reconstruction of the structures responsible for regulating calcium homeostasis, including the internal structure of nanodomains themselves. Extrapolating functional and imaging data from the nanodomain to the whole-heart is non-trivial, yet essential for translational insight into disease mechanisms. Computational modeling has important roles to play in relating structural and functional data at the sub-cellular scale and translating data across the scales. This review covers recent methodological advances that enable image-based modeling of the single nanodomain and whole cardiomyocyte, as well as the development of multi-scale simulation approaches to integrate data from nanometer to whole-heart. Firstly, methods to overcome the computational challenges of simulating spatial calcium dynamics in the nanodomain are discussed, including image-based modeling at this scale. Then, recent whole-cell models, capable of capturing a range of different structures (such as the T-system and mitochondria) and cellular heterogeneity/variability are discussed at two different levels of discretization. Novel methods to integrate the models and data across the scales and simulate stochastic dynamics in tissue-scale models are then discussed, enabling elucidation of the mechanisms by which nanodomain remodeling underlies arrhythmia and contractile dysfunction. Perspectives on model differences and future directions are provided throughout.

Multi-scale approaches for the simulation of cardiac electrophysiology: II - Tissue-level structure and function.

Computational models of the heart, from cell-level models, through one-, two- and three-dimensional tissue-level simplifications, to biophysically-detailed three-dimensional models of the ventricles, atria or whole heart, allow the simulation of excitation and propagation of this excitation, and have provided remarkable insight into the normal and pathological functioning of the heart. In this article we present equations for modelling cellular excitation (i.e. the cell action potential) from both a phenomenological and a biophysical perspective. Hodgkin-Huxley formalism is discussed, along with the current generation of biophysically-detailed cardiac cell models. Alternative Markovian formulations for modelling ionic currents are also presented. Equations describing propagation of this cellular excitation, through one-, two- and three-dimensional idealised or realistic tissues, are then presented. For all types of model, from cell to tissue, methods for discretisation and integration of the underlying equations are discussed. The article finishes with a discussion of two tissue-level experimental imaging techniques - diffusion tensor magnetic resonance imaging and optical imaging - that can be used to provide data for parameterisation and validation of cell- and tissue-level cardiac models.

Multi-scale approaches for the simulation of cardiac electrophysiology: I - Sub-cellular and stochastic calcium dynamics from cell to organ.

Computational models of the heart at multiple spatial scales, from sub-cellular nanodomains to the whole-organ, are a powerful tool for the simulation of cardiac electrophysiology. Application of these models has provided remarkable insight into the normal and pathological functioning of the heart. In these two articles, we present methods for modelling cardiac electrophysiology at all of these spatial scales. In part one, presented here, we discuss methods and approaches for modelling sub-cellular calcium dynamics at the whole-cell and organ scales, valuable for modelling excitation-contraction coupling and mechanisms of arrhythmia triggers.

Investigation of the Role of Myocyte Orientations in Cardiac Arrhythmia Using Image-Based Models.

Cardiac electrical excitation-propagation is influenced by myocyte orientations (cellular organization). Quantitatively understanding this relationship presents a significant research challenge, especially during arrhythmias in which excitation patterns become complex. Tissue-scale simulations of cardiac electrophysiology, incorporating both dynamic action potential behavior and image-based myocardial architecture, provide an approach to investigate three-dimensional (3D) propagation of excitation waves in the heart. In this study, we aimed to assess the importance of natural variation in myocyte orientations on cardiac arrhythmogenesis using 3D tissue electrophysiology simulations. Three anatomical models (i.e., describing myocyte orientations) of healthy rat ventricles-obtained using diffusion tensor imaging at 100 µm resolution-were registered to a single biventricular geometry (i.e., a single cardiac shape), in which the myocyte orientations could be represented by each of the diffusion tensor imaging data sets or by an idealized rule-based description. The Fenton-Karma cellular excitation model was modified to reproduce rat ventricular action potential duration restitution to create reaction-diffusion cardiac electrophysiology models. Over 250 3D simulations were performed to investigate the effects of myocyte orientations on the following: 1) ventricular activation, 2) location-dependent arrhythmia induction via rapid pacing, and 3) dynamics of re-entry averaged over multiple episodes. It was shown that 1) myocyte orientation differences manifested themselves in local activation times, but the influence on total activation time was small; 2) differences in myocyte orientations could critically affect the inducibility and persistence of arrhythmias for specific stimulus-location/cycle-length combinations; and 3) myocyte orientations alone could be an important determinant of scroll wave break, although no significant differences were observed in averaged arrhythmia dynamics between the four myocyte orientation scenarios considered. Our results show that myocyte orientations are an important determinant of arrhythmia inducibility, persistence, and scroll wave break. These findings suggest that where specificity is desired (for example, when predicting location-dependent, patient-specific arrhythmia inducibility), subject-specific myocyte orientations may be important.

Arrhythmia mechanisms and spontaneous calcium release: Bi-directional coupling between re-entrant and focal excitation.

Spontaneous sub-cellular calcium release events (SCRE) are conjectured to promote rapid arrhythmias associated with conditions such as heart failure and atrial fibrillation: they can underlie the emergence of spontaneous action potentials in single cells which can lead to arrhythmogenic triggers in tissue. The multi-scale mechanisms of the development of SCRE into arrhythmia triggers, and their dynamic interaction with the tissue substrate, remain elusive; rigorous and simultaneous study of dynamics from the nanometre to the centimetre scale is a major challenge. The aim of this study was to develop a computational approach to overcome this challenge and study potential bi-directional coupling between sub-cellular and tissue-scale arrhythmia phenomena. A framework comprising a hierarchy of computational models was developed, which includes detailed single-cell models describing spatio-temporal calcium dynamics in 3D, efficient non-spatial cell models, and both idealised and realistic tissue models. A phenomenological approach was implemented to reproduce SCRE morphology and variability in the efficient cell models, comprising the definition of analytical Spontaneous Release Functions (SRF) whose parameters may be randomly sampled from appropriate distributions in order to match either the 3D cell models or experimental data. Pro-arrhythmogenic pacing protocols were applied to initiate re-entry and promote calcium overload, leading to the emergence of SCRE. The SRF accurately reproduced the dynamics of SCRE and its dependence on environment variables under multiple different conditions. Sustained re-entrant excitation promoted calcium overload, and led to the emergence of focal excitations after termination. A purely functional mechanism of re-entry and focal activity localisation was demonstrated, related to the unexcited spiral wave core. In conclusion, a novel approach has been developed to dynamically model SCRE at the tissue scale, which facilitates novel, detailed multi-scale mechanistic analysis. It was revealed that complex re-entrant excitation patterns and SCRE may be bi-directionally coupled, promoting novel mechanisms of arrhythmia perpetuation.

Dynamic clamping human and rabbit atrial calcium current: narrowing ICaL window abolishes early afterdepolarizations.

KEY POINTS: Early-afterdepolarizations (EADs) are abnormal action potential oscillations and a known cause of cardiac arrhythmias. Ventricular EADs involve reactivation of a Ca2+ current (ICaL ) in its 'window region' voltage range. However, electrical mechanisms of atrial EADs, a potential cause of atrial fibrillation, are poorly understood. Atrial cells were obtained from consenting patients undergoing heart surgery, as well as from rabbits. ICaL was blocked with nifedipine and then a hybrid patch clamp/mathematical-modelling technique, 'dynamic clamping', was used to record action potentials at the same time as injecting an artificial, modifiable, ICaL (ICaL,D-C ). Progressively widening the ICaL,D-C window region produced EADs of various types, dependent on window width. EAD production was strongest upon moving the activation (vs. inactivation) side of the window. EADs were then induced by a different method: increasing ICaL,D-C amplitude and/or K+ channel-blockade (4-aminopyridine). Narrowing of the ICaL,D-C window by ∼10 mV abolished these EADs. Atrial ICaL window narrowing is worthy of further testing as a potential anti-atrial fibrillation drug mechanism. ABSTRACT: Atrial early-afterdepolarizations (EADs) may contribute to atrial fibrillation (AF), perhaps involving reactivation of L-type Ca2+ current (ICaL ) in its window region voltage range. The present study aimed (i) to validate the dynamic clamp technique for modifying the ICaL contribution to atrial action potential (AP) waveform; (ii) to investigate the effects of widening the window ICaL on EAD-propensity; and (iii) to test whether EADs from increased ICaL and AP duration are supressed by narrowing the window ICaL . ICaL and APs were recorded from rabbit and human atrial myocytes by whole-cell-patch clamp. During AP recording, ICaL was inhibited (3 µm nifedipine) and replaced by a dynamic clamp model current, ICaL,D-C (tuned to native ICaL characteristics), computed in real-time (every 50 µs) based on myocyte membrane potential. ICaL,D-C -injection restored the nifedipine-suppressed AP plateau. Widening the window ICaL,D-C , symmetrically by stepwise simultaneous equal shifts of half-voltages (V0.5 ) of ICaL,D-C activation (negatively) and inactivation (positively), generated EADs (single, multiple or preceding repolarization failure) in a window width-dependent manner, as well as AP alternans. A stronger EAD-generating effect resulted from independently shifting activation V0.5 (asymmetrical widening) than inactivation V0.5 ; for example, a 15 mV activation shift produced EADs in nine of 17 (53%) human atrial myocytes vs. 0 of 18 from inactivation shift (P < 0.05). In 11 rabbit atrial myocytes in which EADs were generated either by increasing the conductance of normal window width ICaL,D-C or subsequent 4-aminopyridine (2 mm), window ICaL,D-C narrowing (10 mV) abolished EADs of all types (P < 0.05). The present study validated the dynamic clamp for ICaL , which is novel in atrial cardiomyocytes, and showed that EADs of various types are generated by widening (particularly asymmetrically) the window ICaL , as well as abolished by narrowing it. Window ICaL narrowing is a potential therapeutic mechanism worth pursuing in the search for improved anti-AF drugs.

Effects of Heart Rate and Ventricular Wall Thickness on Non-invasive Mapping: An in silico Study.

Background: Non-invasive cardiac mapping-also known as Electrocardiographic imaging (ECGi)-is a novel, painless and relatively economic method to map the electrical activation and repolarization patterns of the heart, providing a valuable tool for early identification and diagnosis of conduction abnormalities and arrhythmias. Moreover, the ability to obtain information on cardiac electrical activity non-invasively using ECGi provides the potential for a priori information to guide invasive surgical procedures, improving success rates, and reducing procedure time. Previous studies have shown the influence of clinical variables, such as heart rate, heart size, endocardial wall, and body composition on surface electrocardiogram (ECG) measurements. The influence of clinical variables on the ECG variability has provided information on cardiovascular control and its abnormalities in various pathologies. However, the effects of such clinical variables on the Body Surface Potential (BSP) and ECGi maps have yet to be systematically investigated. Methods: In this study we investigated the effects of heart size, intracardiac thickness, and heart rate on BSP and ECGi maps using a previously-developed 3D electrophysiologically-detailed ventricles-torso model. The inverse solution was solved using the three different Tikhonov regularization methods. Results: Through comparison of multiple measures of error/accuracy on the ECGi reconstructions, our results showed that using different heart geometries to solve the forward and inverse problems produced a larger estimated focal excitation location. An increase of ~2 mm in the Euclidean distance error was observed for an increase in the heart size. However, the estimation of the location of focal activity was still able to be obtained. Similarly, a Euclidean distance increase was observed when the order of regularization was reduced. For the case of activation maps reconstructed at the same ectopic focus location but different heart rates, an increase in the errors and Euclidean distance was observed when the heart rate was increased. Conclusions: Non-invasive cardiac mapping can still provide useful information about cardiac activation patterns for the cases when a different geometry is used for the inverse problem compared to the one used for the forward solution; rapid pacing rates can induce order-dependent errors in the accuracy of reconstruction.

Three-Dimensional and Chemical Mapping of Intracellular Signaling Nanodomains in Health and Disease with Enhanced Expansion Microscopy.

Nanodomains are intracellular foci which transduce signals between major cellular compartments. One of the most ubiquitous signal transducers, the ryanodine receptor (RyR) calcium channel, is tightly clustered within these nanodomains. Super-resolution microscopy has previously been used to visualize RyR clusters near the cell surface. A majority of nanodomains located deeper within cells have remained unresolved due to limited imaging depths and axial resolution of these modalities. A series of enhancements made to expansion microscopy allowed individual RyRs to be resolved within planar nanodomains at the cell periphery and the curved nanodomains located deeper within the interiors of cardiomyocytes. With a resolution of ∼ 15 nm, we localized both the position of RyRs and their individual phosphorylation for the residue Ser2808. With a three-dimensional imaging protocol, we observed disturbances to the RyR arrays in the nanometer scale which accompanied right-heart failure caused by pulmonary hypertension. The disease coincided with a distinct gradient of RyR hyperphosphorylation from the edge of the nanodomain toward the center, not seen in healthy cells. This spatial profile appeared to contrast distinctly from that sustained by the cells during acute, physiological hyperphosphorylation when they were stimulated with a ß-adrenergic agonist. Simulations of RyR arrays based on the experimentally determined channel positions and phosphorylation signatures showed how the nanoscale dispersal of the RyRs during pathology diminishes its intrinsic likelihood to ignite a calcium signal. It also revealed that the natural topography of RyR phosphorylation could offset potential heterogeneity in nanodomain excitability which may arise from such RyR reorganization.

Human Atrial Arrhythmogenesis and Sinus Bradycardia in KCNQ1-Linked Short QT Syndrome: Insights From Computational Modelling.

Atrial fibrillation (AF) and sinus bradycardia have been reported in patients with short QT syndrome variant 2 (SQT2), which is underlain by gain-of-function mutations in KCNQ1 encoding the α subunit of channels carrying slow delayed rectifier potassium current, I Ks. However, the mechanism(s) underlying the increased atrial arrhythmogenesis and impaired cardiac pacemaking activity arising from increased I Ks remain unclear. Possible pharmacological interventions of AF in the SQT2 condition also remain to be elucidated. Using computational modelling, we assessed the functional impact of SQT2 mutations on human sinoatrial node (SAN) pacemaking, atrial repolarisation and arrhythmogenesis, and efficacy of the anti-arrhythmic drug quinidine. Markov chain formulations of I Ks describing two KCNQ1 mutations - V141M and V307L - were developed from voltage-clamp experimental data and then incorporated into contemporary action potential (AP) models of human atrial and SAN cells, the former of which were integrated into idealised and anatomically detailed tissue models. Both mutations shortened atrial AP duration (APD) through distinct I Ks 'gain-of-function' mechanisms, whereas SAN pacemaking rate was slowed markedly only by the V141M mutation. Differences in APD restitution steepness influenced re-entry dynamics in tissue - the V141M mutation promoted stationary and stable spiral waves whereas the V307L mutation promoted non-stationary and unstable re-entrant waves. Both mutations shortened tissue excitation wavelength through reduced effective refractory period but not conduction velocity, which served to increase the lifespan of re-entrant excitation in a 3D anatomical human atria model, as well as the dominant frequency (DF), which was higher for the V141M mutation. Quinidine was effective at terminating arrhythmic excitation waves associated with the V307L but not V141M mutation, and reduced the DF in a dose-dependent manner under both mutation conditions. This study provides mechanistic insights into different AF/bradycardia phenotypes in SQT2 and the efficacy of quinidine pharmacotherapy.

Description of the Human Atrial Action Potential Derived From a Single, Congruent Data Source: Novel Computational Models for Integrated Experimental-Numerical Study of Atrial Arrhythmia Mechanisms.

Introduction: The development of improved diagnosis, management, and treatment strategies for human atrial fibrillation (AF) is a significant and important challenge in order to improve quality of life for millions and reduce the substantial social-economic costs of the condition. As a complex condition demonstrating high variability and relation to other cardiac conditions, the study of AF requires approaches from multiple disciplines including single-cell experimental electrophysiology and computational modeling. Models of human atrial cells are less well parameterized than those of the human ventricle or other mammal species, largely due to the inherent challenges in patch clamping human atrial cells. Such challenges include, frequently, unphysiologically depolarized resting potentials and thus injection of a compensatory hyperpolarizing current, as well as detecting certain ion currents which may be disrupted by the cell isolation process. The aim of this study was to develop a laboratory specific model of human atrial electrophysiology which reproduces exactly the conditions of isolated-cell experiments, including testing of multiple experimental interventions. Methods: Formulations for the primary ion currents characterized by isolated-cell experiments in the Workman laboratory were fit directly to voltage-clamp data; the fast sodium-current was parameterized based on experiments relating resting membrane potential to maximal action potential upstroke velocity; compensatory hyperpolarizing current was included as a constant applied current. These formulations were integrated with three independent human atrial cell models to provide a family of novel models. Extrapolated intact-cell models were developed through removal of the hyperpolarizing current and introduction of terminal repolarization potassium currents. Results: The isolated-cell models quantitatively reproduced experimentally measured properties of excitation in both control and pharmacological and dynamic-clamp interventions. Comparison of isolated and intact-cell models highlighted the importance of reproducing this cellular environment when comparing experimental and simulation data. Conclusion: We have developed a laboratory specific model of the human atrial cell which directly reproduces the experimental isolated-cell conditions and captures human atrial excitation properties. The model may be particularly useful for directly relating model to experiment, and offers a complementary tool to the available set of human atrial cell models with specific advantages resulting from the congruent input data source.

Trigger vs. Substrate: Multi-Dimensional Modulation of QT-Prolongation Associated Arrhythmic Dynamics by a hERG Channel Activator.

Background: Prolongation of the QT interval of the electrocardiogram (ECG), underlain by prolongation of the action potential duration (APD) at the cellular level, is linked to increased vulnerability to cardiac arrhythmia. Pharmacological management of arrhythmia associated with QT prolongation is typically achieved through attempting to restore APD to control ranges, reversing the enhanced vulnerability to Ca2+-dependent afterdepolarisations (arrhythmia triggers) and increased transmural dispersion of repolarisation (arrhythmia substrate) associated with APD prolongation. However, such pharmacological modulation has been demonstrated to have limited effectiveness. Understanding the integrative functional impact of pharmacological modulation requires simultaneous investigation of both the trigger and substrate. Methods: We implemented a multi-scale (cell and tissue) in silico approach using a model of the human ventricular action potential, integrated with a model of stochastic 3D spatiotemporal Ca2+ dynamics, and parameter modification to mimic prolonged QT conditions. We used these models to examine the efficacy of the hERG activator MC-II-157c in restoring APD to control ranges, examined its effects on arrhythmia triggers and substrates, and the interaction of these arrhythmia triggers and substrates. Results: QT prolongation conditions promoted the development of spontaneous release events underlying afterdepolarisations during rapid pacing. MC-II-157c applied to prolonged QT conditions shortened the APD, inhibited the development of afterdepolarisations and reduced the probability of afterdepolarisations manifesting as triggered activity in single cells. In tissue, QT prolongation resulted in an increased transmural dispersion of repolarisation, which manifested as an increased vulnerable window for uni-directional conduction block. In some cases, MC-II-157c further increased the vulnerable window through its effects on INa. The combination of stochastic release event modulation and transmural dispersion of repolarisation modulation by MC-II-157c resulted in an integrative behavior wherein the arrhythmia trigger is reduced but the arrhythmia substrate is increased, leading to variable and non-linear overall vulnerability to arrhythmia. Conclusion: The relative balance of reduced trigger and increased substrate underlies a multi-dimensional role of MC-II-157c in modulation of cardiac arrhythmia vulnerability associated with prolonged QT interval.

A computational model of spatio-temporal cardiac intracellular calcium handling with realistic structure and spatial flux distribution from sarcoplasmic reticulum and t-tubule reconstructions.

Intracellular calcium cycling is a vital component of cardiac excitation-contraction coupling. The key structures responsible for controlling calcium dynamics are the cell membrane (comprising the surface sarcolemma and transverse-tubules), the intracellular calcium store (the sarcoplasmic reticulum), and the co-localisation of these two structures to form dyads within which calcium-induced-calcium-release occurs. The organisation of these structures tightly controls intracellular calcium dynamics. In this study, we present a computational model of intracellular calcium cycling in three-dimensions (3-D), which incorporates high resolution reconstructions of these key regulatory structures, attained through imaging of tissue taken from the sheep left ventricle using serial block face scanning electron microscopy. An approach was developed to model the sarcoplasmic reticulum structure at the whole-cell scale, by reducing its full 3-D structure to a 3-D network of one-dimensional strands. The model reproduces intracellular calcium dynamics during control pacing and reveals the high-resolution 3-D spatial structure of calcium gradients and intracellular fluxes in both the cytoplasm and sarcoplasmic reticulum. We also demonstrated the capability of the model to reproduce potentially pro-arrhythmic dynamics under perturbed conditions, pertaining to calcium-transient alternans and spontaneous release events. Comparison with idealised cell models emphasised the importance of structure in determining calcium gradients and controlling the spatial dynamics associated with calcium-transient alternans, wherein the probabilistic nature of dyad activation and recruitment was constrained. The model was further used to highlight the criticality in calcium spark propagation in relation to inter-dyad distances. The model presented provides a powerful tool for future investigation of structure-function relationships underlying physiological and pathophysiological intracellular calcium handling phenomena at the whole-cell. The approach allows for the first time direct integration of high-resolution images of 3-D intracellular structures with models of calcium cycling, presenting the possibility to directly assess the functional impact of structural remodelling at the cellular scale.

In silico assessment of genetic variation in KCNA5 reveals multiple mechanisms of human atrial arrhythmogenesis.

A recent experimental study investigating patients with lone atrial fibrillation identified six novel mutations in the KCNA5 gene. The mutants exhibited both gain- and loss-of-function of the atrial specific ultra-rapid delayed rectifier K+ current, IKur. The aim of this study is to elucidate and quantify the functional impact of these KCNA5 mutations on atrial electrical activity. A multi-scale model of the human atria was updated to incorporate detailed experimental data on IKur from both wild-type and mutants. The effects of the mutations on human atrial action potential and rate dependence were investigated at the cellular level. In tissue, we assessed the effects of the mutations on the vulnerability to unidirectional conduction patterns and dynamics of re-entrant excitation waves. Gain-of-function mutations shortened the action potential duration in single cells, and stabilised and accelerated re-entrant excitation in tissue. Loss-of-function mutations had heterogeneous effects on action potential duration and promoted early-after-depolarisations following beta-adrenergic stimulation. In the tissue model, loss-of-function mutations facilitated breakdown of excitation waves at more physiological excitation rates than the wild-type, and the generation of early-after-depolarisations promoted unidirectional patterns of excitation. Gain- and loss-of-function IKur mutations produced multiple mechanisms of atrial arrhythmogenesis, with significant differences between the two groups of mutations. This study provides new insights into understanding the mechanisms by which mutant IKur contributes to atrial arrhythmias. In addition, as IKur is an atrial-specific channel and a number of IKur-selective blockers have been developed as anti-AF agents, this study also helps to understand some contradictory results on both pro- and anti-arrhythmic effects of blocking IKur.