Interactive effect of beta-adrenergic stimulation and mechanical stretch on low-frequency oscillations of ventricular action potential duration in humans

Stretch of the papillary insertion triggers reentrant arrhythmia: an in silico patient study

Background

The electrophysiological mechanism connecting mitral valve prolapse (MVP), premature ventricular complexes and life-threatening ventricular arrhythmia is unknown. A common hypothesis is that stretch activated channels (SACs) play a significant role. SACs can trigger depolarizations or shorten repolarization times in response to myocardial stretch. Through these mechanisms, pathological traction of the papillary muscle (PM), as has been observed in patients with MVP, may induce irregular electrical activity and result in reentrant arrhythmia.

Methods

Based on a patient with MVP and mitral annulus disjunction, we modeled the effect of excessive PM traction in a detailed medical image-derived ventricular model by activating SACs in the PM insertion region. By systematically varying the onset of SAC activation following sinus pacing, we identified vulnerability windows for reentry with 1 ms resolution. We explored how reentry was affected by the SAC reversal potential ( E SAC ) and the size of the region with simulated stretch (SAC region). Finally, the effect of global or focal fibrosis, modeled as reduction in tissue conductivity or mesh splitting (fibrotic microstructure), was investigated.

Results

In models with healthy tissue or fibrosis modeled solely as CV slowing, we observed two vulnerable periods of reentry: ForE SAC of -10 and -30 mV, SAC activated during the T-wave could cause depolarization of the SAC region which lead to reentry. ForE SAC of -40 and -70 mV, SAC activated during the QRS complex could result in early repolarization of the SAC region and subsequent reentry. In models with fibrotic microstructure in the SAC region, we observed micro-reentries and a larger variability in which times of SAC activation triggered reentry. In these models, 86% of reentries were triggered during the QRS complex or T-wave. We only observed reentry for sufficiently large SAC regions ( > = 8 mm radius in models with healthy tissue).

Conclusion

Stretch of the PM insertion region following sinus activation may initiate ventricular reentry in patients with MVP, with or without fibrosis. Depending on the SAC reversal potential and timing of stretch, reentry may be triggered by ectopy due to SAC-induced depolarizations or by early repolarization within the SAC region.

Fiber-Optic-Based System for High-Resolution Monitoring of Stretch in Excised Tissues

Cardiovascular diseases cause a high number of deaths nowadays. To improve these statistics, new strategies to better understand the electrical and mechanical abnormalities underlying them are urgently required. This study focuses on the development of a sensor to measure tissue stretch in excised tissues, enabling improved knowledge of biomechanical properties and allowing greater control in real time. A system made of biocompatible materials is described, which is based on two cantilevered platforms that integrate an optical fiber inside them to quantify the amount of stretch the tissues are exposed to with a precision of μm. The operating principle of the sensor is based on the variation of the optical path with the movement of the platforms onto which the samples are fixed. The conducted tests highlight that this system, based on a simple topology and technology, is capable of achieving the desired purpose (a resolution of ∼1 μm), enabling the tissue to be bathed in any medium within the system.

Steady-state and transient effects of SK channel block and adrenergic stimulation to counteract acetylcholine-induced arrhythmogenic effects in the human atria: A computational study

Hyperactivity of the parasympathetic nervous system has been linked to the development of paroxysmal atrial fibrillation (AF). The parasympathetic neurotransmitter acetylcholine (ACh) causes a reduction in action potential (AP) duration (APD) and an increase in resting membrane potential (RMP), both of which contribute to enhance the risk for reentry. Research suggests that small-conductance calcium activated potassium (SK) channels may be an effective target for treating AF. Therapies targeting the autonomic nervous system, either alone or in combination with other drugs, have been explored and have been shown to decrease the incidence of atrial arrhythmias. This study uses computational modeling and simulation to examine the impact of SK channel block (SKb) and β-adrenergic stimulation through Isoproterenol (Iso) on countering the negative effects of cholinergic activity in human atrial cell and 2D tissue models. The steady-state effects of Iso and/or SKb on AP shape, APD at 90% repolarization (APD₉₀) and RMP were evaluated. The ability to terminate stable rotational activity in cholinergically-stimulated 2D tissue models of AF was also investigated. A range of SKb and Iso application kinetics, which reflect varying drug binding rates, were taken into consideration. The results showed that SKb alone prolonged APD₉₀ and was able to stop sustained rotors in the presence of ACh concentrations up to 0.01 μM. Iso terminated rotors under all tested ACh concentrations, but resulted in highly-variable steady-state outcomes depending on baseline AP morphology. Importantly, the combination of SKb and Iso resulted in greater APD₉₀ prolongation and showed promising anti-arrhythmic potential by stopping stable rotors and preventing re-inducibility.

Insights from an electro-mechanical heart failure cell model: Role of SERCA enhancement on arrhythmogenesis and myocyte contraction

BACKGROUND AND OBJECTIVE

Structural and electrical remodeling in heart failure predisposes the heart to ventricular arrhythmias. Computer modeling approaches, used to complement experimental results, can provide a more mechanistic knowledge of the biophysical phenomena underlying cardiac pathologies. Indeed, previous in-silico studies have improved the understanding of the electrical correlates of heart failure involved in arrhythmogenesis; however, information on the crosstalk between electrical activity, intracellular Ca²⁺ and contraction is still incomplete. This study aims to investigate the electro-mechanical behavior of virtual failing human ventricular myocytes to help in the development of therapies, which should ideally target pump failure and arrhythmias at the same time.

METHODS

We implemented characteristic remodeling of heart failure with reduced ejection fraction by including reported changes in ionic conductances, sarcomere function and cell structure (e.g. T-tubules disarray). Model parametrization was based on published experimental data and the outcome of simulations was validated against experimentally observed patterns. We focused on two aspects of myocardial dysfunction central in heart failure: altered force-frequency relationship and susceptibility to arrhythmogenic early afterdepolarizations. Because biological variability is a major problem in the generalization of in-silico findings based on a unique set of model parameters, we generated and evaluated a population of models.

RESULTS

The population-based approach is crucial in robust identification of parameters at the core of abnormalities and in generalizing the outcome of their correction. As compared to non-failing ones, failing myocytes had prolonged repolarization, a higher incidence of early afterdepolarizations, reduced contraction and a shallower force-frequency relationship, all features peculiar of heart failure. Component analysis applied to the model population identified reduced SERCA function as a relevant contributor to most of these derangements, which were largely reverted or diminished by restoration of SERCA function alone.

CONCLUSIONS

These simulated results encourage the development of strategies comprising SERCA stimulation and highlight the need to evaluate both electrical and mechanical outcomes.

The role of β-adrenergic stimulation in QT interval adaptation to heart rate during stress test

