Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Competitive tuning: Competition's role in setting the frequency-dependence of Ca2+-dependent proteins

  • Daniel R. Romano ,

    Contributed equally to this work with: Daniel R. Romano, Matthew C. Pharris

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Writing – original draft

    Affiliation Weldon School of Biomedical Engineering, Purdue University, West Lafayette, IN, United States of America

    ⨯
  • Matthew C. Pharris ,

    Contributed equally to this work with: Daniel R. Romano, Matthew C. Pharris

    Roles Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Weldon School of Biomedical Engineering, Purdue University, West Lafayette, IN, United States of America

    ⨯
  • Neal M. Patel,

    Roles Data curation, Visualization, Writing – review & editing

    Affiliation Weldon School of Biomedical Engineering, Purdue University, West Lafayette, IN, United States of America

    ⨯
  • Tamara L. Kinzer-Ursem

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing

    tursem@purdue.edu

    Affiliation Weldon School of Biomedical Engineering, Purdue University, West Lafayette, IN, United States of America

    ⨯

Abstract

A number of neurological disorders arise from perturbations in biochemical signaling and protein complex formation within neurons. Normally, proteins form networks that when activated produce persistent changes in a synapse’s molecular composition. In hippocampal neurons, calcium ion (Ca2+) flux through N-methyl-D-aspartate (NMDA) receptors activates Ca2+/calmodulin signal transduction networks that either increase or decrease the strength of the neuronal synapse, phenomena known as long-term potentiation (LTP) or long-term depression (LTD), respectively. The calcium-sensor calmodulin (CaM) acts as a common activator of the networks responsible for both LTP and LTD. This is possible, in part, because CaM binding proteins are “tuned” to different Ca2+ flux signals by their unique binding and activation dynamics. Computational modeling is used to describe the binding and activation dynamics of Ca2+/CaM signal transduction and can be used to guide focused experimental studies. Although CaM binds over 100 proteins, practical limitations cause many models to include only one or two CaM-activated proteins. In this work, we view Ca2+/CaM as a limiting resource in the signal transduction pathway owing to its low abundance relative to its binding partners. With this view, we investigate the effect of competitive binding on the dynamics of CaM binding partner activation. Using an explicit model of Ca2+, CaM, and seven highly-expressed hippocampal CaM binding proteins, we find that competition for CaM binding serves as a tuning mechanism: the presence of competitors shifts and sharpens the Ca2+ frequency-dependence of CaM binding proteins. Notably, we find that simulated competition may be sufficient to recreate the in vivo frequency dependence of the CaM-dependent phosphatase calcineurin. Additionally, competition alone (without feedback mechanisms or spatial parameters) could replicate counter-intuitive experimental observations of decreased activation of Ca2+/CaM-dependent protein kinase II in knockout models of neurogranin. We conclude that competitive tuning could be an important dynamic process underlying synaptic plasticity.

Author summary

Learning and memory formation are likely associated with dynamic fluctuations in the connective strength of neuronal synapses. These fluctuations, called synaptic plasticity, are regulated by calcium ion (Ca2+) influx through ion channels localized to the post-synaptic membrane. Within the post-synapse, the dominant Ca2+ sensor protein, calmodulin (CaM), may activate a variety of downstream binding partners, each contributing to synaptic plasticity outcomes. The conditions at which certain binding partners most strongly activate are increasingly studied using computational models. Nearly all computational studies describe these binding partners in combinations of only one or two CaM binding proteins. In contrast, we combine seven well-studied CaM binding partners into a single model wherein they simultaneously compete for access to CaM. Our dynamic model suggests that competition narrows the window of conditions for optimal activation of some binding partners, mimicking the Ca2+-frequency dependence of some proteins in vivo. Further characterization of CaM-dependent signaling dynamics in neuronal synapses may benefit our understanding of learning and memory formation. Furthermore, we propose that competitive binding may be another framework, alongside feedback and feed-forward loops, signaling motifs, and spatial localization, that can be applied to other signal transduction networks, particularly second messenger cascades, to explain the dynamical behavior of protein activation.

Citation: Romano DR, Pharris MC, Patel NM, Kinzer-Ursem TL (2017) Competitive tuning: Competition's role in setting the frequency-dependence of Ca2+-dependent proteins. PLoS Comput Biol 13(11): e1005820. https://doi.org/10.1371/journal.pcbi.1005820

Editor: Jeffrey J. Saucerman, University of Virginia, UNITED STATES

Received: May 3, 2017; Accepted: October 13, 2017; Published: November 6, 2017

Copyright: © 2017 Romano et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files. However, relevant Mathematica code is available from the Purdue University Research Repository: Romano, D., Pharris, M. C., Patel, N., Kinzer-Ursem, T. L. (2017). Mathematica Files: Competitive tuning: competition’s role in setting the frequency-dependence of Ca2+-dependent proteins. Purdue University Research Repository. doi:10.4231/R7154F7Q.

Funding: We acknowledge support from the Purdue Summer Undergrad Research Fellowship (SURF) program and the National Institute of Neurological Disorders and Stroke (NINDS) of the National Institute of Health (NIH) under award number R21NS095218. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Calcium (Ca2+) is well-recognized as an important second messenger in cellular signaling. One of the most widely expressed Ca2+ binding proteins, calmodulin (CaM), is a highly conserved protein in the EF-hand family [1] (Fig 1A). CaM has over 100 reported downstream binding proteins, including enzymes that regulate a variety of cellular functions, such as neurotransmitter release in presynaptic neuronal axons[2], insulin secretion in the pancreas [3], and contractility in muscle [4]. Ca2+-dependent signaling in postsynaptic dendrites of excitatory neurons has been the frequent subject of computational studies (see a recent review [5]). Indeed, it comprises an ideal system for mathematical modeling. Its parameters (molecular concentrations and kinetic rate constants) have been measured using controlled experiments, and experimental interest has produced an abundance of published values for model parameterization [4, 6–21]. Two highly-studied functions of synaptic Ca2+ signaling are the induction and maintenance of long-term potentiation (LTP) and long-term depression (LTD) [22], which are correlated to learning processes and memory storage in various brain regions [23–26]. Both LTP and LTD are accompanied by persistent changes in postsynaptic gene transcription [27], actin polymerization [28], and AMPA receptor trafficking [29] that adjust cellular excitability and, in turn, synaptic strength. Among the best-studied forms of LTP and LTD are those initiated by transient, localized increases in intracellular Ca2+ through postsynaptic N-methyl-D-aspartate receptors (NMDARs). CaM translates Ca2+ signals into either LTP or LTD by forming Ca2+/CaM complexes that bind and thereby activate downstream proteins (Fig 1C) [30]. Upon activation, these CaM-dependent proteins, which include a variety of enzymes—kinases, phosphatases, cyclases, and synthases—initiate protein signaling cascades that differentially modulate gene transcription, actin polymerization, and AMPA receptor trafficking.

thumbnail
Download:
Fig 1. Schematic of CaM binding.

(A) Structure of CaM (PDB 1CLL), shown in blue, with two Ca2+ ions (gold) at each terminus. (B) Structure of Ca2+/CaM (PDB 2JZI) bound to a calcineurin (CaN) peptide (red). (C) Schematic of CaM interactions with downstream binding partners. CaM may bind Ng in the absence of Ca2+. In the presence of Ca2+, CaM binds to CaN, CaMKII, NOS, MLCK, and AC1 and AC8.

https://doi.org/10.1371/journal.pcbi.1005820.g001

The frequency [31], amplitude, duration, and location [32] of Ca2+ fluxes determine the pattern of activation of CaM-dependent enzymes and, in turn, the fate of the synapse. For example, 1 Hz stimulation for 10–15 minutes both increases activation of the CaM-dependent phosphatase calcineurin (CaN, or PP3) [33] and induces NMDAR-dependent LTD [34]. On the other hand, 100 Hz stimulation for 1 second increases Ca2+/CaM-dependent protein kinase II (CaMKII) activation and induces NMDAR-dependent LTP [35]. These and similar observations have led to the consensus that kinase cascades induce LTP, while phosphatase cascades induce LTD [36]. But more recent studies have found that CaN may also contribute to LTP induction [37], and that activated CaMKII can promote LTD [38]. These results suggest that normal initiation and maintenance of LTP and LTD do not simply depend on the Boolean activation of kinases or phosphatases in response to a given Ca2+ signal, but rather on the precise activation of a variety of often-counteracting proteins. Therefore, elucidation of the mechanisms that regulate NMDAR-dependent long-term plasticity depends on a complete understanding of the endogenous tuning mechanisms that pair precise patterns of enzyme activation to certain Ca2+ signals.

Computational studies have demonstrated the role of binding dynamics [39], feedback loops [40], and spatial effects [41] in regulating enzyme activation during synaptic Ca2+ signaling. In this work, we hypothesize that competition among CaM binding proteins for access to CaM may serve as an additional tuning mechanism. The concentration of CaM binding partners in the cell far exceeds that of CaM itself [42], and in vitro studies have demonstrated competitive inhibition among neuronal CaM binding partners [43–45]. But, despite the implicit presence of competition in many computational models of Ca2+/CaM signaling in neurons [41, 46–51] and cardiac myocytes [52–56], just one study [46] has had the explicit aim of investigating competition among CaM binding partners as a regulator of enzyme activation. Antunes et al. use such a model to investigate competitive binding as a potential facilitator of the frequency-dependence of CaM binding partners at low frequency Ca2+ fluxes (5 mHz to 5 Hz) for generalized sets of CaM binding partners. However, it is worth noting that both He et. al. and Slavov et. al. both mention competition for CaM as a part of their broader studies on the frequency dependent behavior of networks of generalized CaM targets [51] and relative activation of kinase versus phosphatase signaling [50].

In this work we develop models of Ca2+ binding to CaM that explicitly includes Ca2+-binding to each of the two termini (N- and C-termini, Fig 1). Previous experimental work has shown that CaM is able to activate downstream binding proteins at sub-saturating levels of Ca2+[57]. Moreover, a previous computational study explicitly including Ca2+-binding to each of the two binding sites (N- and C-termini) of CaM has shown that Ca2+ bound at the C-terminus likely significantly contributes to activation of downstream binding partners [39]. Our models also include seven experimentally-characterized postsynaptic CaM binding proteins expressed in CA1 hippocampal neurons. These mathematical models are used to investigate competition’s potential role as a regulator of Ca2+-dependent protein activation across a range of Ca2+ flux frequencies (0.1 Hz to 1000 Hz) that spans those found in vivo and oft employed experimentally in vitro. Specifically, we first develop a set of “isolated” models simulating CaM binding to Ca2+ and just one binding protein. We then combine the isolated models into a “competitive” model that simulates Ca2+ binding to CaM and CaM binding to its binding partners. The CaM binding proteins in this study have been chosen because they are known neuronal proteins with relatively well-characterized CaM-binding kinetics: adenylyl cyclase type I (AC1), the adenylyl cyclase type VIII N-terminus (AC8-Nt), the adenylyl cyclase type VIII C2b domain (AC8-Ct), calcineurin (CaN, also known as PP2B and PP3), CaMKII, myosin light chain kinase (MLCK), neurogranin (Ng), and nitric oxide synthase (NOS) (Fig 1C). Because our model is devoid of feedback loops and spatial localization, the differences in CaM-binding between the competitive and isolated models are solely due to competitive effects. We demonstrate the ability of competition to “tune” the binding and activation profiles of CaM-binding proteins at various Ca2+ flux frequencies and use the model to explain the counterintuitive role of neurogranin in CaMKII activation and LTP induction.

