Calmodulin Mediates Differential Sensitivity of CaMKII and Calcineurin to Local Ca2+ in Cardiac Myocytes

AbstractCalmodulin (CaM) mediates Ca-dependent regulation of numerous pathways in the heart, including CaM-dependent kinase (CaMKII) and calcineurin (CaN), yet the local Ca2+ signals responsible for their selective activation are unclear. To assess when and where CaM, CaMKII, and CaN may be activated in the cardiac myocyte, we integrated new mechanistic computational models of CaM, CaMKII, and CaN with the Shannon-Bers model of excitation-contraction coupling in the rabbit ventricular myocyte. These models are validated with independent in vitro data. In the intact myocyte, model simulations predict that CaM is highly activated in the dyadic cleft during each beat, but not appreciably in the cytosol. CaMKII-δC was almost insensitive to cytosolic Ca due to relatively low CaM affinity. Dyadic cleft CaMKII exhibits dynamic frequency-dependent responses to Ca, yet autophosphorylates only when local phosphatases are suppressed. In contrast, dyadic cleft CaN in beating myocytes is predicted to be constitutively active, whereas the extremely high affinity of CaN for CaM allows gradual integration of small cytosolic CaM signals. Reversing CaM affinities for CaMKII and CaN also reverses their characteristic local responses. Deactivation of both CaMKII and CaN seems dominated by Ca dissociation from the complex (versus Ca-CaM dissociation from the target). In summary, the different affinities of CaM for CaMKII and CaN determine their sensitivity to local Ca signals in cardiac myocytes

Graphical Approach to Model Reduction for Nonlinear Biochemical Networks

Model reduction is a central challenge to the development and analysis of multiscale physiology models. Advances in model reduction are needed not only for computational feasibility but also for obtaining conceptual insights from complex systems. Here, we introduce an intuitive graphical approach to model reduction based on phase plane analysis. Timescale separation is identified by the degree of hysteresis observed in phase-loops, which guides a “concentration-clamp” procedure for estimating explicit algebraic relationships between species equilibrating on fast timescales. The primary advantages of this approach over Jacobian-based timescale decomposition are that: 1) it incorporates nonlinear system dynamics, and 2) it can be easily visualized, even directly from experimental data. We tested this graphical model reduction approach using a 25-variable model of cardiac β1-adrenergic signaling, obtaining 6- and 4-variable reduced models that retain good predictive capabilities even in response to new perturbations. These 6 signaling species appear to be optimal “kinetic biomarkers” of the overall β1-adrenergic pathway. The 6-variable reduced model is well suited for integration into multiscale models of heart function, and more generally, this graphical model reduction approach is readily applicable to a variety of other complex biological systems

Transfer Functions for Protein Signal Transduction: Application to a Model of Striatal Neural Plasticity

We present a novel formulation for biochemical reaction networks in the context of signal transduction. The model consists of input-output transfer functions, which are derived from differential equations, using stable equilibria. We select a set of 'source' species, which receive input signals. Signals are transmitted to all other species in the system (the 'target' species) with a specific delay and transmission strength. The delay is computed as the maximal reaction time until a stable equilibrium for the target species is reached, in the context of all other reactions in the system. The transmission strength is the concentration change of the target species. The computed input-output transfer functions can be stored in a matrix, fitted with parameters, and recalled to build discrete dynamical models. By separating reaction time and concentration we can greatly simplify the model, circumventing typical problems of complex dynamical systems. The transfer function transformation can be applied to mass-action kinetic models of signal transduction. The paper shows that this approach yields significant insight, while remaining an executable dynamical model for signal transduction. In particular we can deconstruct the complex system into local transfer functions between individual species. As an example, we examine modularity and signal integration using a published model of striatal neural plasticity. The modules that emerge correspond to a known biological distinction between calcium-dependent and cAMP-dependent pathways. We also found that overall interconnectedness depends on the magnitude of input, with high connectivity at low input and less connectivity at moderate to high input. This general result, which directly follows from the properties of individual transfer functions, contradicts notions of ubiquitous complexity by showing input-dependent signal transmission inactivation.Comment: 13 pages, 5 tables, 15 figure

Perspectives on: Cyclic Nucleotide Microdomains and Signaling Specificity

The author provides an introduction to a special issue on cyclic nucleotides concerning four areas in cAMP compartmentation

Subcellular heterogeneity of ryanodine receptor properties in ventricular myocytes with low T-tubule density