The adaptation lag of the QT interval after heart rate (HR) has been proposed as an arrhythmic risk marker. Most studies have quantified the QT adaptation lag in response to abrupt, step-like changes in HR induced by atrial pacing, in response to tilt test or during ambulatory recordings. Recent studies have introduced novel methods to quantify the QT adaptation lag to gradual, ramp-like HR changes in stress tests by evaluating the differences between the measured QT series and an estimated, memoryless QT series obtained from the instantaneous HR. These studies have observed the QT adaptation lag to progressively reduce when approaching the stress peak, with the underlying mechanisms being still unclear. This study analyzes the contribution of β-adrenergic stimulation to QT interval rate adaptation in response to gradual, ramp-like HR changes. We first quantify the QT adaptation lag in Coronary Artery Disease (CAD) patients undergoing stress test. To uncover the involved mechanisms, we use biophysically detailed computational models coupling descriptions of human ventricular electrophysiology and β-adrenergic signaling, from which we simulate ventricular action potentials and ECG signals. We characterize the adaptation of the simulated QT interval in response to the HR time series measured from each of the analyzed CAD patients. We show that, when the simulated ventricular tissue is subjected to a time-varying β-adrenergic stimulation pattern, with higher stimulation levels close to the stress peak, the simulated QT interval presents adaptation lags during exercise that are more similar to those measured from the patients than when subjected to constant β-adrenergic stimulation. During stress test recovery, constant and time-varying β-adrenergic stimulation patterns render similar adaptation lags, which are generally shorter than during exercise, in agreement with results from the patients. In conclusion, our findings support the role of time-varying β-adrenergic stimulation in contributing to QT interval adaptation to gradually increasing HR changes as those seen during the exercise phase of a stress test.

Periodic repolarization dynamics as predictor of risk for sudden cardiac death in chronic heart failure patients

The two most common modes of death among chronic heart failure (CHF) patients are sudden cardiac death (SCD) and pump failure death (PFD). Periodic repolarization dynamics (PRD) quantifies low-frequency oscillations in the T wave vector of the electrocardiogram (ECG) and has been postulated to reflect sympathetic modulation of ventricular repolarization. This study aims to evaluate the prognostic value of PRD to predict SCD and PFD in a population of CHF patients. 20-min high-resolution (1000 Hz) ECG recordings from 569 CHF patients were analyzed. Patients were divided into two groups, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PRD}^+$$\end{document}PRD+ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PRD}^-$$\end{document}PRD-, corresponding to PRD values above and below the optimum cutoff point of PRD in the study population. Univariate Cox regression analysis showed that SCD risk in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PRD}^+$$\end{document}PRD+ group was double the risk in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {PRD}^-$$\end{document}PRD- group [hazard ratio (95% CI) 2.001 (1.127–3.554), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {p}<0.05$$\end{document}p<0.05]. The combination of PRD with other Holter-based ECG indices, such as turbulence slope (TS) and index of average alternans (IAA), improved SCD prediction by identifying groups of patients at high SCD risk. PFD could be predicted by PRD only when combined with TS [hazard ratio 2.758 (1.572–4.838), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {p}<0.001$$\end{document}p<0.001]. In conclusion, the combination of PRD with IAA and TS can be used to stratify the risk for SCD and PFD, respectively, in CHF patients.

Electrophysiological Remodeling: Cardiac T-Tubules and ß-Adrenoceptors

Beta-adrenoceptors (βAR) are often viewed as archetypal G-protein coupled receptors. Over the past fifteen years, investigations in cardiovascular biology have provided remarkable insights into this receptor family. These studies have shifted pharmacological dogma, from one which centralized the receptor to a new focus on structural micro-domains such as caveolae and t-tubules. Important studies have examined, separately, the structural compartmentation of ion channels and βAR. Despite links being assumed, relatively few studies have specifically examined the direct link between structural remodeling and electrical remodeling with a focus on βAR. In this review, we will examine the nature of receptor and ion channel dysfunction on a substrate of cardiomyocyte microdomain remodeling, as well as the likely ramifications for cardiac electrophysiology. We will then discuss the advances in methodologies in this area with a specific focus on super-resolution microscopy, fluorescent imaging, and new approaches involving microdomain specific, polymer-based agonists. The advent of powerful computational modelling approaches has allowed the science to shift from purely empirical work, and may allow future investigations based on prediction. Issues such as the cross-reactivity of receptors and cellular heterogeneity will also be discussed. Finally, we will speculate as to the potential developments within this field over the next ten years.

ECG Ventricular Repolarization Dynamics during Exercise: Temporal Profile, Relation to Heart Rate Variability and Effects of Age and Physical Health

Periodic repolarization dynamics (PRD) is a novel electrocardiographic marker of cardiac repolarization instability with powerful risk stratification capacity for total mortality and sudden cardiac death. Here, we use a time-frequency analysis approach to continuously quantify PRD at rest and during exercise, assess its dependence on heart rate variability (HRV) and characterize the effects of age (young adults/middle-aged adults/older adults), body mass index (non-overweight/overweight) and cardiorespiratory fitness level (fit/unfit). Sixty-six male volunteers performed an exercise test. RR and dT variabilities (RRV, dTV), as well as the fraction of dT variability unrelated to RR variability, were computed based on time-frequency representations. The instantaneous LF power of dT (PdTV), representing the same concept as PRD, and of its RRV-unrelated component (PdTVuRRV) were quantified. dT angle was found to mostly oscillate in the LF band. Overall, 50-70% of PdTV was linearly unrelated to RRV. The onset of exercise caused a sudden increase in PdTV and PdTVuRRV, which returned to pre-exercise levels during recovery. Clustering analysis identified a group of overweight and unfit individuals with significantly higher PdTV and PdTVuRRV values at rest than the rest of the population. Our findings shed new light on the temporal profile of PRD during exercise, its relationship to HRV and the differences in PRD between subjects according to phenotypic characteristics.

Emerging evidence for a mechanistic link between low-frequency oscillation of ventricular repolarization measured from the electrocardiogram T-wave vector and arrhythmia

Strong recent clinical evidence links the presence of prominent oscillations of ventricular repolarization in the low-frequency range (0.04-0.15 Hz) to the incidence of ventricular arrhythmia and sudden death in post-MI patients and patients with ischaemic and non-ischaemic cardiomyopathy. It has been proposed that these oscillations reflect oscillations of ventricular action potential duration at the sympathetic nerve frequency. Here we review emerging evidence to support that contention and provide insight into possible underlying mechanisms for this association.

Electro-Mechanical Whole-Heart Digital Twins: A Fully Coupled Multi-Physics Approach

Mathematical models of the human heart are evolving to become a cornerstone of precision medicine and support clinical decision making by providing a powerful tool to understand the mechanisms underlying pathophysiological conditions. In this study, we present a detailed mathematical description of a fully coupled multi-scale model of the human heart, including electrophysiology, mechanics, and a closed-loop model of circulation. State-of-the-art models based on human physiology are used to describe membrane kinetics, excitation-contraction coupling and active tension generation in the atria and the ventricles. Furthermore, we highlight ways to adapt this framework to patient specific measurements to build digital twins. The validity of the model is demonstrated through simulations on a personalized whole heart geometry based on magnetic resonance imaging data of a healthy volunteer. Additionally, the fully coupled model was employed to evaluate the effects of a typical atrial ablation scar on the cardiovascular system. With this work, we provide an adaptable multi-scale model that allows a comprehensive personalization from ion channels to the organ level enabling digital twin modeling.

