Lobe specific Ca2+-calmodulin nano-domain in neuronal spines: a single molecule level analysis

SYNAPTIC PLASTICITY IN THE INJURED BRAIN DEPENDS ON THE TEMPORAL PATTERN OF STIMULATION

Neurostimulation protocols are increasingly used as therapeutic interventions, including for brain injury. In addition to the direct activation of neurons, these stimulation protocols are also likely to have downstream effects on those neurons' synaptic outputs. It is well known that alterations in the strength of synaptic connections (long-term potentiation, LTP; long-term depression, LTD) are sensitive to the frequency of stimulation used for induction, however little is known about the contribution of the temporal pattern of stimulation to the downstream synaptic plasticity that may be induced by neurostimulation in the injured brain. We explored interactions of the temporal pattern and frequency of neurostimulation in the normal cerebral cortex and after mild traumatic brain injury (mTBI), to inform therapies to strengthen or weaken neural circuits in injured brains, as well as to better understand the role of these factors in normal brain plasticity. Whole-cell (WC) patch-clamp recordings of evoked postsynaptic potentials (PSPs) in individual neurons, as well as field potential (FP) recordings, were made from layer 2/3 of visual cortex in response to stimulation of layer 4, in acute slices from control (naïve), sham operated, and mTBI rats. We compared synaptic plasticity induced by different stimulation protocols, each consisting of a specific frequency (1 Hz, 10 Hz, or 100 Hz), continuity (continuous or discontinuous), and temporal pattern (perfectly regular, slightly irregular, or highly irregular). At the individual neuron level, dramatic differences in plasticity outcome occurred when the highly irregular stimulation protocol was used at 1 Hz or 10 Hz, producing an overall LTD in controls and shams, but a robust overall LTP after mTBI. Consistent with the individual neuron results, the plasticity outcomes for simultaneous FP recordings were similar, indicative of our results generalizing to a larger scale synaptic network than can be sampled by individual WC recordings alone. In addition to the differences in plasticity outcome between control (naïve or sham) and injured brains, the dynamics of the changes in synaptic responses that developed during stimulation were predictive of the final plasticity outcome. Our results demonstrate that the temporal pattern of stimulation plays a role in the polarity and magnitude of synaptic plasticity induced in the cerebral cortex while highlighting differences between normal and injured brain responses. Moreover, these results may be useful for optimization of neurostimulation therapies to treat mTBI and other brain disorders, in addition to providing new insights into downstream plasticity signaling mechanisms in the normal brain.

Stationary Ca2+ nanodomains in the presence of buffers with two binding sites

We examine closed-form approximations for the equilibrium Ca²⁺ and buffer concentrations near a point Ca²⁺ source representing a Ca²⁺ channel, in the presence of a mobile buffer with two Ca²⁺ binding sites activated sequentially and possessing distinct binding affinities and kinetics. This allows us to model the impact on Ca²⁺ nanodomains of realistic endogenous Ca²⁺ buffers characterized by cooperative Ca²⁺ binding, such as calretinin. The approximations we present involve a combination or rational and exponential functions, whose parameters are constrained using the series interpolation method that we recently introduced for the case of simpler Ca²⁺ buffers with a single Ca²⁺ binding site. We conduct extensive parameter sensitivity analysis and show that the obtained closed-form approximations achieve reasonable qualitative accuracy for a wide range of buffer's Ca²⁺ binding properties and other relevant model parameters. In particular, the accuracy of the derived approximants exceeds that of the rapid buffering approximation in large portions of the relevant parameter space.

Neurogranin stimulates Ca2+/calmodulin-dependent kinase II by suppressing calcineurin activity at specific calcium spike frequencies

Calmodulin sits at the center of molecular mechanisms underlying learning and memory. Its complex and sometimes opposite influences, mediated via the binding to various proteins, are yet to be fully understood. Calcium/calmodulin-dependent protein kinase II (CaMKII) and calcineurin (CaN) both bind open calmodulin, favoring Long-Term Potentiation (LTP) or Depression (LTD) respectively. Neurogranin binds to the closed conformation of calmodulin and its impact on synaptic plasticity is less clear. We set up a mechanistic computational model based on allosteric principles to simulate calmodulin state transitions and its interactions with calcium ions and the three binding partners mentioned above. We simulated calcium spikes at various frequencies and show that neurogranin regulates synaptic plasticity along three modalities. At low spike frequencies, neurogranin inhibits the onset of LTD by limiting CaN activation. At intermediate frequencies, neurogranin facilitates LTD, but limits LTP by precluding binding of CaMKII with calmodulin. Finally, at high spike frequencies, neurogranin promotes LTP by enhancing CaMKII autophosphorylation. While neurogranin might act as a calmodulin buffer, it does not significantly preclude the calmodulin opening by calcium. On the contrary, neurogranin synchronizes the opening of calmodulin's two lobes and promotes their activation at specific frequencies. Neurogranin suppresses basal CaN activity, thus increasing the chance of CaMKII trans-autophosphorylation at high-frequency calcium spikes. Taken together, our study reveals dynamic regulatory roles played by neurogranin on synaptic plasticity, which provide mechanistic explanations for opposing experimental findings.