Results

Model development

Model structure.

The interactions of Ca2+, CaM, and CaM binding partners are quite complex. CaM binds a total of four Ca2+ ions, one pair at each of two EF-hand domains located at its N- and C-termini, respectively (Fig 1A and schematically in Fig 2A) [1]. Ca2+-binding at each terminus is highly cooperative [58], but the Ca2+-binding kinetics between these termini are distinct [21]. Moreover, the binding of Ca2+ to CaM changes its affinity for downstream binding partners. Similarly, the binding of CaM to its binding partners changes its affinity for Ca2+ (Fig 2B) [30]. We develop a mathematical model based on mass action kinetics that uses ordinary differential equations to simulate the dynamics of: Ca2+ ions binding reversibly to CaM, the dynamics of CaM binding reversibly to its binding partners, and the dynamics of Ca2+ ions binding reversibly to CaM when CaM is bound to a binding partner [39, 59] (Fig 2).

thumbnail
Download:
Fig 2. Model of Ca2+/CaM interactions.

(A) Reversible binding of Ca2+ binding to CaM (blue). (B) Reversible binding of Ca2+ to CaM bound to a given binding partner, denoted with ‘B’ (green). (C) Reversible binding of a given binding partner to any state of Ca2+/CaM (orange).

https://doi.org/10.1371/journal.pcbi.1005820.g002

A previous study by Pepke et al. offers two models for describing Ca2+-CaM binding. First, they describe a four-state model in which it is assumed that binding of two Ca2+ at each CaM terminus can be treated as a single event due to the highly cooperative binding of Ca2+ at each terminus. Alternatively, a nine-state model is presented that explicitly accounts for each Ca2+ binding event, for which further details are discussed in S1 Appendix. In the present study, we construct both model types and simulate CaM-binding of seven proteins implicated in hippocampal-dependent memory and long-term plasticity [35–37, 60–74]. We find that differences in the output of the four- and nine-state models are negligible for the purposes of this work (Fig S1 and discussion in S1 Appendix). To reduce computational complexity, all model results are based on a four-state model of Ca2+-CaM binding.

Model parameterization.

Initial concentrations of all proteins were either obtained directly from literature or calculated from published values. Equilibrium dissociation (KD) constants not available in the literature are calculated (see [39, 75]) using the thermodynamic principle of microscopic reversibility. From the dissociation constants, any unmeasured kinetic rates are calculated using the equality: KD = koff/kon, and a pair of experimentally-supported assumptions regarding Ca2+, CaM, and binding partner interactions.

First, it is well known that the affinity of Ca2+ for CaM is increased when CaM is bound to a binding partner (CaM-B) [14, 30]. The change in affinity could be represented by either an increase in the association rate constant of Ca2+ for CaM-B or a decrease in the dissociation rate constant of Ca2+ for CaM-B. Experimental work by Peersen et al. showed that the increased affinity of Ca2+ to target peptide-bound CaM was best explained by a reduction in Ca2+ dissociation rate constant [14]. More recent work has shown in the case of Ng binding to CaM, that the change in affinity for Ca2+ results primarily from a change in the dissociation rate constant [76]. Thus, we assume that the increase in affinity of Ca2+ for CaM when CaM is bound to a binding partner comes from a change in the dissociation, but not association, rate constant of Ca2+ from CaM-B.

Second, the binding of Ca2+ to CaM increases the affinity of CaM for most of its binding proteins [30, 77, 78], with the notable exception of Ng. Like others ([39, 41]), we note experimental observations showing that Ca2+ dissociation from CaM typically precedes CaM dissociation from binding proteins [79], indicating that the increase in affinity of CaM for most of its binding proteins in the presence of Ca2+ may be due to an increase the association (and not the dissociation) rate constant of CaM binding to target proteins. Biophysically, Ca2+ binding to CaM induces a conformational change that exposes hydrophobic patches that then facilitate binding to hydrophobic residues on the target proteins [1, 80, 81]. These Ca2+-binding induced structural changes on CaM could be thought of as increasing the probability of successful binding to a target protein, which would translate to an increased association rate. Thus, we implement an assumption that Ca2+-binding changes the association, but not dissociation, rate of CaM to most of its binding partners (with the exception of Ng). It should be noted that it is likely that the increases in affinity discussed above come from changes in both the association and dissociation rate constant parameters. Current experimental techniques are unable to measure the kinetic rate constants of apo-CaM binding to target proteins (again, with the exception of Ng), and so the exact quantitative values or even relative changes in affinity and dissociate rate constants are unknowable at this time. The assumptions implemented here are our best-educated interpretation of current biophysical understanding of Ca2+, CaM and CaM-target binding.

Published rate constant values that were obtained using full-length proteins are used preferentially over those for oligopeptides, but oligopeptide values are included in setting the physiological ranges for sensitivity analyses. Values for simulations are the geometric means of published values, or derived values (listed in S1 Table). Geometric means were chosen as opposed to arithmetic means so that outlier values less significantly biased the parameter values in the simulations.

Adenylyl Cyclase Type I and VIII (AC1 and AC8).

AC1 and AC8 are 123 kDa and 135 kDa [82] membrane-spanning enzymes expressed in the CA1 pyramidal cells of mammalian hippocampus [83, 84]. The primary function of both AC1 and AC8 is the formation of the second messenger cyclic-AMP from ATP [85]. CaM activation of AC1 is dependent upon the binding of Ca2+/CaM [86] to a single site in its C1b domain [87]. Work by D. Cooper and colleagues has shown that CaM binding to the C1b domain on AC1 requires participation from both the N- and C-lobes of CaM [6, 88]. CaM activation of AC8 is also dependent upon Ca2+/CaM-binding [84], but unlike AC1, each AC8 enzyme contains binding sites at both its N-terminus and C2b domain [89]. C2b-binding is the major contributor to CaM-dependent activation of AC8 and can be substantially activated by binding of the Ca2+-bound N-lobe of CaM [6, 88]. A peptide derived from AC8-C2b was able to pull down a CaM with mutations that prohibited Ca2+-binding at the C-lobe at similar levels to WT CaM [6]; indicating that the N-lobe of CaM mediates most of the binding interaction between CaM and the C-terminus of AC8 (AC8-Ct). In contrast, similar pulldown experiments indicate that CaM binding to the N-terminus of AC8 (AC8-Nt) is mediated by the C-lobe of CaM [6]. CaM binding to AC8-Nt does not activate AC8’s enzymatic activity [6], but has been suggested but CaM-binding at the N-terminus may support activation by increasing the local CaM concentration in a “CaM trapping” mechanism [89]. Each of these binding sites associates to Ca2+/CaM in a 1:1 stoichiometry [90]. Therefore, we model AC8 as a pair of distinct targets, AC8-Nt and AC8-Ct. The concentrations of AC1 and AC8 in CA1 pyramidal cells have been estimated at 42.2 and 41.9 μM, respectively [41]. For all simulations, a concentration of 42 μM is used for AC1, AC8-Nt, and AC8-Ct. All kinetic parameters are either obtained from literature [6, 41] or calculated using previously-described assumptions.

Calcineurin (CaN).

CaN is a 78 kDa [7], PSD-associated [91] enzyme expressed in the CA1 pyramidal cells of mammalian hippocampus [92]. As a heterodimer [93], CaN activation is dependent upon both the association of the catalytic subunit CaNA to the regulatory subunit CaNB as well as the binding of CaNA to Ca2+/CaM [10] in a 1:1 stoichiometry [94]. Although CaNB is a Ca2+-binding protein, Ca2+ binding to CaNB does not affect the affinity of CaNA for either CaNB [93] or Ca2+/CaM [10]. For this reason, the binding of both CaNA and Ca2+ to CaNB are neglected in our model. CaN dephosphorylates the residues of many cellular proteins, including AMPA receptors, NMDARs, protein kinase A, and inhibitor-1 [95]. The concentration of CaN in the hippocampus is 36.4 mg of protein per kg of tissue [96]. Assuming an average protein concentration of 100 mg/mL, or 10% by mass [39], the density of CaN in hippocampus was calculated at 36.4 μg/mL, corresponding to a concentration of 0.47 μM. Here, a concentration of 0.5 μM is used. All kinetic parameters are either obtained from literature [7, 8, 10] or calculated based on previously-described assumptions.

Ca2+/CaM-dependent protein kinase II (CaMKII).

CaMKII is a PSD-associated [97] enzyme expressed in CA1 pyramidal cells of the mammalian hippocampus [98]. As a 650 kDa dodecamer, CaMKII is composed of twelve catalytic subunits [99]. In the hippocampus, the alpha isoform of CaMKII comprises approximately two-thirds of these subunits, while the beta isoform constitutes the remaining one-third [99]. The activation of each of these subunits is dependent upon the binding of Ca2+/CaM [100] in a 1:1 stoichiometry [101], such that the full dodecamer binds Ca2+/CaM in a 1:12 ratio [102]. CaMKII phosphorylates the residues of many cellular proteins, including synapsin I, pyruvate kinase, phenylalanine hydroxylase, tyrosine hydroxylase, phospholamban, MLCK, and MAP-2 [103]. CaMKII monomers can also phosphorylate intramolecular neighbors in an autophosphorylation process [104]. The resulting autophosphorylated CaMKII, termed autonomous CaMKII, remains partially active even after dissociating from Ca2+/CaM [105]. Because our model is non-spatial and generally ignores catalytic processes, CaMKII is modeled in its monomeric form (i.e., as separate, independent subunits). The local concentration of catalytic CaMKII subunits in the dendritic spines of CA1 pyramidal cells has been previously estimated at 74 μM [39]. All kinetic parameters are obtained from literature [39].

Myosin light chain kinase (MLCK).

MLCK is a 146 kDa [106] enzyme expressed in CA1 hippocampal dendrites [69]. Its activation is dependent upon the binding of Ca2+/CaM [107] in a 1:1 stoichiometry [108]. MLCK phosphorylates the regulatory light chain of the molecular motor myosin II [109]. The concentration of MLCK in the hippocampus has not been measured, but it has been observed to be much less than that in smooth muscle [110], where its concentration is about 50 μM [111]. Therefore, its concentration in CA1 pyramidal cells is estimated as one order-of-magnitude less, or 5 μM. Because the amino acid sequence of neuronal MLCK is almost identical to that of smooth muscle MLCK [110], we used the kinetic parameters of smooth muscle MLCK in our model. All kinetic parameters are either obtained from literature [4, 11, 13, 14] or calculated based on previously-described assumptions.

