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

In silico prediction of drug therapy in catecholaminergic polymorphic ventricular tachycardia

Key points

The mechanism of therapeutic efficacy of flecainide for catecholaminergic polymorphic ventricular tachycardia (CPVT) is unclear.

Model predictions suggest that Na⁺ channel effects are insufficient to explain flecainide efficacy in CPVT.

This study represents a first step toward predicting therapeutic mechanisms of drug efficacy in the setting of CPVT and then using these mechanisms to guide modelling and simulation to predict alternative drug therapies.

Abstract

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is an inherited arrhythmia syndrome characterized by fatal ventricular arrhythmias in structurally normal hearts during β‐adrenergic stimulation. Current treatment strategies include β‐blockade, flecainide and ICD implementation – none of which is fully effective and each comes with associated risk. Recently, flecainide has gained considerable interest in CPVT treatment, but its mechanism of action for therapeutic efficacy is unclear. In this study, we performed in silico mutagenesis to construct a CPVT model and then used a computational modelling and simulation approach to make predictions of drug mechanisms and efficacy in the setting of CPVT. Experiments were carried out to validate model results. Our simulations revealed that Na⁺ channel effects are insufficient to explain flecainide efficacy in CPVT. The pure Na⁺ channel blocker lidocaine and the antianginal ranolazine were additionally tested and also found to be ineffective. When we tested lower dose combination therapy with flecainide, β‐blockade and CaMKII inhibition, our model predicted superior therapeutic efficacy than with flecainide monotherapy. Simulations indicate a polytherapeutic approach may mitigate side‐effects and proarrhythmic potential plaguing CPVT pharmacological management today. Importantly, our prediction of a novel polytherapy for CPVT was confirmed experimentally. Our simulations suggest that flecainide therapeutic efficacy in CPVT is unlikely to derive from primary interactions with the Na⁺ channel, and benefit may be gained from an alternative multi‐drug regimen.

The mechanism of therapeutic efficacy of flecainide for catecholaminergic polymorphic ventricular tachycardia (CPVT) is unclear.

Model predictions suggest that Na⁺ channel effects are insufficient to explain flecainide efficacy in CPVT.

This study represents a first step toward predicting therapeutic mechanisms of drug efficacy in the setting of CPVT and then using these mechanisms to guide modelling and simulation to predict alternative drug therapies.

Cardiac voltage-gated calcium channel macromolecular complexes

Over the past 20 years, a new field of research, called channelopathies, investigating diseases caused by ion channel dysfunction has emerged. Cardiac ion channels play an essential role in the generation of the cardiac action potential. Investigators have largely determined the physiological roles of different cardiac ion channels, but little is known about the molecular determinants of their regulation. The voltage-gated calcium channel Ca(v)1.2 shapes the plateau phase of the cardiac action potential and allows the influx of calcium leading to cardiomyocyte contraction. Studies suggest that the regulation of Ca(v)1.2 channels is not uniform in working cardiomyocytes. The notion of micro-domains containing Ca(v)1.2 channels and different calcium channel interacting proteins, called macro-molecular complex, has been proposed to explain these observations. The objective of this review is to summarize the currently known information on the Ca(v)1.2 macromolecular complexes in the cardiac cell and discuss their implication in cardiac function and disorder. This article is part of a Special Issue entitled: Cardiomyocyte Biology: Integration of Developmental and Environmental Cues in the Heart edited by Marcus Schaub and Hughes Abriel.

Modeling calcium regulation of contraction, energetics, signaling, and transcription in the cardiac myocyte

Calcium (Ca(2+)) plays many important regulatory roles in cardiac muscle cells. In the initial phase of the action potential, influx of Ca(2+) through sarcolemmal voltage-gated L-type Ca(2+) channels (LCCs) acts as a feed-forward signal that triggers a large release of Ca(2+) from the junctional sarcoplasmic reticulum (SR). This Ca(2+) drives heart muscle contraction and pumping of blood in a process known as excitation-contraction coupling (ECC). Triggered and released Ca(2+) also feed back to inactivate LCCs, attenuating the triggered Ca(2+) signal once release has been achieved. The process of ECC consumes large amounts of ATP. It is now clear that in a process known as excitation-energetics coupling, Ca(2+) signals exert beat-to-beat regulation of mitochondrial ATP production that closely couples energy production with demand. This occurs through transport of Ca(2+) into mitochondria, where it regulates enzymes of the tricarboxylic acid cycle. In excitation-signaling coupling, Ca(2+) activates a number of signaling pathways in a feed-forward manner. Through effects on their target proteins, these interconnected pathways regulate Ca(2+) signals in complex ways to control electrical excitability and contractility of heart muscle. In a process known as excitation-transcription coupling, Ca(2+) acting primarily through signal transduction pathways also regulates the process of gene transcription. Because of these diverse and complex roles, experimentally based mechanistic computational models are proving to be very useful for understanding Ca(2+) signaling in the cardiac myocyte.

Calmodulin kinase II inhibitor regulates calcium homeostasis changes caused by acute β-adrenergic receptor agonist stimulation in mouse ventricular myocytes

Ca(2+)/calmodulin-dependent kinase II (CaMKII) is an important regulatory molecule under chronic β-adrenergic receptor agonist stimulation but cardiac diseases also occur when β-adrenergic elevated acutely in the circulation, of which the most harmful is lethal arrhythmia. The purpose of this study was to explore the effects of acute isoproterenol (ISO) stimulation on intracellular calcium handling and evaluate whether CaMKII inhibitor may change the effects caused by isoproterenol. Mouse ventricular myocytes were acutely isolated by enzymatic method and divided into four groups: control group, ISO group, KN-93 group, ISO + KN-93 group. The whole-cell patch clamp was used to study the effect of ISO and KN-93 on L-type calcium current (IL-Ca) in isolated ventricular myocytes. The technology of laser scanning confocal microscopy was used to record cardiomyocyte calcium transients after ISO and KN-93 were used. ISO significantly increased current density of IL-Ca (p < 0.01) and decreased the half activation voltage (p < 0.01), half inactivation voltage (p < 0.01), and the recovery time constant (p < 0.01). In the presence of CaMKII inhibitor, KN-93 decreased the increased current density of IL-Ca (p < 0.05), increased the reduced half activation voltage caused by ISO (p < 0.01), and prolonged the shortened recovery time constant caused by ISO (p < 0.01). In addition, KN-93 alone can change the activation, inactivation, and recovery kinetics of L-type calcium channels. Moreover, ISO significantly increased the Ca(2+) transient amplitude during both stimulation frequencies (0.5 Hz: p < 0.01, 1 Hz: p < 0.01) and was easy to induce calcium disorders; in the presence of KN-93, these changes were weakened (0.5 Hz: p < 0.05, 1 Hz: p < 0.05). Therefore, changes of the calcium homeostasis in cardiomyocytes caused by ISO can be adjusted by KN-93, thus KN-93 plays a vital role in regulating calcium homeostasis changes caused by ISO.

