Synergy between CaMKII substrates and β-adrenergic signaling in regulation of cardiac myocyte Ca(2+) handling

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Computational Modeling of Cardiac Electrophysiology

Mathematical modeling and simulation are well-established and powerful tools to integrate experimental data of individual components of cardiac electrophysiology, excitation-contraction coupling, and regulatory signaling pathways, to gain quantitative and mechanistic insight into pathophysiological processes and guide therapeutic strategies. Here, we briefly describe the processes governing cardiac myocyte electrophysiology and Ca2+ handling and their regulation, as well as action potential propagation in tissue. We discuss the models and methods used to describe these phenomena, including procedures for model parameterization and validation, in addition to protocols for model interrogation and analysis and techniques that account for phenotypic variability and parameter uncertainty. Our objective is to provide a summary of basic concepts and approaches as a resource for scientists training in this discipline and for all researchers aiming to gain an understanding of cardiac modeling studies.

Integrative human atrial modelling unravels interactive protein kinase A and Ca2+/calmodulin-dependent protein kinase II signalling as key determinants of atrial arrhythmogenesis

AIMS

Atrial fibrillation (AF), the most prevalent clinical arrhythmia, is associated with atrial remodelling manifesting as acute and chronic alterations in expression, function, and regulation of atrial electrophysiological and Ca2+-handling processes. These AF-induced modifications crosstalk and propagate across spatial scales creating a complex pathophysiological network, which renders AF resistant to existing pharmacotherapies that predominantly target transmembrane ion channels. Developing innovative therapeutic strategies requires a systems approach to disentangle quantitatively the pro-arrhythmic contributions of individual AF-induced alterations.

METHODS AND RESULTS

Here, we built a novel computational framework for simulating electrophysiology and Ca2+-handling in human atrial cardiomyocytes and tissues, and their regulation by key upstream signalling pathways [i.e. protein kinase A (PKA), and Ca2+/calmodulin-dependent protein kinase II (CaMKII)] involved in AF-pathogenesis. Populations of atrial cardiomyocyte models were constructed to determine the influence of subcellular ionic processes, signalling components, and regulatory networks on atrial arrhythmogenesis. Our results reveal a novel synergistic crosstalk between PKA and CaMKII that promotes atrial cardiomyocyte electrical instability and arrhythmogenic triggered activity. Simulations of heterogeneous tissue demonstrate that this cellular triggered activity is further amplified by CaMKII- and PKA-dependent alterations of tissue properties, further exacerbating atrial arrhythmogenesis.

CONCLUSIONS

Our analysis reveals potential mechanisms by which the stress-associated adaptive changes turn into maladaptive pro-arrhythmic triggers at the cellular and tissue levels and identifies potential anti-AF targets. Collectively, our integrative approach is powerful and instrumental to assemble and reconcile existing knowledge into a systems network for identifying novel anti-AF targets and innovative approaches moving beyond the traditional ion channel-based strategy.

A multiscale predictive digital twin for neurocardiac modulation

Cardiac function is tightly regulated by the autonomic nervous system (ANS). Activation of the sympathetic nervous system increases cardiac output by increasing heart rate and stroke volume, while parasympathetic nerve stimulation instantly slows heart rate. Importantly, imbalance in autonomic control of the heart has been implicated in the development of arrhythmias and heart failure. Understanding of the mechanisms and effects of autonomic stimulation is a major challenge because synapses in different regions of the heart result in multiple changes to heart function. For example, nerve synapses on the sinoatrial node (SAN) impact pacemaking, while synapses on contractile cells alter contraction and arrhythmia vulnerability. Here, we present a multiscale neurocardiac modelling and simulator tool that predicts the effect of efferent stimulation of the sympathetic and parasympathetic branches of the ANS on the cardiac SAN and ventricular myocardium. The model includes a layered representation of the ANS and reproduces firing properties measured experimentally. Model parameters are derived from experiments and atomistic simulations. The model is a first prototype of a digital twin that is applied to make predictions across all system scales, from subcellular signalling to pacemaker frequency to tissue level responses. We predict conditions under which autonomic imbalance induces proarrhythmia and can be modified to prevent or inhibit arrhythmia. In summary, the multiscale model constitutes a predictive digital twin framework to test and guide high-throughput prediction of novel neuromodulatory therapy. KEY POINTS: A multi-layered model representation of the autonomic nervous system that includes sympathetic and parasympathetic branches, each with sparse random intralayer connectivity, synaptic dynamics and conductance based integrate-and-fire neurons generates firing patterns in close agreement with experiment. A key feature of the neurocardiac computational model is the connection between the autonomic nervous system and both pacemaker and contractile cells, where modification to pacemaker frequency drives initiation of electrical signals in the contractile cells. We utilized atomic-scale molecular dynamics simulations to predict the association and dissociation rates of noradrenaline with the β-adrenergic receptor. Multiscale predictions demonstrate how autonomic imbalance may increase proclivity to arrhythmias or be used to terminate arrhythmias. The model serves as a first step towards a digital twin for predicting neuromodulation to prevent or reduce disease.

A detailed mathematical model of the human atrial cardiomyocyte: integration of electrophysiology and cardiomechanics

Mechano-electric regulations (MER) play an important role in the maintenance of cardiac performance. Mechano-calcium and mechano-electric feedback (MCF and MEF) pathways adjust the cardiomyocyte contractile force according to mechanical perturbations and affects electro-mechanical coupling. MER integrates all these regulations in one unit resulting in a complex phenomenon. Computational modelling is a useful tool to accelerate the mechanistic understanding of complex experimental phenomena. We have developed a novel model that integrates the MER loop for human atrial cardiomyocytes with proper consideration of feedforward and feedback pathways. The model couples a modified version of the action potential (AP) Koivumäki model with the contraction model by Quarteroni group. The model simulates iso-sarcometric and isometric twitches and the feedback effects on AP and Ca2+ -handling. The model showed a biphasic response of Ca2+ transient (CaT) peak to increasing pacing rates and highlights the possible mechanisms involved. The model has shown a shift of the threshold for AP and CaT alternans from 4.6 to 4 Hz under post-operative atrial fibrillation, induced by depressed SERCA activity. The alternans incidence was dependent on a chain of mechanisms including RyRs availability time, MCF coupling, CaMKII phosphorylation, and the stretch levels. As a result, the model predicted a 10% slowdown of conduction velocity for a 20% stretch, suggesting a role of stretch in creation of substrate formation for atrial fibrillation. Overall, we conclude that the developed model provides a physiological CaT followed by a physiological twitch. This model can open pathways for the future studies of human atrial electromechanics. KEY POINTS: With the availability of human atrial cellular data, interest in atrial-specific model integration has been enhanced. We have developed a detailed mathematical model of human atrial cardiomyocytes including the mechano-electric regulatory loop. The model has gone through calibration and evaluation phases against a wide collection of available human in-vitro data. The usefulness of the model for analysing clinical problems has been preliminaryly tested by simulating the increased incidence of Ca2+ transient and action potential alternans at high rates in post-operative atrial fibrillation condition. The model determines the possible role of mechano-electric feedback in alternans incidence, which can increase vulnerability to atrial arrhythmias by varying stretch levels. We found that our physiologically accurate description of Ca2+ handling can reproduce many experimental phenomena and can help to gain insights into the underlying pathophysiological mechanisms.

Pro-arrhythmic role of adrenergic spatial densities in the human atria: An in-silico study

Chronic stress among young patients (≤ 45 years old) could result in autonomic dysfunction. Autonomic dysfunction could be exhibited via sympathetic hyperactivity, sympathetic nerve sprouting, and diffuse adrenergic stimulation in the atria. Adrenergic spatial densities could alter atrial electrophysiology and increase arrhythmic susceptibility. Therefore, we examined the role of adrenergic spatial densities in creating arrhythmogenic substrates in silico. We simulated three 25 cm2 atrial sheets with varying adrenergic spatial densities (ASD), activation rates, and external transmembrane currents. We measured their effects on spatial and temporal heterogeneity of action potential durations (APD) at 50% and 20%. Increasing ASD shortens overall APD, and maximum spatial heterogeneity (31%) is achieved at 15% ASD. The addition of a few (5% to 10%) adrenergic elements decreases the excitation threshold, below 18 μA/cm2, while ASDs greater than 10% increase their excitation threshold up to 22 μA/cm2. Increase in ASD during rapid activation increases APD50 and APD20 by 21% and 41%, respectively. Activation times of captured beats during rapid activation could change by as much as 120 ms from the baseline cycle length. Rapidly activated atrial sheets with high ASDs significantly increase temporal heterogeneity of APD50 and APD20. Rapidly activated atrial sheets with 10% ASD have a high likelihood (0.7 ± 0.06) of fragmenting otherwise uniform wavefronts due to the transient inexcitability of adrenergically stimulated elements, producing an effective functional block. The likelihood of wave fragmentation due to ASD highly correlates with the spatial variations of APD20 (ρ = 0.90, p = 0.04). Our simulations provide a novel insight into the contributions of ASD to spatial and temporal heterogeneities of APDs, changes in excitation thresholds, and a potential explanation for wave fragmentation in the human atria due to sympathetic hyperactivity. Our work may aid in elucidating an electrophysiological link to arrhythmia initiation due to chronic stress among young patients.