Neurogranin (Ng).

Ng is a 7.8 kDa [112], membrane-associated [113] protein expressed in high quantities in the dendritic spines of CA1 pyramidal cells in the mammalian hippocampus [114]. Ng binds apo-CaM in a 1:1 stoichiometry [115]. It has no enzymatic function [116] but has been found to localize CaM to the cell membrane [117], theoretically resulting in the spatial coupling of CaM to both Ca2+ channels and CaM-dependent enzymes. The concentration of Ng in hippocampus has been estimated at 65 μM [73]. All kinetic parameters are obtained from literature [59].

Nitric oxide synthase (NOS).

NOS is a 155 kDa [118], PSD-associated [119] enzyme expressed in CA1 pyramidal cells of the mammalian hippocampus [120]. Its activation is dependent upon the binding of Ca2+/CaM [121] in a 1:1 stoichiometry[19]. NOS catalyzes the formation of nitric oxide and citrulline from arginine [121]. The active form of NOS is a homodimeric complex [122]. However, because our model is non-spatial and generally ignores catalytic processes, NOS is modeled in its monomeric form. NOS is found in 100x diluted, homogenized rat striatum at a density of 0.7 μg/mL [123], corresponding to a concentration of 0.45 μM. Because the density of NOS is 1.5 times greater in the CA1 region of hippocampus than in striatum [124], and because NOS is localized to dendritic spines, a concentration of 1 μM is used. All kinetic parameters are either obtained from literature [15, 18], or calculated based on previously described assumptions.

Calcium (Ca2+).

In response to a single presynaptic action potential, the transient opening of postsynaptic NMDARs in hippocampal dendritic spines generates a single spike in free Ca2+ concentration that peaks at 12 μM and, as Ca2+ is rapidly buffered, decays with a time constant of 12 milliseconds [125]. Therefore, free Ca2+ fluxes into the system by the equation, [Ca](t) = 12e-t/0.012. In our model, this function is a fixed boundary condition, meaning that the total Ca2+ concentration in the system is not conserved over the course of the simulation. Free Ca2+ is introduced into simulations at frequencies ranging from 0.1 Hz to 1 kHz, which spans one order-of-magnitude past the range of frequencies used in LTD- and LTP-inducing experimental protocols [33, 34]. Before the introduction of free Ca2+ into the system, all simulations are run to steady state for 600 seconds to equilibrate Ca2+-independent binding events.

Model analysis

We use the total concentration of CaM-bound protein as a primary output parameter. This is contrary to most published computational models, which investigate the concentration of Ca2+-saturated CaM (CaM4) bound to each protein. This approach is preferred for three main reasons. First, although most CaM-dependent enzymes are maximally activated by binding CaM4, sub-saturated forms of CaM have also been found to activate these enzymes, albeit at a lower catalytic rate [57]. Therefore, the concentration of CaM4-bound enzyme does not represent the total concentration of active enzyme. Second, not all binding sites in our model increase in catalytic activity upon CaM binding. For these proteins (Ng and AC8-Nt) the CaM4-bound concentration is no more relevant than the concentration bound to apo-CaM or, for that matter, any other sub-saturated form. Third, CaM-binding to non-catalytic sites has been found to influence CaM availability to CaM-dependent enzymes [89, 117], suggesting an important physiological role for minimally-active, yet still CaM-bound, enzymes. Therefore, although the total concentration of CaM bound to each binding site is not a direct measure of its activation, it provides important information about patterns of enzyme activation that cannot be inferred from the concentration bound to CaM4 alone. To obtain a representative measure of total CaM-binding during Ca2+ spiking at a particular frequency, the average value (henceforth designated the average bound concentration, Cb) is calculated by Eq 1: (1) Where the subscript b indexes the binding partners, so the average bound concentration for a given binding partner (Cb) is found by integrating the total concentration of that binding partner (Tb) bound to each CaM state (CaMNiCj, i and j = 0, 1, or 2) over the stimulation period (to until tf) and dividing by the stimulus duration (tf—to). To measure relative levels of CaM-binding across various proteins and experimental conditions, for each binding partner we normalize Cb by its peak value from among all the Ca2+ frequencies simulated.

We observe that for competitive models, the frequency range at which Cb peaks may shift or narrow relative to the isolated case. To quantify this tuning, we define a metric of frequency specificity (Sb), where the subscript b indexes the binding partners. A binding partner with high frequency specificity is one that most significantly binds CaM over a narrow range of frequencies; correspondingly, this binding partner’s frequency-dependence curve would have a tall, narrow peak. First, the frequency-dependence curve is integrated and then normalized by the maximum Cb (Eq 2). We also divide by the total simulated frequency range and subtract from 1 to report Sb as a metric that identifies the most strongly tuned binding partners. In Eq 2, f denotes Ca2+ frequency.

(2)

Sensitivity analysis.

To determine which parameters most greatly impacted our models’ outputs and, therefore, which may benefit most from further characterization in future experiments, we conducted two sets of global sensitivity analyses using Latin Hypercube sampling (LHS) to efficiently sample the input parameter space and partial rank correlation coefficients (PRCC) to quantify the results [75]. In one set, we fixed the kinetic rate constants and investigated the impact of variations in initial concentrations on the average bound concentrations (Cb) of the eight CaM binding partners. In the second, we fixed the initial concentrations and investigated the impact of variations in kinetic rate constants. Each of these analyses was performed at low (1 Hz), moderate (10 Hz), and high (100 Hz) frequency Ca2+ oscillations, allowing us to observe how the impacts of parameter variations change with frequency (see S2 Appendix). To control for total Ca2+ introduced, oscillations were limited to 10 concentration spikes, regardless of frequency. In Tables 1 and 2, we present the results of our 10 Hz sensitivity analysis, listing the parameters that most strongly influence each Cb.

thumbnail
Download:
Table 1. Significant PRCCs for initial protein concentration parameters.

https://doi.org/10.1371/journal.pcbi.1005820.t001

thumbnail
Download:
Table 2. Significant PRCCs for rate parameters.

https://doi.org/10.1371/journal.pcbi.1005820.t002

Although their experimental ranges span several orders of magnitude, the kinetic binding constants of CaM4 binding to CaN (konPPCaM4 and koffPPCaM4), as well as CaM4 binding to NOS (konNOSCaM4 and koffNOSCaM4), only significantly affected the Cb of CaN and NOS, respectively. In contrast, the rate constant of CaM4 binding to CaMKII (konKCaM4), despite having an experimental range that varies only four-fold, significantly impacted almost all outputs of CaM binding partners at each of the three frequencies. Future competitive computational models may benefit from more accurate measurement of konKCaM4 than from the more accurate measurement of konPPCaM4, koffPPCaM4, konNOSCaM4, or koffNOSCaM4, despite the clear experimental uncertainty in measurements of the latter four.

Competition regulates CaM-binding dynamics

To investigate how competition alters the CaM-binding dynamics of each of the eight binding partners, we plotted the normalized concentrations of individual partners bound to different CaM states: apo-CaM (CaM0), CaM bound to two Ca2+ ions at its N-terminus (CaM2N), CaM bound to two Ca2+ ions at its C-terminus (CaM2C), and CaM4 (Fig 3). In each simulation, 10 Ca2+ fluxes (not plotted) were introduced at 10 Hz, corresponding to the logarithmic midpoint of our chosen frequency range. In Fig 3, the different colors of the plotted traces correspond to the concentration of binding partner bound to each of the four CaM states normalized to the total concentration of all CaM-bound binding partner (CaMtot). The time-course of CaM binding partners bound to various states of CaM in micromolar for 1 second of 10 Hz Ca2+ flux is plotted in Fig S3 in S1 Appendix.

thumbnail
Download:
Fig 3. Competition for CaM alters binding dynamics.

Time-course of CaM binding partners bound to various states of CaM for 1 second of 10 Hz Ca2+ flux: CaM0 (blue), CaM2N (red), CaM2C (green), CaM4 (purple), and CaMtot (orange). CaMtot is the sum of all CaM-bound states for a given protein. The concentration of each species is normalized against its maximum value of CaMtot. Solid lines denote the isolated model. Dotted lines denote the competitive model. The differences between isolated and competitive behavior are more significant for some CaM binding partners than others.

https://doi.org/10.1371/journal.pcbi.1005820.g003

As expected, the presence of competitors decreases the concentration of CaM bound to each binding partner. Because the relative contributions of the various CaM states to each binding partner’s CaMtot in the competitive model were similar to those in the isolated model, competition did not appear to have a disproportionately large effect on the binding of any one CaM state. This suggests that CaM, and not Ca2+, is the major limiting factor in the activation of CaM-dependent enzymes in hippocampal dendritic spines. Furthermore, competition appears to change not just the concentration of CaM bound to each partner, but also the CaM-binding dynamics. To paraphrase, concentrations in the competitive model are not simply scaled versions of their counterparts in the isolated model. Instead, competition seems to change how each binding partner responds to rapid Ca2+ transients, including how CaM-binding changes with each subsequent Ca2+ flux. For example, after just three Ca2+ fluxes, the concentration of CaM-bound MLCK no longer changes in the isolated model, while it continues to increase in the competitive model. Conversely, while the CaM-binding of Ng decreases with each subsequent Ca2+ spike in the competitive model, it does not change in the isolated model. Therefore, the dynamic behavior of CaM targets in cellular environments cannot necessarily be inferred from computational studies that model them in isolation.

Finally, although competition attenuates the CaM-binding of all binding partners, the magnitude of their attenuation varies considerably in our model. For example, while NOS experiences virtually no change in CaM-binding in the presence of competitors, CaN experiences a more than 20-fold reduction in CaM-binding in the competitive model. Therefore, the binding partners are unequally competitive under the simulated conditions. From these observations, we hypothesize that the competitiveness of each binding partner (i.e., the ability of a binding partner to bind CaM in the presence of other binding partners) might not be absolute and, instead, that the competitiveness of each protein may change across environmental conditions. In this case, competition for CaM is well-positioned to serve as a tuning mechanism, suppressing the CaM-binding of each binding partner for all but a small range of internal conditions and external stimuli and allowing for the tight control of enzyme activation needed for the precise regulation of LTP, LTD, and other neurological processes. Therefore, we investigate how competition may tune the CaM-binding of each neuronal protein to certain Ca2+ frequencies.

Competition tunes CaM binding to certain Ca2+ frequencies

To investigate our hypothesis that competition affects the frequency-dependence of CaM-binding, we construct frequency-dependence curves for all eight CaM binding sites (distinguishing between each AC8 terminus) using both the isolated and competitive models (Fig S2 in S1 Appendix). The frequency dependence of Cb is then projected onto heat maps (Fig 4A and 4B). For all simulations, Ca2+ oscillations consisted of 100 concentration spikes ranging from 0.1 Hz to 1 kHz.