Lessons learned from multi-scale modeling of the failing heart

Heart failure constitutes a major public health problem worldwide. Affected patients experience a number of changes in the electrical function of the heart that predispose to potentially lethal cardiac arrhythmias. Due to the multitude of electrophysiological changes that may occur during heart failure, the scientific literature is complex and sometimes ambiguous, perhaps because these findings are highly dependent on the etiology, the stage of heart failure, and the experimental model used to study these changes. Nevertheless, a number of common features of failing hearts have been documented. Prolongation of the action potential (AP) involving ion channel remodeling and alterations in calcium handling have been established as the hallmark characteristics of myocytes isolated from failing hearts. Intercellular uncoupling and fibrosis are identified as major arrhythmogenic factors. Multi-scale computational simulations are a powerful tool that complements experimental and clinical research. The development of biophysically detailed computer models of single myocytes and cardiac tissues has contributed greatly to our understanding of processes underlying excitation and repolarization in the heart. The electrical, structural, and metabolic remodeling that arises in cardiac tissues during heart failure has been addressed from different computational perspectives to further understand the arrhythmogenic substrate. This review summarizes the contributions from computational modeling and simulation to predict the underlying mechanisms of heart failure phenotypes and their implications for arrhythmogenesis, ranging from the cellular level to whole-heart simulations. The main aspects of heart failure are presented in several related sections. An overview of the main electrophysiological and structural changes that have been observed experimentally in failing hearts is followed by the description and discussion of the simulation work in this field at the cellular level, and then in 2D and 3D cardiac structures. The implications for arrhythmogenesis in heart failure are also discussed including therapeutic measures, such as drug effects and cardiac resynchronization therapy. Finally, the future challenges in heart failure modeling and simulation will be discussed.

Real-time relationship between PKA biochemical signal network dynamics and increased action potential firing rate in heart pacemaker cells: Kinetics of PKA activation in heart pacemaker cells

cAMP-PKA protein kinase is a key nodal signaling pathway that regulates a wide range of heart pacemaker cell functions. These functions are predicted to be involved in regulation of spontaneous action potential (AP) generation of these cells. Here we investigate if the kinetics and stoichiometry of increase in PKA activity match the increase in AP firing rate in response to β-adrenergic receptor (β-AR) stimulation or phosphodiesterase (PDE) inhibition, that alters the AP firing rate of heart sinoatrial pacemaker cells. In cultured adult rabbit pacemaker cells infected with an adenovirus expressing the FRET sensor AKAR3, the EC50 in response to graded increases in the intensity of β-AR stimulation (by Isoproterenol) the magnitude of the increases in PKA activity and the spontaneous AP firing rate were similar (0.4±0.1nM vs. 0.6±0.15nM, respectively). Moreover, the kinetics (t1/2) of the increases in PKA activity and spontaneous AP firing rate in response to β-AR stimulation or PDE inhibition were tightly linked. We characterized the system rate-limiting biochemical reactions by integrating these experimentally derived data into a mechanistic-computational model. Model simulations predicted that phospholamban phosphorylation is a potent target of the increase in PKA activity that links to increase in spontaneous AP firing rate. In summary, the kinetics and stoichiometry of increases in PKA activity in response to a physiological (β-AR stimulation) or pharmacological (PDE inhibitor) stimuli match those of changes in the AP firing rate. Thus Ca(2+)-cAMP/PKA-dependent phosphorylation limits the rate and magnitude of increase in spontaneous AP firing rate.

β-adrenergic effects on cardiac myofilaments and contraction in an integrated rabbit ventricular myocyte model

A five-state model of myofilament contraction was integrated into a well-established rabbit ventricular myocyte model of ion channels, Ca(2+) transporters and kinase signaling to analyze the relative contribution of different phosphorylation targets to the overall mechanical response driven by β-adrenergic stimulation (β-AS). β-AS effect on sarcoplasmic reticulum Ca(2+) handling, Ca(2+), K(+) and Cl(-) currents, and Na(+)/K(+)-ATPase properties was included based on experimental data. The inotropic effect on the myofilaments was represented as reduced myofilament Ca(2+) sensitivity (XBCa) and titin stiffness, and increased cross-bridge (XB) cycling rate (XBcy). Assuming independent roles of XBCa and XBcy, the model reproduced experimental β-AS responses on action potentials and Ca(2+) transient amplitude and kinetics. It also replicated the behavior of force-Ca(2+), release-restretch, length-step, stiffness-frequency and force-velocity relationships, and increased force and shortening in isometric and isotonic twitch contractions. The β-AS effect was then switched off from individual targets to analyze their relative impact on contractility. Preventing β-AS effects on L-type Ca(2+) channels or phospholamban limited Ca(2+) transients and contractile responses in parallel, while blocking phospholemman and K(+) channel (IKs) effects enhanced Ca(2+) and inotropy. Removal of β-AS effects from XBCa enhanced contractile force while decreasing peak Ca(2+) (due to greater Ca(2+) buffering), but had less effect on shortening. Conversely, preventing β-AS effects on XBcy preserved Ca(2+) transient effects, but blunted inotropy (both isometric force and especially shortening). Removal of titin effects had little impact on contraction. Finally, exclusion of β-AS from XBCa and XBcy while preserving effects on other targets resulted in preserved peak isometric force response (with slower kinetics) but nearly abolished enhanced shortening. β-AS effects on XBCa and XBcy have greater impact on isometric and isotonic contraction, respectively.

Role of CaMKII in cardiac arrhythmias

Protein phosphorylation is a central mechanism in vertebrates for the regulation of signaling. With regard to the cardiovascular system, phosphorylation of myocyte targets is critical for the regulation of excitation contraction coupling, metabolism, intracellular calcium regulation, mitochondrial activity, transcriptional regulation, and cytoskeletal dynamics. In fact, pathways that tune protein kinase signaling have been a mainstay for cardiovascular therapies for the past 60 years. The calcium/calmodulin-dependent protein kinase II (CaMKII) is a multifunctional serine/threonine kinase with numerous roles in human physiology. Dysfunction in CaMKII-based signaling has been linked with a host of cardiovascular phenotypes including heart failure and arrhythmia, and CaMKII levels are elevated in human and animal disease models of heart disease. While nearly a decade has been invested in targeting CaMKII for the treatment of heart failure and arrhythmia phenotypes, to date, approaches to target the molecule for antiarrhythmic benefit have been unsuccessful for reasons that are still not entirely clear, although (1) lack of compound specificity and (2) the multitude of downstream targets are likely contributing factors. This review will provide an update on current pathways regulated by CaMKII with the goal of illustrating potential upstream regulatory mechanisms and downstream targets that may be modulated for the prevention of cardiac electrical defects. While the review will cover multiple aspects of CaMKII dysfunction in cardiovascular disease, we have given special attention to the potential of CaMKII-associated late Na(+) current as a novel therapeutic target for cardiac arrhythmia.