Characterization of Spatio-Temporal Cardiac Action Potential Variability at Baseline and Under β-Adrenergic Stimulation by Combined Unscented Kalman Filter and Double Greedy Dimension Reduction

OBJECTIVE

Elevated spatio-temporal variability of human ventricular repolarization has been related to increased risk for ventricular arrhythmias and sudden cardiac death, particularly under β-adrenergic stimulation ( β-AS). This work presents a methodology for theoretical characterization of temporal and spatial repolarization variability at baseline conditions and in response to β-AS. For any measured voltage trace, the proposed methodology estimates the parameters and state variables of an underlying human ventricular action potential (AP) model by combining Double Greedy Dimension Reduction (DGDR) with automatic selection of biomarkers and the Unscented Kalman Filter (UKF). Such theoretical characterization can facilitate subsequent characterization of underlying variability mechanisms.

MATERIAL AND METHODS

Given an AP trace, initial estimates for the ionic conductances in a stochastic version of the baseline human ventricular O'Hara et al. model were obtained by DGDR. Those estimates served to initialize and update model parameter estimates by the UKF method based on formulation of an associated nonlinear state-space representation and joint estimation of model parameters and state variables. Similarly, β-AS-induced phosphorylation levels of cellular substrates were estimated by the DGDR-UKF methodology. Performance was tested by building an experimentally-calibrated population of virtual cells, from which synthetic AP traces were generated for baseline and β-AS conditions.

RESULTS

The combined DGDR-UKF methodology led to 25% reduction in the error associated with estimation of ionic current conductances at baseline conditions and phosphorylation levels under β-AS with respect to individual DGDR and UKF methods. This improvement was not at the expense of higher computational load, which was diminished by 90% with respect to the individual UKF method. Both temporal and spatial AP variability of repolarization were accurately characterized by the DGDR-UKF methodology.

CONCLUSIONS

A combined DGDR-UKF methodology is proposed for parameter and state variable estimation of human ventricular cell models from available AP traces at baseline and under β-AS. This methodology improves the estimation performance and reduces the convergence time with respect to individual DGDR and UKF methods and renders a suitable approach for computational characterization of spatio-temporal repolarization variability to be used for ascertainment of variability mechanisms and its relation to arrhythmogenesis.

Minimally invasive system to reliably characterize ventricular electrophysiology from living donors

Cardiac tissue slices preserve the heterogeneous structure and multicellularity of the myocardium and allow its functional characterization. However, access to human ventricular samples is scarce. We aim to demonstrate that slices from small transmural core biopsies collected from living donors during routine cardiac surgery preserve structural and functional properties of larger myocardial specimens, allowing accurate electrophysiological characterization. In pigs, we compared left ventricular transmural core biopsies with transmural tissue blocks from the same ventricular region. In humans, we analyzed transmural biopsies and papillary muscles from living donors. All tissues were vibratome-sliced. By histological analysis of the transmural biopsies, we showed that tissue architecture and cellular organization were preserved. Enzymatic and vital staining methods verified viability. Optically mapped transmembrane potentials confirmed that action potential duration and morphology were similar in pig biopsies and tissue blocks. Action potential morphology and duration in human biopsies and papillary muscles agreed with published ranges. In both pigs and humans, responses to increasing pacing frequencies and β-adrenergic stimulation were similar in transmural biopsies and larger tissues. We show that it is possible to successfully collect and characterize tissue slices from human myocardial biopsies routinely extracted from living donors, whose behavior mimics that of larger myocardial preparations both structurally and electrophysiologically.

Estimation of the second ventilatory threshold through ventricular repolarization profile analysis

Under the hypothesis that sympathetic control of ventricular repolarization may change once the second ventilatory threshold (VT2) has been reached, a novel methodology for non-invasive VT2 estimation based on the analysis of the T wave from the electrocardiogram (ECG) is proposed, and potential underlying physiological mechanisms are suggested. 25 volunteers (33.4 ± 5.2 years) underwent an incremental power cycle ergometer test (25 W/minute). During the test, respiratory gas exchange and multi-lead ECG were acquired. The former was employed to determine VT2, used here as a reference, whereas the latter was used to compute the temporal profiles of an index of ventricular repolarization instability (dT) and its low-frequency (LF) oscillations (LFdT). The sudden increases observed in dT and LFdT profiles above an established heart rate threshold were employed to derive VT2 estimates, referred to as VT2_{d T} and VT2_{LF d T} , respectively. Estimation errors of -4.7 ± 25.2 W were obtained when considering VT2_{d T} . Errors were lower than the one-minute power increment of 25 W in 68% of the subjects and lower than 50 W in 89.5% of them. When using VT2_{LF d T} , estimation error was of 15.3 ± 32.4 W. Most of the subjects shared common characteristic dT and LFdT profiles, which could be reflecting changes in the autonomic control of ventricular repolarization before and after reaching VT2. The analysis of ventricular repolarization dynamics during exercise allows non-invasive ECG-based estimation of VT2, possibly in relation to changes in the autonomic control of ventricular electrical activity when VT2 is reached.

Cardiac Mechano-Electric Coupling: Acute Effects of Mechanical Stimulation on Heart Rate and Rhythm

The heart is vital for biological function in almost all chordates, including humans. It beats continually throughout our life, supplying the body with oxygen and nutrients while removing waste products. If it stops, so does life. The heartbeat involves precise coordination of the activity of billions of individual cells, as well as their swift and well-coordinated adaption to changes in physiological demand. Much of the vital control of cardiac function occurs at the level of individual cardiac muscle cells, including acute beat-by-beat feedback from the local mechanical environment to electrical activity (as opposed to longer term changes in gene expression and functional or structural remodeling). This process is known as mechano-electric coupling (MEC). In the current review, we present evidence for, and implications of, MEC in health and disease in human; summarize our understanding of MEC effects gained from whole animal, organ, tissue, and cell studies; identify potential molecular mediators of MEC responses; and demonstrate the power of computational modeling in developing a more comprehensive understanding of ‟what makes the heart tick.ˮ.

Low-Frequency Oscillations in Cardiac Sympathetic Neuronal Activity