Rationale: In ventricular myocytes of large mammals, not all ryanodine receptor (RyR) clusters are associated with T-tubules (TTs); this fraction increases with cellular remodeling after myocardial infarction (MI). Objective: To characterize RyR functional properties in relation to TT proximity, at baseline and after MI. Methods: Myocytes were isolated from left ventricle of healthy pigs (CTRL) or from the area adjacent to a myocardial infarction (MI). Ca2+ transients were measured under whole-cell voltage clamp during confocal linescan imaging (fluo-3) and segmented according to proximity of TTs (sites of early Ca2+ release, F>F50 within 20 ms) or their absence (delayed areas). Spontaneous Ca2+ release events during diastole, Ca2+ sparks, reflecting RyR activity and properties, were subsequently assigned to either category. Results: In CTRL, spark frequency was higher in proximity of TTs, but spark duration was significantly shorter. Block of Na+/Ca2+ exchanger (NCX) prolonged spark duration selectively near TTs, while block of Ca2+ influx via Ca2+ channels did not affect sparks properties. In MI, total spark mass was increased in line with higher SR Ca2+ content. Extremely long sparks (>47.6 ms) occurred more frequently. The fraction of near-TT sparks was reduced; frequency increased mainly in delayed sites. Increased duration was seen in near-TT sparks only; Ca2+ removal by NCX at the membrane was significantly lower in MI. Conclusion: TT proximity modulates RyR cluster properties resulting in intracellular heterogeneity of diastolic spark activity. Remodeling in the area adjacent to MI differentially affects these RyR subpopulations. Reduction of the number of sparks near TTs and reduced local NCX removal limit cellular Ca2+ loss and raise SR Ca2+ content, but may promote Ca2+ waves

Mechanistic Systems Modeling to Improve Understanding and Prediction of Cardiotoxicity Caused by Targeted Cancer Therapeutics

Tyrosine kinase inhibitors (TKIs) are highly potent cancer therapeutics that have been linked with serious cardiotoxicity, including left ventricular dysfunction, heart failure, and QT prolongation. TKI-induced cardiotoxicity is thought to result from interference with tyrosine kinase activity in cardiomyocytes, where these signaling pathways help to control critical processes such as survival signaling, energy homeostasis, and excitation–contraction coupling. However, mechanistic understanding is limited at present due to the complexities of tyrosine kinase signaling, and the wide range of targets inhibited by TKIs. Here, we review the use of TKIs in cancer and the cardiotoxicities that have been reported, discuss potential mechanisms underlying cardiotoxicity, and describe recent progress in achieving a more systematic understanding of cardiotoxicity via the use of mechanistic models. In particular, we argue that future advances are likely to be enabled by studies that combine large-scale experimental measurements with Quantitative Systems Pharmacology (QSP) models describing biological mechanisms and dynamics. As such approaches have proven extremely valuable for understanding and predicting other drug toxicities, it is likely that QSP modeling can be successfully applied to cardiotoxicity induced by TKIs. We conclude by discussing a potential strategy for integrating genome-wide expression measurements with models, illustrate initial advances in applying this approach to cardiotoxicity, and describe challenges that must be overcome to truly develop a mechanistic and systematic understanding of cardiotoxicity caused by TKIs

Development of FRET Assay into Quantitative and High-throughput Screening Technology Platforms for Protein–Protein Interactions

Förster resonance energy transfer (FRET) technology has been widely used in biological and biomedical research and is a very powerful tool in elucidating protein interactions in many cellular processes. Ubiquitination and SUMOylation are multi-step cascade reactions, involving multiple enzymes and protein–protein interactions. Here we report the development of dissociation constant (Kd) determination for protein–protein interaction and cell-based high-throughput screening (HTS) assay in SUMOylation cascade using FRET technology. These developments are based on steady state and high efficiency of fluorescent energy transfer between CyPet and YPet fused with SUMO1 and Ubc9, respectively. The developments in theoretical and experimental procedures for protein interaction Kd determination and cell-based HTS provide novel tools in affinity measurement and protein interaction inhibitor screening. The Kd determined by FRET between SUMO1 and Ubc9 is compatible with those determined with other traditional approaches, such as isothermal titration calorimetry (ITC) and surface plasmon resonance (SPR). The FRET-based HTS is pioneer in cell-based HTS. Both Kd determination and cell-based HTS, carried out in 384-well plate format, provide powerful tools for large-scale and high-throughput applications

Cross-Talk between Signaling Pathways Can Generate Robust Oscillations in Calcium and cAMP

BACKGROUND:To control and manipulate cellular signaling, we need to understand cellular strategies for information transfer, integration, and decision-making. A key feature of signal transduction is the generation of only a few intracellular messengers by many extracellular stimuli. METHODOLOGY/PRINCIPAL FINDINGS:Here we model molecular cross-talk between two classic second messengers, cyclic AMP (cAMP) and calcium, and show that the dynamical complexity of the response of both messengers increases substantially through their interaction. In our model of a non-excitable cell, both cAMP and calcium concentrations can oscillate. If mutually inhibitory, cross-talk between the two second messengers can increase the range of agonist concentrations for which oscillations occur. If mutually activating, cross-talk decreases the oscillation range, but can generate 'bursting' oscillations of calcium and may enable better filtering of noise. CONCLUSION:We postulate that this increased dynamical complexity allows the cell to encode more information, particularly if both second messengers encode signals. In their native environments, it is unlikely that cells are exposed to one stimulus at a time, and cross-talk may help generate sufficiently complex responses to allow the cell to discriminate between different combinations and concentrations of extracellular agonists