Arrhythmogenic transient dynamics in cardiac myocytes

Cardiac action potential alternans and early afterdepolarizations (EADs) are linked to cardiac arrhythmias. Periodic action potentials (period 1) in healthy conditions bifurcate to other states such as period 2 or chaos when alternans or EADs occur in pathological conditions. The mechanisms of alternans and EADs have been extensively studied under steady-state conditions, but lethal arrhythmias often occur during the transition between steady states. Why arrhythmias tend to develop during the transition is unclear. We used low-dimensional mathematical models to analyze dynamical mechanisms of transient alternans and EADs. We show that depending on the route from one state to another, action potential alternans and EADs may occur during the transition between two periodic steady states. The route taken depends on the time course of external perturbations or intrinsic signaling, such as β-adrenergic stimulation, which regulate cardiac calcium and potassium currents with differential kinetics.

How does β-adrenergic signalling affect the transitions from ventricular tachycardia to ventricular fibrillation?

AIMS

Ventricular tachycardia (VT) and fibrillation (VF) are the most lethal cardiac arrhythmias. The degeneration of VT into VF is associated with the breakup of a spiral wave of the action potential in cardiac tissue. β-Adrenergic (βAR) signalling potentiates the L-type Ca current (ICaL) faster than the slow delayed rectifier potassium current (IKs), which transiently prolongs the action potential duration (APD) and promotes early after depolarizations. In this study, we aimed at investigating how βAR signalling affects the transition from VT to VF.

METHODS AND RESULTS

We used a physiologically detailed computer model of the rabbit ventricular myocyte in a two-dimensional tissue to determine how spiral waves respond to βAR activation following administration of isoproterenol. A simplified mathematical model was also used to investigate the underlying dynamics. We found that the spatiotemporal behaviour of spiral waves strongly depends on the kinetics of βAR activation. When βAR activation is rapid, a stable spiral wave turns into small fragments and its electrocardiogram reveals the transition from VT to VF. This is due to the transiently steepened APD restitution induced by the faster activation of ICaL vs. IKs upon sudden βAR activation. The spiral wave may also disappear if its transient wavelength is too large to be supported by the tissue size upon sudden strong βAR activation that prolongs APD transiently. When βAR activation is gradual, a stable spiral wave remains such, because of more limited increase in both APD and slope of APD restitution due to more contemporaneous ICaL and IKs activation.

CONCLUSION

Changes in APD restitution during βAR activation revealed a novel transient spiral wave dynamics; this spatiotemporal characteristic strongly depends on the protocol of isoproterenol application.

Ca(2+) current facilitation is CaMKII-dependent and has arrhythmogenic consequences

The cardiac voltage gated Ca²⁺ current (I_Ca) is critical to the electrophysiological properties, excitation-contraction coupling, mitochondrial energetics, and transcriptional regulation in heart. Thus, it is not surprising that cardiac I_Ca is regulated by numerous pathways. This review will focus on changes in I_Ca that occur during the cardiac action potential (AP), with particular attention to Ca²⁺-dependent inactivation (CDI), Ca²⁺-dependent facilitation (CDF) and how calmodulin (CaM) and Ca²⁺-CaM dependent protein kinase (CaMKII) participate in the regulation of Ca²⁺ current during the cardiac AP. CDI depends on CaM pre-bound to the C-terminal of the L-type Ca²⁺ channel, such that Ca²⁺ influx and Ca²⁺ released from the sarcoplasmic reticulum bind to that CaM and cause CDI. In cardiac myocytes CDI normally pre-dominates over voltage-dependent inactivation. The decrease in I_Ca via CDI provides direct negative feedback on the overall Ca²⁺ influx during a single beat, when myocyte Ca²⁺ loading is high. CDF builds up over several beats, depends on CaMKII-dependent Ca²⁺ channel phosphorylation, and results in a staircase of increasing I_Ca peak, with progressively slower inactivation. CDF and CDI co-exist and in combination may fine-tune the I_Ca waveform during the cardiac AP. CDF may partially compensate for the tendency for Ca²⁺ channel availability to decrease at higher heart rates because of accumulating inactivation. CDF may also allow some reactivation of I_Ca during long duration cardiac APs, and contribute to early afterdepolarizations, a form of triggered arrhythmias.

Modeling CaMKII-mediated regulation of L-type Ca(2+) channels and ryanodine receptors in the heart

Excitation-contraction coupling (ECC) in the cardiac myocyte is mediated by a number of highly integrated mechanisms of intracellular Ca²⁺ transport. Voltage- and Ca²⁺-dependent L-type Ca²⁺ channels (LCCs) allow for Ca²⁺ entry into the myocyte, which then binds to nearby ryanodine receptors (RyRs) and triggers Ca²⁺ release from the sarcoplasmic reticulum in a process known as Ca²⁺-induced Ca²⁺ release. The highly coordinated Ca²⁺-mediated interaction between LCCs and RyRs is further regulated by the cardiac isoform of the Ca²⁺/calmodulin-dependent protein kinase (CaMKII). Because CaMKII targets and modulates the function of many ECC proteins, elucidation of its role in ECC and integrative cellular function is challenging and much insight has been gained through the use of detailed computational models. Multiscale models that can both reconstruct the detailed nature of local signaling events within the cardiac dyad and predict their functional consequences at the level of the whole cell have played an important role in advancing our understanding of CaMKII function in ECC. Here, we review experimentally based models of CaMKII function with a focus on LCC and RyR regulation, and the mechanistic insights that have been gained through their application.

Calcium and IP3 dynamics in cardiac myocytes: experimental and computational perspectives and approaches