Lysosomal calcium loading promotes spontaneous calcium release by potentiating ryanodine receptors

Spontaneous calcium release by ryanodine receptors (RyRs) due to intracellular calcium overload results in delayed afterdepolarizations, closely associated with life-threatening arrhythmias. In this regard, inhibiting lysosomal calcium release by two-pore channel 2 (TPC2) knockout has been shown to reduce the incidence of ventricular arrhythmias under β-adrenergic stimulation. However, mechanistic investigations into the role of lysosomal function on RyR spontaneous release remain missing. We investigate the calcium handling mechanisms by which lysosome function modulates RyR spontaneous release, and determine how lysosomes are able to mediate arrhythmias by its influence on calcium loading. Mechanistic studies were conducted using a population of biophysically detailed mouse ventricular models including for the first time modeling of lysosomal function, and calibrated by experimental calcium transients modulated by TPC2. We demonstrate that lysosomal calcium uptake and release can synergistically provide a pathway for fast calcium transport, by which lysosomal calcium release primarily modulates sarcoplasmic reticulum calcium reuptake and RyR release. Enhancement of this lysosomal transport pathway promoted RyR spontaneous release by elevating RyR open probability. In contrast, blocking either lysosomal calcium uptake or release revealed an antiarrhythmic impact. Under conditions of calcium overload, our results indicate that these responses are strongly modulated by intercellular variability in L-type calcium current, RyR release, and sarcoplasmic reticulum calcium-ATPase reuptake. Altogether, our investigations identify that lysosomal calcium handling directly influences RyR spontaneous release by regulating RyR open probability, suggesting antiarrhythmic strategies and identifying key modulators of lysosomal proarrhythmic action.

Diversity of cells and signals in the cardiovascular system

This white paper is the outcome of the seventh UC Davis Cardiovascular Research Symposium on Systems Approach to Understanding Cardiovascular Disease and Arrhythmia. This biannual meeting aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2022 Symposium was 'Cell Diversity in the Cardiovascular System, cell-autonomous and cell-cell signalling'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies, and challenges in examining cell and signal diversity, co-ordination and interrelationships involved in cardiovascular function. This paper originates from the topics of formal presentations and informal discussions from the Symposium, which aimed to develop a holistic view of how the multiple cell types in the cardiovascular system integrate to influence cardiovascular function, disease progression and therapeutic strategies. The first section describes the major cell types (e.g. cardiomyocytes, vascular smooth muscle and endothelial cells, fibroblasts, neurons, immune cells, etc.) and the signals involved in cardiovascular function. The second section emphasizes the complexity at the subcellular, cellular and system levels in the context of cardiovascular development, ageing and disease. Finally, the third section surveys the technological innovations that allow the interrogation of this diversity and advancing our understanding of the integrated cardiovascular function and dysfunction.

Mechanisms of spontaneous Ca2+ release-mediated arrhythmia in a novel 3D human atrial myocyte model: I. Transverse-axial tubule variation

Intracellular calcium (Ca2+ ) cycling is tightly regulated in the healthy heart ensuring effective contraction. This is achieved by transverse (t)-tubule membrane invaginations that facilitate close coupling of key Ca2+ -handling proteins such as the L-type Ca2+ channel and Na+ -Ca2+ exchanger (NCX) on the cell surface with ryanodine receptors (RyRs) on the intracellular Ca2+ store. Although less abundant and regular than in the ventricle, t-tubules also exist in atrial myocytes as a network of transverse invaginations with axial extensions known as the transverse-axial tubule system (TATS). In heart failure and atrial fibrillation, there is TATS remodelling that is associated with aberrant Ca2+ -handling and Ca2+ -induced arrhythmic activity; however, the mechanism underlying this is not fully understood. To address this, we developed a novel 3D human atrial myocyte model that couples electrophysiology and Ca2+ -handling with variable TATS organization and density. We extensively parameterized and validated our model against experimental data to build a robust tool examining TATS regulation of subcellular Ca2+ release. We found that varying TATS density and thus the localization of key Ca2+ -handling proteins has profound effects on Ca2+ handling. Following TATS loss, there is reduced NCX that results in increased cleft Ca2+ concentration through decreased Ca2+ extrusion. This elevated Ca2+ increases RyR open probability causing spontaneous Ca2+ releases and the promotion of arrhythmogenic waves (especially in the cell interior) leading to voltage instabilities through delayed afterdepolarizations. In summary, the present study demonstrates a mechanistic link between TATS remodelling and Ca2+ -driven proarrhythmic behaviour that probably reflects the arrhythmogenic state observed in disease. KEY POINTS: Transverse-axial tubule systems (TATS) modulate Ca2+ handling and excitation-contraction coupling in atrial myocytes, with TATS remodelling in heart failure and atrial fibrillation being associated with altered Ca2+ cycling and subsequent arrhythmogenesis. To investigate the poorly understood mechanisms linking TATS variation and spontaneous Ca2+ release, we built, parameterized and validated a 3D human atrial myocyte model coupling electrophysiology and spatially-detailed subcellular Ca2+ handling governed by the TATS. Simulated TATS loss causes diastolic Ca2+ and voltage instabilities through reduced Na+ -Ca2+ exchanger-mediated Ca2+ removal, cleft Ca2+ accumulation and increased ryanodine receptor open probability, resulting in spontaneous Ca2+ release and promotion of arrhythmogenic waves and delayed afterdepolarizations. At fast electrical rates typical of atrial tachycardia/fibrillation, spontaneous Ca2+ releases are larger and more frequent in the cell interior than at the periphery. Our work provides mechanistic insight into how atrial TATS remodelling can lead to Ca2+ -driven instabilities that may ultimately contribute to the arrhythmogenic state in disease.

Models of the cardiac L-type calcium current: A quantitative review

The L-type calcium current (I CaL ) plays a critical role in cardiac electrophysiology, and models ofI CaL are vital tools to predict arrhythmogenicity of drugs and mutations. Five decades of measuring and modelingI CaL have resulted in several competing theories (encoded in mathematical equations). However, the introduction of new models has not typically been accompanied by a data-driven critical comparison with previous work, so that it is unclear which model is best suited for any particular application. In this review, we describe and compare 73 published mammalianI CaL models and use simulated experiments to show that there is a large variability in their predictions, which is not substantially diminished when grouping by species or other categories. We provide model code for 60 models, list major data sources, and discuss experimental and modeling work that will be required to reduce this huge list of competing theories and ultimately develop a community consensus model ofI CaL . This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Molecular and Cellular Physiology.

Mutation-Specific Differences in Kv7.1 (KCNQ1) and Kv11.1 (KCNH2) Channel Dysfunction and Long QT Syndrome Phenotypes

The electrocardiogram (ECG) empowered clinician scientists to measure the electrical activity of the heart noninvasively to identify arrhythmias and heart disease. Shortly after the standardization of the 12-lead ECG for the diagnosis of heart disease, several families with autosomal recessive (Jervell and Lange-Nielsen Syndrome) and dominant (Romano-Ward Syndrome) forms of long QT syndrome (LQTS) were identified. An abnormally long heart rate-corrected QT-interval was established as a biomarker for the risk of sudden cardiac death. Since then, the International LQTS Registry was established; a phenotypic scoring system to identify LQTS patients was developed; the major genes that associate with typical forms of LQTS were identified; and guidelines for the successful management of patients advanced. In this review, we discuss the molecular and cellular mechanisms for LQTS associated with missense variants in KCNQ1 (LQT1) and KCNH2 (LQT2). We move beyond the "benign" to a "pathogenic" binary classification scheme for different KCNQ1 and KCNH2 missense variants and discuss gene- and mutation-specific differences in K+ channel dysfunction, which can predispose people to distinct clinical phenotypes (e.g., concealed, pleiotropic, severe, etc.). We conclude by discussing the emerging computational structural modeling strategies that will distinguish between dysfunctional subtypes of KCNQ1 and KCNH2 variants, with the goal of realizing a layered precision medicine approach focused on individuals.

Association Study of Genetic Variants in Calcium Signaling-Related Genes With Cardiovascular Diseases