Intracellular calcium stores mediate metaplasticity at hippocampal dendritic spines

Key points

Calcium (Ca²⁺) entry mediated by NMDA receptors is considered central to the induction of activity‐dependent synaptic plasticity in hippocampal area CA1; this description does not, however, take into account the potential contribution of endoplasmic reticulum (ER) Ca²⁺ stores.

The ER has a heterogeneous distribution in CA1 dendritic spines, and may introduce localized functional differences in Ca²⁺ signalling between synapses, as suggested by experiments on metabotropic receptor‐dependent long‐term depression.

A physiologically detailed computational model of Ca²⁺ dynamics at a CA3–CA1 excitatory synapse characterizes the contribution of spine ER via metabotropic signalling during plasticity induction protocols.

ER Ca²⁺ release via IP₃ receptors modulates NMDA receptor‐dependent plasticity in a graded manner, to selectively promote synaptic depression with relatively diminished effect on LTP induction; this may temper further strengthening at the stronger synapses which are preferentially associated with ER‐containing spines.

Acquisition of spine ER may thus represent a local, biophysically plausible ‘metaplastic switch’ at potentiated CA1 synapses, contributing to the plasticity–stability balance in neural circuits.

Abstract

Long‐term plasticity mediated by NMDA receptors supports input‐specific, Hebbian forms of learning at excitatory CA3–CA1 connections in the hippocampus. There exists an additional layer of stabilizing mechanisms that act globally as well as locally over multiple time scales to ensure that plasticity occurs in a constrained manner. Here, we investigated the role of calcium (Ca²⁺) stores associated with the endoplasmic reticulum (ER) in the local regulation of plasticity at individual CA1 synapses. Our study was spurred by (1) the curious observation that ER is sparsely distributed in dendritic spines, but over‐represented in larger spines that are likely to have undergone activity‐dependent strengthening, and (2) evidence suggesting that ER motility at synapses can be rapid, and accompany activity‐regulated spine remodelling. We constructed a physiologically realistic computational model of an ER‐bearing CA1 spine, and examined how IP₃‐sensitive Ca²⁺ stores affect spine Ca²⁺ dynamics during activity patterns mimicking the induction of long‐term potentiation and long‐term depression (LTD). Our results suggest that the presence of ER modulates NMDA receptor‐dependent plasticity in a graded manner that selectively enhances LTD induction. We propose that ER may locally tune Ca²⁺‐based plasticity, providing a braking mechanism to mitigate runaway strengthening at potentiated synapses. Our study provides a biophysically accurate description of postsynaptic Ca²⁺ regulation, and suggests that ER in the spine may promote the re‐use of hippocampal synapses with saturated strengths.

Calcium (Ca²⁺) entry mediated by NMDA receptors is considered central to the induction of activity‐dependent synaptic plasticity in hippocampal area CA1; this description does not, however, take into account the potential contribution of endoplasmic reticulum (ER) Ca²⁺ stores.

The ER has a heterogeneous distribution in CA1 dendritic spines, and may introduce localized functional differences in Ca²⁺ signalling between synapses, as suggested by experiments on metabotropic receptor‐dependent long‐term depression.

A physiologically detailed computational model of Ca²⁺ dynamics at a CA3–CA1 excitatory synapse characterizes the contribution of spine ER via metabotropic signalling during plasticity induction protocols.

ER Ca²⁺ release via IP₃ receptors modulates NMDA receptor‐dependent plasticity in a graded manner, to selectively promote synaptic depression with relatively diminished effect on LTP induction; this may temper further strengthening at the stronger synapses which are preferentially associated with ER‐containing spines.

Acquisition of spine ER may thus represent a local, biophysically plausible ‘metaplastic switch’ at potentiated CA1 synapses, contributing to the plasticity–stability balance in neural circuits.

Extension of Rapid Buffering Approximation to Ca2+ Buffers with Two Binding Sites