Calcium plays a crucial role in excitation-contraction coupling (ECC), but it is also a pivotal second messenger activating Ca²⁺-dependent transcription factors in a process termed excitation-transcription coupling (ETC). Evidence accumulated over the past decade indicates a pivotal role of inositol 1,4,5-trisphosphate receptor (IP₃R)-mediated Ca²⁺ release in the regulation of cytosolic and nuclear Ca²⁺ signals. IP₃ is generated by stimulation of plasma membrane receptors that couple to phospholipase C (PLC), liberating IP₃ from phosphatidylinositol 4,5-bisphosphate (PIP₂). An intriguing aspect of IP₃ signaling is the presence of the entire PIP₂-PLC-IP₃ signaling cascade as well as the presence of IP₃Rs at the inner and outer membranes of the nuclear envelope (NE) which functions as a Ca²⁺ store. The observation that the nucleus is surrounded by its own putative Ca²⁺ store raises the possibility that nuclear IP₃-dependent Ca²⁺ release plays a critical role in ETC. This provides a potential mechanism of regulation that acts locally and autonomously from the global cytosolic Ca²⁺ signal underlying ECC. Moreover, there is evidence that: (i) the sarcoplasmic reticulum (SR) and NE are a single contiguous Ca²⁺ store; (ii) the nuclear pore complex is the major gateway for Ca²⁺ and macromolecules to pass between the cytosol and the nucleoplasm; (iii) the inner membrane of the NE hosts key Ca²⁺ handling proteins including the Na⁺/Ca²⁺ exchanger (NCX)/GM1 complex, ryanodine receptors (RyRs), nicotinic acid adenine dinucleotide phosphate receptors (NAADPRs), Na⁺/K⁺ ATPase, and Na⁺/H⁺ exchanger. Thus, it appears that the nucleus represents a Ca²⁺ signaling domain equipped with its own ion channels and transporters that allow for complex local Ca²⁺ signals. Many experimental and modeling approaches have been used for the study of intracellular Ca²⁺ signaling but the key to the understanding of the dual role of Ca²⁺ mediating ECC and ECT lays in quantitative differences of local [Ca²⁺] in the nuclear and cytosolic compartment. In this review, we discuss the state of knowledge regarding the origin and the physiological implications of nuclear Ca²⁺ transients in different cardiac cell types (adult atrial and ventricular myocytes) as well as experimental and mathematical approaches to study Ca²⁺ and IP₃ signaling in the cytosol and nucleus. In particular, we focus on the concept that highly localized Ca²⁺ signals are required to translocate and activate Ca²⁺-dependent transcription factors (e.g., nuclear factor of activated T-cells, NFAT; histone deacetylase, HDAC) through phosphorylation/dephosphorylation processes.

Modeling CaMKII in cardiac physiology: from molecule to tissue

Post-translational modification of membrane proteins (e.g., ion channels, receptors) by protein kinases is an essential mechanism for control of excitable cell function. Importantly, loss of temporal and/or spatial control of ion channel post-translational modification is common in congenital and acquired forms of cardiac disease and arrhythmia. The multifunctional Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) regulates a number of diverse cellular functions in heart, including excitation-contraction coupling, gene transcription, and apoptosis. Dysregulation of CaMKII signaling has been implicated in human and animal models of disease. Understanding of CaMKII function has been advanced by mathematical modeling approaches well-suited to the study of complex biological systems. Early kinetic models of CaMKII function in the brain characterized this holoenzyme as a bistable molecular switch capable of storing information over a long period of time. Models of CaMKII activity have been incorporated into models of the cell and tissue (particularly in the heart) to predict the role of CaMKII in regulating organ function. Disease models that incorporate CaMKII overexpression clearly demonstrate a link between its excessive activity and arrhythmias associated with congenital and acquired heart disease. This review aims at discussing systems biology approaches that have been applied to analyze CaMKII signaling from the single molecule to intact cardiac tissue. In particular, efforts to use computational biology to provide new insight into cardiac disease mechanisms are emphasized.

A novel computational model of mouse myocyte electrophysiology to assess the synergy between Na+ loading and CaMKII

Ca(2+)-calmodulin-dependent protein kinase II (CaMKII) hyperactivity in heart failure causes intracellular Na(+) ([Na(+)]i) loading (at least in part by enhancing the late Na(+) current). This [Na(+)]i gain promotes intracellular Ca(2+) ([Ca(2+)]i) overload by altering the equilibrium of the Na(+)-Ca(2+) exchanger to impair forward-mode (Ca(2+) extrusion), and favour reverse-mode (Ca(2+) influx) exchange. In turn, this Ca(2+) overload would be expected to further activate CaMKII and thereby form a pathological positive feedback loop of ever-increasing CaMKII activity, [Na(+)]i, and [Ca(2+)]i. We developed an ionic model of the mouse ventricular myocyte to interrogate this potentially arrhythmogenic positive feedback in both control conditions and when CaMKIIδC is overexpressed as in genetically engineered mice. In control conditions, simulation of increased [Na(+)]i causes the expected increases in [Ca(2+)]i, CaMKII activity, and target phosphorylation, which degenerate into unstable Ca(2+) handling and electrophysiology at high [Na(+)]i gain. Notably, clamping CaMKII activity to basal levels ameliorates but does not completely offset this outcome, suggesting that the increase in [Ca(2+)]i per se plays an important role. The effect of this CaMKII-Na(+)-Ca(2+)-CaMKII feedback is more striking in CaMKIIδC overexpression, where high [Na(+)]i causes delayed afterdepolarizations, which can be prevented by imposing low [Na(+)]i, or clamping CaMKII phosphorylation of L-type Ca(2+) channels, ryanodine receptors and phospholamban to basal levels. In this setting, Na(+) loading fuels a vicious loop whereby increased CaMKII activation perturbs Ca(2+) and membrane potential homeostasis. High [Na(+)]i is also required to produce instability when CaMKII is further activated by increased Ca(2+) loading due to β-adrenergic activation. Our results support recent experimental findings of a synergistic interaction between perturbed Na(+) fluxes and CaMKII, and suggest that pharmacological inhibition of intracellular Na(+) loading can contribute to normalizing Ca(2+) and membrane potential dynamics in heart failure.

A KCNQ1 mutation contributes to the concealed type 1 long QT phenotype by limiting the Kv7.1 channel conformational changes associated with protein kinase A phosphorylation

BACKGROUND

Type 1 long QT syndrome (LQT1) is caused by loss-of-function mutations in the KCNQ1-encoded Kv7.1 channel that conducts the slowly activating component of the delayed rectifier K(+) current (IKs). Clinically, the diagnosis of LQT1 is complicated by variable phenotypic expressivity, whereby approximately 25% of genotype-positive individuals present with concealed LQT1 (resting corrected QT [QTc] interval ≤460 ms).

OBJECTIVE

To determine whether a specific molecular mechanism contributes to concealed LQT1.

METHODS

We identified a multigenerational LQT1 family whereby 79% of the patients genotype-positive for p.Ile235Asn-KCNQ1 (I235N-Kv7.1) have concealed LQT1. We assessed the effect I235N-Kv7.1 has on IKs and the ventricular action potential (AP) by using in vitro analysis and computational simulations.

RESULTS