Sudden cardiac death caused by ventricular arrhythmias is among the leading causes of mortality, with approximately half of all deaths attributed to heart disease worldwide. Periodic repolarization dynamics (PRD) is a novel marker of repolarization instability and strong predictor of death in patients post-myocardial infarction that is believed to occur in association with low-frequency oscillations in sympathetic nerve activity. However, this hypothesis is based on associations of PRD with indices of sympathetic activity that are not directly linked to cardiac function, such as muscle vasoconstrictor activity and the variability of cardiovascular autospectra. In this review article, we critically evaluate existing scientific evidence obtained primarily in experimental animal models, with the aim of identifying the neuronal networks responsible for the generation of low-frequency sympathetic rhythms along the neurocardiac axis. We discuss the functional significance of rhythmic sympathetic activity on neurotransmission efficacy and explore its role in the pathogenesis of ventricular repolarization instability. Most importantly, we discuss important gaps in our knowledge that require further investigation in order to confirm the hypothesis that low frequency cardiac sympathetic oscillations play a causative role in the generation of PRD.

Sensitivity analysis of a strongly-coupled human-based electromechanical cardiac model: Effect of mechanical parameters on physiologically relevant biomarkers

The human heart beats as a result of multiscale nonlinear dynamics coupling subcellular to whole organ processes, achieving electrophysiologically-driven mechanical contraction. Computational cardiac modelling and simulation have achieved a great degree of maturity, both in terms of mathematical models of underlying biophysical processes and the development of simulation software. In this study, we present the detailed description of a human-based physiologically-based, and fully-coupled ventricular electromechanical modelling and simulation framework, and a sensitivity analysis focused on its mechanical properties. The biophysical detail of the model, from ionic to whole-organ, is crucial to enable future simulations of disease and drug action. Key novelties include the coupling of state-of-the-art human-based electrophysiology membrane kinetics, excitation-contraction and active contraction models, and the incorporation of a pre-stress model to allow for pre-stressing and pre-loading the ventricles in a dynamical regime. Through high performance computing simulations, we demonstrate that 50% to 200% - 1000% variations in key parameters result in changes in clinically-relevant mechanical biomarkers ranging from diseased to healthy values in clinical studies. Furthermore mechanical biomarkers are primarily affected by only one or two parameters. Specifically, ejection fraction is dominated by the scaling parameter of the active tension model and its scaling parameter in the normal direction ( k ort 2 ); the end systolic pressure is dominated by the pressure at which the ejection phase is triggered ( P ej ) and the compliance of the Windkessel fluid model ( C ); and the longitudinal fractional shortening is dominated by the fibre angle ( ϕ ) and k ort 2 . The wall thickening does not seem to be clearly dominated by any of the considered input parameters. In summary, this study presents in detail the description and implementation of a human-based coupled electromechanical modelling and simulation framework, and a high performance computing study on the sensitivity of mechanical biomarkers to key model parameters. The tools and knowledge generated enable future investigations into disease and drug action on human ventricles.

Cardiomyocyte calcium handling in health and disease: Insights from in vitro and in silico studies

Calcium (Ca²⁺) plays a central role in cardiomyocyte excitation-contraction coupling. To ensure an optimal electrical impulse propagation and cardiac contraction, Ca²⁺ levels are regulated by a variety of Ca²⁺-handling proteins. In turn, Ca²⁺ modulates numerous electrophysiological processes. Accordingly, Ca²⁺-handling abnormalities can promote cardiac arrhythmias via various mechanisms, including the promotion of afterdepolarizations, ion-channel modulation and structural remodeling. In the last 30 years, significant improvements have been made in the computational modeling of cardiomyocyte Ca²⁺ handling under physiological and pathological conditions. However, numerous questions involving the Ca²⁺-dependent regulation of different macromolecular complexes, cross-talk between Ca²⁺-dependent regulatory pathways operating over a wide range of time scales, and bidirectional interactions between electrophysiology and mechanics remain to be addressed by in vitro and in silico studies. A better understanding of disease-specific Ca²⁺-dependent proarrhythmic mechanisms may facilitate the development of improved therapeutic strategies. In this review, we describe the fundamental mechanisms of cardiomyocyte Ca²⁺ handling in health and disease, and provide an overview of currently available computational models for cardiomyocyte Ca²⁺ handling. Finally, we discuss important uncertainties and open questions about cardiomyocyte Ca²⁺ handling and highlight how synergy between in vitro and in silico studies may help to answer several of these issues.

Complex Interaction Between Low-Frequency APD Oscillations and Beat-to-Beat APD Variability in Humans Is Governed by the Sympathetic Nervous System

BACKGROUND

Recent clinical, experimental and modeling studies link oscillations of ventricular repolarization in the low frequency (LF) (approx. 0.1 Hz) to arrhythmogenesis. Sympathetic provocation has been shown to enhance both LF oscillations of action potential duration (APD) and beat-to-beat variability (BVR) in humans. We hypothesized that beta-adrenergic blockade would reduce LF oscillations of APD and BVR of APD in humans and that the two processes might be linked.

METHODS AND RESULTS

Twelve patients with normal ventricles were studied during routine electrophysiological procedures. Activation-recovery intervals (ARI) as a conventional surrogate for APD were recorded from 10 left and 10 right ventricular endocardial sites before and after acute beta-adrenergic adrenergic blockade. Cycle length was maintained constant with right ventricular pacing. Oscillatory behavior of ARI was quantified by spectral analysis and BVR as the short-term variability. Beta-adrenergic blockade reduced LF ARI oscillations (8.6 ± 4.5 ms2 vs. 5.5 ± 3.5 ms2, p = 0.027). A significant correlation was present between the initial control values and reduction seen following beta-adrenergic blockade in LF ARI (r s = 0.62, p = 0.037) such that when initial values are high the effect is greater. A similar relationship was also seen in the beat-to beat variability of ARI (r s = 0.74, p = 0.008). There was a significant correlation between the beta-adrenergic blockade induced reduction in LF power of ARI and the witnessed reduction of beat-to-beat variability of ARI (r s = 0.74, p = 0.01). These clinical results accord with recent computational modeling studies which provide mechanistic insight into the interactions of LF oscillations and beat-to-beat variability of APD at the cellular level.

CONCLUSION

Beta-adrenergic blockade reduces LF oscillatory behavior of APD (ARI) in humans in vivo. Our results support the importance of LF oscillations in modulating the response of BVR to beta-adrenergic blockers, suggesting that LF oscillations may play role in modulating beta-adrenergic mechanisms underlying BVR.

Time Course of Low-Frequency Oscillatory Behavior in Human Ventricular Repolarization Following Enhanced Sympathetic Activity and Relation to Arrhythmogenesis

Background and Objectives: Recent studies in humans and dogs have shown that ventricular repolarization exhibits a low-frequency (LF) oscillatory pattern following enhanced sympathetic activity, which has been related to arrhythmic risk. The appearance of LF oscillations in ventricular repolarization is, however, not immediate, but it may take up to some minutes. This study seeks to characterize the time course of the action potential (AP) duration (APD) oscillatory behavior in response to sympathetic provocations, unveil its underlying mechanisms and establish a potential link to arrhythmogenesis under disease conditions.

Materials and Methods: A representative set of human ventricular computational models coupling cellular electrophysiology, calcium dynamics, β-adrenergic signaling, and mechanics was built. Sympathetic provocation was modeled via phasic changes in β-adrenergic stimulation (β-AS) and mechanical stretch at Mayer wave frequencies within the 0.03–0.15 Hz band.