Fundamental cell processes such as synaptic neurotransmitter release, endocrine hormone secretion, and myocyte contraction are controlled by highly localized calcium (Ca²⁺) signals resulting from brief openings of trans-membrane Ca²⁺ channels. On short temporal and spatial scales, the corresponding local Ca²⁺ nanodomains formed in the vicinity of a single or several open Ca²⁺ channels can be effectively approximated by quasi-stationary solutions. The rapid buffering approximation (RBA) is one of the most powerful of such approximations, and is based on the assumption of instantaneous equilibration of the bimolecular Ca²⁺ buffering reaction, combined with the conservation condition for the total Ca²⁺ and buffer molecule numbers. Previously, RBA has been generalized to an arbitrary arrangement of Ca²⁺ channels on a flat membrane, in the presence of any number of simple Ca²⁺ buffers with one-to-one Ca²⁺ binding stoichiometry. However, many biological buffers have multiple binding sites. For example, buffers and sensors phylogenetically related to calmodulin consist of two Ca²⁺-binding domains (lobes), with each domain binding two Ca²⁺ ions in a cooperative manner. Here we consider an extension of RBA to such buffers with two interdependent Ca²⁺ binding sites. We show that in the presence of such buffers, RBA solution is given by the solution to a cubic equation, analogous to the quadratic equation describing RBA in the case of a simple, one-to-one Ca²⁺ buffer. We examine in detail the dependence of RBA accuracy on buffering parameters, to reveal conditions under which RBA provides sufficient precision.

Control of Ca2+ Influx and Calmodulin Activation by SK-Channels in Dendritic Spines

The key trigger for Hebbian synaptic plasticity is influx of Ca²⁺ into postsynaptic dendritic spines. The magnitude of [Ca²⁺] increase caused by NMDA-receptor (NMDAR) and voltage-gated Ca²⁺ -channel (VGCC) activation is thought to determine both the amplitude and direction of synaptic plasticity by differential activation of Ca²⁺ -sensitive enzymes such as calmodulin. Ca²⁺ influx is negatively regulated by Ca²⁺ -activated K⁺ channels (SK-channels) which are in turn inhibited by neuromodulators such as acetylcholine. However, the precise mechanisms by which SK-channels control the induction of synaptic plasticity remain unclear. Using a 3-dimensional model of Ca²⁺ and calmodulin dynamics within an idealised, but biophysically-plausible, dendritic spine, we show that SK-channels regulate calmodulin activation specifically during neuron-firing patterns associated with induction of spike timing-dependent plasticity. SK-channel activation and the subsequent reduction in Ca²⁺ influx through NMDARs and L-type VGCCs results in an order of magnitude decrease in calmodulin (CaM) activation, providing a mechanism for the effective gating of synaptic plasticity induction. This provides a common mechanism for the regulation of synaptic plasticity by neuromodulators.

Hebbian or associative plasticity is triggered by postsynaptic Ca²⁺ influx which activates calmodulin and CaMKII. The influx of Ca²⁺ through voltage-dependent NMDA receptors and Ca²⁺ channels is regulated by Ca²⁺ -activated K⁺ channels (SK-channels) providing negative feedback regulation of postsynaptic [Ca²⁺]. Using 3-dimensional modeling of Ca²⁺ and calmodulin dynamics within dendritic spines we show that the non-linear relationship between Ca²⁺ influx and calmodulin activation endows SK-channels with the ability to “gate” calmodulin activation and therefore the induction of Hebbian synaptic plasticity. Since SK-channels are inhibited by several neuromodulator receptors including acetylcholine and noradrenaline, the gating of synaptic plasticity by SK-channels could represent a common mechanism by which neuromodulators control the induction of synaptic plasticity.

Calcium: amplitude, duration, or location?

Calcium plays a role in long-term plasticity by triggering postsynaptic signaling pathways for both the strengthening (LTP) and weakening (LTD) of synapses. Since these are opposing processes, several hypotheses have been developed to explain how calcium can trigger LTP in some situations and LTD in others. These hypotheses fall broadly into three categories, based on the amplitude of calcium concentration, the duration of the calcium elevation, and the location of the calcium influx. Here we review the experimental evidence for and against each of these hypotheses and the recent computational models utilizing each. We argue that with new experimental techniques for the precise visualization of calcium and new computational techniques for the modeling of calcium diffusion, it is time to take a new look at the location hypothesis.

Neurogranin alters the structure and calcium binding properties of calmodulin