Background: Calcium ions (Ca2+) play an essential role in excitation-contraction coupling in the heart. The association between cardiovascular diseases (CVDs) and genetic polymorphisms in key regulators of Ca2+ homeostasis is well established but still inadequately understood. Methods: The associations of 11,274 genetic variants located in nine calcium signaling-related genes with 118 diseases of the circulatory system were explored using a large sample from the United Kingdom Biobank (N = 308,366). The clinical outcomes in electronic health records were mapped to the phecode system. Survival analyses were employed to study the role of variants in CVDs incidence and mortality. Phenome-wide association studies (PheWAS) were performed to investigate the effect of variants on cardiovascular risk factors. Results: The reported association between rs1801253 in β1-adrenergic receptor (ADRB1) and hypertension was successfully replicated, and we additionally found the blood pressure-lowering G allele of this variant was associated with a delayed onset of hypertension and a decreased level of apolipoprotein A. The association of rs4484922 in calsequestrin 2 (CASQ2) with atrial fibrillation/flutter was identified, and this variant also displayed nominal evidence of association with QRS duration and carotid intima-medial thickness. Moreover, our results indicated suggestive associations of rs79613429 in ryanodine receptor 2 (RYR2) with precordial pain. Conclusion: Multiple novel associations established in our study highlight genetic testing as a useful method for CVDs diagnosis and prevention.

Quantitative cross-species translators of cardiac myocyte electrophysiology: Model training, experimental validation, and applications

Animal experimentation is key in the evaluation of cardiac efficacy and safety of novel therapeutic compounds. However, interspecies differences in the mechanisms regulating excitation-contraction coupling can limit the translation of experimental findings from animal models to human physiology and undermine the assessment of drugs’ efficacy and safety. Here, we built a suite of translators for quantitatively mapping electrophysiological responses in ventricular myocytes across species. We trained these statistical operators using a broad dataset obtained by simulating populations of our biophysically detailed computational models of action potential and Ca2+ transient in mouse, rabbit, and human. We then tested our translators against experimental data describing the response to stimuli, such as ion channel block, change in beating rate, and β-adrenergic challenge. We demonstrate that this approach is well suited to predicting the effects of perturbations across different species or experimental conditions and suggest its integration into mechanistic studies and drug development pipelines.

Extracorporeal Shock Wave Enhanced Exogenous Mitochondria into Adipose-Derived Mesenchymal Stem Cells and Further Preserved Heart Function in Rat Dilated Cardiomyopathy

This study tested whether extracorporeal shock wave (ECSW) supported-exogenous mitochondria (Mito) into adipose-derived mesenchymal stem cells (ADMSCs) would preserve left-ventricular-ejection-fraction (LVEF) in doxorubicin/12 mg/kg-induced dilated cardiomyopathy (DCM) rat. Adult-male-SD rats were equally categorized into group 1 (sham-control), group 2 (DCM), group 3 (DCM + ECSW/1.5 mJ/mm² for 140 shots/week × 3 times/since day 14 after DCM induction), group 4 (DCM + ECSW/1.5 mJ/mm²/100 shots-assisted mito delivery (500 μg) into ADMSCs/1.2 × 10⁶ cells, then implanted into LV myocardium day 14 after DCM induction) and group 5 (DCM + ECSW-assisted mito delivery into ADMSCs/1.2 × 10⁶ cells, then implanted into LV, followed by ECSW/1.5 mJ/mm² for 140 shots/week × 3 times/since day 14 after DCM induction) and euthanized by day 49. Microscopic findings showed mitochondria were abundantly enhanced by ECSW into H9C2 cells. The q-PCR showed a significant increase in relative number of mitDNA in mitochondrial-transferred H9C2 cells than in control group (p < 0.01). The angiogenesis/angiogenesis factors (VEGF/SDF-1α/IG-F1) in HUVECs were significantly progressively increased by a stepwise-increased amount of ECSW energy (0.1/0.25/0.35 mJ/mm²) (all p < 0.001). The 49-day LVEF was highest in group 1 and significantly progressively increased from groups 2 to 5 (all p < 0.0001). Cardiomyocyte size/fibrosis exhibited an opposite pattern of LVEF, whereas cellular/protein levels of angiogenesis factors (VEGF/SDF-1α) in myocardium were significantly progressively increased from groups 1 to 5 (all p < 0.0001). The protein expressions of apoptotic/mitochondrial (cleaved-caspase-3/cleaved-PARP/mitochondrial-Bax/cytosolic-cytochrome-C), fibrotic (p-Smad3/TGF-ß), oxidative-stress (NOX-1/NOX-2) and pressure-overload/heart failure (BNP/ß-MHC) biomarkers exhibited an opposite pattern of LVEF among the five groups (all p < 0.0001). ECSW-assisted mitochondrial-delivery into ADMSCs plus ECSW offered an additional benefit for preserving LVEF in DCM rat.

Molecular noise filtering in the β-adrenergic signaling network by phospholamban pentamers

Phospholamban (PLN) is an important regulator of cardiac calcium handling due to its ability to inhibit the calcium ATPase SERCA. β-Adrenergic stimulation reverses SERCA inhibition via PLN phosphorylation and facilitates fast calcium reuptake. PLN also forms pentamers whose physiological significance has remained elusive. Using mathematical modeling combined with biochemical and cell biological experiments, we show that pentamers regulate both the dynamics and steady-state levels of monomer phosphorylation. Substrate competition by pentamers and a feed-forward loop involving inhibitor-1 can delay monomer phosphorylation by protein kinase A (PKA), whereas cooperative pentamer dephosphorylation enables bistable PLN steady-state phosphorylation. Simulations show that phosphorylation delay and bistability act as complementary filters that reduce the effect of random fluctuations in PKA activity, thereby ensuring consistent monomer phosphorylation and SERCA activity despite noisy upstream signals. Preliminary analyses suggest that the PLN mutation R14del could impair noise filtering, offering a new perspective on how this mutation causes cardiac arrhythmias.

Dynamics of adrenergic signaling in cardiac myocytes and implications for pharmacological treatment

Dense innervation of the heart by the sympathetic nervous system (SNS) allows cardiac output to respond appropriately to the needs of the body under varying conditions, but occasionally the abrupt onset of SNS activity can trigger cardiac arrhythmias. Sympathetic activity leads to the release of norepinephrine (NE) onto cardiomyocytes, activating β₁-adrenergic receptors (β₁-ARs) and leading to the production of the second messenger cyclic AMP (cAMP). Upon sudden activation of β₁-ARs in experiments, intracellular cAMP can transiently rise to a high concentration before converging to a steady state level. Although changes to cellular cAMP concentration are important in modulating the overall cardiovascular response to sympathetic tone, the underlying mechanisms of the cAMP transients and the parameters that control their magnitude are unclear.

We reduce a detailed computational model of the β₁-adrenergic signaling cascade to a system of two differential equations by eliminating extraneous variables and applying quasi-steady state approximation. The structure of the reduced model reveals that the large cAMP transients associated with abrupt β₁-AR activation are generated by the interplay of production/degradation of cAMP and desensitization/resensitization of β₁-ARs. The reduced model is used to predict how the dynamics of intracellular cAMP depend on the concentrations of norepinephrine (NE), phosphodiesterases 3 and 4 (PDE3,4), G-protein coupled receptor kinase 2 (GRK2), and β₁-AR, in healthy conditions and a simple model of early stages of heart failure.

The key findings of the study are as follows: 1) Applying a reduced model of the dynamics of cardiac sympathetic signaling we show that the concentrations of two variables, cAMP and non-desensitized β₁-AR, capture the overall dynamics of sympathetic signaling; 2) The key factors influencing cAMP production are AC activity and PDE3,4 activity, while those that directly impact β₁-AR phosphorylation are GRK2 and PKA₁. Thus, disease states that affect sympathetic control of the heart can be thoroughly assessed by studying AC activity, PDE3,4, GRK2 and PKA activity, as these factors directly impact cAMP production/degradation and β₁-AR (de) phosphorylation and are therefore predicted to comprise the most effective pharmaceutical targets in diseases affecting cardiac β₁-adrenergic signaling.

Protective Effects of 3,4-Seco-Lupane Triterpenes from Food Raw Materials of the Leaves of Eleutherococcus Senticosus and Eleutherococcus Sessiliflorus on Arrhythmia Induced by Barium Chloride

As a traditional wild vegetable and food raw material, the leaves of Eleutherococcus senticosus and Eleutherococcus sessiliflorus are rich in 3,4-seco-lupane triterpenes, including chiisanoside (CSS), divaroside (DVS), sessiloside-A (SSA), and chiisanogenin (CSG). This study was conducted to evaluate the anti-arrhythmic effects of these 3,4-seco-lupane triterpenes. Evaluation of the cytotoxicity of compounds was performed by measuring cell viability and apoptosis with the CCK-8 assay. In vivo, arrhythmia was induced by rapid injection of BaCl₂ via rat caudal vein. The occurrence time and duration of arrhythmias in rats were studied. The levels of SOD and MDA in serum, and Na⁺ -K⁺ -ATPase and Ca²⁺ -Mg²⁺ -ATPase in myocardial homogenate were detected by ELISA. The histopathological changes of rats myocardial were observed by HE staining. Changes in the expression of PKA and related proteins were detected by Western blot. The 3,4-seco-lupane triterpenes interactions with protein kinase A were analyzed by molecular docking. In the present study, we found that 3,4-seco-lupane triterpenes exhibited powerful anti-arrhythmic activity, especially DVS completely relieved the ventricular arrhythmia induced by BaCl₂ . This study suggests that the leaves of E. senticosus and E. sessiliflorus might be used as functional food materials to prevent arrhythmia, and DVS can potentially be further developed as an anti-arrhythmic drug.