Results: Our results show that there are large inter-individual differences in the time lapse for the development of LF oscillations in APD following sympathetic provocation, with some cells requiring just a few seconds and other cells needing more than 3 min. Whereas, the oscillatory response to phasic mechanical stretch is almost immediate, the response to β-AS is much more prolonged, in line with experimentally reported evidences, thus being this component the one driving the slow development of APD oscillations following enhanced sympathetic activity. If β-adrenoceptors are priorly stimulated, the time for APD oscillations to become apparent is remarkably reduced, with the oscillation time lapse being an exponential function of the pre-stimulation level. The major mechanism underlying the delay in APD oscillations appearance is related to the slow I_Ks phosphorylation kinetics, with its relevance being modulated by the I_Ks conductance of each individual cell. Cells presenting short oscillation time lapses are commonly associated with large APD oscillation magnitudes, which facilitate the occurrence of pro-arrhythmic events under disease conditions involving calcium overload and reduced repolarization reserve.

Conclusions: The time course of LF oscillatory behavior of APD in response to increased sympathetic activity presents high inter-individual variability, which is associated with different expression and PKA phosphorylation kinetics of the I_Ks current. Short time lapses in the development of APD oscillations are associated with large oscillatory magnitudes and pro-arrhythmic risk under disease conditions.

Dynamic Changes of Cardiac Repolarization Instability during Exercise Testing

INTRODUCTION

Physical exercise triggers efferent cardiac sympathetic activation. Here, we tracked the spatiotemporal properties of cardiac repolarization on a beat-to-beat basis throughout a standardized exercise test and hypothesized a detectable change at the point of the anaerobic threshold (AT).

METHODS

The study included 20 healthy adults (age 35.3 ± 6.7 yr) undergoing a standardized incremental exercise test on a cycle ergometer. During exercise testing, high-resolution (2000 Hz) ECG monitoring in Frank lead configuration was performed. Three-dimensional beat-to-beat repolarization instability (dT°) was assessed by a novel vector-based method according to a previously published technology. In parallel, the lactate threshold (LT) was detected according to Dickhuth and Mader.

RESULTS

We could identify a characteristic pattern of dT° signal during exercise testing. With increasing physical activity, dT° increased concordantly to heart rate. At an average of 164 ± 38 W, dT° and heart rate abruptly showed a discordant behavior, characterized by a transient drop of dT°. The maximal discordance between dT° and heart rate was defined as ATdT° and highly significantly correlated with LTDickhuth (r = 0.841, P < 0.001) and LTMader (r = 0.819, P < 0.001), which were at 156 ± 39 and 172 ± 46 W, respectively. The characteristic of dT° could not be provoked by fast atrial pacing in the absence of exercise.

CONCLUSIONS

Repolarization instability shows a characteristic pattern during standardized exercise in healthy individuals that allows for a noninvasive estimation of AT.

Long-Term Microgravity Exposure Increases ECG Repolarization Instability Manifested by Low-Frequency Oscillations of T-Wave Vector

Ventricular arrhythmias and sudden cardiac death during long-term space missions are a major concern for space agencies. Long-duration spaceflight and its ground-based analog head-down bed rest (HDBR) have been reported to markedly alter autonomic and cardiac functioning, particularly affecting ventricular repolarization of the electrocardiogram (ECG). In this study, novel methods are developed, departing from previously published methodologies, to quantify the index of Periodic Repolarization Dynamics (PRD), an arrhythmic risk marker that characterizes sympathetically-mediated low-frequency oscillations in the T-wave vector. PRD is evaluated in ECGs from 42 volunteers at rest and during an orthostatic tilt table test recorded before and after 60-day –6° HDBR. Our results indicate that tilt test, on top of enhancing sympathetic regulation of heart rate, notably increases PRD, both before and after HDBR, thus supporting previous evidence on PRD being an indicator of sympathetic modulation of ventricular repolarization. Importantly, long-term microgravity exposure is shown to lead to significant increases in PRD, both when evaluated at rest and, even more notably, in response to tilt test. The extent of microgravity-induced changes in PRD has been associated with arrhythmic risk in prior studies. An exercise-based, but not a nutrition-based, countermeasure is able to partially reverse microgravity-induced effects on PRD. In conclusion, long-term exposure to microgravity conditions leads to elevated low-frequency oscillations of ventricular repolarization, which are potentiated following sympathetic stimulation and are related to increased risk for repolarization instabilities and arrhythmias. Tested countermeasures are only partially effective in counteracting microgravity effects.

The Effect of Emotional Valence on Ventricular Repolarization Dynamics Is Mediated by Heart Rate Variability: A Study of QT Variability and Music-Induced Emotions

Background

Emotions can affect cardiac activity, but their impact on ventricular repolarization variability, an important parameter providing information about cardiac risk and autonomic nervous system activity, is unknown. The beat-to-beat variability of the QT interval (QTV) from the body surface ECG is a non-invasive marker of repolarization variability, which can be decomposed into QTV related to RR variability (QTVrRRV) and QTV unrelated to RRV (QTVuRRV), with the latter thought to be a marker of intrinsic repolarization variability.

Aim

To determine the effect of emotional valence (pleasant and unpleasant) on repolarization variability in healthy volunteers by means of QTV analysis.

Methods

75 individuals (24.5 ± 3.2 years, 36 females) without a history of cardiovascular disease listened to music-excerpts that were either felt as pleasant (n = 6) or unpleasant (n = 6). Excerpts lasted about 90 s and were presented in a random order along with silent intervals (n = 6). QTV and RRV were derived from the ECG and the time-frequency spectrum of RRV, QTV, QTVuRRV and QTVrRRV as well as time-frequency coherence between QTV and RRV were estimated. Analysis was performed in low-frequency (LF), high frequency (HF) and total spectral bands.

Results

The heart rate-corrected QTV showed a small but significant increase from silence (median 347/interquartile range 31 ms) to listening to music felt as unpleasant (351/30 ms) and pleasant (355/32 ms). The dynamic response of QTV to emotional valence showed a transient phase lasting about 20 s after the onset of each musical excerpt. QTV and RRV were highly correlated in both HF and LF (mean coherence ranging 0.76–0.85). QTV and QTVrRRV decreased during listening to music felt as pleasant and unpleasant with respect to silence and further decreased during listening to music felt as pleasant. QTVuRRV was small and not affected by emotional valence.

Conclusion

Emotional valence, as evoked by music, has a small but significant effect on QTV and QTVrRRV, but not on QTVuRRV. This suggests that the interaction between emotional valence and ventricular repolarization variability is mediated by cycle length dynamics and not due to intrinsic repolarization variability.

Prediction of mortality benefit based on periodic repolarisation dynamics in patients undergoing prophylactic implantation of a defibrillator: a prospective, controlled, multicentre cohort study

BACKGROUND