Neurogranin (Ng) is a member of the IQ motif class of calmodulin (CaM)-binding proteins, and interactions with CaM are its only known biological function. In this report we demonstrate that the binding affinity of Ng for CaM is weakened by Ca(2+) but to a lesser extent (2-3-fold) than that previously suggested from qualitative observations. We also show that Ng induced a >10-fold decrease in the affinity of Ca(2+) binding to the C-terminal domain of CaM with an associated increase in the Ca(2+) dissociation rate. We also discovered a modest, but potentially important, increase in the cooperativity in Ca(2+) binding to the C-lobe of CaM in the presence of Ng, thus sharpening the threshold for the C-domain to become Ca(2+)-saturated. Domain mapping using synthetic peptides indicated that the IQ motif of Ng is a poor mimetic of the intact protein and that the acidic sequence just N-terminal to the IQ motif plays an important role in reproducing Ng-mediated decreases in the Ca(2+) binding affinity of CaM. Using NMR, full-length Ng was shown to make contacts largely with residues in the C-domain of CaM, although contacts were also detected in residues in the N-terminal domain. Together, our results can be consolidated into a model where Ng contacts residues in the N- and C-lobes of both apo- and Ca(2+)-bound CaM and that although Ca(2+) binding weakens Ng interactions with CaM, the most dramatic biochemical effect is the impact of Ng on Ca(2+) binding to the C-terminal lobe of CaM.

Dynamic modulation of spike timing-dependent calcium influx during corticostriatal upstates

The striatum of the basal ganglia demonstrates distinctive upstate and downstate membrane potential oscillations during slow-wave sleep and under anesthetic. The upstates generate calcium transients in the dendrites, and the amplitude of these calcium transients depends strongly on the timing of the action potential (AP) within the upstate. Calcium is essential for synaptic plasticity in the striatum, and these large calcium transients during the upstates may control which synapses undergo plastic changes. To investigate the mechanisms that underlie the relationship between calcium and AP timing, we have developed a realistic biophysical model of a medium spiny neuron (MSN). We have implemented sophisticated calcium dynamics including calcium diffusion, buffering, and pump extrusion, which accurately replicate published data. Using this model, we found that either the slow inactivation of dendritic sodium channels (NaSI) or the calcium inactivation of voltage-gated calcium channels (CDI) can cause high calcium corresponding to early APs and lower calcium corresponding to later APs. We found that only CDI can account for the experimental observation that sensitivity to AP timing is dependent on NMDA receptors. Additional simulations demonstrated a mechanism by which MSNs can dynamically modulate their sensitivity to AP timing and show that sensitivity to specifically timed pre- and postsynaptic pairings (as in spike timing-dependent plasticity protocols) is altered by the timing of the pairing within the upstate. These findings have implications for synaptic plasticity in vivo during sleep when the upstate-downstate pattern is prominent in the striatum.

On the precision of quasi steady state assumptions in stochastic dynamics

Many biochemical networks have complex multidimensional dynamics and there is a long history of methods that have been used for dimensionality reduction for such reaction networks. Usually a deterministic mass action approach is used; however, in small volumes, there are significant fluctuations from the mean which the mass action approach cannot capture. In such cases stochastic simulation methods should be used. In this paper, we evaluate the applicability of one such dimensionality reduction method, the quasi-steady state approximation (QSSA) [L. Menten and M. Michaelis, "Die kinetik der invertinwirkung," Biochem. Z 49, 333369 (1913)] for dimensionality reduction in case of stochastic dynamics. First, the applicability of QSSA approach is evaluated for a canonical system of enzyme reactions. Application of QSSA to such a reaction system in a deterministic setting leads to Michaelis-Menten reduced kinetics which can be used to derive the equilibrium concentrations of the reaction species. In the case of stochastic simulations, however, the steady state is characterized by fluctuations around the mean equilibrium concentration. Our analysis shows that a QSSA based approach for dimensionality reduction captures well the mean of the distribution as obtained from a full dimensional simulation but fails to accurately capture the distribution around that mean. Moreover, the QSSA approximation is not unique. We have then extended the analysis to a simple bistable biochemical network model proposed to account for the stability of synaptic efficacies; the substrate of learning and memory [J. E. Lisman, "A mechanism of memory storage insensitive to molecular turnover: A bistable autophosphorylating kinase," Proc. Natl. Acad. Sci. U.S.A. 82, 3055-3057 (1985)]. Our analysis shows that a QSSA based dimensionality reduction method results in errors as big as two orders of magnitude in predicting the residence times in the two stable states.

AKAP79/150 interacts with the neuronal calcium-binding protein caldendrin