thumbnail
Download:
Fig 4. Competition tunes activation frequencies.

(A) and (B) show normalized activation of CaM as a function of frequency for the isolated and competitive models, respectively. Red denotes peak activation; blue denotes minimal activation. Frequency windows of peak activation tend to narrow and shift for many of the binding partners in the competitive case. Indeed, (C) indicates a sharpening of activation frequency windows as an increase in specificity in the competitive model, at least for most proteins. Specificity is Sb multiplied by 100 percent.

https://doi.org/10.1371/journal.pcbi.1005820.g004

The introduction of competition shifts the frequency-dependence curves of almost all binding partners. For some, such as AC1, AC8-Ct, CaMKII, and MLCK, this shift is slight, but apparent. For other partners, such as AC8-Nt and CaN, this shift is dramatic. In the competitive model, maximal CaM-binding occurs at frequencies almost one order of magnitude lower for AC8-Nt (10 Hz in the competitive model, as compared to 60 Hz in the isolated model). For CaN, maximal CaM-binding occurs at frequencies over two orders of magnitude lower (0.3 Hz in the competitive model, as compared to 80 Hz in the isolated model). For NOS, a frequency shift is present but not visible in Fig 4B.

Although, as stated earlier, total CaM-binding and enzymatic activation are not the same (particularly for CaN, which is subject to dual regulation by Ca2+/CaM and CaNB), it is worth noting that CaN is activated by low, but not high, frequency stimulation in vivo [33]. Therefore, it would be expected that maximal CaM-binding of CaN occurs at a similarly low frequency. The fact that this held true in the competitive, but not in the isolated, model suggests that the in vivo frequency-dependence of CaN may be reliant upon the presence of cellular competitors. Because of both the established role of CaN in LTD induction [36, 64] and the demonstrated ability of low frequency stimulation to induce LTD [33], our results further suggest that competition for CaM may be essential to normal LTD induction. Furthermore, because activated CaN downregulates LTP induction [63], competitive suppression of CaM-binding to CaN at high frequencies may be equally essential to normal LTP induction.

To investigate the effects of competition on each CaM binding partner’s level of preference for a certain frequency range, we used the frequency-dependence curves to calculate the frequency specificity of each binding partner in both the isolated and competitive models as defined in Eq 2. If a binding partner were only active at one frequency, it would have a frequency specificity of 100 percent.

The introduction of competitors sharpens the frequency-dependence curves of almost all binding partners, as also indicated by increased frequency specificity values in the competitive models relative to the isolated models (Fig 4C); frequency specificities increased for AC1 (42.80%, as compared to 38.73%), AC8-Nt (58.55%, as compared to 24.79%), CaN (39.13%, as compared to 27.19%), CaMKII (54.17%, as compared to 37.89%), MLCK (27.45%, as compared to 5.70%), Ng (58.23%, as compared to 16.64%), and NOS (1.77%, as compared to 0.08%). The sole decrease, AC8-Ct, was small (37.86%, as compared to 39.50%). Therefore, competition for CaM not only regulates CaM-binding by changing the frequencies of maximal CaM binding, but also by narrowing the range over which appreciable CaM binding occurs.

Competition for CaM mediates Ng/CaMKII crosstalk

Two studies have reported decreased CaMKII autophosphorylation and CaMKII activity in CA1 hippocampal slices harvested from Ng genetic knockout (Ng-/-) mice [126, 127]. Although both studies reported about a 30% decrease in CaMKII autophosphorylation and CaMKII activity, they were in disagreement concerning the effect of the genetic knockout (Ng-/-) on LTP induction. Pak et al. (2000) found that wild type (Ng+/+) mice required a single tetanus to achieve potentiation, while Ng-/- mice required multiple tetanic stimulations [126]. In direct contrast, Krucker et al. (2002) found that Ng-/- mice required only a single tetanus to induce LTP [127]. Despite these inconsistent results, both sets of authors suggested that this phenomenon may be caused by abnormal regulation of local Ca2+ and CaM concentrations, a proposal that has since been supported by several studies.

For example, Huang et al. (2004) attributed diminished LTP in Ng-/- mice to lower levels of free Ca2+ following high frequency stimulation [73]. And using two sets of Ng mutants which, respectively cannot bind, and constitutively bind, CaM, Zhong et al. (2009) provided evidence that abnormal regulation of local CaM concentrations may also be responsible. Using a model of the interactions of Ca2+, CaM, CaMKII, CaN, and AMPARs, Zhabotinsky et al. (2006) reproduced the effects of Ng knockout on LTP induction reported by Huang et al., but did not address the diminished CaMKII activity reported by both Pak et al. and Krucker et al. To date, no mathematical model has replicated the paradoxical effect of Ng genetic knockout on autonomous CaMKII activity.

We hypothesize that these phenomena could be explained by competitive tuning. We simulate autonomous CaMKII activation by extending our model according to a previously-published model of CaMKII autophosphorylation by Pepke et al. (see Fig 6 in [39]). In that work, two CaM-bound (active) CaMKII monomers form a complex that enzymatically catalyzes the phosphorylation of one of the monomers. We stimulate this extended model according to an LTP induction protocol followed by Krucker et al., in which hippocampal slices were subjected to two tetanic stimuli of 100 pulses at 100 Hz, 20 seconds apart. Using this protocol, we assess our isolated (Fig 5A) and competitive (Fig 5B) models’ responses to simulated Ng knockout at 600 seconds after the last stimulus. Normalized results from the same experimental stimulation protocol by Krucker et al. are shown in Fig 5C (see activity data in Fig 1F in [127]).

thumbnail
Download:
Fig 5. Competitive tuning explains intermolecular crosstalk.

(A) Simulations of CaMKII phosphorylation in our isolated model with and without inclusion of Ng. (B) Simulations of CaMKII phosphorylation in our competitive model with and without Ng. (C) CaMKII activity in WT and Ng-/- knockout mice from Krucker et al. Simulations were performed to replicate the experimental method of Krucker et al. as closely as possible. (D) The average bound concentration (Cb) of each CaM binding protein in semi-isolated models as a function of Ng concentration. AC8-Ct and AC8-Nt exhibit the greatest relative change in CaM-binding (Cb, Eq 1) as Ng concentration decreases. (E) The average bound concentration (Cb) of each CaM binding protein in the competitive model as a function of Ng concentration. For a decreasing Ng concentration, AC8-Ct and AC8-Nt again exhibit the greatest relative change in CaM-binding. (F) Comparing the semi-isolated (dotted traces) to the competitive (solid traces) model shows that only in the competitive model does summed AC8 (AC8-Nt + AC8-Ct, dark red) mirror the loss in CaM-CaMKII binding as Ng concentration decreases.

https://doi.org/10.1371/journal.pcbi.1005820.g005

In the absence of other competitors, the isolated model elicits similar levels of CaMKII autophosphorylation (pCaMKII) whether in the presence or absence of Ng. That is, the complete removal of Ng, which competes with CaMKII for CaM, results in only a slight increase in pCaMKII (Fig 5A). In contrast, in the presence of competitors for CaM, simulated Ng knockout decreases pCaMKII levels by 44% compared to WT (Fig 5C). Notably, this decrease in pCaMKII is quantitatively similar to the roughly 33% loss of Ca2+-independent CaMKII activity indicated by Krucker et al. [127]. Further, our competitive model results are also consistent with Pak et al., who report a 40% decrease in pCaMKII in KO Ng-/- mice compared to WT (Ng+/+) mice [126].

Because our model does not allow for either spatial effects or variations in free Ca2+ concentration, these results suggest that competition for CaM alone could explain the paradoxical effect of Ng genetic knockout on CaMKII autophosphorylation and activity. pCaMKII levels seem to be regulated, at least in-part, by the competition for CaM established by Ng. With Ng, an abundance of the CaM not bound to Ng preferentially binds CaMKII (at moderate Ca2+ levels) because CaMKII can out-compete the other candidate binding partners. Without Ng, as in Ng-/- knockout mice, this competitive advantage of CaMKII to bind CaM is reduced, likely because the CaM that would normally bind Ng instead binds other proteins that do not dissociate as readily when high levels of Ca2+ are introduced.

This interpretation predicts that the decreased CaMKII autophosphorylation and activity seen in the Ng-/- knockouts occurs as a result of increased CaM-binding to other partners. To identify which other partners most preferentially bind CaM upon decreasing Ng, we first employ “semi-isolated” models containing only Ng and one of the seven other CaM binding partners (Fig 5D). Semi-isolated models are utilized in Fig 5D to help ensure that shifts in binding partner activation with decreasing Ng are in fact due to decreasing Ng. The partners that experience the greatest relative increase in CaM-binding as Ng concentration is decreased are AC8-Ct and AC8-Nt (calculated according the average CaM bound concentration, Cb, in Eq 1). A more pronounced increase in the average bound concentration of AC8-Ct and AC8-Nt is seen in full competitive model simulations at decreasing Ng concentrations (Fig 5E). This could indicate that the decrease in CaMKII autophosphorylation and activity in Ng-/- mice is due to the shift in availability of (that is, the competition for) CaM due to its increased binding to AC8 during high frequency stimulation. To investigate this, the average bound concentrations of AC8-Ct and AC8-Nt are summed together into AC8 in Fig 5E and plotted along with the average bound concentration of CaMKII as a function of initial Ng concentration for both isolated and competitive model simulations. CaM-binding to AC8 appears sufficient to explain these changes, with the amount of increase in the average bound AC8 concentration at decreasing Ng concentration closely mirroring the decrease in the average bound CaMKII concentration.

Discussion

In the present study, we use a system of ordinary differential equations to model the dynamic interactions of Ca2+, CaM, and seven CaM target proteins implicated in LTP and LTD of hippocampal synapses. By developing both “isolated” and “competitive” models of this system, we observe competition among these target proteins for CaM-binding and investigate competition’s role in regulating the frequency-dependent activation of downstream CaM binding proteins. The dynamic behavior of our model is largely determined by kinetic rate constants that describe the binding of CaM to Ca2+ and CaM binding to downstream binding to CaM binding proteins. Our models are parameterized using published values where available, and are calculated by applying experimentally supported assumptions and the thermodynamic principle of microscopic reversibility. Global sensitivity analyses are performed to determine the impact of these assumptions on our conclusions, and we find that very few of the parameters that significantly impacted our results are derived from these assumptions.

One of the major results of this work is that competitive binding could be among the mechanisms by which protein activation is dynamically tuned and regulated. We find that the presence of competitors affects not only the concentration of all respective CaM-bound proteins, but also the CaM-binding dynamics of these targets. Based on the results of the present work, we recommend at least the inclusion of Ng into models simulating the activation of CaM-dependent proteins in response to low frequency Ca2+ transients and the inclusion of CaMKII into models simulating the activation of CaM-dependent proteins in response to high frequency Ca2+ transients. Based on the results of our global sensitivity analyses, these two proteins appear to have the most significant impact on the CaM-binding of other CaM targets at these frequency ranges.