Evolution of mathematical models of cardiomyocyte electrophysiology

Advanced computational techniques and mathematical modeling have become more and more important to the study of cardiac electrophysiology. In this review, we provide a brief history of the evolution of cardiomyocyte electrophysiology models and highlight some of the most important ones that had a major impact on our understanding of the electrical activity of the myocardium and associated transmembrane ion fluxes in normal and pathological states. We also present the use of these models in the study of various arrhythmogenesis mechanisms, particularly the integration of experimental pharmacology data into advanced humanized models for in silico proarrhythmogenic risk prediction as an essential component of the Comprehensive in vitro Proarrhythmia Assay (CiPA) drug safety paradigm.

Computational modeling approaches to cAMP/PKA signaling in cardiomyocytes

The cAMP/PKA pathway is a fundamental regulator of excitation-contraction coupling in cardiomyocytes. Activation of cAMP has a variety of downstream effects on cardiac function including enhanced contraction, accelerated relaxation, adaptive stress response, mitochondrial regulation, and gene transcription. Experimental advances have shed light on the compartmentation of cAMP and PKA, which allow for control over the varied targets of these second messengers and is disrupted in heart failure conditions. Computational modeling is an important tool for understanding the spatial and temporal complexities of this system. In this review article, we outline the advances in computational modeling that have allowed for deeper understanding of cAMP/PKA dynamics in the cardiomyocyte in health and disease, and explore new modeling frameworks that may bring us closer to a more complete understanding of this system. We outline various compartmental and spatial signaling models that have been used to understand how β-adrenergic signaling pathways function in a variety of simulation conditions. We also discuss newer subcellular models of cardiovascular function that may be used as templates for the next phase of computational study of cAMP and PKA in the heart, and outline open challenges which are important to consider in future models.

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.

Context-specific network modeling identifies new crosstalk in β-adrenergic cardiac hypertrophy

Cardiac hypertrophy is a context-dependent phenomenon wherein a myriad of biochemical and biomechanical factors regulate myocardial growth through a complex large-scale signaling network. Although numerous studies have investigated hypertrophic signaling pathways, less is known about hypertrophy signaling as a whole network and how this network acts in a context-dependent manner. Here, we developed a systematic approach, CLASSED (Context-specific Logic-bASed Signaling nEtwork Development), to revise a large-scale signaling model based on context-specific data and identify main reactions and new crosstalks regulating context-specific response. CLASSED involves four sequential stages with an automated validation module as a core which builds a logic-based ODE model from the interaction graph and outputs the model validation percent. The context-specific model is developed by estimation of default parameters, classified qualitative validation, hybrid Morris-Sobol global sensitivity analysis, and discovery of missing context-dependent crosstalks. Applying this pipeline to our prior-knowledge hypertrophy network with context-specific data revealed key signaling reactions which distinctly regulate cell response to isoproterenol, phenylephrine, angiotensin II and stretch. Furthermore, with CLASSED we developed a context-specific model of β-adrenergic cardiac hypertrophy. The model predicted new crosstalks between calcium/calmodulin-dependent pathways and upstream signaling of Ras in the ISO-specific context. Experiments in cardiomyocytes validated the model’s predictions on the role of CaMKII-Gβγ and CaN-Gβγ interactions in mediating hypertrophic signals in ISO-specific context and revealed a difference in the phosphorylation magnitude and translocation of ERK1/2 between cardiac myocytes and fibroblasts. CLASSED is a systematic approach for developing context-specific large-scale signaling networks, yielding insights into new-found crosstalks in β-adrenergic cardiac hypertrophy.

Pathological cardiac hypertrophy is a disease in which the heart grows abnormally in response to different motivators such as high blood pressure or variations in hormones and growth factors. The shape of the heart after its growth depends on the context in which it grows. Since cell signaling in the cardiac cells plays a key role in the determination of heart shape, a thorough understanding of cardiac cells signaling in each context enlightens the mechanisms which control response of cardiac cells. However, cell signaling in cardiac hypertrophy comprises a complex web of pathways with numerous interactions, and predicting how these interactions control the hypertrophic signal in each context is not achievable by only experiments or general computational models. To address this need, we developed an approach to bring together the experimental data of each context with a signaling network curated from literature to identify the main players of cardiac cells response in each context and attain the context-specific models of cardiac hypertrophy. By utilizing our approach, we identified the main regulators of cardiac hypertrophy in four important contexts. We developed a network model of β-adrenergic cardiac hypertrophy, and predicted and validated new interactions that regulate cardiac cells response in this context.

Ranolazine as an Alternative Therapy to Flecainide for SCN5A V411M Long QT Syndrome Type 3 Patients

The prolongation of the QT interval represents the main feature of the long QT syndrome (LQTS), a life-threatening genetic disease. The heterozygous SCN5A V411M mutation of the human sodium channel leads to a LQTS type 3 with severe proarrhythmic effects due to an increase in the late component of the sodium current (INaL). The two sodium blockers flecainide and ranolazine are equally recommended by the current 2015 ESC guidelines to treat patients with LQTS type 3 and persistently prolonged QT intervals. However, awareness of pro-arrhythmic effects of flecainide in LQTS type 3 patients arose upon the study of the SCN5A E1784K mutation. Regarding SCN5A V411M individuals, flecainide showed good results albeit in a reduced number of patients and no evidence supporting the use of ranolazine has ever been released. Therefore, we ought to compare the effect of ranolazine and flecainide in a SCN5A V411M model using an in-silico modeling and simulation approach. We collected clinical data of four patients. Then, we fitted four Markovian models of the human sodium current (INa) to experimental and clinical data. Two of them correspond to the wild type and the heterozygous SCN5A V411M scenarios, and the other two mimic the effects of flecainide and ranolazine on INa. Next, we inserted them into three isolated cell action potential (AP) models for endocardial, midmyocardial and epicardial cells and in a one-dimensional tissue model. The SCN5A V411M mutation produced a 15.9% APD90 prolongation in the isolated endocardial cell model, which corresponded to a 14.3% of the QT interval prolongation in a one-dimensional strand model, in keeping with clinical observations. Although with different underlying mechanisms, flecainide and ranolazine partially countered this prolongation at the isolated endocardial model by reducing the APD90 by 8.7 and 4.3%, and the QT interval by 7.2 and 3.2%, respectively. While flecainide specifically targeted the mutation-induced increase in peak INaL, ranolazine reduced it during the entire AP. Our simulations also suggest that ranolazine could prevent early afterdepolarizations triggered by the SCN5A V411M mutation during bradycardia, as flecainide. We conclude that ranolazine could be used to treat SCN5A V411M patients, specifically when flecainide is contraindicated.

The role of β-adrenergic system remodeling in human heart failure: A mechanistic investigation

β-adrenergic receptor antagonists (β-blockers) are extensively used to improve cardiac performance in heart failure (HF), but the electrical improvements with these clinical treatments are not fully understood. The aim of this study was to analyze the electrophysiological effects of β-adrenergic system remodeling in heart failure with reduced ejection fraction and the underlying mechanisms. We used a combined mathematical model that integrated β-adrenergic signaling with electrophysiology and calcium cycling in human ventricular myocytes. HF remodeling, both in the electrophysiological and signaling systems, was introduced to quantitatively analyze changes in electrophysiological properties due to the stimulation of β-adrenergic receptors in failing myocytes. We found that the inotropic effect of β-adrenergic stimulation was reduced in HF due to the altered Ca²⁺ dynamics resulting from the combination of structural, electrophysiological and signaling remodeling. Isolated cells showed proarrhythmic risk after sympathetic stimulation because early afterdepolarizations appeared, and the vulnerability was greater in failing myocytes. When analyzing coupled cells, β-adrenergic stimulation reduced transmural repolarization gradients between endocardium and epicardium in normal tissue, but was less effective at reducing these gradients after HF remodeling. The comparison of the selective activation of β-adrenergic isoforms revealed that the response to β₂-adrenergic receptors stimulation was blunted in HF while β₁-adrenergic receptors downstream effectors regulated most of the changes observed after sympathetic stimulation. In conclusion, this study was able to reproduce an altered β-adrenergic activity on failing myocytes and to explain the mechanisms involved. The derived predictions could help in the treatment of HF and guide in the design of future experiments.

Role of Oxidation-Dependent CaMKII Activation in the Genesis of Abnormal Action Potentials in Atrial Cardiomyocytes: A Simulation Study