Clinical data showed that all 10 patients with I235N-Kv7.1 have normal resting QTc intervals but abnormal QTc interval prolongation during the recovery phase of an electrocardiographic treadmill stress test. Voltage-clamping HEK293 cells coexpressing wild-type Kv7.1 and I235N-Kv7.1 (to mimic the patients' genotypes) showed that I235N-Kv7.1 generated relatively normal functioning Kv7.1 channels but were insensitive to protein kinase A (PKA) activation. Phosphomimetic and quinidine sensitivity studies suggest that I235N-Kv7.1 limits the conformational changes in Kv7.1 channels, which are necessary to upregulate IKs after PKA phosphorylation. Computational ventricular AP simulations predicted that the PKA insensitivity of I235N-Kv7.1 is primarily responsible for prolonging the AP with β-adrenergic stimulation, especially at slower cycle lengths.

CONCLUSIONS

KCNQ1 mutations that generate relatively normal Kv7.1 channels, but limit the upregulation of IKs by PKA activation, likely contribute to concealed LQT1.

Ranolazine for congenital and acquired late INa-linked arrhythmias: in silico pharmacological screening

RATIONALE

The antianginal ranolazine blocks the human ether-a-go-go-related gene-based current IKr at therapeutic concentrations and causes QT interval prolongation. Thus, ranolazine is contraindicated for patients with preexisting long-QT and those with repolarization abnormalities. However, with its preferential targeting of late INa (INaL), patients with disease resulting from increased INaL from inherited defects (eg, long-QT syndrome type 3 or disease-induced electric remodeling (eg, ischemic heart failure) might be exactly the ones to benefit most from the presumed antiarrhythmic properties of ranolazine.

OBJECTIVE

We developed a computational model to predict if therapeutic effects of pharmacological targeting of INaL by ranolazine prevailed over the off-target block of IKr in the setting of inherited long-QT syndrome type 3 and heart failure.

METHODS AND RESULTS

We developed computational models describing the kinetics and the interaction of ranolazine with cardiac Na(+) channels in the setting of normal physiology, long-QT syndrome type 3-linked ΔKPQ mutation, and heart failure. We then simulated clinically relevant concentrations of ranolazine and predicted the combined effects of Na(+) channel and IKr blockade by both the parent compound ranolazine and its active metabolites, which have shown potent blocking effects in the therapeutically relevant range. Our simulations suggest that ranolazine is effective at normalizing arrhythmia triggers in bradycardia-dependent arrhythmias in long-QT syndrome type 3 as well tachyarrhythmogenic triggers arising from heart failure-induced remodeling.

CONCLUSIONS

Our model predictions suggest that acute targeting of INaL with ranolazine may be an effective therapeutic strategy in diverse arrhythmia-provoking situations that arise from a common pathway of increased pathological INaL.

The promise of CaMKII inhibition for heart disease: preventing heart failure and arrhythmias

INTRODUCTION

Calcium-calmodulin-dependent protein kinase II (CaMKII) has emerged as a central mediator of cardiac stress responses which may serve several critical roles in the regulation of cardiac rhythm, cardiac contractility and growth. Sustained and excessive activation of CaMKII during cardiac disease has, however, been linked to arrhythmias, and maladaptive cardiac remodeling, eventually leading to heart failure (HF) and sudden cardiac death.

AREAS COVERED

In the current review, the authors describe the unique structural and biochemical properties of CaMKII and focus on its physiological effects in cardiomyocytes. Furthermore, they provide evidence for a role of CaMKII in cardiac pathologies, including arrhythmogenesis, myocardial ischemia and HF development. The authors conclude by discussing the potential for CaMKII as a target for inhibition in heart disease.

EXPERT OPINION

CaMKII provides a promising nodal point for intervention that may allow simultaneous prevention of HF progression and development of arrhythmias. For future studies and drug development there is a strong rationale for the development of more specific CaMKII inhibitors. In addition, an improved understanding of the differential roles of CaMKII subtypes is required.

Post-translational modifications of the cardiac Na channel: contribution of CaMKII-dependent phosphorylation to acquired arrhythmias

The voltage-gated Na channel isoform 1.5 (NaV1.5) is the pore forming α-subunit of the voltage-gated cardiac Na channel, which is responsible for the initiation and propagation of cardiac action potentials. Mutations in the SCN5A gene encoding NaV1.5 have been linked to changes in the Na current leading to a variety of arrhythmogenic phenotypes, and alterations in the NaV1.5 expression level, Na current density, and/or gating have been observed in acquired cardiac disorders, including heart failure. The precise mechanisms underlying these abnormalities have not been fully elucidated. However, several recent studies have made it clear that NaV1.5 forms a macromolecular complex with a number of proteins that modulate its expression levels, localization, and gating and is the target of extensive post-translational modifications, which may also influence all these properties. We review here the molecular aspects of cardiac Na channel regulation and their functional consequences. In particular, we focus on the molecular and functional aspects of Na channel phosphorylation by the Ca/calmodulin-dependent protein kinase II, which is hyperactive in heart failure and has been causally linked to cardiac arrhythmia. Understanding the mechanisms of altered NaV1.5 expression and function is crucial for gaining insight into arrhythmogenesis and developing novel therapeutic strategies.

Role of CaMKII in post acidosis arrhythmias: a simulation study using a human myocyte model

Postacidotic arrhythmias have been associated to increased sarcoplasmic reticulum (SR) Ca(2+) load and Ca(2+)/calmodulin-dependent protein kinase II (CaMKII) activation. However, the molecular mechanisms underlying these arrhythmias are still unclear. To better understand this process, acidosis produced by CO2 increase from 5% to 30%, resulting in intracellular pH (pHi) change from 7.15 to 6.7, was incorporated into a myocyte model of excitation-contraction coupling and contractility, including acidotic inhibition of L-type Ca(2+) channel (I(CaL)), Na(+)-Ca(2+) exchanger, Ca(2+) release through the SR ryanodine receptor (RyR2) (I(rel)), Ca(2+) reuptake by the SR Ca(2+) ATPase2a (I(up)), Na(+)-K(+) pump, K(+) efflux through the inward rectifier K(+) channel and the transient outward K(+) flow (I(to)) together with increased activity of the Na(+)-H(+) exchanger (I(NHE)). Simulated CaMKII regulation affecting I(rel), I(up), I(CaL), I(NHE) and I(to) was introduced in the model to partially compensate the acidosis outcome. Late Na(+) current increase by CaMKII was also incorporated. Using this scheme and assuming that diastolic Ca(2+) leak through the RyR2 was modulated by the resting state of this channel and the difference between SR and dyadic cleft [Ca(2+)], postacidotic delayed after depolarizations (DADs) were triggered upon returning to normal pHi after 6 min acidosis. The model showed that DADs depend on SR Ca(2+) load and on increased Ca(2+) leak through RyR2. This postacidotic arrhythmogenic pattern relies mainly on CaMKII effect on I(CaL) and I(up), since its individual elimination produced the highest DAD reduction. The model further revealed that during the return to normal pHi, DADs are fully determined by SR Ca(2+) load at the end of acidosis. Thereafter, DADs are maintained by SR Ca(2+) reloading by Ca(2+) influx through the reverse NCX mode during the time period in which [Na(+)]i is elevated.

Compartmentation of cAMP signalling in cardiomyocytes in health and disease

3',5'-cyclic adenosine monophosphate (cAMP) is a ubiquitous second messenger critically involved in the regulation of heart function. It has been shown to act in discrete subcellular signalling compartments formed by differentially localized receptors, phosphodiesterases and protein kinases. Cardiac diseases such as hypertrophy or heart failure are associated with structural and functional remodelling of these microdomains which leads to changes in cAMP compartmentation. In this review, we will discuss recent key findings which provided new insights into cAMP compartmentation in cardiomyocytes with a particular focus on its alterations in heart disease.

Compartmentalization of β-adrenergic signals in cardiomyocytes

Under β-adrenergic stimulation, the distribution of cAMP is highly restricted at distinct intracellular domains for compartmentalized activation of protein kinase A, which promotes selective phosphorylation of proteins for contractile responses in cardiomyocytes. This is primarily due to a concerted effort between restrictions of cAMP distribution by a family of phosphodiesterases and locally anchored protein kinase A by a family of scaffold A kinase-anchoring proteins. Moreover, these regulatory mechanisms underlie the cross talk between β-adrenergic signals and other receptor-stimulated signaling cascades, which alters the compartmentalized β-adrenergic signals for proper contractility in myocardium. Maintaining integrity of compartmentalized β-adrenergic signals is critical for physiological cardiac function and for preventing development of cardiac diseases.

β-adrenergic stimulation activates early afterdepolarizations transiently via kinetic mismatch of PKA targets

Sympathetic stimulation regulates cardiac excitation-contraction coupling in hearts but can also trigger ventricular arrhythmias caused by early afterdepolarizations (EADs) in pathological conditions. Isoproterenol (ISO) stimulation can transiently cause EADs which could result from differential kinetics of L-type Ca current (ICaL) vs. delayed rectifier potassium current (IKs) effects, but multiple PKA targets complicate mechanistic analysis. Utilizing a biophysically detailed model integrating Ca and β-adrenergic signaling, we investigate how different phosphorylation kinetics and targets influence β-adrenergic-induced transient EADs. We found that: 1) The faster time course of ICaL vs. IKs increases recapitulates experimentally observed ISO-induced transient EADs (which are due to ICaL reactivation). These EADs disappear at steady state ISO and do not occur during more gradual ISO application. 2) This ICaL vs. IKs kinetic mismatch with ISO can also induce transient EADs due to spontaneous sarcoplasmic reticulum (SR) Ca release and Na/Ca exchange current. The increased ICaL, SR Ca uptake and action potential duration (APD) raise SR Ca to cause spontaneous SR Ca release, but eventual IKs activation and APD shortening abolish these EADs. 3) Phospholemman (PLM) phosphorylation decreases both types of EADs by increasing outward Na/K-ATPase current (INaK) for ICaL-mediated EADs, and reducing intracellular Na and Ca loading for SR Ca-release-mediated EADs. Slowing PLM phosphorylation kinetics abolishes this protective effect. 4) Blocking phospholamban (PLB) phosphorylation has little effect on ICaL-mediated transient EADs, but abolishes SR Ca-release-mediated transient EADs by limiting SR Ca loading. 5) RyR phosphorylation has little effect on either transient EAD type. Our study emphasizes the importance of understanding non-steady state kinetics of several systems in mediating β-adrenergic-induced EADs and arrhythmias.