A small proportion of patients undergoing primary prophylactic implantation of implantable cardioverter defibrillators (ICDs) experiences malignant arrhythmias. We postulated that periodic repolarisation dynamics, a novel marker of sympathetic-activity-associated repolarisation instability, could be used to identify electrically vulnerable patients who would benefit from prophylactic implantation of ICDs by way of a reduction in mortality.

METHODS

We did a prespecified substudy of EUropean Comparative Effectiveness Research to Assess the Use of Primary ProphylacTic Implantable Cardioverter Defibrillators (EU-CERT-ICD), a prospective, investigator-initiated, non-randomised, controlled cohort study done at 44 centres in 15 EU countries. Patients aged 18 years or older with ischaemic or non-ischaemic cardiomyopathy and reduced left ventricular ejection fraction (≤35%) were eligible for inclusion if they met guideline-based criteria for primary prophylactic implantation of ICDs. Periodic repolarisation dynamics from 24-h Holter recordings were assessed blindly in patients the day before ICD implantation or on the day of study enrolment in patients who were conservatively managed. The primary endpoint was all-cause mortality. Propensity scoring and multivariable models were used to assess the interaction between periodic repolarisation dynamics and the treatment effect of ICDs on mortality.

FINDINGS

Between May 12, 2014, and Sept 7, 2018, 1371 patients were enrolled in our study. 968 of these patients underwent ICD implantation, and 403 were treated conservatively. During follow-up (median 2·7 years [IQR 2·0-3·3] in the ICD group and 1·2 years [0·8-2·7] in the control group), 138 (14%) patients died in the ICD group and 64 (16%) patients died in the control group. We noted a 43% reduction in mortality in the ICD group compared with the control group (adjusted hazard ratio [HR] 0·57 [95% CI 0·41-0·79]; p=0·0008). Periodic repolarisation dynamics significantly predicted the treatment effect of ICDs on mortality (adjusted p=0·0307). The mortality benefits associated with ICD implantation were greater in patients with periodic repolarisation dynamics of 7·5 deg or higher (n=199; adjusted HR 0·25 [95% CI 0·13-0·47] for the ICD group vs the control group; p<0·0001) than in those with periodic repolarisation dynamics less than 7·5 deg (n=1166; adjusted HR 0·69 [95% CI 0·47-1·00]; p=0·0492; p_interaction=0·0056). The number needed to treat was 18·3 (95% CI 10·6-4895·3) in patients with periodic repolarisation dynamics less than 7·5 deg and 3·1 (2·6-4·8) in those with periodic repolarisation dynamics of 7·5 deg or higher.

INTERPRETATION

Periodic repolarisation dynamics predict mortality reductions associated with prophylactic implantation of ICDs in contemporarily treated patients with ischaemic or non-ischaemic cardiomyopathy. Periodic repolarisation dynamics could help to guide treatment decisions about prophylactic ICD implantation.

FUNDING

The European Community's 7th Framework Programme.

Pro-Arrhythmic Ventricular Remodeling Is Associated With Increased Respiratory and Low-Frequency Oscillations of Monophasic Action Potential Duration in the Chronic Atrioventricular Block Dog Model

In addition to beat-to-beat fluctuations, action potential duration (APD) oscillates at (1) a respiratory frequency and (2) a low frequency (LF) (<0.1 Hz), probably caused by bursts of sympathetic nervous system discharge. This study investigates whether ventricular remodeling in the chronic AV block (CAVB) dog alters these oscillations of APD and whether this has consequences for arrhythmogenesis. We performed a retrospective analysis of 39 dog experiments in sinus rhythm (SR), acute AV block (AAVB), and after 2 weeks of chronic AV block. Spectral analysis of left ventricular monophasic action potential duration (LV MAPD) was done to quantify respiratory frequency (RF) power and LF power. Dofetilide (0.025 mg/kg in 5 min) was infused to test for inducibility of Torsade de Pointes (TdP) arrhythmias. RF power was significantly increased at CAVB compared to AAVB and SR (log[RF] of -1.13 ± 1.62 at CAVB vs. log[RF] of -2.82 ± 1.24 and -3.29 ± 1.29 at SR and AAVB, respectively, p < 0.001). LF power was already significantly increased at AAVB and increased even further at CAVB (-3.91 ± 0.70 at SR vs. -2.52 ± 0.85 at AAVB and -1.14 ± 1.62 at CAVB, p < 0.001). In addition, LF power was significantly larger in inducible CAVB dogs (log[LF] -0.6 ± 1.54 in inducible dogs vs. -2.56 ± 0.43 in non-inducible dogs, p < 0.001). In conclusion, ventricular remodeling in the CAVB dog results in augmentation of respiratory and low-frequency (LF) oscillations of LV MAPD. Furthermore, TdP-inducible CAVB dogs show increased LF power.

Mechanisms Underlying Interactions Between Low-Frequency Oscillations and Beat-to-Beat Variability of Celullar Ventricular Repolarization in Response to Sympathetic Stimulation: Implications for Arrhythmogenesis

Background and Objectives: Enhanced beat-to-beat variability of ventricular repolarization (BVR) has been linked to arrhythmias and sudden cardiac death. Recent experimental studies on human left ventricular epicardial electrograms have shown that BVR closely interacts with low-frequency (LF) oscillations of activation recovery interval during sympathetic provocation. In this work human ventricular computational cell models are developed to reproduce the experimentally observed interactions between BVR and its LF oscillations, to assess underlying mechanisms and to establish a relationship with arrhythmic risk.

Materials and Methods: A set of human ventricular action potential (AP) models covering a range of experimental electrophysiological characteristics was constructed. These models incorporated stochasticity in major ionic currents as well as descriptions of β-adrenergic stimulation and mechanical effects to investigate the AP response to enhanced sympathetic activity. Statistical methods based on Automatic Relevance Determination and Canonical Correlation Analysis were developed to unravel individual and common factors contributing to BVR and LF patterning of APD in response to sympathetic provocation.

Results: Simulated results reproduced experimental evidences on the interactions between BVR and LF oscillations of AP duration (APD), with replication of the high inter-individual variability observed in both phenomena. I_CaL, I_Kr and I_K1 currents were identified as common ionic modulators of the inter-individual differences in BVR and LF oscillatory behavior and were shown to be crucial in determining susceptibility to arrhythmogenic events.

Conclusions: The calibrated family of human ventricular cell models proposed in this study allows reproducing experimentally reported interactions between BVR and LF oscillations of APD. Ionic factors involving I_CaL, I_Kr and I_K1 currents are found to underlie correlated increments in both phenomena in response to sympathetic provocation. A link to arrhythmogenesis is established for concomitantly elevated levels of BVR and its LF oscillations.

Periodic repolarization dynamics as a risk predictor after myocardial infarction: Prospective validation study

BACKGROUND

Periodic repolarization dynamics (PRD) is a novel electrocardiographic phenomenon that refers to sympathetic activity-associated low-frequency modulations of cardiac repolarization. Retrospective post-myocardial infarction (MI) studies revealed that increased PRD indicates an increased risk of subsequent death.