Atrial fibrillation is a common cardiac arrhythmia with an increasing incidence rate. Particularly for the aging population, understanding the underlying mechanisms of atrial arrhythmia is important in designing clinical treatment. Recently, experiments have shown that atrial arrhythmia is associated with oxidative stress. In this study, an atrial cell model including oxidative-dependent Ca²⁺/calmodulin- (CaM-) dependent protein kinase II (CaMKII) activation was developed to explore the intrinsic mechanisms of atrial arrhythmia induced by oxidative stress. The simulation results showed that oxidative stress caused early afterdepolarizations (EADs) of action potentials by altering the dynamics of transmembrane currents and intracellular calcium cycling. Oxidative stress gradually elevated the concentration of calcium ions in the cytoplasm by enhancing the L-type Ca²⁺ current and sarcoplasmic reticulum (SR) calcium release. Owing to increased intracellular calcium concentration, the inward Na⁺/Ca²⁺ exchange current was elevated which slowed down the repolarization of the action potential. Thus, the action potential was prolonged and the L-type Ca²⁺ current was reactivated, resulting in the genesis of EAD. Furthermore, based on the atrial single-cell model, a two-dimensional (2D) ideal tissue model was developed to explore the effect of oxidative stress on the electrical excitation wave conduction in 2D tissue. Simulation results demonstrated that, under oxidative stress conditions, EAD hindered the conduction of electrical excitation and caused an unstable spiral wave, which could disrupt normal cardiac rhythm and cause atrial arrhythmia. This study showed the effects of excess reactive oxygen species on calcium cycling and action potential in atrial myocytes and provided insights regarding atrial arrhythmia induced by oxidative stress.

Quantitative analysis of variability in an integrated model of human ventricular electrophysiology and β-adrenergic signaling

In ventricular myocytes, stimulation of β-adrenergic receptors activates critical cardiac signaling pathways, leading to shorter action potentials and increased contraction strength during the "fight-or-flight" response. These changes primarily result, at the cellular level, from the coordinated phosphorylation of multiple targets by protein kinase A. Although mathematical models of the intracellular signaling downstream of β-adrenergic receptor activation have previously been described, only a limited number of studies have explored quantitative interactions between intracellular signaling and electrophysiology in human ventricular myocytes. Accordingly, our objective was to develop an integrative mathematical model of β-adrenergic receptor signaling, electrophysiology, and intracellular calcium (Ca²⁺) handling in the healthy human ventricular myocyte. We combined published mathematical models of intracellular signaling and electrophysiology, then calibrated the model results against voltage clamp data and physiological changes occurring after stimulation of β-adrenergic receptors with isoproterenol. We subsequently: (1) explored how molecular variability in different categories of model parameters translated into phenotypic variability; (2) identified the most important parameters determining physiological cellular outputs in the model before and after β-adrenergic receptor stimulation; and (3) investigated which molecular level alterations can produce a phenotype indicative of heart failure with preserved ejection fraction (HFpEF). Major results included: (1) variability in parameters that controlled intracellular signaling caused qualitatively different behavior than variability in parameters controlling ion transport pathways; (2) the most important model parameters determining action potential duration and intracellular Ca²⁺ transient amplitude were generally consistent before and after β-adrenergic receptor stimulation, except for a shift in the importance of K⁺ currents in determining action potential duration; and (3) decreased Ca²⁺ uptake into the sarcoplasmic reticulum, increased Ca²⁺ extrusion through Na⁺/Ca²⁺ exchanger and decreased Ca²⁺ leak from the sarcoplasmic reticulum may contribute to HFpEF. Overall, this study provided novel insight into the phenotypic consequences of molecular variability, and our integrated model may be useful in the design and interpretation of future experimental studies of interactions between β-adrenergic signaling and cellular physiology in human ventricular myocytes.

Mathematical model for β1-adrenergic regulation of the mouse ventricular myocyte contraction

The β1-adrenergic regulation of cardiac myocyte contraction plays an important role in regulating heart function. Activation of this system leads to an increased heart rate and stronger myocyte contraction. However, chronic stimulation of the β1-adrenergic signaling system can lead to cardiac hypertrophy and heart failure. To understand the mechanisms of action of β1-adrenoceptors, a mathematical model of cardiac myocyte contraction that includes the β1-adrenergic system was developed and studied. The model was able to simulate major experimental protocols for measurements of steady-state force-calcium relationships, cross-bridge release rate and force development rate, force-velocity relationship, and force redevelopment rate. It also reproduced quite well frequency and isoproterenol dependencies for intracellular Ca2+ concentration ([Ca2+]i) transients, total contraction force, and sarcomere shortening. The mathematical model suggested the mechanisms of increased contraction force and myocyte shortening on stimulation of β1-adrenergic receptors is due to phosphorylation of troponin I and myosin-binding protein C and increased [Ca2+]i transient resulting from activation of the β1-adrenergic signaling system. The model was used to simulate work-loop contractions and estimate the power during the cardiac cycle as well as the effects of 4-aminopyridine and tedisamil on the myocyte contraction. The developed mathematical model can be used further for simulations of contraction of ventricular myocytes from genetically modified mice and myocytes from mice with chronic cardiac diseases.

NEW & NOTEWORTHY A new mathematical model of mouse ventricular myocyte contraction that includes the β1-adrenergic system was developed. The model simulated major experimental protocols for myocyte contraction and predicted the effects of 4-aminopyridine and tedisamil on the myocyte contraction. The model also allowed for simulations of work-loop contractions and estimation of the power during the cardiac cycle.

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.

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.

Different paths, same destination: divergent action potential responses produce conserved cardiac fight-or-flight response in mouse and rabbit hearts

KEY POINTS

Cardiac electrophysiology and Ca²⁺ handling change rapidly during the fight-or-flight response to meet physiological demands. Despite dramatic differences in cardiac electrophysiology, the cardiac fight-or-flight response is highly conserved across species. In this study, we performed physiological sympathetic nerve stimulation (SNS) while optically mapping cardiac action potentials and intracellular Ca²⁺ transients in innervated mouse and rabbit hearts. Despite similar heart rate and Ca²⁺ handling responses between mouse and rabbit hearts, we found notable species differences in spatio-temporal repolarization dynamics during SNS. Species-specific computational models revealed that these electrophysiological differences allowed for enhanced Ca²⁺ handling (i.e. enhanced inotropy) in each species, suggesting that electrophysiological responses are fine-tuned across species to produce optimal cardiac fight-or-flight responses.

ABSTRACT

Sympathetic activation of the heart results in positive chronotropy and inotropy, which together rapidly increase cardiac output. The precise mechanisms that produce the electrophysiological and Ca²⁺ handling changes underlying chronotropic and inotropic responses have been studied in detail in isolated cardiac myocytes. However, few studies have examined the dynamic effects of physiological sympathetic nerve activation on cardiac action potentials (APs) and intracellular Ca²⁺ transients (CaTs) in the intact heart. Here, we performed bilateral sympathetic nerve stimulation (SNS) in fully innervated, Langendorff-perfused rabbit and mouse hearts. Dual optical mapping with voltage- and Ca²⁺ -sensitive dyes allowed for analysis of spatio-temporal AP and CaT dynamics. The rabbit heart responded to SNS with a monotonic increase in heart rate (HR), monotonic decreases in AP and CaT duration (APD, CaTD), and a monotonic increase in CaT amplitude. The mouse heart had similar HR and CaT responses; however, a pronounced biphasic APD response occurred, with initial prolongation (50.9 ± 5.1 ms at t = 0 s vs. 60.6 ± 4.1 ms at t = 15 s, P < 0.05) followed by shortening (46.5 ± 9.1 ms at t = 60 s, P = NS vs. t = 0). We determined the biphasic APD response in mouse was partly due to dynamic changes in HR during SNS and was exacerbated by β-adrenergic activation. Simulations with species-specific cardiac models revealed that transient APD prolongation in mouse allowed for greater and more rapid CaT responses, suggesting more rapid increases in contractility; conversely, the rabbit heart requires APD shortening to produce optimal inotropic responses. Thus, while the cardiac fight-or-flight response is highly conserved between species, the underlying mechanisms orchestrating these effects differ significantly.

A demonstration of modularity, reuse, reproducibility, portability and scalability for modeling and simulation of cardiac electrophysiology using Kepler Workflows

Multi-scale computational modeling is a major branch of computational biology as evidenced by the US federal interagency Multi-Scale Modeling Consortium and major international projects. It invariably involves specific and detailed sequences of data analysis and simulation, often with multiple tools and datasets, and the community recognizes improved modularity, reuse, reproducibility, portability and scalability as critical unmet needs in this area. Scientific workflows are a well-recognized strategy for addressing these needs in scientific computing. While there are good examples if the use of scientific workflows in bioinformatics, medical informatics, biomedical imaging and data analysis, there are fewer examples in multi-scale computational modeling in general and cardiac electrophysiology in particular. Cardiac electrophysiology simulation is a mature area of multi-scale computational biology that serves as an excellent use case for developing and testing new scientific workflows. In this article, we develop, describe and test a computational workflow that serves as a proof of concept of a platform for the robust integration and implementation of a reusable and reproducible multi-scale cardiac cell and tissue model that is expandable, modular and portable. The workflow described leverages Python and Kepler-Python actor for plotting and pre/post-processing. During all stages of the workflow design, we rely on freely available open-source tools, to make our workflow freely usable by scientists.