Crucial role for Ca2(+)/calmodulin-dependent protein kinase-II in regulating diastolic stress of normal and failing hearts via titin phosphorylation

RATIONALE

Myocardial diastolic stiffness and cardiomyocyte passive force (F(passive)) depend in part on titin isoform composition and phosphorylation. Ca(2+)/calmodulin-dependent protein kinase-II (CaMKII) phosphorylates ion channels, Ca(2+)-handling proteins, and chromatin-modifying enzymes in the heart, but has not been known to target titin.

OBJECTIVE

To elucidate whether CaMKII phosphorylates titin and modulates F(passive) in normal and failing myocardium.

METHODS AND RESULTS

Titin phosphorylation was assessed in CaMKIIδ/γ double-knockout (DKO) mouse, transgenic CaMKIIδC-overexpressing mouse, and human hearts, by Pro-Q-Diamond/Sypro-Ruby staining, autoradiography, and immunoblotting using phosphoserine-specific titin-antibodies. CaMKII-dependent site-specific titin phosphorylation was quantified in vivo by mass spectrometry using stable isotope labeling by amino acids in cell culture mouse heart mixed with wild-type (WT) or DKO heart. F(passive) of single permeabilized cardiomyocytes was recorded before and after CaMKII-administration. All-titin phosphorylation was reduced by >50% in DKO but increased by up to ≈100% in transgenic versus WT hearts. Conserved CaMKII-dependent phosphosites were identified within the PEVK-domain of titin by quantitative mass spectrometry and confirmed in recombinant human PEVK-fragments. CaMKII also phosphorylated the cardiac titin N2B-unique sequence. Phosphorylation at specific PEVK/titin N2B-unique sequence sites was decreased in DKO and amplified in transgenic versus WT hearts. F(passive) was elevated in DKO and reduced in transgenic compared with WT cardiomyocytes. CaMKII-administration lowered F(passive) of WT and DKO cardiomyocytes, an effect blunted by titin antibody pretreatment. Human end-stage failing hearts revealed higher CaMKII expression/activity and phosphorylation at PEVK/titin N2B-unique sequence sites than nonfailing donor hearts.

CONCLUSIONS

CaMKII phosphorylates the titin springs at conserved serines/threonines, thereby lowering F(passive). Deranged CaMKII-dependent titin phosphorylation occurs in heart failure and contributes to altered diastolic stress.

Human atrial fibrillation: insights from computational electrophysiological models

Computational electrophysiology has proven useful to investigate the mechanisms of cardiac arrhythmias at various spatial scales, from isolated myocytes to the whole heart. This article reviews how mathematical modeling has aided our understanding of human atrial myocyte electrophysiology to study the contribution of structural and electrical remodeling to human atrial fibrillation. Potential new avenues of investigation and model development are suggested.

Decoding myocardial Ca²⁺ signals across multiple spatial scales: a role for sensitivity analysis

Numerous studies have employed mathematical modeling to quantitatively understand release of Ca(2+) from the sarcoplasmic reticulum (SR) in the heart. Models have been used to investigate physiologically important phenomena such as triggering of SR Ca(2+) release by Ca(2+) entry across the cell membrane and spontaneous leak of Ca(2+) from the SR in quiescent heart cells. In this review we summarize studies that have modeled myocardial Ca(2+) at different spatial scales: the sub-cellular level, the cellular level, and the multicellular level. We discuss each category of models from the standpoint of parameter sensitivity analysis, a common simulation procedure that can generate quantitative, comprehensive predictions about how changes in conditions influence model output. We propose that this is a useful perspective for conceptualizing models, in part because a sensitivity analysis requires the investigator to define the relevant parameters and model outputs. This procedure therefore helps to illustrate the capabilities and limitations of each model. We further suggest that in future studies, sensitivity analyses will aid in simplifying complex models and in suggesting experiments to differentiate between competing models built with different assumptions. We conclude with a discussion of unresolved questions that are likely to be addressed over the next several years.