The A kinase-anchoring protein AKAP79/150 is a postsynaptic scaffold molecule and a key regulator of signaling events. At the postsynapse it coordinates phosphorylation and dephosphorylation of receptors via anchoring kinases and phosphatases near their substrates. Interactions between AKAP79 and two Ca(2+) -binding proteins caldendrin and calmodulin have been investigated here. Calmodulin is a known interaction partner of AKAP79/150 that has been shown to regulate activity of the kinase PKC in a Ca(2+) -dependent manner. Pull-down experiments and surface plasmon resonance biosensor analyses have been used here to demonstrate that AKAP79 can also interact with caldendrin, a neuronal calcium-binding protein implicated in regulation of Ca(2+) -influx and release. We demonstrate that calmodulin and caldendrin compete for a partially overlapping binding site on AKAP79 and that their binding is differentially dependent on calcium. Therefore, this competition is regulated by calcium levels. Moreover, both proteins have different binding characteristics suggesting that the two proteins might play complementary roles. The postsynaptic enrichment, the complex binding mechanism, and the competition with calmodulin, makes caldendrin an interesting novel player in the signaling toolkit of the AKAP interactome.

The C-terminus of human Ca(v)2.3 voltage-gated calcium channel interacts with alternatively spliced calmodulin-2 expressed in two human cell lines

Ca(v)2.3 containing voltage-activated Ca(2+) channels are expressed in excitable cells and trigger neurotransmitter and peptide-hormone release. Their expression remote from the fast release sites leads to the accumulation of presynaptic Ca(2+) which can both, facilitate and inhibit the influx of Ca(2+) ions through Ca(v)2.3. The facilitated Ca(2+) influx was recently related to hippocampal postsynaptic facilitation and long term potentiation. To analyze Ca(2+) mediated modulation of cellular processes more in detail, protein partners of the carboxy terminal tail of Ca(v)2.3 were identified by yeast-2-hybrid screening, leading in two human cell lines to the detection of a novel, extended and rarely occurring splice variant of calmodulin-2 (CaM-2), called CaM-2-extended (CaM-2-ext). CaM-2-ext interacts biochemically with the C-terminus of Ca(v)2.3 similar to the classical CaM-2 as shown by co-immunoprecipitation. Functionally, only CaM-2-ext reduces whole cell inward currents significantly. The insertion of the novel 46 nts long exon and the consecutive expression of CaM-2-ext must be dependent on a new upstream translation initiation site which is only rarely used in the tested human cell lines. The structure of the N-terminal extension is predicted to be more hydrophobic than the remaining CaM-2-ext protein, suggesting that it may help to dock it to the lipophilic membrane surrounding.

Ca(2+) sensor proteins in dendritic spines: a race for Ca(2+)

Dendritic spines are believed to be micro-compartments of Ca²⁺ regulation. In a recent study, it was suggested that the ubiquitous and evolutionarily conserved Ca²⁺ sensor, calmodulin (CaM), is the first to intercept Ca²⁺ entering the spine and might be responsible for the fast decay of Ca²⁺ transients in spines. Neuronal calcium sensor (NCS) and neuronal calcium-binding protein (nCaBP) families consist of Ca²⁺ sensors with largely unknown synaptic functions despite an increasing number of interaction partners. Particularly how these sensors operate in spines in the presence of CaM has not been discussed in detail before. The limited Ca²⁺ resources and the existence of common targets create a highly competitive environment where Ca²⁺ sensors compete with each other for Ca²⁺ and target binding. In this review, we take a simple numerical approach to put forth possible scenarios and their impact on signaling via Ca²⁺ sensors of the NCS and nCaBP families. We also discuss the ways in which spine geometry and properties of ion channels, their kinetics and distribution, alter the spatio-temporal aspects of Ca²⁺ transients in dendritic spines, whose interplay with Ca²⁺ sensors in turn influences the race for Ca²⁺.

Between promiscuity and specificity: novel roles of EF-hand calcium sensors in neuronal Ca2+ signalling

In recent years, substantial progress has been made towards an understanding of the physiological function of EF-hand calcium sensor proteins of the Calmodulin (CaM) superfamily in neurons. This deeper appreciation is based on the identification of novel target interactions, structural studies and the discovery of novel signalling mechanisms in protein trafficking and synaptic plasticity, in which CaM-like sensor proteins appear to play a role. However, not all interactions are of plausible physiological relevance and in many cases it is not yet clear how the CaM signaling network relates to the proposed function of other EF-hand sensors. In this review, we will summarize these findings and address some of the open questions on the functional role of EF-hand calcium binding proteins in neurons.