We present a computational workflow as a proof of concept for integration and implementation of a reusable and reproducible cardiac multi-scale electrophysiology model that is expandable, modular and portable. This framework enables scientists to create intuitive, user-friendly and flexible end-to-end automated scientific workflows using a graphical user interface. Kepler is an advanced open-source platform that supports multiple models of computation. The underlying workflow engine handles scalability, provenance, reproducibility aspects of the code, performs orchestration of data flow, and automates execution on heterogeneous computing resources. One of the main advantages of workflow utilization is the integration of code written in multiple languages Standardization occurs at the interfaces of the workflow elements and allows for general applications and easy comparison and integration of code from different research groups or even multiple programmers coding in different languages for various purposes from the same group. A workflow driven problem-solving approach enables domain scientists to focus on resolving the core science questions, and delegates the computational and process management burden to the underlying Workflow. The workflow driven approach allows scaling the computational experiment with distributed data-parallel execution on multiple computing platforms, such as, HPC resources, GPU clusters, Cloud etc. The workflow framework tracks software version information along with hardware information to allow users an opportunity to trace any variation in workflow outcome to the system configurations.

Enhanced Depolarization Drive in Failing Rabbit Ventricular Myocytes: Calcium-Dependent and β-Adrenergic Effects on Late Sodium, L-Type Calcium, and Sodium-Calcium Exchange Currents

BACKGROUND

Heart failure (HF) is characterized by electrophysiological remodeling resulting in increased risk of cardiac arrhythmias. Previous reports suggest that elevated inward ionic currents in HF promote action potential (AP) prolongation, increased short-term variability of AP repolarization, and delayed afterdepolarizations. However, the underlying changes in late Na+ current (INaL), L-type Ca2+ current, and NCX (Na+/Ca2+ exchanger) current are often measured in nonphysiological conditions (square-pulse voltage clamp, slow pacing rates, exogenous Ca2+ buffers).

METHODS

We measured the major inward currents and their Ca2+- and β-adrenergic dependence under physiological AP clamp in rabbit ventricular myocytes in chronic pressure/volume overload-induced HF (versus age-matched control).

RESULTS

AP duration and short-term variability of AP repolarization were increased in HF, and importantly, inhibition of INaL decreased both parameters to the control level. INaL was slightly increased in HF versus control even when intracellular Ca2+ was strongly buffered. But under physiological AP clamp with normal Ca2+ cycling, INaL was markedly upregulated in HF versus control (dependent largely on CaMKII [Ca2+/calmodulin-dependent protein kinase II] activity). β-Adrenergic stimulation (often elevated in HF) further enhanced INaL. L-type Ca2+ current was decreased in HF when Ca2+ was buffered, but CaMKII-mediated Ca2+-dependent facilitation upregulated physiological L-type Ca2+ current to the control level. Furthermore, L-type Ca2+ current response to β-adrenergic stimulation was significantly attenuated in HF. Inward NCX current was upregulated at phase 3 of AP in HF when assessed by combining experimental data and computational modeling.

CONCLUSIONS

Our results suggest that CaMKII-dependent upregulation of INaL in HF significantly contributes to AP prolongation and increased short-term variability of AP repolarization, which may lead to increased arrhythmia propensity, and is further exacerbated by adrenergic stress.

Methyl-β-Cyclodextrin Impairs the Phosphorylation of the β₂ Subunit of L-Type Calcium Channels and Cytosolic Calcium Homeostasis in Mature Cerebellar Granule Neurons

The activation of L-type calcium channels (LTCCs) prevents cerebellar granule neurons (CGNs) from entering low-K⁺-induced apoptosis. In previous works, we showed that LTCCs are largely associated with caveolin-1-rich lipid rafts in the CGN plasma membrane. In this work, we show that protein kinase A (PKA) and calmodulin-dependent protein kinase II (CaMK-II) are associated with caveolin-1-rich lipid rafts of mature CGNs, and we further show that treatment with the cholesterol-trapping and lipid raft-disrupting agent methyl-β-cyclodextrin decreases the phosphorylation level of the LTCC β₂ subunit and the steady-state calcium concentration in neuronal somas ([Ca²⁺]_i) to values close to those measured in 5 mM KCl proapoptotic conditions. These effects correlate with the effects produced by a short (15 min) treatment of CGNs with H-89 and KN-93—inhibitors of PKA and CaMK-II, respectively—in 25 mM KCl medium. Moreover, only a 15 min incubation of CGNs with H-89 produces about a 90% inhibition of the calcium entry that would normally occur through LTCCs to increase [Ca²⁺]_i upon raising the extracellular K⁺ from 5 to 25 mM, i.e., from proapoptotic to survival conditions. In conclusion, the results of this work suggest that caveolin-1-rich lipid rafts play a major role in the control of the PKA- and CaMK-II-induced phosphorylation level of the LTCC β₂ subunit, thus preventing CGNs from entering apoptosis.

Mechanical Strain Regulates Osteoblast Proliferation Through Ca2+-CaMK-CREB Signal Pathway

Objective To investigate the effects of mechanical strain on Ca²⁺-calmodulin dependent kinase (CaMK)-cAMP response element binding protein (CREB) signal pathway and proliferation of osteoblasts.Methods Using a four-point bending device, MC3T3-E1 cells were exposed to mechanical tensile strains of 2500 µs and 5000 µs at 0.5 Hz respectively. The intracellular free Ca²⁺ ([Ca²⁺]i) concentration and calmodulin activity were assayed by fluorospectrophotometry, CaMK II β, CREB, and phosphorylated (activated) CREB (p-CREB) were assessed by Western blot, and cells proliferation was assayed with MTT. Pretreatment with verapamil was carried out to block Ca²⁺ channel, and inhibitor U73122 was used to inhibit phospholipase C (PLC).Results Mechanical strains of 2500 µs and 5000 µs for 1 to 10 minutes both increased [Ca²⁺]i level of the cells. The 2500 µs strain, a periodicity of 1 h/d for 3 days, activated calmodulin, elevated protein levels of CaMK II β and p-CREB, and promoted cells proliferation, which were attenuated by pretreatment of verapamil or U73122. The effects of 5000 µs strain on calmodulin, CaMK II β, p-CREB and proliferation were contrary to 2500 µs strain.Conclusion The mechanical strain regulates osteoblasts proliferation through Ca²⁺-CaMK-CREB signal pathway via Ca²⁺ channel and PLC/IP₃ transduction cascades.

A new myofilament contraction model with ATP consumption for ventricular cell model

A new contraction model of cardiac muscle was developed by combining previously described biochemical and biophysical models. The biochemical component of the new contraction model represents events in the presence of Ca2+-crossbridge attachment and power stroke following inorganic phosphate release, detachment evoked by the replacement of ADP by ATP, ATP hydrolysis, and recovery stroke. The biophysical component focuses on Ca2+ activation and force (F b) development assuming an equivalent crossbridge. The new model faithfully incorporates the major characteristics of the biochemical and biophysical models, such as F b activation by transient Ca2+ ([Ca2+]-F b), [Ca2+]-ATP hydrolysis relations, sarcomere length-F b, and F b recovery after jumps in length under the isometric mode and upon sarcomere shortening after a rapid release of mechanical load under the isotonic mode together with the load-velocity relationship. ATP consumption was obtained for all responses. When incorporated in a ventricular cell model, the contraction model was found to share approximately 60% of the total ATP usage in the cell model.

A Heart for Diversity: Simulating Variability in Cardiac Arrhythmia Research

In cardiac electrophysiology, there exist many sources of inter- and intra-personal variability. These include variability in conditions and environment, and genotypic and molecular diversity, including differences in expression and behavior of ion channels and transporters, which lead to phenotypic diversity (e.g., variable integrated responses at the cell, tissue, and organ levels). These variabilities play an important role in progression of heart disease and arrhythmia syndromes and outcomes of therapeutic interventions. Yet, the traditional in silico framework for investigating cardiac arrhythmias is built upon a parameter/property-averaging approach that typically overlooks the physiological diversity. Inspired by work done in genetics and neuroscience, new modeling frameworks of cardiac electrophysiology have been recently developed that take advantage of modern computational capabilities and approaches, and account for the variance in the biological data they are intended to illuminate. In this review, we outline the recent advances in statistical and computational techniques that take into account physiological variability, and move beyond the traditional cardiac model-building scheme that involves averaging over samples from many individuals in the construction of a highly tuned composite model. We discuss how these advanced methods have harnessed the power of big (simulated) data to study the mechanisms of cardiac arrhythmias, with a special emphasis on atrial fibrillation, and improve the assessment of proarrhythmic risk and drug response. The challenges of using in silico approaches with variability are also addressed and future directions are proposed.

Quantitative systems models illuminate arrhythmia mechanisms in heart failure: Role of the Na+ -Ca2+ -Ca2+ /calmodulin-dependent protein kinase II-reactive oxygen species feedback