OBJECTIVE

This is the first prospective study to validate PRD in patients after MI receiving up-to-date treatment.

METHODS

Four hundred fifty-five survivors of MI (age ≤80 years) in sinus rhythm were enrolled. PRD was assessed from 20-minute electrocardiographic recordings (2048 Hz) and prospectively dichotomized at 5.75 deg². Primary and secondary end points were total mortality and cardiovascular mortality, respectively. Multivariable analyses additionally included Global Registry of Acute Coronary Events score (dichotomized at >140), left ventricular ejection fraction (dichotomized at ≤35%), diabetes mellitus, and deceleration capacity of heart rate (dichotomized at ≤2.5 ms). The prognostic power of PRD was evaluated using receiver operating characteristic curve analysis, Cox regression analysis, and the integrated discrimination improvement index.

RESULTS

During a median follow-up period of 27 months, 47 patients died. Twenty-three of these deaths were classified as cardiovascular. Increased PRD was significantly associated with both end points, yielding areas under receiver operating characteristic curves of 69.3% (60.2%-77.8%) and 79.1% (69.7%-86.7%) for total mortality and cardiovascular mortality, respectively (P < .001 for both). In multivariable analysis, increased PRD indicated a 2.2- and 9.5-fold risk of total mortality and cardiovascular mortality (P = .024 and P = .003, respectively). Addition of PRD to the models significantly improved the integrated discrimination improvement index for total (P = .047) and cardiovascular mortality (P = .007).

CONCLUSION

PRD is a strong and independent predictor of total mortality and cardiovascular mortality in patients after MI treated with contemporary therapy.

Risk prediction in post-infarction patients with moderately reduced left ventricular ejection fraction by combined assessment of the sympathetic and vagal cardiac autonomic nervous system

AIM

Most deaths after myocardial infarction (MI) occur in patients with normal or moderately reduced left ventricular ejection fraction (LVEF >35%). Periodic repolarization dynamics (PRD) and deceleration capacity (DC) are novel ECG-based markers related to sympathetic and vagal cardiac autonomic nervous system activity. Here, we test the combination of PRD and DC to predict risk in post-infarction patients with LVEF >35%.

METHODS AND RESULTS

We included 823 survivors of acute MI with LVEF >35%, aged ≤80years and in sinus rhythm. PRD and DC were obtained from 30-min ECG-recordings within the second week after index infarction and dichotomized at established cut-off values of ≥5.75deg² and ≤2.5ms, respectively. Patients were classified as having normal (CAF 0), partly abnormal (DC or PRD abnormal; CAF 1) or abnormal cardiac autonomic function (DC and PRD abnormal; CAF 2). Primary endpoint was 5-year all-cause mortality. Within the first 5years of follow-up, 51 patients died (6.2%). PRD and DC effectively stratified patients into low-risk (CAF 0; n=562), intermediate-risk (CAF 1; n=193) and high-risk patients (CAF 2; n=68) with cumulative 5-year mortality rates of 2.9%, 9.4% and 25.2%, respectively (p<0.001). On multivariable analyses, CAF was independent from established risk factors (GRACE-score, diabetes mellitus, mean heart rate, heart rate variability). Addition of CAF significantly improved the model (increase of C-statistics from 0.732 (0.651-0.812) to 0.777 (0.703-0.850), p=0.047; continuous NRI (0.400, 95% CI 0.230-0.560, p<0.001); IDI (0.056, 95% CI 0.022-0.122, p<0.001)).

CONCLUSION

CAF identifies new high-risk post-MI patients with LVEF >35% which might benefit from prophylactic strategies.

Prediction of sudden and non-sudden cardiac death in post-infarction patients with reduced left ventricular ejection fraction by periodic repolarization dynamics: MADIT-II substudy

Aims

To test the value of Periodic Repolarization Dynamics (PRD), a recently validated electrocardiographic marker of sympathetic activity, as a novel approach to predict sudden cardiac death (SCD) and non-sudden cardiac death (N-SCD) and to improve identification of patients that profit from ICD-implantation.

Methods and results

We included 856 post-infarction patients with left-ventricular ejection fraction (LVEF) ≤30% of the MADIT-II trial in sinus rhythm. Of these, 507 and 348 patients were randomized to ICD or conventional treatment. PRD was assessed from multipolar 10-min baseline ECGs. Primary and secondary endpoints were total mortality, SCD and N-SCD. Multivariable analyses included treatment group, QRS-duration, New York Heart Association classification, blood-urea nitrogen, diabetes mellitus, beta-blocker therapy and LVEF. During follow-up of 20.4 months, 119 patients died (53 SCD and 36 N-SCD). On multivariable analyses, increased PRD was a significant predictor of mortality (standardized coefficient 1.37[1.19–1.59]; P < 0.001) and SCD (1.40 [1.13–1.75]; P = 0.003) but also predicted N-SCD (1.41[1.10–1.81]; P = 0.006). While increased PRD predicted SCD in conventionally treated patients (1.61[1.23–2.11]; P < 0.001), it was predictive of N-SCD (1.63[1.28–2.09]; P < 0.001) and adequate ICD-therapies (1.20[1.03–1.39]; P = 0.017) in ICD-treated patients. ICD-treatment substantially reduced mortality in the lowest three PRD-quartiles by 53% (P = 0.001). However, there was no effect in the highest PRD-quartile (mortality increase by 29%; P = 0.412; P < 0.001 for difference) as the reduction of SCD was compensated by an increase of N-SCD.

Conclusion

In post-infarction patients with impaired LVEF, PRD is a significant predictor of SCD and N-SCD. Assessment of PRD is a promising tool to identify post-MI patients with reduced LVEF who might benefit from intensified treatment.

Predictive model identifies key network regulators of cardiomyocyte mechano-signaling

Mechanical strain is a potent stimulus for growth and remodeling in cells. Although many pathways have been implicated in stretch-induced remodeling, the control structures by which signals from distinct mechano-sensors are integrated to modulate hypertrophy and gene expression in cardiomyocytes remain unclear. Here, we constructed and validated a predictive computational model of the cardiac mechano-signaling network in order to elucidate the mechanisms underlying signal integration. The model identifies calcium, actin, Ras, Raf1, PI3K, and JAK as key regulators of cardiac mechano-signaling and characterizes crosstalk logic imparting differential control of transcription by AT1R, integrins, and calcium channels. We find that while these regulators maintain mostly independent control over distinct groups of transcription factors, synergy between multiple pathways is necessary to activate all the transcription factors necessary for gene transcription and hypertrophy. We also identify a PKG-dependent mechanism by which valsartan/sacubitril, a combination drug recently approved for treating heart failure, inhibits stretch-induced hypertrophy, and predict further efficacious pairs of drug targets in the network through a network-wide combinatorial search.