Calmodulin-dependent protein kinase II: linking heart failure and arrhythmias

Understanding relationships between heart failure and arrhythmias, important causes of suffering and sudden death, remains an unmet goal for biomedical researchers and physicians. Evidence assembled over the past decade supports a view that activation of the multifunctional Ca(2+) and calmodulin-dependent protein kinase II (CaMKII) favors myocardial dysfunction and cell membrane electrical instability. CaMKII activation follows increases in intracellular Ca(2+) or oxidation, upstream signals with the capacity to transition CaMKII into a Ca(2+) and calmodulin-independent constitutively active enzyme. Constitutively active CaMKII appears poised to participate in disease pathways by catalyzing the phosphorylation of classes of protein targets important for excitation-contraction coupling and cell survival, including ion channels and Ca(2+) homeostatic proteins, and transcription factors that drive hypertrophic and inflammatory gene expression. This rich diversity of downstream targets helps to explain the potential for CaMKII to simultaneously affect mechanical and electrical properties of heart muscle cells. Proof-of-concept studies from a growing number of investigators show that CaMKII inhibition is beneficial for improving myocardial performance and for reducing arrhythmias. We review the molecular physiology of CaMKII and discuss CaMKII actions at key cellular targets and results of animal models of myocardial hypertrophy, dysfunction, and arrhythmias that suggest CaMKII inhibition may benefit myocardial function while reducing arrhythmias.

Computational approaches to understand cardiac electrophysiology and arrhythmias

Cardiac rhythms arise from electrical activity generated by precisely timed opening and closing of ion channels in individual cardiac myocytes. These impulses spread throughout the cardiac muscle to manifest as electrical waves in the whole heart. Regularity of electrical waves is critically important since they signal the heart muscle to contract, driving the primary function of the heart to act as a pump and deliver blood to the brain and vital organs. When electrical activity goes awry during a cardiac arrhythmia, the pump does not function, the brain does not receive oxygenated blood, and death ensues. For more than 50 years, mathematically based models of cardiac electrical activity have been used to improve understanding of basic mechanisms of normal and abnormal cardiac electrical function. Computer-based modeling approaches to understand cardiac activity are uniquely helpful because they allow for distillation of complex emergent behaviors into the key contributing components underlying them. Here we review the latest advances and novel concepts in the field as they relate to understanding the complex interplay between electrical, mechanical, structural, and genetic mechanisms during arrhythmia development at the level of ion channels, cells, and tissues. We also discuss the latest computational approaches to guiding arrhythmia therapy.

Regulation of sarco(endo)plasmic reticulum Ca2+-ATPase and calsequestrin gene expression in the heart

The precise control of Ca2+levels during the contraction–relaxation cycle in cardiac myocytes is extremely important for normal beat-to-beat contractile activity. The sarcoplasmic reticulum (SR) plays a key role controlling calcium concentration in the cytosol. The SR Ca2+-ATPase (SERCA2) transports Ca2+inside the SR lumen during relaxation of the cardiac myocyte. Calsequestrin (Casq2) is the main protein in the SR lumen, functioning as a Ca2+buffer and participating in Ca2+release by interacting with the ryanodine receptor 2 (RyR2) Ca2+-release channel. Alterations in normal Ca2+handling significantly contribute to the contractile dysfunction observed in cardiac hypertrophy and in heart failure. Transcriptional regulation of the SERCA2 gene has been extensively studied and some of the mechanisms regulating its expression have been elucidated. Overexpression of Sp1 factor in cardiac hypertrophy downregulates SERCA2 gene expression and increased levels of thyroid hormone up-regulates its transcription. Other hormones such norepinephrine, angiotensin II, endothelin-1, parathyroid hormone, prostaglandin-F2α, as well the cytokines tumor necrosis factor-α and interleukin-6 also downregulate SERCA2 expression. Calcium acting through the calcineurin–NFAT (nuclear factor of activated T cells) pathway has been suggested to regulate SERCA2 and CASQ2 gene expression. This review focuses on the current knowledge regarding transcriptional regulation of SERCA2 and CASQ2 genes in the normal and pathologic heart.

Theoretical study of L-type Ca(2+) current inactivation kinetics during action potential repolarization and early afterdepolarizations

Sarcoplasmic reticulum (SR) Ca(2+) release mediates excitation–contraction coupling (ECC) in cardiac myocytes. It is triggered upon membrane depolarization by entry of Ca(2+) via L-type Ca(2+) channels (LTCCs), which undergo both voltage- and Ca(2+)-dependent inactivation (VDI and CDI, respectively). We developed improved models of L-type Ca(2+) current and SR Ca(2+) release within the framework of the Shannon-Bers rabbit ventricular action potential (AP) model. The formulation of SR Ca(2+) release was modified to reproduce high ECC gain at negative membrane voltages. An existing LTCC model was extended to reflect more faithfully contributions of CDI and VDI to total inactivation. Ba(2+) current inactivation included an ion-dependent component (albeit small compared with CDI), in addition to pure VDI. Under physiological conditions (during an AP) LTCC inactivates predominantly via CDI, which is controlled mostly by SR Ca(2+) release during the initial AP phase, but by Ca(2+) through LTCCs for the remaining part. Simulations of decreased CDI or K(+) channel block predicted the occurrence of early and delayed after depolarizations. Our model accurately describes ECC and allows dissection of the relative contributions of different Ca(2+) sources to total CDI, and the relative roles of CDI and VDI, during normal and abnormal repolarization.

Phospholemman is a negative feed-forward regulator of Ca2+ in β-adrenergic signaling, accelerating β-adrenergic inotropy

Sympathetic stimulation enhances cardiac contractility by stimulating β-adrenergic signaling and protein kinase A (PKA). Recently, phospholemman (PLM) has emerged as an important PKA substrate capable of regulating cytosolic Ca(2+) transients. However, it remains unclear how PLM contributes to β-adrenergic inotropy. Here we developed a computational model to clarify PLM's role in the β-adrenergic signaling response. Simulating Na(+) and sarcoplasmic reticulum (SR) Ca(2+) clamps, we identify an effect of PLM phosphorylation on SR unloading as the key mechanism by which PLM confers cytosolic Ca(2+) adaptation to long-term β-adrenergic receptor (β-AR) stimulation. Moreover, we show that phospholamban (PLB) opposes and overtakes these actions on SR load, forming a negative feed-forward loop in the β-adrenergic signaling cascade. This network motif dominates the negative feedback conferred by β-AR desensitization and accelerates β-AR-induced inotropy. Model analysis therefore unmasks key actions of PLM phosphorylation during β-adrenergic signaling, indicating that PLM is a critical component of the fight-or-flight response.