Quantitative systems modeling aims to integrate knowledge in different research areas with models describing biological mechanisms and dynamics to gain a better understanding of complex clinical syndromes. Heart failure (HF) is a chronic complex cardiac disease that results from structural or functional disorders impairing the ability of the ventricle to fill with or eject blood. Highly interactive and dynamic changes in mechanical, structural, neurohumoral, metabolic, and electrophysiological properties collectively predispose the failing heart to cardiac arrhythmias, which are responsible for about a half of HF deaths. Multiscale cardiac modeling and simulation integrate structural and functional data from HF experimental models and patients to improve our mechanistic understanding of this complex arrhythmia syndrome. In particular, they allow investigating how disease-induced remodeling alters the coupling of electrophysiology, Ca²⁺ and Na⁺ handling, contraction, and energetics that lead to rhythm derangements. The Ca²⁺ /calmodulin-dependent protein kinase II, which expression and activity are enhanced in HF, emerges as a critical hub that modulates the feedbacks between these various subsystems and promotes arrhythmogenesis. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Mechanistic Models Models of Systems Properties and Processes > Cellular Models Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

Cardiophysiology Illustrated by Comparing Ventricular Volumes in Healthy Adult Males and Females

Recent advances in cardiac imaging techniques have substantially contributed to a growing interest in the analysis of global cardiac chamber dimensions and regional myocardial deformation. During the cardiac cycle, ventricular luminal volume varies due to the contraction process, which also confers a shape change including substantial alteration of long axis length, as well as rotation of the base compared to the apex. Local deformation can be assessed by strain (rate) analysis. Reviewing the present literature, it must be concluded that there is no single metric available to comprehensively characterize ventricular function. Every candidate advanced thus far has been found to incompletely reflect ventricular performance. This observation is not surprising in view of the complexity of the cardiac pump system. Additionally, sex-specific modifiers may play a role. More than three decades ago, it was shown that on average the ventricular volume is smaller in healthy women compared to matched males. Therefore, the present contribution concerns the interpretation of data derived from the healthy heart in both men and women. Starting from the classical Starling concept, we apply a simple mathematical transformation which permits an insightful representation of ventricular mechanics. Relating end-systolic volume (ESV) to end-diastolic volume creates the ventricular volume regulation graph which features the pertinent working point of an individual heart. This fundamental approach illustrates why certain proposed performance indexes cannot individually reveal the essence of ventricular systolic function. We demonstrate that particular metrics are highly interconnected and just tell us the same story in a different disguise. It is imperative to understand which associations exist and if they expectedly are (nearly) linear or frankly nonlinear. Notably, ejection fraction (EF) is primarily determined by ESV, while in turn EF is not much different from ventriculo-arterial coupling (VAC). Insight into cardiac function is promoted by identification of the paramount/essential components involved. The smaller ESV (p < 0.0001) implies that EF is higher in women and may also have consequences for VAC.

Estimating the probabilities of rare arrhythmic events in multiscale computational models of cardiac cells and tissue

Ectopic heartbeats can trigger reentrant arrhythmias, leading to ventricular fibrillation and sudden cardiac death. Such events have been attributed to perturbed Ca2+ handling in cardiac myocytes leading to spontaneous Ca2+ release and delayed afterdepolarizations (DADs). However, the ways in which perturbation of specific molecular mechanisms alters the probability of ectopic beats is not understood. We present a multiscale model of cardiac tissue incorporating a biophysically detailed three-dimensional model of the ventricular myocyte. This model reproduces realistic Ca2+ waves and DADs driven by stochastic Ca2+ release channel (RyR) gating and is used to study mechanisms of DAD variability. In agreement with previous experimental and modeling studies, key factors influencing the distribution of DAD amplitude and timing include cytosolic and sarcoplasmic reticulum Ca2+ concentrations, inwardly rectifying potassium current (IK1) density, and gap junction conductance. The cardiac tissue model is used to investigate how random RyR gating gives rise to probabilistic triggered activity in a one-dimensional myocyte tissue model. A novel spatial-average filtering method for estimating the probability of extreme (i.e. rare, high-amplitude) stochastic events from a limited set of spontaneous Ca2+ release profiles is presented. These events occur when randomly organized clusters of cells exhibit synchronized, high amplitude Ca2+ release flux. It is shown how reduced IK1 density and gap junction coupling, as observed in heart failure, increase the probability of extreme DADs by multiple orders of magnitude. This method enables prediction of arrhythmia likelihood and its modulation by alterations of other cellular mechanisms.

Ca2+/calmodulin-dependent kinase II-dependent regulation of atrial myocyte late Na+ current, Ca2+ cycling, and excitability: a mathematical modeling study

Atrial fibrillation (AF) affects more than three million people per year in the United States and is associated with high morbidity and mortality. Both electrical and structural remodeling contribute to AF, but the molecular pathways underlying AF pathogenesis are not well understood. Recently, a role for Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) in the regulation of persistent "late" Na⁺ current ( I_Na,L) has been identified. Although I_Na,L inhibition is emerging as a potential antiarrhythmic strategy in patients with AF, little is known about the mechanism linking I_Na,L to atrial arrhythmogenesis. A computational approach was used to test the hypothesis that increased CaMKII-activated I_Na,L in atrial myocytes disrupts Ca²⁺ homeostasis, promoting arrhythmogenic afterdepolarizations. Dynamic CaMKII activity and regulation of multiple downstream targets [ I_Na,L, L-type Ca²⁺ current, phospholamban, and the ryanodine receptor sarcoplasmic reticulum Ca²⁺-release channel (RyR2)] were incorporated into an existing well-validated computational model of the human atrial action potential. Model simulations showed that constitutive CaMKII-dependent phosphorylation of Na_v1.5 and the subsequent increase in I_Na,L effectively disrupt intracellular atrial myocyte ion homeostasis and CaMKII signaling. Specifically, increased I_Na,L promotes intracellular Ca²⁺ overload via forward-mode Na⁺/Ca²⁺ exchange activity, which greatly increases RyR2 open probability beyond that observed for CaMKII-dependent phosphorylation of RyR2 alone. Increased I_Na,L promotes atrial myocyte repolarization defects (afterdepolarizations and alternans) in the setting of acute β-adrenergic stimulation. We anticipate that our modeling efforts will help identify new mechanisms for atrial Na_V1.5 regulation with direct relevance for human AF. NEW & NOTEWORTHY Here, we present a novel computational model to study the effects of late Na⁺ current ( I_Na,L) in human atrial myocytes. Simulations predict that I_Na,L promotes intracellular accumulation of Ca²⁺, with subsequent dysregulation of Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) signaling and ryanodine receptor 2-mediated Ca²⁺ release. Although I_Na,L plays a small role in regulating atrial myocyte excitability at baseline, CaMKII-dependent enhancement of the current promoted arrhythmogenic dynamics. Listen to this article's corresponding podcast at http://ajpheart.podbean.com/e/camkii-dependent-regulation-of-atrial-late-sodium-current-and-excitability/ .

Mechanisms of the cyclic nucleotide cross-talk signaling network in cardiac L-type calcium channel regulation

Regulation of L-type Calcium (Ca2+) Channel (LCC) gating is critical to shaping the cardiac action potential (AP) and triggering the initiation of excitation-contraction (EC) coupling in cardiac myocytes. The cyclic nucleotide (cN) cross-talk signaling network, which encompasses the β-adrenergic and the Nitric Oxide (NO)/cGMP/Protein Kinase G (PKG) pathways and their interaction (cross-talk) through distinctively-regulated phosphodiesterase isoenzymes (PDEs), regulates LCC current via Protein Kinase A- (PKA) and PKG-mediated phosphorylation. Due to the tightly-coupled and intertwined biochemical reactions involved, it remains to be clarified how LCC gating is regulated by the signaling network from receptor to end target. In addition, the large number of EC coupling-related phosphorylation targets of PKA and PKG makes it difficult to quantify and isolate changes in L-type Ca2+ current (ICaL) responses regulated by the signaling network. We have developed a multi-scale, biophysically-detailed computational model of LCC regulation by the cN signaling network that is supported by experimental data. LCCs are modeled with functionally distinct PKA- and PKG-phosphorylation dependent gating modes. The model exhibits experimentally observed single channel characteristics, as well as whole-cell LCC currents upon activation of the cross-talk signaling network. Simulations show 1) redistribution of LCC gating modes explains changes in whole-cell current under various stimulation scenarios of the cN cross-talk network; 2) NO regulation occurs via potentiation of a gating mode characterized by prolonged closed times; and 3) due to compensatory actions of cross-talk and antagonizing functions of PKA- and PKG-mediated phosphorylation of LCCs, the effects of individual inhibitions of PDEs 2, 3, and 4 on ICaL are most pronounced at low levels of β-adrenergic stimulation. Simulations also delineate the contribution of the following two mechanisms to overall LCC regulation, which have otherwise been challenging to distinguish: 1) regulation of PKA and PKG activation via cN cross-talk (Mechanism 1); and 2) LCC interaction with activated PKA and PKG (Mechanism 2). These results provide insights into how cN signals transduced via the cN cross-talk signaling network are integrated via LCC regulation in the heart.