Common stresses such as high blood pressure or heart attack can lead to heart failure, which afflicts over 25 million people worldwide. These stresses cause cardiomyocytes to grow and remodel, which may initially be beneficial but ultimately worsen heart function. Current heart failure drugs such as beta-blockers counteract biochemical cues prompting cardiomyocyte growth, yet mechanical cues to cardiomyocytes such as stretch are just as important in driving cardiac dysfunction. However, no pharmacological treatments have yet been approved that specifically target mechano-signaling, in part because it is not clear how cardiomyocytes integrate signals from multiple mechano-responsive sensors and pathways into their decision to grow. To address this challenge, we built a systems-level computational model that represents 125 interactions between 94 stretch-responsive signaling molecules. The model correctly predicts 134 of 172 previous independent experimental observations, and identifies the key regulators of stretch-induced cardiomyocyte remodeling. Although cardiomyocytes have many mechano-signaling pathways that function largely independently, we find that cooperation between them is necessary to cause growth and remodeling. We identify mechanisms by which a recently approved heart failure drug pair affects mechano-signaling, and we further predict additional pairs of drug targets that could be used to help reverse heart failure.

A Time-Varying Nonparametric Methodology for Assessing Changes in QT Variability Unrelated to Heart Rate Variability

OBJECTIVE

To propose and test a novel methodology to measure changes in QT interval variability (QTV) unrelated to RR interval variability (RRV) in nonstationary conditions.

METHODS

Time-frequency coherent and residual spectra representing QTV related (QTVrRRV) and unrelated (QTVuRRV) to RRV, respectively, are estimated using time-frequency Cohen's class distributions. The proposed approach decomposes the nonstationary output spectrum of any two-input one-output model with uncorrelated inputs into two spectra representing the information related and unrelated to one of the two inputs, respectively. An algorithm to correct for the bias of the time-frequency coherence function between QTV and RRV is proposed to provide accurate estimates of both QTVuRRV and QTVrRRV. Two simulation studies were conducted to assess the methodology in challenging nonstationary conditions and data recorded during head-up tilt in 16 healthy volunteers were analyzed.

RESULTS

In the simulation studies, QTVuRRV changes were tracked with only a minor delay due to the filtering necessary to estimate the nonstationary spectra. The correlation coefficient between theoretical and estimated patterns was even for extremely noisy recordings (signal to noise ratio (SNR) in QTV dB). During head-up tilt, QTVrRRV explained the largest proportion of QTV, whereas QTVuRRV showed higher relative increase than QTV or QTVrRRV in all spectral bands ( for most pairwise comparisons).

CONCLUSION

The proposed approach accurately tracks changes in QTVuRRV. Head-up tilt induced a slightly greater increase in QTVuRRV than in QTVrRRV.

SIGNIFICANCE

The proposed index QTVuRRV may represent an indirect measure of intrinsic ventricular repolarization variability, a marker of cardiac instability associated with sympathetic ventricular modulation and sudden cardiac death.

Identifying high-risk post-infarction patients by autonomic testing - Below the tip of the iceberg

Despite major advances in medical therapies late mortality after myocardial infarction (MI) is still high. A substantial proportion of post-MI patients die from sudden cardiac death. Prophylactic implantable-cardioverter defibrillator (ICD) therapy has been established for post-MI patients with reduced left ventricular ejection fraction (LVEF ≤35%). However, most patients who die after MI have an LVEF >35%. For this large group of patients, no specific prophylactic strategies exist. There is strong evidence that measures of cardiac autonomic dysfunction after MI provide important prognostic information in post-MI patients with preserved LVEF. Combinations of autonomic markers can identify high-risk patients after MI with LVEF >35% whose prognosis is equally worse than that of patients with LVEF ≤35%. The ongoing REFINE-ICD (NCT00673842) and SMART-MI trials (NCT02594488) test different preventive strategies in high-risk post-MI patients with cardiac autonomic dysfunction and LVEF 36-50%. While REFINE-ICD follows the traditional concept of ICD-implantation, SMART-MI uses implantable cardiac monitors with remote monitoring capabilities to sensitively detect asymptomatic, but prognostically relevant arrhythmias that could trigger specific diagnostic and therapeutic interventions.

The control of cardiac ventricular excitability by autonomic pathways

Central to the genesis of ventricular cardiac arrhythmia are variations in determinants of excitability. These involve individual ionic channels and transporters in cardiac myocytes but also tissue factors such as variable conduction of the excitation wave, fibrosis and source-sink mismatch. It is also known that in certain diseases and particularly the channelopathies critical events occur with specific stressors. For example, in hereditary long QT syndrome due to mutations in KCNQ1 arrhythmic episodes are provoked by exercise and in particular swimming. Thus not only is the static substrate important but also how this is modified by dynamic signalling events associated with common physiological responses. In this review, we examine the regulation of ventricular excitability by signalling pathways from a cellular and tissue perspective in an effort to identify key processes, effectors and potential therapeutic approaches. We specifically focus on the autonomic nervous system and related signalling pathways.

Interactions between Activation and Repolarization Restitution Properties in the Intact Human Heart: In-Vivo Whole-Heart Data and Mathematical Description

BACKGROUND

The restitution of the action potential duration (APDR) and conduction velocity (CVR) are mechanisms whereby cardiac excitation and repolarization adapt to changes in heart rate. They modulate the vulnerability to dangerous arrhythmia, but the mechanistic link between restitution and arrhythmogenesis remains only partially understood.

METHODS

This paper provides an experimental and theoretical study of repolarization and excitation restitution properties and their interactions in the intact human epicardium. The interdependence between excitation and repolarization dynamic is studied in 8 patients (14 restitution protocols, 1722 restitution curves) undergoing global epicardial mapping with multi-electrode socks before open heart surgery. A mathematical description of the contribution of both repolarization and conduction dynamics to the steepness of the APDR slope is proposed.

RESULTS

This study demonstrates that the APDR slope is a function of both activation and repolarization dynamics. At short cycle length, conduction delay significantly increases the APDR slope by interacting with the diastolic interval. As predicted by the proposed mathematical formulation, the APDR slope was more sensitive to activation time prolongation than to the simultaneous shortening of repolarization time. A steep APDR slope was frequently identified, with 61% of all cardiac sites exhibiting an APDR slope > 1, suggesting that a slope > 1 may not necessarily promote electrical instability in the human epicardium. APDR slope did not change for different activation or repolarization times, and it was not a function of local baseline APD. However, it was affected by the spatial organization of electrical excitation, suggesting that in tissue APDR is not a unique function of local electrophysiological properties. Spatial heterogeneity in both activation and repolarization restitution contributed to the increase in the modulated dispersion of repolarization, which for short cycle length was as high as 250 ms. Heterogeneity in conduction velocity restitution can translate into both activation and repolarization dispersion and increase cardiac instability. The proposed mathematical formulation shows an excellent agreement with the experimental data (correlation coefficient r = 0.94) and provides a useful tool for the understanding of the complex interactions between activation and repolarization restitution properties as well as between their measurements.