Another major result of this work is that competitive tuning may be able to explain the counter-intuitive results from studies of Ng knockouts in mice (Ng-/-) in which CaMKII autophosphorylation and activity levels were seen to decrease in the Ng-/- compared to WT. Our results suggest that under tetanic stimulation and normal initial Ng concentration, Ng buffers CaM from AC8 but not CaMKII. At low concentrations or in the absence of Ng, AC and particularly AC-Ct, is able to bind more CaM, while CaMKII binds less CaM (Fig 5F). Although the KD value of CaM4 binding to CaMKII and AC-Ct are only within 2-fold of each other (1.7 μM and 0.8 μM, respectively), they exhibit very different binding dynamics based on their binding of sub-saturated CaM (CaM2C and CaM2N). This is best seen in Fig 3. For AC-Ct, the dominant species of CaM that binds is CaM2N, making up greater than 50% of the total CaM species bound to AC-Ct. In contrast, for CaMKII there is no dominant species of CaM that binds; CaM2N and CaM4 are major contributors to the total CaMKII-CaM bound species. The binding dynamics of CaM-CaMKII interactions that are seen in the competitive model suggest, as previous work has suggested [39], that CaMKII binds to CaM2C and this CaM2C is then converted to CaM4 while still bound to CaMKII, as noted by the coincident decline in CaMKII-CaM2 and increase in CaMKII-CaM4 in Fig 3 and S3 Fig in S1 Appendix. The binding dynamics of AC8-Ct seem to indicate that AC8-Ct binds CaM2N and stays bound until the next Ca2+ spike. Thus, we hypothesize that AC8-Nt is able to out compete CaMKII for CaM binding in absence of Ng because of its relatively high affinity for CaM2N. Since the dynamic behavior that we see in the competitive model is so dependent on the rate parameters it would be ideal if more of them could be experimentally determined in the future. To test the hypothesis that AC8 activity would be increased in a Ng-/- model, we suggest an experiment in which cAMP production is measured in CA1 hippocampal slices from Ng+/+ and Ng-/- mice while employing forskolin and specific AC1 blockers to control for cAMP production by AC1 and G protein activation, respectively. If our proposed model is accurate, then increased cAMP production will be observed in Ng-/- mice.

Protein networks for which the initiating ligand is a limiting resource, such as the Ca2+/CaM network studied here, are common in biology. As in vivo ligand concentrations often approach the dissociation constants of their binding partners, the concentration of bound ligand could exceed that of free ligand, resulting in the phenomenon of ligand depletion [128]. Ligand depletion, as described by Edelstein et al., reduces cooperative interactions and broadens the range of signals to which the ligand is most responsive. It may be that we observe ligand depletion phenomena in our isolated models (Fig 4A), given the broad range of Ca2+ frequencies at which many binding partners are activated, especially for AC8 and MLCK. However, if ligand depletion really were the predominant regulatory phenomenon, we would expect that by introducing more binding partners (Fig 4B), the broadening effect of ligand depletion would become more conspicuous. Instead, we see a shift and narrowing of the Ca2+ frequencies over which the binding partners are activated. Thus, we are confident that it is competition among the CaM-binding proteins that is the mechanism underlying this tuning behavior.

Because competition seems to be important in our neuron-based model, we sought to compare our results to a different biological system with Ca2+/CaM-dependent signaling. The 2008 publication by Saucerman and Bers examines activation of CaMKII and CaN in a compartmentalized model of cardiomyocytes, stimulated at Ca2+ frequencies ranging from 0–4 Hz [52, 53]. Although this frequency range is much narrower than that used in our competitive model, we can still compare trends of frequency-dependent protein activation. For example, CaMKII activation increases with frequency for both models. Additionally, our isolated model agrees with the Saucerman-Bers model without CaM buffers, in which CaN activation dramatically increases over 0-4Hz. In our competitive model CaN activation is attenuated, in agreement with the Saucerman-Bers model with CaM buffers. This agreement lends further confidence to our model, as the Saucerman-Bers results were subsequently verified experimentally [129]. It appears our model using explicitly-defined CaM buffers (binding proteins) is consistent with the Saucerman-Bers implementation of generalized, unidentified CaM buffers.

The 2008 model by Saucerman and Bers, though not explicitly spatial, highlights how protein localization may affect model output. In Saucerman’s model, Ca2+ frequency-dependent activation levels are different for the cytosolic and membrane-localized (dyadic) CaN sub-populations. Our current model excludes spatial effects in order to scrutinize competitive binding in the absence of confounding factors. However, we acknowledge that spatial effects likely alter competition for CaM, especially in the PSD. Future work would investigate the effect of spatial localization on competition for CaM binding; in particular instantiating membrane-localized proteins such as AC1, AC8, NOS and especially Ng at or near the membrane. Sub-populations of CaN may also be localized to the PSD through binding with scaffolding proteins such as AKAP79 [130–132]. Indeed, because we describe Ng as freely diffusing, it is possible our model exaggerates the ability of Ng to compete for CaM relative to other proteins in our model. Therefore, it would be interesting to assess whether a competitive model accounting for membrane localization can still explain the paradoxical effect of Ng-/- on CaMKII autophosphorylation.

Together, our results suggest that the frequency-dependence of CaM targets observed in vivo is not an inherent property of these proteins, but rather may be an emergent property of their competitive environment. This competitive tuning may provide a mechanism by which otherwise-independent protein pathways can engage in crosstalk through the limited availability of CaM. We propose that competitive tuning, alongside binding dynamics, feedback loops, and spatial localization, may serve as a major regulator of CaM target protein activation. Furthermore, we have attempted to explain the paradoxical decrease in CaMKII activity seen in Ng-/- mice as a result of the dysregulation of this competitive tuning mechanism. In the absence of spatial effects or aperiodic variations in free Ca2+ concentration, competitive tuning is able to offer an explanation for this phenomenon. It is important to note that other proteins, mechanisms, or pathways not included in this model likely lend robustness and further regulatory mechanisms of this phenomenon. Further, it is unlikely that seven CaM-target proteins studied here are the only CaM target proteins that engage in this type of crosstalk through limiting CaM. If competitive tuning facilitates crosstalk among CaM binding proteins, then genetic disorders, neurological diseases, normal aging processes, and therapeutics that disrupt any one CaM target protein may have non-intuitive effects that extend into other signaling pathways. Computational modeling and analysis will continue to play a large role deciphering these oft counter-intuitive regulatory mechanisms that when disrupted, give rise to complex neurological disorders and other important diseases.

Methods

Simulation methods

All numerical integration and data manipulation were performed in Mathematica as described in Model Analysis. Reaction equations were implemented using Mathematica [133] with the XCellerator package [134]. XCellerator uses the Law of Mass Action to create ordinary differential equations describing the time rate of change in concentration for each binding partner and their respective CaM-bound states. In Eq 3 we monitor the concentration of a generalized Ca2+/CaM state complexed with an arbitrary binding partner, Tb: (3) where i and j = 0, 1, or 2.

For simulations involving autophosphorylation of CaMKII, we extend the system of differential equations generalized in Eq 3 to describe formation of a complex between two active (CaM-bound) CaMKII monomers (Eq 4). Finally, complexes of CaMKII monomers react such that one monomer behaves as an enzyme and the other becomes the phosphorylated substrate (Eq 5). As stated, we refer directly to the previously-published model of CaMKII autophosphorylation by Pepke et al. (see Fig 6 in [39]). (4) (5) Where i, j, m, and n = 0, 1, or 2. Phosphorylated CaMKII monomers may also be one of the two participating species in Eq 4.

All the equations for this model can be found in S3 Appendix. Mathematica files for the complete models can be found on the Purdue PURR database: Romano, D.; Pharris, M. C.; Patel, N.; Kinzer-Ursem, T. L. (2017), "Mathematica Files: Competitive tuning: competition’s role in setting the frequency-dependence of Ca2+-dependent proteins." (DOI: 10.4231/R7154F7Q). Our model code is also being uploaded to the BioModels Database [135–137].

Sensitivity analysis

Despite our best efforts to constrain our models’ parameter values to those that have been experimentally-measured or those which can be calculated by the principle of thermodynamic equilibrium, it was still a valuable exercise to investigate the effects of the previously-described calculations and assumptions on model conclusions. Therefore, a global sensitivity analysis was used to investigate how uncertainty in parameter values impacted model outputs. Latin hypercube sampling (LHS) was used to simultaneously sample input parameter spaces, and partial rank correlation coefficients (PRCC) were calculated to measure the correlation between variation in parameter values and variation in model outputs. These methods have been previously described (see [39, 75]). In short, for each CaM target, a uniform probability distribution of input parameter values was assumed to either span the experimental range specified in S1 Table or, if a range of experimental values is not present, 50–200% of experimental, calculated, or assumed values. A perfect positive correlation gave a PRCC of 1, whereas a perfect negative correlation gives a PRCC of -1. A threshold of 0.5 was used to select for only the parameters that significantly impacted (either positively or negatively) the average bound concentration of each binding partner, and parameters were then ranked by the absolute value of their PRCCs. For the sake of completeness, the sensitivity analysis was done for the nine-state model of Ca2+-CaM binding. Further discussion and enumeration of our sensitivity analysis is in S2 Appendix.

Supporting information

S1 Table. Model parameter description and values.

Refer to the Model Parameterization sub-section for explanation of parameters and justification of our model assumptions.

https://doi.org/10.1371/journal.pcbi.1005820.s001

(PDF)

S1 Appendix. Includes a discussion and comparison of 4-state and 9-state models of CaM-protein binding.

Detailed explanation of the differences between these two model types.

https://doi.org/10.1371/journal.pcbi.1005820.s002

(PDF)

S2 Appendix. Discussion and results of sensitivity analysis.

Partial rank correlation coefficients for all initial concentrations and kinetic parameters and their influence on the average bound concentration of CaM to each of its downstream binding partners included in this study.

https://doi.org/10.1371/journal.pcbi.1005820.s003

(PDF)

S3 Appendix. Equations for four-state model, generated by the Mathematica module XCellerator.

https://doi.org/10.1371/journal.pcbi.1005820.s004

(PDF)

Acknowledgments

The authors thank David Umulis, Katherine Clayton, and Elsje Pienaar for helpful conversations and critical feedback on the development of this manuscript.