Role of CaMKII and PKA in Early Afterdepolarization of Human Ventricular Myocardium Cell: A Computational Model Study

Early afterdepolarization (EAD) plays an important role in arrhythmogenesis. Many experimental studies have reported that Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) and β-adrenergic signaling pathway are two important regulators. In this study, we developed a modified computational model of human ventricular myocyte to investigate the combined role of CaMKII and β-adrenergic signaling pathway on the occurrence of EADs. Our simulation results showed that (1) CaMKII overexpression facilitates EADs through the prolongation of late sodium current's (I_NaL) deactivation progress; (2) the combined effect of CaMKII overexpression and activation of β-adrenergic signaling pathway further increases the risk of EADs, where EADs could occur at shorter cycle length (2000 ms versus 4000 ms) and lower rapid delayed rectifier K⁺ current (I_Kr) blockage (77% versus 85%). In summary, this study computationally demonstrated the combined role of CaMKII and β-adrenergic signaling pathway on the occurrence of EADs, which could be useful for searching for therapy strategies to treat EADs related arrhythmogenesis.

Anti-arrhythmic strategies for atrial fibrillation: The role of computational modeling in discovery, development, and optimization

Atrial fibrillation (AF), the most common cardiac arrhythmia, is associated with increased risk of cerebrovascular stroke, and with several other pathologies, including heart failure. Current therapies for AF are targeted at reducing risk of stroke (anticoagulation) and tachycardia-induced cardiomyopathy (rate or rhythm control). Rate control, typically achieved by atrioventricular nodal blocking drugs, is often insufficient to alleviate symptoms. Rhythm control approaches include antiarrhythmic drugs, electrical cardioversion, and ablation strategies. Here, we offer several examples of how computational modeling can provide a quantitative framework for integrating multiscale data to: (a) gain insight into multiscale mechanisms of AF; (b) identify and test pharmacological and electrical therapy and interventions; and (c) support clinical decisions. We review how modeling approaches have evolved and contributed to the research pipeline and preclinical development and discuss future directions and challenges in the field.

A computational model predicts adjunctive pharmacotherapy for cardiac safety via selective inhibition of the late cardiac Na current

Background

The QT interval is a phase of the cardiac cycle that corresponds to action potential duration (APD) including cellular repolarization (T-wave). In both clinical and experimental settings, prolongation of the QT interval of the electrocardiogram (ECG) and related proarrhythmia have been so strongly associated that a prolonged QT interval is largely accepted as surrogate marker for proarrhythmia. Accordingly, drugs that prolong the QT interval are not considered for further preclinical development resulting in removal of many promising drugs from development. While reduction of drug interactions with hERG is an important goal, there are promising means to mitigate hERG block. Here, we examine one possibility and test the hypothesis that selective inhibition of the cardiac late Na current (I_NaL) by the novel compound GS-458967 can suppress proarrhythmic markers.

Methods and results

New experimental data has been used to calibrate I_NaL in the Soltis-Saucerman computationally based model of the rabbit ventricular action potential to study effects of GS-458967 on I_NaL during the rabbit ventricular AP. We have also carried out systematic in silico tests to determine if targeted block of I_NaL would suppress proarrhythmia markers in ventricular myocytes described by TRIaD: Triangulation, Reverse use dependence, beat-to-beat Instability of action potential duration, and temporal and spatial action potential duration Dispersion.

Conclusions

Our computer modeling approach based on experimental data, yields results that suggest that selective inhibition of I_NaL modifies all TRIaD related parameters arising from acquired Long-QT Syndrome, and thereby reduced arrhythmia risk. This study reveals the potential for adjunctive pharmacotherapy via targeted block of I_NaL to mitigate proarrhythmia risk for drugs with significant but unintended off-target hERG blocking effects.

Mechanistic Investigation of the Arrhythmogenic Role of Oxidized CaMKII in the Heart

Oxidative stress and calcium (Ca(2+))/calmodulin (CaM)-dependent protein kinase II (CaMKII) both play important roles in the pathogenesis of cardiac disease. Although the pathophysiological relevance of reactive oxygen species (ROS) and CaMKII has been appreciated for some time, recent work has shown that ROS can directly oxidize CaMKII, leading to its persistent activity and an increase of the likelihood of cellular arrhythmias such as early afterdepolarizations (EADs). Because CaMKII modulates the function of many proteins involved in excitation-contraction coupling, elucidation of its role in cardiac function, in both healthy and oxidative stress conditions, is challenging. To investigate this role, we have developed a model of CaMKII activation that includes both the phosphorylation-dependent and the newly identified oxidation-dependent activation pathways. This model is incorporated into our previous local-control model of the cardiac myocyte that describes excitation-contraction coupling via stochastic simulation of individual Ca(2+) release units and CaMKII-mediated phosphorylation of L-type Ca(2+) channels (LCCs), ryanodine receptors and sodium (Na(+)) channels. The model predicts the experimentally measured slow-rate dependence of H2O2-induced EADs. Upon increased H2O2, simulations suggest that selective activation of late Na(+) current (INaL), although it prolongs action potential duration, is not by itself sufficient to produce EADs. Similar results are obtained if CaMKII effects on LCCs and ryanodine receptors are considered separately. However, EADs emerge upon simultaneous activation of both LCCs and Na(+) channels. Further modeling results implicate activation of the Na(+)-Ca(2+) exchanger (NCX) as an important player in the generation of EADs. During bradycardia, the emergence of H2O2-induced EADs was correlated with a shift in the timing of NCX current reversal toward the plateau phase earlier in the action potential. Using the timing of NCX current reversal as an indicator event for EADs, the model identified counterintuitive ionic changes-difficult to experimentally dissect-that have the greatest influence on ROS-related arrhythmia propensity.

A novel heterozygous mutation in cardiac calsequestrin causes autosomal dominant catecholaminergic polymorphic ventricular tachycardia

BACKGROUND

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is a lethal inherited arrhythmia syndrome characterized by adrenergically stimulated ventricular tachycardia. Mutations in the cardiac ryanodine receptor gene (RYR2) cause an autosomal dominant form of CPVT, while mutations in the cardiac calsequestrin 2 gene (CASQ2) cause an autosomal recessive form.

OBJECTIVE

The aim of this study was to clinically and genetically evaluate a large family with severe autosomal dominant CPVT.

METHODS

Clinical evaluation of family members was performed, including detailed history, physical examination, electrocardiogram, exercise stress test, and autopsy review of decedents. We performed genome-wide linkage analysis in 12 family members and exome sequencing in 2 affected family members. In silico models of mouse and rabbit myocyte electrophysiology were used to predict potential disease mechanisms.

RESULTS

Severe CPVT with dominant inheritance in 6 members was diagnosed in a large family with 2 sudden deaths, 2 resuscitated cardiac arrests, and multiple appropriate implantable cardioverter-defibrillator shocks. A comprehensive analysis of cardiac arrhythmia genes did not reveal a pathogenic variant. Exome sequencing identified a novel heterozygous missense variant in CASQ2 (Lys180Arg) affecting a highly conserved residue, which cosegregated with disease and was absent in unaffected family members. Genome-wide linkage analysis confirmed a single linkage peak at the CASQ2 locus (logarithm of odds ratio score 3.01; θ = 0). Computer simulations predicted that haploinsufficiency was unlikely to cause the severe CPVT phenotype and suggested a dominant negative mechanism.

CONCLUSION

We show for the first time that a variant in CASQ2 causes autosomal dominant CPVT. Genetic testing in dominant CPVT should include screening for heterozygous CASQ2 variants.

Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention

Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.

Molecular pathogenesis of long QT syndrome type 2

The molecular mechanisms underlying congenital long QT syndrome (LQTS) are now beginning to be understood. New insights into the etiology and therapeutic strategies are emerging from heterologous expression studies of LQTS-linked mutant proteins, as well as inducible pluripotent stem cell derived cardiomyocytes (iPSC-CMs) from LQTS patients. This review focuses on the major molecular mechanism that underlies LQTS type 2 (LQT2). LQT2 is caused by loss of function (LOF) mutations in KCNH2 (also known as the human Ether-à-go-go-Related Gene or hERG). Most LQT2-linked mutations are missense mutations and functional studies suggest that ~90% of them disrupt the intracellular transport (trafficking) of KCNH2-encoded Kv11.1 proteins to the cell membrane. Trafficking deficient LQT2 mutations disrupt Kv11.1 protein folding and misfolded Kv11.1 proteins are retained in the endoplasmic reticulum (ER) until they are degraded in the ER associated degradation pathway (ERAD). This review focuses on the quality control mechanisms in the ER that contribute to the folding and ERAD of Kv11.1 proteins; the mechanism for ER export of Kv11.1 proteins in the secretory pathway; different subclasses of trafficking deficient LQT2 mutations; and strategies being developed to mitigate or correct trafficking deficient LQT2-related phenotypes.