Human atrial action potential and Ca2+ model: sinus rhythm and chronic atrial fibrillation

RATIONALE

Understanding atrial fibrillation (AF) requires integrated understanding of ionic currents and Ca2+ transport in remodeled human atrium, but appropriate models are limited.

OBJECTIVE

To study AF, we developed a new human atrial action potential (AP) model, derived from atrial experimental results and our human ventricular myocyte model.

METHODS AND RESULTS

Atria versus ventricles have lower I(K1), resulting in more depolarized resting membrane potential (≈7 mV). We used higher I(to,fast) density in atrium, removed I(to,slow), and included an atrial-specific I(Kur). I(NCX) and I(NaK) densities were reduced in atrial versus ventricular myocytes according to experimental results. SERCA function was altered to reproduce human atrial myocyte Ca2+ transients. To simulate chronic AF, we reduced I(CaL), I(to), I(Kur) and SERCA, and increased I(K1),I(Ks) and I(NCX). We also investigated the link between Kv1.5 channelopathy, [Ca2+]i, and AF. The sinus rhythm model showed a typical human atrial AP morphology. Consistent with experiments, the model showed shorter APs and reduced AP duration shortening at increasing pacing frequencies in AF or when I(CaL) was partially blocked, suggesting a crucial role of Ca2+ and Na+ in this effect. This also explained blunted Ca2+ transient and rate-adaptation of [Ca2+]i and [Na+]i in chronic AF. Moreover, increasing [Na+]i and altered I(NaK) and I(NCX) causes rate-dependent atrial AP shortening. Blocking I(Kur) to mimic Kv1.5 loss-of-function increased [Ca2+]i and caused early afterdepolarizations under adrenergic stress, as observed experimentally.

CONCLUSIONS

Our study provides a novel tool and insights into ionic bases of atrioventricular AP differences, and shows how Na+ and Ca2+ homeostases critically mediate abnormal repolarization in AF.

Calmodulin binding proteins provide domains of local Ca2+ signaling in cardiac myocytes

Calmodulin (CaM) acts as a common Ca(2+) sensor for many signaling pathways, transducing local Ca(2+) signals into specific cellular outcomes. Many of CaM's signaling functions can be explained by its unique biochemical properties, including high and low affinity Ca(2+)-binding sites with slow and fast kinetics, respectively. CaM is expected to have a limited spatial range of action, emphasizing its role in local Ca(2+) signaling. Interactions with target proteins further fine-tune CaM signal transduction. Here, we focus on only three specific cellular targets for CaM signaling in cardiac myocytes: the L-type Ca(2+) channel, the ryanodine receptor, and the IP(3) receptor. We elaborate a working hypothesis that each channel is regulated by two distinct functional populations of CaM: dedicated CaM and promiscuous CaM. Dedicated CaM is typically tethered to each channel and directly regulates channel activity. In addition, a local pool of promiscuous CaM appears poised to sense local Ca(2+) signals and trigger downstream pathways such as Ca(2+)/CaM dependent-protein kinase II and calcineurin. Understanding how promiscuous CaM coordinates multiple distinct signaling pathways remains a challenge, but is aided by the use of mathematical modeling and a new generation of fluorescent biosensors. This article is part of a special issue entitled "Local Signaling in Myocytes."

Local control of β-adrenergic stimulation: Effects on ventricular myocyte electrophysiology and Ca(2+)-transient

Local signaling domains and numerous interacting molecular pathways and substrates contribute to the whole-cell response of myocytes during β-adrenergic stimulation (βARS). We aimed to elucidate the quantitative contribution of substrates and their local signaling environments during βARS to the canine epicardial ventricular myocyte electrophysiology and calcium transient (CaT). We present a computational compartmental model of βARS and its electrophysiological effects. Novel aspects of the model include localized signaling domains, incorporation of β1 and β2 receptor isoforms, a detailed population-based approach to integrate the βAR and Ca(2+)/Calmodulin kinase (CaMKII) signaling pathways and their effects on a wide range of substrates that affect whole-cell electrophysiology and CaT. The model identifies major roles for phosphodiesterases, adenylyl cyclases, PKA and restricted diffusion in the control of local cAMP levels and shows that activation of specific cAMP domains by different receptor isoforms allows for specific control of action potential and CaT properties. In addition, the model predicts increased CaMKII activity during βARS due to rate-dependent accumulation and increased Ca(2+) cycling. CaMKII inhibition, reduced compartmentation, and selective blockade of β1AR is predicted to reduce the occurrence of delayed afterdepolarizations during βARS. Finally, the relative contribution of each PKA substrate to whole-cell electrophysiology is quantified by comparing simulations with and without phosphorylation of each target. In conclusion, this model enhances our understanding of localized βAR signaling and its whole-cell effects in ventricular myocytes by incorporating receptor isoforms, multiple pathways and a detailed representation of multiple-target phosphorylation; it provides a basis for further studies of βARS under pathological conditions.

Computational models reduce complexity and accelerate insight into cardiac signaling networks

Cardiac signaling networks exhibit considerable complexity in size and connectivity. The intrinsic complexity of these networks complicates the interpretation of experimental findings. This motivates new methods for investigating the mechanisms regulating cardiac signaling networks and the consequences these networks have on cardiac physiology and disease. Next-generation experimental techniques are also generating a wealth of genomic and proteomic data that can be difficult to analyze or interpret. Computational models are poised to play a key role in addressing these challenges. Computational models have a long history in contributing to the understanding of cardiac physiology and are useful for identifying biological mechanisms, inferring multiscale consequences to cell signaling activities and reducing the complexity of large data sets. Models also integrate well with experimental studies to explain experimental observations and generate new hypotheses. Here, we review the contributions computational modeling approaches have made to the analysis of cardiac signaling networks and forecast opportunities for computational models to accelerate cardiac signaling research.

Cardiac models in drug discovery and development: a review

Cardiovascular diseases are among the leading causes of death in the developed world. Developing novel therapies for diseases like heart failure is crucial, but this is hampered by the high attrition rate in drug development. The withdrawal of drugs at the final hurdle of approval is mostly because of their unpredictable effects on normal cardiac rhythm. The advent of cardiac computational modeling in the last 5 decades has aided the understanding of heart function significantly. Recently, these models increasingly have been applied toward designing and understanding therapies for cardiac disease. This article will discuss how cellular models of electrophysiology, cell signaling, and metabolism have been used to investigate pharmacologic therapies for cardiac diseases including arrhythmia, ischemia, and heart failure.