References

  1. 1. Babu YS, Bugg CE, Cook WJ. Structure of calmodulin refined at 2.2 A resolution. J Mol Biol. 1988;204(1):191–204. pmid:3145979
  2. 2. Ando K. Calmodulin-dependent regulation of neurotransmitter release differs in subsets of neuronal cells. Brain Res. 2013;1535:1–13. pmid:23973605
  3. 3. Tabuchi H, Yamamoto H, Matsumoto K, Ebihara K, Takeuchi Y, Fukunaga K, et al. Regulation of insulin secretion by overexpression of Ca2+/calmodulin-dependent protein kinase II in insulinoma MIN6 cells. Endocrinology. 2000;141(7):2350–60. Epub JUL 2000. pmid:10875234
  4. 4. Torok K, Trentham DR. Mechanism of 2-chloro-(epsilon-Lys75)-(6-(N,N-diethylamino)phenyl)-(1,3,5-triazin-4-yl)calmodulin interactions with smooth muscle myosin light chain kinase and derived peptides. Biochemistry. 1994;33(43):12807–14. pmid:7947686
  5. 5. He Y, Kulasiri D, Samarasinghe S. Systems biology of synaptic plasticity: A review on N-methyl-D-aspartate receptor mediated biochemical pathways and related mathematical models. BioSyst. 2014;122:7–18.
  6. 6. Masada N, Schaks S, Jackson SE, Sinz A, Cooper DM. Distinct mechanisms of calmodulin binding and regulation of adenylyl cyclases 1 and 8. Biochemistry. 2012;51(40):7917–29. pmid:22971080
  7. 7. Quintana AR, Wang D, Forbes JE, Waxham MN. Kinetics of calmodulin binding to calcineurin. Biochem Biophys Res Commun. 2005;334(2):674–80. pmid:16009337
  8. 8. Hubbard MJ, Klee CB. Calmodulin binding by calcineurin: Ligand-induced renaturation of protein immobilized on nitrocellulose. J Biol Chem. 1987;262(31):15062–70. pmid:2822712
  9. 9. Takano E, Hatanaka M, Maki M. Real-time analysis of the calcium-dependent interaction between calmodulin and a synthetic oligopeptide of calcineurin by a surface plasmon resonance biosensor. FEBS Lett. 1994;352(2):247–50. pmid:7925982
  10. 10. Stemmer PM, Klee CB. Dual calcium ion regulation of calcineurin by calmodulin and calcineurin B. Biochemistry. 1994;33(22):6859–66. pmid:8204620
  11. 11. Kasturi R, Vasulka C, Johnson JD. Ca2+, caldesmon, and myosin light chain kinase exchange with calmodulin. J Biol Chem. 1993;268(11):7958–64. pmid:8463316
  12. 12. Tsvetkov PO, Protasevich II, Gilli R, Lafitte D, Lobachov VM, Haiech J, et al. Apocalmodulin binds to the myosin light chain kinase calmodulin target site. J Biol Chem. 1999;274(26):18161–4. pmid:10373414
  13. 13. Johnson JD, Snyder C, Walsh M, Flynn M. Effects of myosin light chain kinase and peptides on Ca2+ exchange with the N- and C-terminal Ca2+ binding sites of calmodulin. J Biol Chem. 1996;271(2):761–7. pmid:8557684
  14. 14. Peersen OB, Madsen TS, Falke JJ. Intermolecular tuning of calmodulin by target peptides and proteins: Differential effects on Ca2+ binding and implications for kinase activation. Protein Sci. 1997;6(4):794–807. pmid:9098889
  15. 15. McMurry JL, Chrestensen CA, Scott IM, Lee EW, Rahn AM, Johansen AM, et al. Rate, affinity and calcium dependence of nitric oxide synthase isoform binding to the primary physiological regulator calmodulin. FEBS J. 2011;278(24):4943–54. pmid:22004458
  16. 16. Wu G, Berka V, Tsai A-L. Binding kinetics of calmodulin with target peptides of three nitric oxide synthase isozymes. J Inorg Biochem. 2011;105(9):1226–37. pmid:21763233
  17. 17. Zoche M, Bienert M, Beyermann M, Koch K-W. Distinct molecular recognition of calmodulin-binding sites in the neuronal and macrophage nitric oxide synthases: a surface plasmon resonance study. Biochemistry. 1996;35(26):8742–6. pmid:8679637
  18. 18. Sheta EA, McMillan K, Masters BS. Evidence for a bidomain structure of constitutive cerebellar nitric oxide synthase. J Biol Chem. 1994;269(21):15147–53. pmid:7515050
  19. 19. Weissman BA, Jones CL, Liu Q, Gross SS. Activation and inactivation of neuronal nitric oxide synthase: Characterization of Ca2+-dependent [125I]calmodulin binding. Eur J Pharmacol. 2002;435(1):9–18. pmid:11790373
  20. 20. Censarek P, Beyermann M, Koch K-W. Target recognition of apocalmodulin by nitric oxide synthase I peptides. Biochemistry. 2002;41(27):8598–7. pmid:12093276
  21. 21. Persechini A, White HD, Gansz KJ. Different mechanisms for Ca2+ dissociation from complexes of calmodulin with nitric oxide synthase or myosin light chain kinase. J Biol Chem. 1996;271(1):62–7. pmid:8550626
  22. 22. Zucker RS. Calcium- and activity-dependent synaptic plasticity. Curr Opin Neurobiol. 1999;9(3):305–13. pmid:10395573
  23. 23. Whitlock JR, Heynen AJ, Shuler MG, Bear MF. Learning induces long-term potentiation in the hippocampus. Science. 2006;313(5790):1093–7. pmid:16931756
  24. 24. Maren S. Long-term potentiation in the amygdala: A mechanism for emotional learning and memory. Trends Neurosci. 1999;22(12):561–7. pmid:10542437
  25. 25. Ge Y, Dong Z, Bagot RC, Howland JG, Phillips AG, Wong TP, et al. Hippocampal long-term depression is required for the consolidation of spatial memory. Proc Natl Acad Sci U S A. 2010;107(38):16697–702. pmid:20823230
  26. 26. Kemp A, Manahan-Vaughan D. Hippocampal long-term depression: Master or minion in declarative memory processes? Trends Neurosci. 2007;30(3):111–8. pmid:17234277
  27. 27. Platenik J, Kuramoto N, Yoneda Y. Molecular mechanisms associated with long-term consolidation of the NMDA signals. Life Sci. 2000;67(4):335–64. pmid:11003045
  28. 28. Carlisle HJ, Kennedy MB. Spine architecture and synaptic plasticity. Trends Neurosci. 2005;28(4):182–7. pmid:15808352
  29. 29. Anggono V, Huganir RL. Regulation of AMPA receptor trafficking and synaptic plasticity. Curr Opin Neurobiol. 2012;22(3):461–9. pmid:22217700
  30. 30. Xia Z, Storm DR. The role of calmodulin as a signal integrator for synaptic plasticity. Nature Reviews Neuroscience. 2005;6(4):267–76. pmid:15803158
  31. 31. Mayford M, Wang J, Kandel ER, O'Dell TJ. CaMKII regulates the frequency-response function of hippocampal synapses for the production of both LTD and LTP. Cell. 1995;81(6):891–904. pmid:7781066
  32. 32. Evans RC, Blackwell KT. Calcium: Amplitude, duration, or location? Biol Bull. 2015;228(1):75–83. pmid:25745102
  33. 33. Yasuda H, Higashi H, Kudo Y, Inoue T, Hata Y, Mikoshiba K, et al. Imaging of calcineurin activated by long-term depression-inducing synaptic inputs in living neurons of rat visual cortex. Eur J Neurosci. 2003;17(2):287–97. pmid:12542665
  34. 34. Dudek SM, Bear MF. Homosynaptic long-term depression in area CA1 of hippocampus and effects of N-methyl-D-aspartate receptor blockade. Proc Natl Acad Sci U S A. 1992;89(10):4363–7. pmid:1350090
  35. 35. Giese KP, Fedorov NB, Filipkowski RK, Silva AJ. Autophosphorylation at Thr286 of the alpha calcium-calmodulin kinase II in LTP and learning. Science. 1998;279(5352):870–3. pmid:9452388
  36. 36. Lee H-K. Synaptic plasticity and phosphorylation. Pharmacol Ther. 2006;112(3):810–32. pmid:16904750
  37. 37. Lu Y-F, Hayashi Y, Moriwaki A, Tomizawa K, Matsui H. FK506, a Ca2+/calmodulin-dependent phosphatase inhibitor, inhibits the induction of long-term potentiation in the rat hippocampus. Neurosci Lett. 1996;205(2):103–6. pmid:8907327
  38. 38. Pi HJ, Otmakhov N, Lemelin D, De Koninck P, Lisman J. Autonomous CaMKII can promote either long-term potentiation or long-term depression, depending on the state of T305/T306 phosphorylation. J Neurosci. 2010;30(26):8704–9. pmid:20592192
  39. 39. Pepke S, Kinzer-Ursem T, Mihalas S, Kennedy MB. A dynamic model of interactions of Ca2+, calmodulin, and catalytic subunits of Ca2+/calmodulin-dependent protein kinase II. PLoS Comput Biol. 2010;6(2):e1000675. pmid:20168991
  40. 40. Bhalla US, Iyengar R. Emergent properties of networks of biological signaling pathways. Science. 1999;283(5400):381–7. pmid:9888852
  41. 41. Kim M, Park AJ, Havekes R, Chay A, Guercio LA, Oliveira RF, et al. Colocalization of protein kinase A with adenylyl cyclase enhances protein kinase A activity during induction of long-lasting long-term-potentiation. PLoS Comput Biol. 2011;7(6):e1002084. pmid:21738458
  42. 42. Sanabria H, Digman MA, Gratton E, Waxham NM. Spatial diffusivity and availability of intracellular calmodulin. Biophys J. 2008;95(12):6002–15. pmid:18820232
  43. 43. Martzen MR, Slemmon JR. The dendritic peptide neurogranin can regulate a calmodulin-dependent target. J Neurochem. 1995;64(1):92–100. pmid:7528268
  44. 44. Tran Q, Black DJ, Persechini A. Intracellular coupling via limiting calmodulin. J Biol Chem. 2003;278(27):24247–50. pmid:12738782
  45. 45. Blumenthal DK, Charbonneau H, Edelman AM, Hinds TR, Rosenberg GB, Storm DR, et al. Synthetic peptides based on the calmodulin-binding domain of myosin light chain kinase inhibit activation of other calmodulin-dependent enzymes. Biochem Biophys Res Commun. 1988;156(2):860–5. pmid:2903735
  46. 46. Antunes G, Roque AC, Simoes de Souza FM. Modelling intracellular competition for calcium: kinetic and thermodynamic control of different molecular modes of signal decoding. Sci Rep. 2016;6:23730. pmid:27033299
  47. 47. Li L, Stefan MI, Le Novere N. Calcium input frequency, duration and amplitude differentially modulate the relative activation of calcineurin and CaMKII. PLoS One. 2012;7(9):e43810. pmid:22962589
  48. 48. Hayer A, Bhalla US, Friston KJ. Molecular switches at the synapse emerge from receptor and kinase traffic. PLoS Comput Biol. 2005;1(2):e20.
  49. 49. Zhabotinsky AM, Camp RN, Epstein IR, Lisman JE. Role of the neurogranin concentrated in spines in the induction of long-term potentiation. J Neurosci. 2006;26(28):7337–47. pmid:16837580
  50. 50. He Y, Kulasiri D, Samarasinghe S. Modelling bidirectional modulations in synaptic plasticity: A biochemical pathway model to understand the emergence of long term potentiation (LTP) and long term depression (LTD). J Theor Biol. 2016;403:159–77. pmid:27185535
  51. 51. Slavov N, Carey J, Linse S. Calmodulin transduces Ca2+ oscillations into differential regulation of its target proteins. ACS Chem Neurosci. 2013;4(4):601–12. pmid:23384199
  52. 52. Saucerman JJ, Bers DM. Calmodulin Mediates Differential Sensitivity of CaMKII and Calcineurin to Local Ca2+ in Cardiac Myocytes. Biophys J. 2008;95(10):4597–612. pmid:18689454
  53. 53. Saucerman JJ, Bers DM. Calmodulin binding proteins provide domains of local Ca2+ signaling in cardiac myocytes. J Mol Cell Cardiol. 2012;52(2):312–6. pmid:21708171
  54. 54. Song Q, Saucerman JJ, Bossuyt J, Bers DM. Differential Integration of Ca2+-Calmodulin Signal in Intact Ventricular Mocytes at Low and High Affinity Ca2+-Calmodulin Targets. Biol Chem. 2008;283(46).
  55. 55. Onal B, Unudurthi S, Hund TJ. Modeling CaMKII in cardiac physiology: from molecule to tissue. Front Pharmacol. 2014;5(1).
  56. 56. Lakatta E, Maltsev V, Vinogradova T. A coupled system of intracellular Ca2+ clocks and surface membrane voltage controls the timekeeping mechanism of the heart's pacemaker. Circ Res. 2010;106(4).
  57. 57. Shifman JM, Choi MH, Mihalas S, Mayo SL, Kennedy MB. Ca2+/calmodulin-dependent protein kinase II (CaMKII) is activated by calmodulin with two bound calciums. Proc Natl Acad Sci U S A. 2006;103(38):13968–73. Epub 2006/09/13. pmid:16966599; PubMed Central PMCID: PMCPMC1599897.
  58. 58. Linse S, Helmersson A, Forsen S. Calcium binding to calmodulin and its globular domains. J Biol Chem. 1991;266(13):8050–4. pmid:1902469
  59. 59. Kubota Y, Putkey J, Waxham M. Neurogranin controls the spatiotemporal pattern of postsynaptic Ca2+/CaM signaling. Biophys J. 2007;93(11):3848–59. pmid:17704141
  60. 60. Wu Z-L, Thomas SA, Villacres EC, Xia Z, Simmons ML, Chavkin C, et al. Altered behavior and long-term potentiation in type I adenylyl cyclase mutant mice. Proc Natl Acad Sci U S A. 1995;92(1):220–4. pmid:7816821
  61. 61. Wong ST, Athos J, Figueroa XA, Pineda VV. Calcium-stimulated adenylyl cyclase activity is critical for hippocampus-dependent long-term memory and late phase LTP. Neuron. 1999;23(4):787–98. pmid:10482244
  62. 62. Schaefer ML, Wong ST, Wozniak DF, Muglia LM, Liauw JA, Zhuo M, et al. Altered stress-induced anxiety in adenylyl cyclase type VIII-deficient mice. J Neurosci. 2000;20(13):4809–20. pmid:10864938
  63. 63. Wang JH, Kelly PT. Postsynaptic calcineurin activity downregulates synaptic transmission by weakening intracellular Ca2+ signaling mechanisms in hippocampal CA1 neurons. J Neurosci. 1997;17(12):4600–11. pmid:9169521
  64. 64. Mulkey RM, Endo S, Shenolikar S, Malenka RC. Involvement of a calcineurin/inhibitor-1 phosphatase cascade in hippocampal long-term depression. Nature. 1994;369(6480):486–8. pmid:7515479
  65. 65. Malleret G, Haditsch U, Genoux DJ, M. W., Bliss TV, Vanhoose AM, Weitlauf C, et al. Inducible and reversible enhancement of learning, memory, and long-term potentiation by genetic inhibition of calcineurin. Cell. 2001;104(5):675–86. pmid:11257222
  66. 66. Otmakhov N, Griffith LC, Lisman JE. Postsynaptic inhibitors of calcium/calmodulin-dependent protein kinase type II block induction but not maintenance of pairing-induced long-term potentiation. J Neurosci. 1997;17(14):5357–65. pmid:9204920
  67. 67. Silva AJ, Paylor R, Wehner JM, Tonegawa S. Impaired spatial learning in alpha-calcium-calmodulin kinase II mutant mice. Science. 1992;257(5067):206–11. pmid:1321493
  68. 68. Silva AJ, Stevens CF, Tonegawa S, Wang Y. Deficient hippocampal long-term potentiation in alpha-calcium-calmodulin kinase II mutant mice. Science. 1992;257(5067):201–6. pmid:1378648
  69. 69. Lei S, Czerwinska E, Czerwinski W, Walsh MP, MacDonald JF. Regulation of NMDA receptor activity by F-actin and myosin light chain kinase. J Neurosci. 2001;21(21):8464–72. pmid:11606635
  70. 70. Rex CS, Gavin CF, Rubio MD, Kramar EA, Chen LY, Jia Y, et al. Myosin IIb regulates actin dynamics during synaptic plasticity and memory formation. Neuron. 2010;67(4):603–17. pmid:20797537
  71. 71. Fedorov NB, Pasinelli P, Oestreicher AB, Degraan PN, Reymann KG. Antibodies to postsynaptic PKC substrate neurogranin prevent long-term potentiation in hippocampal CA1 neurons. Eur J Neurosci. 1995;7(4):819–22. pmid:7620629
  72. 72. Zhong L, Cherry T, Bies CE, Florence MA, Gerges NZ. Neurogranin enhances synaptic strength through its interaction with calmodulin. EMBO J. 2009;28(19):3027–39. pmid:19713936
  73. 73. Huang K-P, Huang FL, Jager T, Li J, Reymann KG, Balschun D. Neurogranin/RC3 enhances long-term potentiation and learning by promoting calcium-mediated signaling. J Neurosci. 2004;24(47):10660–9. pmid:15564582
  74. 74. Bohme GA, Bon C, Lemaire M, Reibaud M, Piot O, Stutzmann JM, et al. Altered synaptic plasticity and memory formation in nitric oxide synthase inhibitor-treated rats. Proc Natl Acad Sci U S A. 1993;90(19):9191–4. pmid:7692445
  75. 75. Kinzer-Ursem TL, Linderman JJ. Both ligand- and cell-specific parameters control ligand agonism in a kinetic model of G protein-coupled receptor signaling. PLoS Comput Biol. 2007;3(1):e6. pmid:17222056
  76. 76. Hoffman L, Chandrasekar A, Wang X, Putkey JA, Waxham MN. Neurogranin Alters the Structure and Calcium Binding Properties of Calmodulin. J Biol Chem. 2014;289(21):14644–55. pmid:24713697
  77. 77. Olwin BB, Edelman AM, Krebs EG, Storm DR. Quantitation of energy coupling between Ca2+, calmodulin, skeletal-muscle myosin light chain kinase, and kinase substrates. J Biol Chem. 1984;259(17):949–55.
  78. 78. Olwin B, Storm D. Calcium binding to complexes of calmodulin and calmodulin binding proteins. Biochemistry. 1985;24(27).
  79. 79. Brown SE, Martin SR, Bayley PM. Kinetic control of the dissociation pathway of calmodulin-peptide complexes. J Biol Chem. 1997;272(6):3389–97. pmid:9013581
  80. 80. Chattopadhyaya R, Meador W, Means A, Quiocho F. Calmodulin structure refined at 1.7 A resolution. J Mol Biol. 1992;228(4).
  81. 81. Ikura M, Clore G, Gronenborn A, Zhu G, Klee C, Bax A. Solution structure of a calmodulin-target peptide complex by multidimensional NMR. Science. 1992;256(5057).
  82. 82. Conti AC, Maas JW, Muglia LM, Dave BA, Vogt SK, Tran TT, et al. Distinct regional and subcellular localization of adenylyl cyclases type 1 and 8 in mouse brain. Neuroscience. 2007;146(2):713–29. pmid:17335981
  83. 83. Mons N, Yoshimura M, Cooper DM. Discrete expression of Ca2+/calmodulin-sensitive and Ca2+-insensitive adenylyl cyclases in the rat brain. Synapse. 1993;14(1):51–9. pmid:8511718
  84. 84. Cali JJ, Zwaagstra JC, Mons N, Cooper DM, Krupinski J. Type VIII adenylyl cyclase: A Ca2+/calmodulin-stimulated enzyme expressed in discrete regions of rat brain. J Biol Chem. 1994;269(16):12190–5. pmid:8163524
  85. 85. Swislocki NI, Tierney J. Solubilization, stabilization, and partial purification of brain adenylate cyclase from rat. Biochemistry. 1973;12(10):1862–6.
  86. 86. Brostrom CO, Huang Y-C, Wolff DJ. Identification of a calcium-binding protein as a calcium-dependent regulator of brain adenylate cyclase. Proc Natl Acad Sci U S A. 1975;72(1):64–8. pmid:164029
  87. 87. Levin LR, Reed RR. Identification of functional domains of adenylyl cyclase using in vivo chimeras. J Biol Chem. 1995;270(13):7573–9. pmid:7706305
  88. 88. Masada N, Ciruela A, MacDougall DA, Cooper DMF. Distinct Mechanisms of Regulation by Ca(2+)/Calmodulin of Type 1 and 8 Adenylyl Cyclases Support Their Different Physiological Roles. J Biol Chem. 2009;284(7):4451–63. pmid:19029295
  89. 89. Gu C, Cooper DM. Calmodulin-binding sites on adenylyl cyclase type VIII. J Biol Chem. 1999;274(12):8012–21. pmid:10075700
  90. 90. Herbst S, Masada N, Pfennig S, Ihling C, Cooper D, Sinz A. Structural insights into calmodulin/adenylyl cyclase 8 interaction. Anal Bioanal Chem. 2013;405(429):9333–42.
  91. 91. Lilly SM, Alvarez FJ, Tietz EI. Synaptic and subcellular localization of A-kinase anchoring protein 150 in rat hippocampal CA1 pyramidal cells: Co-localization with excitatory synaptic markers. Neuroscience. 2005;134(1):155–63. pmid:15951119
  92. 92. Morioka M, Nagahiro S, Fukunaga K, Miyamoto E, Ushio Y. Calcineurin in the adult rat hippocampus: Different distribution in CA1 and CA3 subfields. Neuroscience. 1997;78(3):673–84. pmid:9153650
  93. 93. Klee CB, Crouch TH, Krinks MH. Calcineurin: A calcium- and calmodulin-binding protein of the nervous system. Proc Natl Acad Sci U S A. 1979;76(12):6270–3. pmid:293720
  94. 94. Dunlap TB, Guo H-F, Cook EC, Holbrook E, Rumi-Masante J, Lester TE, et al. Stoichiometry of the calcineurin regulatory domain-calmodulin complex. Biochemistry. 2014;53(36):5779–90. pmid:25144868
  95. 95. Huiming L, Rao A, Hogan PG. Interaction of calcineurin with substrates and targeting proteins. Trends Cell Biol. 2011;21(2):91–103. pmid:21115349
  96. 96. Wallace RW, Tallant EA, Cheung WY. High levels of a heat-labile calmodulin-binding protein (CaM-BP80) in bovine neostriatum. Biochemistry. 1980;19(9):1831–7. pmid:6246935
  97. 97. Miller SG, Kennedy MB. Distinct forebrain and cerebellar isozymes of type II calcium/calmodulin-dependent protein kinase associate differently with the postsynaptic density fraction. J Biol Chem. 1985;260(15):9039–46. pmid:4019461
  98. 98. Benson DL, Gall CM, Isackson PJ. Dendritic localization of type II calcium calmodulin-dependent protein kinase mRNA in normal and reinnervated rat hippocampus. Neuroscience. 1992;46(4):851–7. pmid:1311815
  99. 99. Bennett MK, Erondu NE, Kennedy MB. Purification and characterization of a calmoduli-dependent protein kinase that is highly concentrated in brain. J Biol Chem. 1983;258(20):12735–44. pmid:6313675
  100. 100. Kuret J, Schulman H. Mechanism of autophosphorylation of the multifunctional Ca2+/calmodulin-dependent protein kinase. J Biol Chem. 1985;260(10):6427–33. pmid:3997831
  101. 101. Katoh T, Fujisawa H. Calmodulin-dependent protein kinase II: Kinetic studies on the interaction with substrates and calmodulin. Biochim Biophys Acta. 1991;1091(2):205–12. pmid:1847304
  102. 102. Kim SA, Heinze KG, Bacia K, Waxham MN, Schwille P. Two-photon cross-correlation analysis of intracellular reactions with variable stoichiometry. Biophys J. 2005;88(6):4319–36. pmid:15792970
  103. 103. Hanson PI, Schulman H. Neuronal Ca2+/calmodulin-dependent protein kinases. Annu Rev Biochem. 1992;61:559–601. pmid:1323238
  104. 104. Irvine EE, von Hertzen LS, Plattner F, Peter K. alphaCaMKII autophosphorylation: A fast track to memory. Trends Neurosci. 2006;29(8):459. pmid:16806507
  105. 105. Meyer T, Hanson P. Calmodulin trapping by calcium-calmodulin-dependent protein kinase. Science. 1992;256(5060):1199–202. pmid:1317063
  106. 106. Edelman AM, Higgins DM, Bowman CL, Haber SN, Rabin RA, Cho-Lee J. Myosin light chain kinase is expressed in neurons and glia: Immunoblotting and immunocytochemical studies. Mol Brain Res. 1992;14(1):27–34.
  107. 107. Dabrowska R, Sherry JM, Aromatorio DK, Hartshorne DJ. Modulator protein as a component of the myosin light chain kinase from chicken gizzard. Biochemistry. 1978;17(2):253–8. pmid:202300
  108. 108. Lukas TJ, Burgess WH, Prendergast FG, Lau W, Watterson DM. Calmodulin binding domains: Characterization of a phosphorylation and calmodulin binding site from myosin light chain kinase. Biochemistry. 1986;25(6):1458–64. pmid:3754463
  109. 109. Conti MA, Adelstein RS. The relationship between calmodulin binding and phosphorylation of smooth muscle myosin kinase by the catalytic subunit of 3':5' cAMP-dependent protein kinase. J Biol Chem. 1981;256(7):3178–81. pmid:6259152
  110. 110. Potier M-C, Chelot E, Pekarsky Y, Gardiner K, Rossier J. The human myosin light chain kinase (MLCK) from hippocampus: Cloning, sequencing, expression, and localization to 3qcen-q21. Genomics. 1995;29(3):562–70. pmid:8575746
  111. 111. Szymanski P, Szymanska G, Goyal R. Differences in calmodulin and calmodulin-binding proteins in phasic and tonic smooth muscle. Am J Physiol. 2002;51(1):C94–C104.
  112. 112. Baudier J, Deloulme JC, Van Dorsselaer A, Black D, Matthes HW. Purification and characterization of a brain-specific protein kinase C substrate, neurogranin (p17). J Biol Chem. 1991;266(1):229–37. pmid:1824695
  113. 113. Dominguez-Gonzalez I, Vazquez-Cuesta SN, Algaba A, Diez-Guerra FJ. Neurogranin binds to phosphatidic acid and associates to cellular membranes. Biochem J. 2007;404(1):31–43. pmid:17295609
  114. 114. Represa A, Deloulme JC, Sensenbrenner M, Ben-Ari Y, Baudier J. Neurogranin: Immunocytochemical localization of a brain-specific protein kinase C substrate. J Neurosci. 1990;10(12):3782–92. pmid:2269883
  115. 115. Gerendasy DD, Herron SR, Watson JB, Sutcliffe JG. Mutational and biophysical studies suggest RC3/neurogranin regulates calmodulin availability. J Biol Chem. 1994;269(35):22420–6. pmid:8071370
  116. 116. Huang K-P, Huang FL, Chen HC. Characterization of a 7.5-kDa protein kinase C substrate (RC3 protein, neurogranin) from rat brain. Arch Biochem Biophys. 1993;305(2):570–80. pmid:8080473
  117. 117. Zhong L, Gerges NZ. Neurogranin targets calmodulin and lowers the threshold for the induction of long-term potentiation. PLoS One. 2012;7(7):e41275. pmid:22848456
  118. 118. Herbison AE, Simonian SX, Norris PJ, Emson PC. Relationship of neuronal nitric oxide synthase immunoreactivity to GnRH neurons in the ovariectomized and intact female rat. J Neuroendocrinol. 1996;8(1):73–82. pmid:8932739
  119. 119. Burette A, Zabel U, Weinberg RJ, Schmidt HH, Valtschanoff JG. Synaptic localization of nitric oxide synthase and soluble guanylyl cyclase in the hippocampus. J Neurosci. 2002;22(20):8961–70. pmid:12388603
  120. 120. Doyle CA, Slater P. Localization of neuronal and endothelial nitric oxide synthase isoforms in human hippocampus. Neuroscience. 1997;76(2):387–95. pmid:9015324
  121. 121. Bredt DS, Snyder SH. Isolation of nitric oxide synthetase, a calmodulin-requiring enzyme. Proc Natl Acad Sci U S A. 1990;87(2):682–5. pmid:1689048
  122. 122. Panda K, Ghosh S, Stuehr DJ. Calmodulin activates intersubunit electron transfer in the neuronal nitric-oxide synthase dimer. J Biol Chem. 2001;276(26):23349–56. pmid:11325964
  123. 123. Koh WC, Choe ES, Lee DK, Chang S-C, Shim Y-B. Monitoring the activation of neuronal nitric oxide synthase in brain tissue and cells with a potentiometric immunosensor. Biosensors Bioelectron. 2009;25(1):211–7.
  124. 124. Lourenco CF, Ferreira NR, Santos RM, Lukacova N, Barbosa RM, Laranjinha J. The pattern of glutamate-induced nitric oxide dynamics in vivo and its correlation with nNOS expression in rat hippocampus, cerebral cortext and striatum. Brain Res. 2014;1554:1–11. pmid:24495843
  125. 125. Sabatini BL, Oertner TG, Svoboda K. The life cycle of Ca2+ ions in dendritic spines. Neuron. 2002;33(3):439–52. pmid:11832230
  126. 126. Pak JH, Huang FL, Li J, Balschun D, Reymann K, Chiang C, et al. Involvement of neurogranin in the modulation of calcium/calmodulin-dependent protein kinase II, synaptic plasticity, and spatial learning: A study with knockout mice. Proc Natl Acad Sci U S A. 2000;97(21):11232–7. pmid:11016969
  127. 127. Krucker T, Siggins GR, NcNamara RK, Lindsley KA, Dao A, Allison DW, et al. Targeted disruption of RC3 reveals a calmodulin-based mechanism for regulating metaplasticity in the hippocampus. J Neurosci. 2002;22(13):5525–35. pmid:12097504
  128. 128. Edelstein SJ, Stefan MI, Le Novere N. Ligand Depletion in vivo Modulates the Dynamic Range and Cooperativity of Signal Transduction. PLoS One. 2010;5(1):10. pmid:20052284
  129. 129. Bazzazi H, Sang L, Dick IE, Joshi-Mukherjee R, Yang W, Yue DT. Novel fluorescence resonance energy transfer-based reporter reveals differential calcineurin activation in neonatal and adult cardiomyocytes. J Physiol. 2015;593(17):3865–84. PubMed pmid:26096996.
  130. 130. Sanderson JL, Gorski JA, Gibson ES, Lam P, Freund RK, Chick WS, et al. AKAP150-anchored calcineurin regulates synaptic plasticity by limiting synaptic incorporation of Ca2+-permeable AMPA receptors. J Neurosci. 2012;32(43):15036–52. Epub 2012/10/27. pmid:23100425; PubMed Central PMCID: PMCPMC3504485.
  131. 131. Sanderson JL, Dell'Acqua ML. AKAP Signaling Complexes in Regulation of Excitatory Synaptic Plasticity. Neuroscientist. 2011;17(3):321–36. PubMed pmid:21498812.
  132. 132. Dell'Acqua ML, Dodge KL, Tavalin SJ, Scott JD. Mapping the protein phosphatase-2B anchoring site on AKAP79. Binding and inhibition of phosphatase activity are mediated by residues 315–360. J Biol Chem. 2002;277(50):48796–802. Epub 2002/10/02. pmid:12354762; PubMed Central PMCID: PMCPMC3923414.
  133. 133. Research W. Mathematica Edition 10.0 Champaign, IL: 2016.
    • 134. Shapiro B, Levchenko A, Meyerowitz E, Wold B, Mjolsness E. Cellerator: Extending a computer algebra system to include biochemical arrows for signal transduction simulations. Bioinformatics. 2002;19:677–8.
    • 135. Chelliah V, Juty N, Ajmera I, Ali R, Dumousseau M, Glont M, et al. BioModels: ten-year anniversary. Nucleic Acids Res. 2015;43(D1):D542–D8. pmid:25414348
    • 136. Li C, Donizelli M, Rodriguez N, Dharuri H, Endler L, Chelliah V, et al. BioModels Database: An enhanced, curated and annotated resource for published quantitative kinetic models. BMC Syst Biol. 2010;4:14.
    • 137. Le Novere N, Bornstein B, Broicher A, Courtot M, Donizelli M, Dharuri H, et al. BioModels Database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems. Nucleic Acids Res. 2006;34:D689–D91. pmid:16381960