Computational analysis of calcium signaling and membrane electrophysiology in cerebellar Purkinje neurons associated with ataxia
© Brown and Loew; licensee BioMed Central Ltd. 2012
Received: 11 March 2012
Accepted: 16 May 2012
Published: 15 June 2012
Mutations in the smooth endoplasmic reticulum (sER) calcium channel Inositol Trisphosphate Receptor type 1 (IP3R1) in humans with the motor function coordination disorders Spinocerebellar Ataxia Types 15 and 16 (SCA15/16) and in a corresponding mouse model, the IP3R1delta18/delta18 mice, lead to reduced IP3R1 levels. We posit that increasing IP3R1 sensitivity to IP3 in ataxias with reduced IP3R1 could restore normal calcium response. On the other hand, in mouse models of the human polyglutamine (polyQ) ataxias, SCA2, and SCA3, the primary finding appears to be hyperactive IP3R1-mediated calcium release. It has been suggested that the polyQ SCA1 mice may also show hyperactive IP3R1. Yet, SCA1 mice show downregulated gene expression of IP3R1, Homer, metabotropic glutamate receptor (mGluR), smooth endoplasmic reticulum Ca-ATP-ase (SERCA), calbindin, parvalbumin, and other calcium signaling proteins.
We create a computational model of pathological alterations in calcium signaling in cerebellar Purkinje neurons to investigate several forms of spinocerebellar ataxia associated with changes in the abundance, sensitivity, or activity of the calcium channel IP3R1. We find that increasing IP3R1 sensitivity to IP3 in computational models of SCA15/16 can restore normal calcium response if IP3R1 abundance is not too low. The studied range in IP3R1 levels reflects variability found in human and mouse ataxic models. Further, the required fold increases in sensitivity are within experimental ranges from experiments that use IP3R1 phosphorylation status to adjust its sensitivity to IP3. Results from our simulations of polyglutamine SCAs suggest that downregulation of some calcium signaling proteins may be partially compensatory. However, the downregulation of calcium buffer proteins observed in the SCA1 mice may contribute to pathology. Finally, our model suggests that the calcium-activated voltage-gated potassium channels may provide an important link between calcium metabolism and membrane potential in Purkinje cell function.
Thus, we have established an initial platform for computational evaluation and prediction of ataxia pathophysiology. Specifically, the model has been used to investigate SCA15/16, SCA1, SCA2, and SCA3. Results suggest that experimental studies treating mouse models of any of these ataxias with appropriately chosen peptides resembling the C-terminal of IP3R1 could adjust receptor sensitivity, and thereby modulate calcium release and normalize IP3 response. In addition, the model supports the hypothesis of IP3R1 supersensitivity in SCA1.
KeywordsModel IP3R Homer CAG repeat instability KCa channels Ataxia Computational BK channel Spinocerebellar Virtual cell
Ataxic mouse models relevant to this study
Mutated gene product
Compensation *Observed/**Suggested Therapeutic
Large conductance calcium-activated voltage-gated potassium channel
Altered Precision Timing
Inositol triphosphate receptor 1
Inositol triphosphate receptor 1
Inositol triphosphate receptor 1
Inositol triphosphate receptor 1
Inositol triphosphate receptor 1
Inositol triphosphate receptor 1
Plasma Membrane Calcium ATPase
The mechanism of pathology is different in polyglutamine ataxias. Humans with spinocerebellar Ataxia 1 (SCA1), spinocerebellar ataxia 2 (SCA2), and spinocerebellar ataxia 3 (SCA3), and their corresponding mouse models, have large numbers of unstable polyglutamine (polyQ, CAG) repeats in the Ataxin-1 [23, 24], Ataxin-2 , and Ataxin-3 [26–28] proteins, respectively. All three mouse models show increased calcium release in response to IP3, relative to wild type mice [6, 8, 10], with results from SCA2 and SCA3 suggesting a mechanism involving supersensitive IP3R1 [6, 10, 18]. SCA1 and SCA3 mice have reduced expression of IP3R1 and other glutamatergic signaling proteins , in addition to increased IP3R1-mediated calcium release, relative to wild type mice . This paradox is yet to be explained . Gardner suggests that SCA15/16 in humans found to have actual mutations in ITPR1, the gene that encodes the IP3R1 protein, be termed ITPR1-associated ataxia. However, we have termed a subset of ataxias, not just with reduced IP3R1 but also with supersensitive IP3R1, and regardless of the mutated gene (e.g., ITPR1, Ataxin-1, Ataxin-2, Ataxin-3) as ‘IP3R1-associated ataxias’ .
A prior modeling study by Hernjak et al. explores the significance of experimental observations in normal cerebellar Purkinje neurons, which show markedly lower sensitivity and higher abundance of IP3R1 relative to other cell types [33–35]. This is the opposite of what is seen in SCA3 and suggested in SCA1. Hernjak et al. demonstrated that the normally high levels of IP3R1 ensure robust calcium signals in an individual spine, and that the normally low sensitivity of IP3R1 to IP3 confines the calcium transients to the stimulated spine, without spillover to adjacent spines . Thus, they argued that IP3R1 abundance and sensitivity have independent effects on calcium signals and that changes in one do not appear to appropriately compensate for abnormalities in the other. In this work, however, we explore this argument by computationally adjusting IP3R1 sensitivity to IP3 in SCA15/16 (in the context of reduced IP3R1 levels) to see if interventions of this type can actually restore normal calcium response indistinguishable from wild type. We also suggest that, in polyglutamine SCAs with hyperactive IP3R1, experimentally observed downregulation of IP3R1 and other calcium signaling proteins (but not calcium buffers) could serve to partially restore normal calcium transients in ataxias.
We test our hypotheses in a novel version of a recently developed computational model that links biochemistry with electrophysiology . We simulate the application of different experimentally available peptides that resemble various portions of the C-terminal of IP3R1, in order to model competitive inhibition of the binding of mutant ataxin proteins. Improvements are made to the model, resulting in an ability to look at changes in membrane potential in response to IP3R1-mediated calcium release or supralinear calcium transients. We show that adjusting IP3R1 sensitivity to IP3, using peptides resembling the C-terminal of IP3R1, could potentially stabilize calcium, without altering fine-tuning of supralinear calcium release required for normal Purkinje neuron function and long term-depression. The model also demonstrates that decreased levels of IP3R1 and various other calcium signaling proteins (except for calcium buffers) experimentally observed in SCA1 mice may partially compensate for the IP3R1 supersensitivity suggested by Liu et al. . Finally, we link IP3R1 activity with the influence of calcium on membrane electrophysiology. Results from our study may provide further insight into the roles that abundance and sensitivity of the calcium channel IP3R1 play in cerebellar function and coincident detection at IP3R1. Results may also provide a foundation for future studies on SCA1, SCA2, SCA3, SCA15, SCA16, and other IP3R1-associated ataxias, which all involve aberrant calcium signaling.
Merged biochemical and electrophysiological modeling
The cerebellar Purkinje neuron geometry is complex. In order to combine detailed biochemistry with electrophysiology in tractable Virtual Cell models (VCell; http://www.vcell.org[37–41], we developed a geometry reduction algorithm preserved path reduction or PPR. This algorithm simplifies the Purkinje neuron geometry to just the essential dendritic segments and branches in a path between a prototypical dendritic spine and the soma, while maintaining the dynamic electrical features of the full geometry in these respective compartments. The resulting model has only 17 compartments and allows us to readily couple the biochemistry in the spine compartment with the electrophysiology of the entire system. Reference  provides full details of the algorithm and the resultant geometry.
The baseline VCell model for normal physiology of the system was also previously published by our laboratory . This compartmental model focused on IP3R1-mediated calcium signaling in the Purkinje cell, taking into account the increased abundance and decreased sensitivity of the Purkinje cell IP3R1 relative to other cell types. The input to the system is provided as a train of IP3 pulses (representing IP3 production from multiple parallel fiber stimuli) with variable parameters that control the number, amplitude and duration of the IP3 pulses. In addition, the model included influx of calcium representing the opening of voltage-gated calcium channels when the Purkinje neuron is stimulated by a climbing fiber. It therefore allowed the study of coincident activation of the Purkinje spine by multiple parallel fibers and a single climbing fiber. The model also included the downstream effects leading to IP3R1-mediated calcium release, as well as calcium buffering and SERCA, which resides on the sER. All of the spine cytosol species are free to diffuse through the spine neck into the adjacent dendrite. IP3 also diffuses in the cytosol to bind its calcium channel-coupled receptor on the sER membrane. Calcium is released from sER stores into the cytoplasm and also co-activates IP3R1, as well as binds to various buffers. Some of this calcium is also pumped back into the sER by SERCA. The model definition in VCell includes specifying reactions, membrane transport mechanisms, stoichiometry, and kinetic rate expressions that correspond with these processes.
We subsequently created a 3D spatial model that explored  the membrane signaling mechanisms upstream of IP3R. We distinguished between candidate sources of sufficient phosphatidylinositol-4,5-bisphosphate (PIP2) in the spine for the production of requisite IP3. Local PIP2 sequestration provided an efficient means of fine-tuning coincident activation of parallel fiber and climbing fiber stimuli in cerebellar Purkinje neuron spines. The results reported in this study concluded that local high PIP2 achieved through a sequestration mechanism could provide sufficient IP3 for supralinear calcium release.
We used optimized parameters  to study physiological effects of varying IP3R1 sensitivity and abundance with a physiologically appropriate IP3 signal amplitude and duration. To manipulate IP3R1 abundance and sensitivity, we adjusted the parameters Jmax and dIP3, respectively, in our comprehensive model (see a table of parameters in ). Small Jmax values represent reduced IP3R1 abundance. The default value is from Hernjak et al. , and is 10x larger than the IP3R1 determined for neuroblastoma cells, to account for the increased levels of IP3R1 in Purkinje neurons. The dIP3 parameter represents the dissociation constant of IP3 from IP3R1, so that high dIP3 values in our model represent low affinity of IP3R for IP3. We use the apparent affinity of IP3R1 for IP3 as a measure of sensitivity, as in the experimental study by Tang et al. . Manipulation of these IP3R1parameters was done in the context of a train of 4 parallel fiber (PF) stimuli followed by a single climbing fiber (CF) stimulus, following the optimal protocol designed by Wang et al. . This allowed us to examine the effects of IP3R1 insufficiency and supersensitivity in the context of Purkinje spine stimulation patterns similar to those in published experimental studies [35, 42, 44–46].
“Brown et al. 2011 Analysis of SCA15-SCA16 with IC4 Peptide Application”,
“Brown et al. 2011 SCA1 Compensation Analysis”,
“Brown et al. 2011 Analysis of SCA1-SCA2-SCA3 with IC-G2736X Peptide”,
“Brown et al. 2011 Purkinje Biochem-Electrophysiol SCA”,
“Brown et al. 2011 Purkinje 3D Spine – Spinocerebellar Ataxia study”, and
“Brown et al. 2011 Combined Purkinje - Several Spines - Current Injection”.
Adapting the PPR model to the SCA model
In order to investigate how biochemical signaling at the spine could affect electrophysiology, we needed to allow calcium release to influence the voltage-gated KCa channels. In the PPR model, KCa channels were only influenced by calcium influx into a phenomenological submembrane shell. This allowed for very high transient concentrations of calcium (500 μM), as well as a high concentration of resting calcium (~ 4 μM) near the inner surface of the plasma membrane. Since such a high amount of basal calcium is not physiologically observed or observable (in the absence of Purkinje neuron stimulation by nerve endings from regulatory cells) and since calcium concentration beyond the submembrane shell was not explicitly considered in the PPR model, we felt it was important to model calcium less phenomenologically. Therefore, in this work, we removed submembrane shell calcium from the model only in the spine and the adjacent dendrite, where calcium signaling biochemistry is modeled explicitly. The adjacent dendrite was adjusted in tandem with the spine, because calcium and other molecules involved in biochemical processes are allowed to diffuse freely between the spine and the adjacent dendrite, restricted only by the geometry of the spine neck. In a previous study, Hernjak et al.  used a one-dimensional model to determine a dendritic distance, 28 μm from a stimulated spine, at which biochemical concentrations remained at steady state during synaptic activity. The model therefore assumes that there is no net molecular flow between compartments, except for the adjacent dendrite (given a length of 28 μm) and the spine.
In other compartments, calcium influx feeds only into submembrane shell calcium. In the spine and adjacent dendrite, calcium influx is allowed to directly alter concentration of calcium involved in biochemical signaling; we term this calcium ‘global calcium’. The concept of the submembrane shell in electrophysiological models [47–49] assumes a characteristic relaxation time of 2 ms. This relaxation time represents the time it could take for submembrane shell calcium to relax to the ‘global calcium’ steady state, 40nM in the model. The submembrane shell relaxation time conceptually encompasses: (i) diffusion from beneath the membrane into the rest of the compartment, for example, the spine head, (ii) diffusion into adjacent compartments, for example, from the spine to the adjacent dendrite, (iii) buffering by calcium binding proteins, for example, calbindin and parvalbumin, (iv) uptake into the smooth endoplasmic reticulum by SERCA, and (v) extrusion through the plasma membrane by the Na-Ca exchanger. In addition, calcium extrusion from the cell by the plasma membrane Ca-ATPase (PMCA, calcium pump) is modeled explicitly as a Michaelis-Menten expression.
In order to properly replace submembrane shell calcium with global calcium in the spine and the adjacent dendrite, (i) – (v) are all modeled explicitly. The model by Hernjak et al.  assumed instantaneous diffusion in the spine, and included diffusion between the spine and adjacent dendrite. It also included explicit buffering by calbindin and parvalbumin, in addition to calcium uptake by SERCA. That model did not include calcium extrusion from the spine or adjacent dendrite. Here, we add in calcium efflux by including the Michaelis-Menten Ca-ATPase from the Purkinje neuron electrophysiology models [48, 49] and a Na-Ca exchanger from a study on cerebellar granule cells  (a similar Na-Ca exchanger in a Purkinje neuron model was not available). The Na-Ca exchanger is not electrogenic and simply provides an additional means of calcium extrusion; its kinetics are unchanged in our model. The total concentration of the Ca-ATPase is increased to allow for calcium efflux appropriate for maintaining a physiological steady state calcium concentration in the Purkinje neuron spine. The kinetics for the voltage-gated calcium-activated potassium channels were adjusted to be influenced by global calcium. This reflects results from experimental studies which suggest that BK channels can be activated by global calcium concentrations less than 10 μM .
The ICpeptide study
Parameters for new features added to model
Forward rate of IC4 binding PP1α
Reverse rate of IC4 binding PP1α
Forward rate of IP3R1 C-terminal binding PP1α
Reverse rate of IP3R1 C-terminal binding PP1α
Forward rate of IC4 binding mAtaxin
Reverse rate of IC4 binding mAtaxin
Forward rate of IP3R1 C-terminal binding mAtaxin
Reverse rate of IP3R1 C-terminal binding mAtaxin
Mutant Ataxin PolyQ Protein
Directly binds tip of IP3R1 C-terminal
Tip of IP3R1 C-terminal
Competitively inhibits PP1α binding
Associates with polyQ region
Competitively inhibits polyQ binding
The C-terminal of IP3R1 can also bind mutant Ataxin. The affinity of the IP3R1 C-terminal tip was taken to be greater for the mutant protein than for PP1α. This was done to incorporate two phenomena into the binding reaction. First, planar lipid bilayer reconstitution experiments from Chen et al. and Liu et al. [6, 10] suggest that direct association of the mutant ataxin protein with the C terminal of IP3R1 leads to supersensitivity of the receptor. Second, the binding domain of the mutant protein  on IP3R1 overlaps the PP1α binding motif  on IP3R1. Mutant protein binding IP3R1 in situ (as in cerebellar slices from mouse models that can be used to test simulation predictions; see Discussion) may decrease PP1α interaction with the receptor. These two factors are effectively modeled by assigning a larger affinity for mutant protein binding to IP3R1 than PP1a binding the receptor. Affinities were captured by the rate constants, all available in Table 2.
Two main categories of ataxias in humans and mouse models are considered in this study: those with reduced IP3R1 and those with supersensitive IP3R1 (Table 1). However, some ataxias involve both phenomena, adding a layer of complexity. Further, a common component of all these ataxias is abnormal calcium metabolism.
SCA15/16: pathophysiology of IP3R1 insufficiency
Application of ICpeptide to adjust sensitivity in SCA15/16
Suggested use of ICpeptides relevant to this study
Minimal binding region of target
Use suggested by this study
in vivo competitive inhibition of PP1α
(Tang et al., 2003b,
with reduced levels of IP3R1
Tu et al., 2004)
in vitro/vivo competitive inhibition of mutant Ataxin
(Tang et al., 2009)
with supersensitive IP3R1
possible alternative to IC-G2726X
(Tang et al., 2009)
The simulations were repeated with the climbing fiber stimulus applied 150 ms after the start of the train of the PF stimuli, unlike the previous results that have a timing difference of 50 ms. A 150 ms timing between the CF and PF stimuli falls outside the optimal time window suggested by both experiments [44, 56] and modeling . Figure 3b shows that with this 150 ms time difference, supralinear calcium release is lost in the wild type case, as well as in the heterozygous case. This suggests that restoration of normal calcium release in the model does not alter fine-tuning of coincidence detection suggested by Brown et al. . Thus, increasing IP3R1 sensitivity could be explored as a therapeutic tactic in SCA15/16 mouse models, e.g. IP3R1delta18/delta18 or IP3R1+/− mice.
SCA2: application of ICpeptide in ataxias with pathologically increased sensitivity
Tang et al. also found that the huntingtin binding site on IP3R1 is immediately upstream of the PP1α-binding motif on the tip of the IP3R1 C-terminal . The predicted minimal region on IP3R1 for huntingtin polyglutamine repeat binding (F2627-G2736) actually overlaps a small portion (5 aa) of the minimal PP1-binding motif (R2731-A2749). The huntingtin protein binding IP3R1 may therefore physically occlude or partially hinder binding by PP1αlpha in vivo and in situ. This would reduce PP1αlpha dephosphorylation of the receptor, leaving it in a hyperphosphorylated and therefore supersensitive state. A similar mechanism should occur with mutant Ataxins in SCA2  and SCA3  and also SCA1 [8, 10, 57], given that the mutant Ataxin and huntingtin proteins all have similarly expanded CAG repeats.
Figure 4 shows that adding 2 μM mutant Ataxin, e.g. mutant Ataxin-2 in SCA2, to the wild type model triples the increase in simulated fluorescence of the calcium indicator. Furthermore, the IP3R1 supersensitivity allows calcium concentration to start rising even before coincidence with calcium influx. In addition, we simulated the application of a different ICpeptide in our model. This ICpeptide, IC-G2736X (D2590-G2736)(Table 3), was also created by Tang et al., and encompasses the part of IP3R1 that encodes the minimal binding region for the expanded polyglutamine repeats in huntingtin, and likely the mutant ataxin proteins . Further, Tang et al. used a β–galactosidase assay to show that IC-G2736X strongly associates with the mutant polyQ huntingtin . In fact, they also used viral infection to insert IC10 into striatum of ataxic HD mice, and observed consequent stabilized calcium signaling, as well as neuroprotective effects . In our model, IC-G2736X competitively binds the mutant Ataxin protein and blocks mAtaxin from binding endogenous IP3R1. The normal sensitivity of IP3R1 to activation by IP3 is therefore restored, and IP3R1 is no longer supersensitive. Figure 4 shows the restoration of normal supralinear calcium release. This result allows us to predict that experimental application of such an ICpeptide in cerebellar slices from mice with ataxias associated with supersensitive IP3R1 and increased calcium release may similarly restore normal calcium response. We expect that various concentrations of ICpeptide could be needed to normalize IP3R1 sensitivity and consequently calcium release in SCA1, SCA2, or SCA3, based on the degrees to which IP3R1 sensitivity to IP3 is pathologically increased. This could depend on the concentration of the mutant protein in the spine, as well as other unknown factors.
SCA1 (& SCA3): compensation versus pathology
Although expression of IP3R1 and other calcium signaling proteins is reduced in SCA3 mice  and in SCA1 mice and humans [9, 11], reduced calcium release is not experimentally observed in either disease. It is therefore not understood what role reduced IP3R1 plays in the pathology of these polyQ diseases. Figure 3 suggested that in SCA15/16 and other ataxias with reduced IP3R1 abundance and consequently reduced calcium release, therapeutic strategies to increase IP3R1 sensitivity to activation by IP3 could restore normal calcium release. This begs the question of whether, conversely, the downregulation of IP3R1 in ataxias with supersensitive IP3R1 could partially or wholly restore normal calcium response. A glimpse into a compensatory effect could lie in examining actual experimentally observed relative expression levels in SCA1 [9, 11], and using those observed levels in the model with supersensitive IP3R1. This allows the modeling and simulations to be driven by experimental observations. Accordingly, a 1.8-fold decrease (approximated from Lin et al. and Serra et al.) [9, 11] was simulated in the model for SCA1 pathology.
Expression of calcium-binding proteins, parvalbumin and calbindin, are also experimentally observed to be reduced [9, 11]. It is not known whether this contributes to pathology or compensation. Calbindin knockout mice crossed with SCA1 mice worsens the phenotype . This suggests that low levels of these buffer proteins could accelerate pathology. It is possible that normal or high levels of buffer proteins may help to maintain homeostasis in the presence of supersensitive IP3R1 early on in the SCA1 disease, then, when their expression is reduced, homeostasis is lost. Modeling results support this suggestion. Tuning down the calcium buffers parvalbumin and calbindin in the model makes supralinear calcium more robust than with down regulation of IP3R1, mGluR and SERCA alone (Figure 5a,b). Accordingly, the downward trend in the peak concentration of calcium is reversed with down regulation of calcium buffers. This is consistent with studies in Purkinje cells from parvalbumin and calbindin knockout mice that are ataxic [60, 61], where peak amplitudes of calcium transients evoked by climbing fiber stimulation were double those in wild type mice . Notably, it has been suggested that because reduced expression of buffer proteins is observed in SCA1 mice before onset of symptoms, the timing of buffer protein downregulation may help determine the onset of symptoms [63, 64].
Notably, Liu et al. report that mutant SCA1 pathologically associates with IP3R1 , as do SCA2 and SCA3 [6, 10], suggesting a common supersensitive IP3R1 pathology, but they do not actually measure or report changes in IP3R1 sensitivity. Thus, we thought it was important to consider changes in calcium resulting from signaling protein downregulation, even if hyperactive IP3R1 is not assumed. Figures 5c and 5d illustrate that calcium responses resulting from IP3R1 and other calcium signaling protein downregulation, but with normal IP3R1 sensitivity, are less robust than wild type. This indicates that downregulation of these proteins in the absence IP3R1 supersensitivity does not appear consistent with the elevated calcium response relative to wild type, as observed in SCA1 mice . The model therefore supports the hypothesis of IP3R1 supersensitivity in SCA1.
Influence of calcium on voltage-gated calcium-activated potassium channels
To better understand the interplay between calcium signaling and electrophysiology, we investigated the relationships between calcium influx, IP3R1-induced calcium release, and the voltage-gated calcium-activated potassium (KCa) channels. Our model includes the BK channel (‘big’ conductance calcium-activated voltage-gated potassium channel) and IK (intermediate conductance calcium-activated voltage-gated potassium channel) channels. Since more is known about the BK channels, we focused our study design and interpretation on these channels (see schematic in Figure 1b). The BK channel in the cerebellum is thought to contribute to repolarization of membrane potential transients in the dendrites, particularly calcium spikes (transient depolarization of membrane potential in spines and dendrites) resulting from the opening of voltage-dependent calcium channels in the dendrites . In addition, experimental results by Sausbier et al. using ataxic BK knockout mice suggest that these channels are also important for the afterhyperpolarization of action potentials at the soma . BK channels contribute to precision timing, and are thus important for normal firing of the cerebellar Purkinje neurons [65–70].
IP3R1 interacts closely with the BK channel, in glioma cells . The BK channels seem to be in lipid rafts in the plasma membrane apposed to the sER. It is thought that in other cell types, including some neurons, the BK channels may form physical complexes with various plasma membrane calcium channels, bringing them within a few nanometers of the calcium channel pores . We suspect that in the cerebellar Purkinje neuron spines, BK activity is also closely linked with the IP3R1-mediated calcium release. This idea is supported by direct activation of the BK channel by IP3R1-mediated calcium release in pyramidal neurons  and cereberal artery smooth muscle cells . We therefore set out to investigate possible activation of BK channels (and IK channels) by IP3R1-mediated calcium release.
Taken together, these results suggest that in Purkinje neurons, coincident activation of the Purkinje spine by current injection at the soma and IP3R1-mediated calcium release with or without concurrent CF stimulus can alter membrane excitability. The results indicate that KCa activity increased by calcium release led to earlier action potential spikes than in the absence of biochemically induced calcium transients. This is reminiscent of a study by Sausbier et al.  in which control mice with wild type levels of BK activation showed more spontaneous discharge relative to ataxic BK−/− mice. This is presumably due to another important result from that study: BK activity generates afterhyperpolarization of sodium (action potential) spikes. Thus, augmented BK channel activity polarizes the plasma membrane and increases excitability by more quickly resetting the sodium channels.
Computational models that combine the cell biology of calcium signaling with membrane electrophysiology involving calcium channels and calcium-activated ion channels can reveal emergent properties that are not observed in models of purely biochemical or purely electrophysiological models. The models that are presented in this work can be useful for studying these emergent properties in normal physiological Purkinje neurons, as well as in ataxia. Future studies can be developed to study spontaneous Purkinje neuron firing rate or precision of spike timing and interspike intervals in various classes of ataxias with mutations in genes that code for calcium channels [9, 76–78] and potassium channels [65, 79], as well as IP3R1 or related proteins (including SCA1, SCA2, SCA3, SCA15/16) [3, 4, 6, 7, 9–16].
We determined that the Virtual Cell modeling environment  can be used to study cellular dynamics associated with disease states, particularly spinocerebellar ataxias with altered IP3R1 abundance or sensitivity or both. We modified a quantitative multicompartmental model recently developed in Virtual Cell to investigate the coupling of detailed biochemistry of IP3 and calcium signaling with membrane potential. Results show that computational modeling allows us to simulate the behavior of calcium dynamics and membrane potentials in neurons at the level of the spine and of the soma, in response to various stimuli. Our models include a myriad of ion channels and molecules important for normal Purkinje neuron function. They merge features of IP3 production, phosphoinositide and calcium signaling with conceptual representations of parallel fiber and climbing fiber stimulation, local PIP2 synthesis and sequestration, PIP2 hydrolysis, calcium influx, coincident detection at IP3R1, modulation of IP3R1 abundance and sensitivity, IP3R1-mediated supralinear calcium release, and components of membrane electrophysiology into comprehensive models that can have compartmental and spatial counterparts.
Results suggest that ICpeptides may be used to modulate calcium release in various IP3R1-associated ataxias (ataxia s in which IP3R1 is less abundant, e.g., SCA15/SCA16 as well as SCA1 and SCA3, or more sensitive than in wild type mice or unaffected individuals, e.g., SCA2 and SCA3 and likely SCA1). Simulations predict that pathophysiology associated with decreases in IP3R1 abundance can be overcome by adjusting IP3R1 sensitivity to activation by IP3, and that targeting PP1α can restore normal calcium transients. Various concentrations of ICpeptide may be needed to increase IP3R1 sensitivity in SCA15/SCA16, based on the levels of reduced abundance. Similarly, various concentrations of ICpeptide could be needed to normalize IP3R1 sensitivity and consequently calcium release in SCA1, SCA2, and SCA3, based on the degrees to which IP3R1 sensitivity to IP3 is pathologically increased (dependent on the concentration of mutant protein and the relative balance with PP1α).
This project complements results by Hernjak et al. from a study on Purkinje spine calcium signaling that considered the importance of low sensitivity and high abundance of IP3R1 suggested by experiments . Findings from that study suggested that increasing IP3R1 sensitivity to IP3 could not rescue calcium release in the context of low IP3. However, our results show that increased sensitivity can restore normal IP3 response if the abundance is not too low. This provides further insight into the roles that IP3R1 abundance and sensitivity play in normal cerebellar functioning and coincident detection at IP3R1. IP3R1 abundance and sensitivity can each function as the primary pathology in IP3R1-associated ataxias or can serve to partially or wholly compensate the effects of pathology. Other important conclusions from the study by Hernjak et al. were that the low IP3R1 sensitivity observed in wild type Purkinje neurons relative to other cell types may be important to restrict biochemical signals and synaptic plasticity to one spine, and that the high IP3R1 abundance in Purkinje neurons is important for ensuring generation of a robust calcium signal in an individual spine. The current study confirms that the balance of IP3R1 abundance and sensitivity is critical for obtaining robust, but not hyperactive, calcium transients. It is possible that in ataxias with reduced abundance of IP3R1, some spines will not generate robust calcium signals, while in ataxias with supersensitive IP3R1, stimulation signals may not be confined to single spines. This could affect both compartmentalization of biochemical signals, and the influence of the biochemistry on membrane electrophysiology. Similarly, therapeutically increased IP3R1 sensitivity (by application of an appropriate ICpeptide) could potentially lead to spillover of calcium into adjacent spines, if sensitivity is not finely tuned. Modeling could thus be useful for helping to determine the adequate levels of ICpeptide to be administered experimentally.
Experimental observations in SCA1 mice appear to be paradoxical  (Table 1). The expression of various molecules (for example, IP3R1, Homer, and SERCA) involved in glutamatergic calcium signaling and IP3R1-mediated calcium release is reduced in mice and humans with SCA1 [9, 11]. However, IP3-induced calcium release is hyperactive in SCA1 mice . Insight may be gained by examining the pathology of SCA2 and SCA3. Liu et al.  and Chen et al.  found that mutant Ataxin-2 and Ataxin-3, respectively, bind the C-terminal of IP3R1 and increase the receptor’s sensitivity to activation by IP3. This corresponds to hyperactive IP3-induced calcium release in both SCA2 and SCA3. Above normal numbers of CAG repeats in SCA1, SCA2, and SCA3 give each respective mutant ataxin protein a toxic gain of function that disrupts calcium homeostasis in neurons. Thus, there may be a common mechanism underlying some of the pathology of the three polyQ ataxias. This mechanism could be binding of the polyQ-expanded protein to the C-terminal of IP3R1, thereby increasing the receptor’s sensitivity to IP3. This mechanism may be extendable to a large variety of polyglutamine diseases. In fact, Liu et al. reported that the mutant Ataxin-1 protein and polyQ-expanded atrophin-1, the protein mutated in Dentatorubral-pallidoluysian atrophy (DRPLA), both associate with the C-terminal of IP3R1 , though IP3R1 supersensitivity has not yet been assessed experimentally in either disease. Of note, a prominent feature of DRPLA is cerebellar ataxia[80–83].
The pathological gain of function of IP3R1 sensitivity that destabilizes calcium signaling  is also observed in a mouse model of Huntington’s disease (HD), another polyQ disease [18, 52, 53, 85–87]. Lentiviral and adenoviral infection of affected cultured medial striatal neurons or in the striatum of the ataxic HD mice themselves with the IC10 peptide lead to reduced brain atrophy and improved motor coordination, respectively . Perhaps similar genetic targeting of Purkinje neurons in the cerebellum of existing mouse models with IC-G2736X could mitigate the pathology of SCA1, SCA2, and SCA3. This could first be tested in in situ experiments using cerebellar slices from existing mouse models, in which mutant ataxin1 binding IP3R1 directly upstream of the PP1αlpha binding site, would be close enough to preclude binding of PP1αlpha.
Reduced gene expression of IP3R1, Homer/MyoVa, mGluR and other molecules in SCA1 may be due to direct pathology due to the effect of the mutant proteins on transcription [7, 9]; it could also be due to the tight regulation of calcium homeostasis in the cerebellar Purkinje neuron . Lin et al. suggest that reduced gene expression is a part of the pathology , but our model indicates that with exception of the calcium buffer proteins the effect on IP3R1-mediated calcium release is partially compensatory. Using the reduced fold expression observed in SCA1, normal calcium release was restored to varying degrees, depending on how much mutant Ataxin protein was placed in the model. Reduced expression of several calcium signaling genes is also observed in some Purkinje neuron subtypes in plasma membrane calcium ATP-ase (PMCA) knockout mice . This could also have the effect of partial compensation in that mouse model, since PMCA contributes to expelling calcium from the Purkinje neuron cytoplasm. ‘Compensatory pathology’ is thus an interesting characteristic of various ataxias, and provides a level of complexity to the study of these neurological disorders. Disruption of calcium signaling and homeostasis, whether due to reduced IP3R1 levels or supersensitive IP3R1, whether in IP3R1-associated ataxias or in other ataxias, can lead to dysfunction of Purkinje cells, and impaired long-term depression and synaptic plasticity that are involved in learning and memory.
Calcium homeostasis is critical for normal function of Purkinje neurons and is thought to be tightly regulated . Downregulation of parvalbumin and calbindin could therefore lead to a loss of compensation, by destabilizing supralinear calcium release and disrupting any homeostasis achieved by downregulation of other molecules. Furthermore, it has been postulated that normal or high concentration of calcium buffer proteins have a protective role in certain neurons . Vig et al. suspect that the decrease in parvalbumin expression in Purkinje cells from SCA1 patients may reflect alterations in a regulatory biochemical pathway that may be important for neuronal survival . Further, our model suggests that these changes (downregulation of various calcium channels and buffers) alone may not be sufficient to reproduce elevated calcium response as observed experimentally in SCA1 mice, without assuming concurrent supersensitive IP3R1.
Purkinje neurons are largely spared in SCA3 , while nerve cells in the pons and substantia nigra are substantially damaged. Chou et al. suggest that although prominent neuronal loss was not found in the cerebellum, the SCA3 mice displayed pronounced ataxic symptoms, suggesting that instead of neuronal demise, mutant Ataxin3 causes neuronal dysfunction and resulting ataxia. While Ataxin1 is found in the nucleus, Ataxin 2 and Ataxin3 are cytoplasmic proteins, under normal conditions. However, in brains from patients with SCA3, mutant Ataxin3 accumulates in the nucleus [28, 89], as does mutant Ataxin1 . Hence, results from our study of SCA1 could be extended to SCA3 modeling.
Long term feeding of SCA2 and SCA3 mice with dantrolene improved motor coordination and slowed brain atrophy [6, 10]. Dantrolene is thought of as a ‘calcium stabilizer’ [6, 10, 91] and has been shown to inhibit the ryanodine receptor (RYR), which is another calcium channel on the smooth endoplasmic reticulum [92, 93]. However, all the targets of dantrolene are not known . The details of the mechanism of action of dantrolene are incomplete , though it has been proposed as a possible therapeutic drug for SCA2  and SCA3 . The drug is currently approved to treat malignant hypothermia as a one-time application in response to adverse reaction due to anesthesia [94–96]. It is also used to reduce muscle spasticity in patients with neurological incidents or disorders. However, dantrolene leads to fluid buildup in the lungs, among other adverse effects . Yet, considering the disturbed neuronal calcium signaling observed in polyQ ataxias, it is likely that inhibiting or downregulating either or both of the sER intracellular calcium release channels (IP3R1 and RYR) should attenuate IP3R1-mediated supralinear calcium release into the cytosol and thereby alleviate SCA2 and SCA3. This is particularly so, since both the RYR and IP3R1 channels increase in function when directly bound by cytosolic calcium. Therefore the supralinear calcium release cascade initiated by IP3 binding IP3R1 may involve cytosolic calcium also binding RYR to increase calcium release from the sER. Adding RYR to the model in the future will facilitate further study of the signaling interactions among RYR, IP3R1, other calcium channels, and the calcium-activated potassium channels. However, it should be noted that there is evidence that RYR is localized to the dendritic shaft and is excluded from spines [97, 98]. Yet, there is also evidence that these physically separate calcium release sites functionally interact .
Predictions from our computer model involving SCA15/16 could be compared with novel experiments in an existing mouse model using GST-IC4 to dissociate PP1α from IP3R1 in Purkinje neurons in cerebellar slices from ataxic mice. This could restore normal calcium and membrane potential response in ataxic mice. For such experiments, it would be useful to select a mouse model that: (i) does not completely knock out IP3R1 expression as in the IP3R1−/− knockout mice, but has reduced expression of IP3R1 protein, (ii) shows motor discoordination, (iii) does not possess a mutation in neuronal IP3R1 at the preferred site of PKA phosphorylation (Ser-1755) , and (iv) does not possess an IP3R1 mutation in the carboxyl terminal PP1α-binding site (2731–2749) . For assessing IP3R1 haploinsufficiency in particular, as found in SCA15 and SCA16 patients [3, 12, 15, 16], it would also be useful to select a mouse model with a large heterozygous deletion mutation in IP3R1. The IP3R1 +/− mice [2, 5] fulfill all of these desired features. The IP3R1delta18/delta18 mice have also been suggested as a mouse model for SCA15/16 . The ITPR1opt/opt mouse model (see Figure 1b) shows reduced IP3R1 protein levels on Western blot, but has a mutation that deletes the preferred PKA phosphorylation site (Ser-1755) [4, 100] in neuronal IP3R1, and therefore violates (iv). Given that a very small region of the gene is deleted and that it occurs at such an important regulation site, we suspect that IP3R1 is dysregulated in these mice. As such, we expect that these mice possess cellular pathophysiology resulting from dysregulation in addition to effects of IP3R1 insufficiency. Consequently, it is yet unclear whether results from this study will be easily extendable to calcium signaling and membrane excitability in ITPR1opt/opt mice.
Results from our models suggest that IP3R1-mediated calcium release can activate voltage-gated KCa channels and thereby alter membrane excitability in the Purkinje neuron. This can have implications for ataxias that involve disruption of intracellular calcium homeostasis. It is possible that pathological alterations in calcium transients result in pathological activation of BK channels. This could lead to variations in the timing of action potentials and other electrophysiological events in the Purkinje neuron, which controls modulation of neurons downstream of the cerebellum.
In summary, we created a series of models incorporating biochemistry and electrophysiology that unify observations in various SCAs. These models employ several novel concepts and approaches and provide a framework for the study not only of IP3R1-associated ataxias, but of various SCAs involving mutations of other molecules in the model, such as potassium channels [65, 79, 101, 102] and calcium channels [29, 76, 77, 103].
Model results indicate that, in mouse models of various ataxias associated with activity of the calcium channel IP3R1, ICpeptides may be used to stabilize intracellular calcium concentration. Further, restoration of normal calcium release in the model does not alter fine-tuning of coincidence detection suggested by Brown et al. . The hypothesis of IP3R1 supersensitivity in SCA1 is supported by simulation results. Even more, IP3R1 downregulation experimentally observed in SCA1 mice may partially compensate for the receptor’s supersensitivity. Homer and MyoVa downregulation are further compensatory. However, downregulation of calcium buffer proteins accelerates pathology. The model demonstrates that IP3-mediated calcium release in the Purkinje neuron could activate voltage-gated KCa channels, namely BK and IK, and provides insight into the interplay between IP3R1 sensitivity and abundance in the function and dysfunction of the Purkinje cell. Results help to explain experimental findings in mice, and can be used to make predictions for further experiments, which may ultimately be translated to ataxic individuals with reduced IP3R1 protein levels or increased sensitivity. IP3R1 abundance and sensitivity are components involved in calcium signaling, but by no means the only factors involved in the signaling systems of these SCAs.
We thank Dr. James Schaff, Dr. Ion Moraru, Fei Gao, and Anuradha Lakshminarayana for support in various aspects of this work. We are grateful to Dr. Arthur Hand for help with analysis of published electron microscope images, and to Dr. Louise McCullough for reading the manuscript. We are also pleased to acknowledge Dr. Ann Cowan, Dr. Raquell Holmes, and other members of the R. D. Berlin Center for Cell Analysis & Modeling at the University of Connecticut Health Center for helpful discussions.
This project was supported by grants from the National Center for Research Resources (2P41RR013186-15) and the National Institute of General Medical Sciences (9 P41 GM103313-15) from the National Institutes of Health.
- Goetz CG: Textbook of Clinical Neurology. St. Louis, MO: WB Saunders; 2003.Google Scholar
- Matsumoto M, Nakagawa T, Inoue T, Nagata E, Tanaka K, Takano H, Minowa O, Kuno J, Sakakibara S, Yamada M, et al: Ataxia and epileptic seizures in mice lacking type 1 inositol 1,4,5-trisphosphate receptor. Nature 1996, 379: 168-171. 10.1038/379168a0Google Scholar
- van de Leemput J, Chandran J, Knight M, Holtzclaw L, Scholz S, Cookson M, Houlden H, Gwinn-Hardy K, Fung H, Lin X, et al: Deletion at ITPR1 underlies ataxia in mice and spinocerebellar ataxia 15 in humans. PLoS Genet 2007, 3: e108. 10.1371/journal.pgen.0030108Google Scholar
- Street V, Bosma M, Demas V, Regan M, Lin D, Robinson L, Agnew W, Tempel B: The type 1 inositol 1,4,5-trisphosphate receptor gene is altered in the opisthotonos mouse. J Neurosci 1997, 17: 635-645.Google Scholar
- Ogura H, Matsumoto M, Mikoshiba K: Motor discoordination in mutant mice heterozygous for the type 1 inositol 1,4,5-trisphosphate receptor. Behav Brain Res 2001, 122: 215-219. 10.1016/S0166-4328(01)00187-5Google Scholar
- Chen X, Tang T, Tu H, Nelson O, Pook M, Hammer R, Nukina N, Bezprozvanny I: Deranged calcium signaling and neurodegeneration in spinocerebellar ataxia type 3. J Neurosci 2008, 28: 12713-12724. 10.1523/JNEUROSCI.3909-08.2008Google Scholar
- Chou A, Yeh T, Ouyang P, Chen Y, Chen S, Wang H: Polyglutamine-expanded ataxin-3 causes cerebellar dysfunction of SCA3 transgenic mice by inducing transcriptional dysregulation. Neurobiol Dis 2008, 31: 89-101. 10.1016/j.nbd.2008.03.011Google Scholar
- Inoue T, Lin X, Kohlmeier K, Orr H, Zoghbi H, Ross W: Calcium dynamics and electrophysiological properties of cerebellar Purkinje cells in SCA1 transgenic mice. J Neurophysiol 2001, 85: 1750-1760.Google Scholar
- Lin X, Antalffy B, Kang D, Orr H, Zoghbi H: Polyglutamine expansion down-regulates specific neuronal genes before pathologic changes in SCA1. Nat Neurosci 2000, 3: 157-163. 10.1038/72101Google Scholar
- Liu J, Tang T, Tu H, Nelson O, Herndon E, Huynh D, Pulst S, Bezprozvanny I: Deranged calcium signaling and neurodegeneration in spinocerebellar ataxia type 2. J Neurosci 2009, 29: 9148-9162. 10.1523/JNEUROSCI.0660-09.2009Google Scholar
- Serra H, Byam C, Lande J, Tousey S, Zoghbi H, Orr H: Gene profiling links SCA1 pathophysiology to glutamate signaling in Purkinje cells of transgenic mice. Hum Mol Genet 2004, 13: 2535-2543. 10.1093/hmg/ddh268Google Scholar
- Iwaki A, Kawano Y, Miura S, Shibata H, Matsuse D, Li W, Furuya H, Ohyagi Y, Taniwaki T, Kira J, Fukumaki Y: Heterozygous deletion of ITPR1, but not SUMF1, in spinocerebellar ataxia type 16. J Med Genet 2008, 45: 32-35.Google Scholar
- Desaiah D, Vig P, Subramony S, Currier R: Inositol 1,4,5-trisphosphate receptors and protein kinase C in olivopontocerebellar atrophy. Brain Res 1991, 552: 36-40. 10.1016/0006-8993(91)90656-GGoogle Scholar
- Zecevic N, Milosevic A, Ehrlich B: Calcium signaling molecules in human cerebellum at midgestation and in ataxia. Early Hum Dev 1999, 54: 103-116. 10.1016/S0378-3782(98)00090-5Google Scholar
- Di Gregorio E, Orsi L, Godani M, Vaula G, Jensen S, Salmon E, Ferrari G, Squadrone S, Abete M, Cagnoli C, et al: Two Italian Families with ITPR1 Gene Deletion Presenting a Broader Phenotype of SCA15. Cerebellum 2010, 9: 115-123. 10.1007/s12311-009-0154-0Google Scholar
- Hara K, Shiga A, Nozaki H, Mitsui J, Takahashi Y, Ishiguro H, Yomono H, Kurisaki H, Goto J, Ikeuchi T, et al: Total deletion and a missense mutation of ITPR1 in Japanese SCA15 families. Neurology 2008, 71: 547-551. 10.1212/01.wnl.0000311277.71046.a0Google Scholar
- Novak M, Davis M, Li A, Goold R, Tabrizi S, Sweeney M, Houlden H, Treacy C, Giunti P: PAW32 ITPR1 gene deletion causes spinocerebellar ataxia 15/16: a genetic, clinical and radiological description of a novel kindred. J Neurol Neurosurg Psychiatry 2010, 81: e32.Google Scholar
- Bezprozvanny I: Role of Inositol 1,4,5-Trishosphate Receptors in Pathogenesis of Huntington's Disease and Spinocerebellar Ataxias. Neurochem Res 2011, 36: 1186-1197. 10.1007/s11064-010-0393-yGoogle Scholar
- Schorge S, van de Leemput J, Singleton A, Houlden H, Hardy J: Human ataxias: a genetic dissection of inositol triphosphate receptor (ITPR1)-dependent signaling. Trends Neurosci 2010, 33: 211-219. 10.1016/j.tins.2010.02.005Google Scholar
- Kurnellas M, Lee A, Li H, Deng L, Ehrlich D, Elkabes S: Molecular alterations in the cerebellum of the plasma membrane calcium ATPase 2 (PMCA2)-null mouse indicate abnormalities in Purkinje neurons. Mol Cell Neurosci 2007, 34: 178-188. 10.1016/j.mcn.2006.10.010Google Scholar
- Tu H, Tang T, Wang Z, Bezprozvanny I: Association of type 1 inositol 1,4,5-trisphosphate receptor with AKAP9 (Yotiao) and protein kinase A. J Biol Chem 2004, 279: 19375-19382. 10.1074/jbc.M313476200Google Scholar
- Bezprozvanny I: The inositol 1,4,5-trisphosphate receptors. Cell Calcium 2005, 38: 261-272. 10.1016/j.ceca.2005.06.030Google Scholar
- Burright E, Clark H, Servadio A, Matilla T, Feddersen R, Yunis W, Duvick L, Zoghbi H, Orr H: SCA1 transgenic mice: a model for neurodegeneration caused by an expanded CAG trinucleotide repeat. Cell 1995, 82: 937-948. 10.1016/0092-8674(95)90273-2Google Scholar
- Orr H, Chung M, Banfi S, Kwiatkowski TJ, Servadio A, Beaudet A, McCall A, Duvick L, Ranum L, Zoghbi H: Expansion of an unstable trinucleotide CAG repeat in spinocerebellar ataxia type 1. Nat Genet 1993, 4: 221-226. 10.1038/ng0793-221Google Scholar
- Pulst S, Nechiporuk A, Nechiporuk T, Gispert S, Chen X, Lopes-Cendes I, Pearlman S, Starkman S, Orozco-Diaz G, Lunkes A, et al: Moderate expansion of a normally biallelic trinucleotide repeat in spinocerebellar ataxia type 2. Nat Genet 1996, 14: 269-276.Google Scholar
- Boy J, Schmidt T, Schumann U, Grasshoff U, Unser S, Holzmann C, Schmitt I, Karl T, Laccone F, Wolburg H, et al: A transgenic mouse model of spinocerebellar ataxia type 3 resembling late disease onset and gender-specific instability of CAG repeats. Neurobiol Dis 2010, 37: 284-293. 10.1016/j.nbd.2009.08.002Google Scholar
- Dürr A, Stevanin G, Cancel G, Duyckaerts C, Abbas N, Didierjean O, Chneiweiss H, Benomar A, Lyon-Caen O, Julien J, et al: Spinocerebellar ataxia 3 and Machado-Joseph disease: clinical, molecular, and neuropathological features. Ann Neurol 1996, 39: 490-499. 10.1002/ana.410390411Google Scholar
- Paulson H, Perez M, Trottier Y, Trojanowski J, Subramony S, Das S, Vig P, Mandel J, Fischbeck K, Pittman R: Intranuclear inclusions of expanded polyglutamine protein in spinocerebellar ataxia type 3. Neuron 1997, 19: 333-344. 10.1016/S0896-6273(00)80943-5Google Scholar
- Mantuano E, Veneziano L, Spadaro M, Giunti P, Guida S, Leggio M, Verriello L, Wood N, Jodice C, Frontali M: Clusters of non-truncating mutations of P/Q type Ca2+ channel subunit Ca(v)2.1 causing episodic ataxia 2. J Med Genet 2004, 41: e82. 10.1136/jmg.2003.015396Google Scholar
- Inoue T: Dynamics of calcium and its roles in the dendrite of the cerebellar Purkinje cell. Keio J Med 2003, 52: 244-249. 10.2302/kjm.52.244Google Scholar
- Gardner R: "SCA16" is really SCA15. J Med Genet 2008, 45: 192.Google Scholar
- Brown S-A, Loew LM: Toward A Computational Model Of IP3R1-associated Ataxia. Biophysical Journal 2009,96(3:S1):96a.Google Scholar
- Khodakhah K, Ogden D: Functional heterogeneity of calcium release by inositol trisphosphate in single Purkinje neurones, cultured cerebellar astrocytes, and peripheral tissues. Proc Natl Acad Sci U S A 1993, 90: 4976-4980. 10.1073/pnas.90.11.4976Google Scholar
- Fujiwara A, Hirose K, Yamazawa T, Iino M: Reduced IP3 sensitivity of IP3 receptor in Purkinje neurons. Neuroreport 2001, 12: 2647-2651. 10.1097/00001756-200108280-00012Google Scholar
- Hernjak N, Slepchenko B, Fernald K, Fink C, Fortin D, Moraru I, Watras J, Loew L: Modeling and analysis of calcium signaling events leading to long-term depression in cerebellar Purkinje cells. Biophys J 2005, 89: 3790-3806. 10.1529/biophysj.105.065771Google Scholar
- Brown SA, Moraru II, Schaff JC, Loew LM: Virtual NEURON: a strategy for merged biochemical and electrophysiological modeling. J Comput Neurosci 2011, 31: 385-400. 10.1007/s10827-011-0317-0Google Scholar
- Loew L, Schaff J: The Virtual Cell: a software environment for computational cell biology. Trends Biotechnol 2001, 19: 401-406. 10.1016/S0167-7799(01)01740-1Google Scholar
- Schaff J, Fink C, Slepchenko B, Carson J, Loew L: A general computational framework for modeling cellular structure and function. Biophys J 1997, 73: 1135-1146. 10.1016/S0006-3495(97)78146-3Google Scholar
- Schaff J, Slepchenko B, Loew L: Physiological modeling with virtual cell framework. Methods Enzymol 2000, 321: 1-23.Google Scholar
- Cowan AE, Moraru II, Schaff JC, Slepchenko BM, Loew LM: Spatial modeling of cell signaling networks. Methods Cell Biol 2012, 110: 195-221.Google Scholar
- Moraru I, Schaff J, Slepchenko B, Loew L: The virtual cell: an integrated modeling environment for experimental and computational cell biology. Ann N Y Acad Sci 2002, 971: 595-596. 10.1111/j.1749-6632.2002.tb04535.xGoogle Scholar
- Brown S, Morgan F, Watras J, Loew L: Analysis of phosphatidylinositol-4,5-bisphosphate signaling in cerebellar Purkinje spines. Biophys J 2008, 95: 1795-1812. 10.1529/biophysj.108.130195Google Scholar
- Tang T, Tu H, Wang Z, Bezprozvanny I: Modulation of type 1 inositol (1,4,5)-trisphosphate receptor function by protein kinase a and protein phosphatase 1alpha. J Neurosci 2003, 23: 403-415.Google Scholar
- Wang S, Denk W, Häusser M: Coincidence detection in single dendritic spines mediated by calcium release. Nat Neurosci 2000, 3: 1266-1273. 10.1038/81792Google Scholar
- Finch E, Augustine G: Local calcium signalling by inositol-1,4,5-trisphosphate in Purkinje cell dendrites. Nature 1998, 396: 753-756. 10.1038/25541Google Scholar
- Takechi H, Eilers J, Konnerth A: A new class of synaptic response involving calcium release in dendritic spines. Nature 1998, 396: 757-760. 10.1038/25547Google Scholar
- Destexhe A, Babloyantz A, Sejnowski T: Ionic mechanisms for intrinsic slow oscillations in thalamic relay neurons. Biophys J 1993, 65: 1538-1552. 10.1016/S0006-3495(93)81190-1Google Scholar
- De Schutter E, Bower J: An active membrane model of the cerebellar Purkinje cell. I. Simulation of current clamps in slice. J Neurophysiol 1994, 71: 375-400.Google Scholar
- Miyasho T, Takagi H, Suzuki H, Watanabe S, Inoue M, Kudo Y, Miyakawa H: Low-threshold potassium channels and a low-threshold calcium channel regulate Ca2+ spike firing in the dendrites of cerebellar Purkinje neurons: a modeling study. Brain Res 2001, 891: 106-115. 10.1016/S0006-8993(00)03206-6Google Scholar
- Gabbiani F, Midtgaard J, Knöpfel T: Synaptic integration in a model of cerebellar granule cells. J Neurophysiol 1994, 72: 999-1009.Google Scholar
- Cui J, Yang H, Lee U: Molecular mechanisms of BK channel activation. Cell Mol Life Sci 2009, 66: 852-875. 10.1007/s00018-008-8609-xGoogle Scholar
- Tang T, Tu H, Chan E, Maximov A, Wang Z, Wellington C, Hayden M, Bezprozvanny I: Huntingtin and huntingtin-associated protein 1 influence neuronal calcium signaling mediated by inositol-(1,4,5) triphosphate receptor type 1. Neuron 2003, 39: 227-239. 10.1016/S0896-6273(03)00366-0Google Scholar
- Tang T, Guo C, Wang H, Chen X, Bezprozvanny I: Neuroprotective effects of inositol 1,4,5-trisphosphate receptor C-terminal fragment in a Huntington's disease mouse model. J Neurosci 2009, 29: 1257-1266. 10.1523/JNEUROSCI.4411-08.2009Google Scholar
- Brown SA, Holmes RM, Loew LM: Spatial Organization and Diffusion in Neuronal Signaling. In Computational Systems Neurobiology. Edited by: Novère N. Dordrecht: Springer; 2012. Chapter 5, in pressGoogle Scholar
- Tu H, Wang Z, Bezprozvanny I: Modulation of mammalian inositol 1,4,5-trisphosphate receptor isoforms by calcium: a role of calcium sensor region. Biophys J 2005, 88: 1056-1069. 10.1529/biophysj.104.049601Google Scholar
- Sarkisov D, Wang S: Order-dependent coincidence detection in cerebellar Purkinje neurons at the inositol trisphosphate receptor. J Neurosci 2008, 28: 133-142. 10.1523/JNEUROSCI.1729-07.2008Google Scholar
- Kasumu A, Bezprozvanny I: Deranged Calcium Signaling in Purkinje Cells and Pathogenesis in Spinocerebellar Ataxia 2 (SCA2) and Other Ataxias. Cerebellum 2010. 10.1007/s12311-010-0182-9Google Scholar
- Miyata M, Finch E, Khiroug L, Hashimoto K, Hayasaka S, Oda S, Inouye M, Takagishi Y, Augustine G, Kano M: Local calcium release in dendritic spines required for long-term synaptic depression. Neuron 2000, 28: 233-244. 10.1016/S0896-6273(00)00099-4Google Scholar
- Wagner W, Hammer Jr: Myosin V and the endoplasmic reticulum: the connection grows. J Cell Biol 2003, 163: 1193-1196. 10.1083/jcb.200311077Google Scholar
- Vecellio M, Schwaller B, Meyer M, Hunziker W, Celio M: Alterations in Purkinje cell spines of calbindin D-28 k and parvalbumin knock-out mice. Eur J Neurosci 2000, 12: 945-954. 10.1046/j.1460-9568.2000.00986.xGoogle Scholar
- Vig PJ, Wei J, Shao Q, Lopez ME, Halperin R, Gerber J: Suppression of Calbindin-D28k Expression Exacerbates SCA1 Phenotype in a Disease Mouse Model. Cerebellum 2011. 10.1007/s12311-011-0323-9Google Scholar
- Schmidt H, Stiefel K, Racay P, Schwaller B, Eilers J: Mutational analysis of dendritic Ca2+ kinetics in rodent Purkinje cells: role of parvalbumin and calbindin D28k. J Physiol 2003, 551: 13-32. 10.1113/jphysiol.2002.035824Google Scholar
- Vig P, Subramony S, Burright E, Fratkin J, McDaniel D, Desaiah D, Qin Z: Reduced immunoreactivity to calcium-binding proteins in Purkinje cells precedes onset of ataxia in spinocerebellar ataxia-1 transgenic mice. Neurology 1998, 50: 106-113. 10.1212/WNL.50.1.106Google Scholar
- Vig P, Subramony S, McDaniel D: Calcium homeostasis and spinocerebellar ataxia-1 (SCA-1). Brain Res Bull 2001, 56: 221-225. 10.1016/S0361-9230(01)00595-0Google Scholar
- Sausbier M, Hu H, Arntz C, Feil S, Kamm S, Adelsberger H, Sausbier U, Sailer C, Feil R, Hofmann F, et al: Cerebellar ataxia and Purkinje cell dysfunction caused by Ca2 + −activated K + channel deficiency. Proc Natl Acad Sci U S A 2004, 101: 9474-9478. 10.1073/pnas.0401702101Google Scholar
- Cheron G, Sausbier M, Sausbier U, Neuhuber W, Ruth P, Dan B, Servais L: BK channels control cerebellar Purkinje and Golgi cell rhythmicity in vivo. PLoS One 2009, 4: e7991. 10.1371/journal.pone.0007991Google Scholar
- Kent J, Meredith A: BK channels regulate spontaneous action potential rhythmicity in the suprachiasmatic nucleus. PLoS One 2008, 3: e3884. 10.1371/journal.pone.0003884Google Scholar
- Womack M, Khodakhah K: Dendritic control of spontaneous bursting in cerebellar Purkinje cells. J Neurosci 2004, 24: 3511-3521. 10.1523/JNEUROSCI.0290-04.2004Google Scholar
- Walter J, Alviña K, Womack M, Chevez C, Khodakhah K: Decreases in the precision of Purkinje cell pacemaking cause cerebellar dysfunction and ataxia. Nat Neurosci 2006, 9: 389-397. 10.1038/nn1648Google Scholar
- Oliver D, Taberner A, Thurm H, Sausbier M, Arntz C, Ruth P, Fakler B, Liberman M: The role of BKCa channels in electrical signal encoding in the mammalian auditory periphery. J Neurosci 2006, 26: 6181-6189. 10.1523/JNEUROSCI.1047-06.2006Google Scholar
- Weaver A, Olsen M, McFerrin M, Sontheimer H: BK channels are linked to inositol 1,4,5-triphosphate receptors via lipid rafts: a novel mechanism for coupling [Ca(2+)](i) to ion channel activation. J Biol Chem 2007, 282: 31558-31568. 10.1074/jbc.M702866200Google Scholar
- Dai S, Hall D, Hell J: Supramolecular assemblies and localized regulation of voltage-gated ion channels. Physiol Rev 2009, 89: 411-452. 10.1152/physrev.00029.2007Google Scholar
- Sakagami Y, Yamamoto K, Sugiura S, Inokuchi K, Hayashi T, Kato N: Essential roles of Homer-1a in homeostatic regulation of pyramidal cell excitability: a possible link to clinical benefits of electroconvulsive shock. Eur J Neurosci 2005, 21: 3229-3239. 10.1111/j.1460-9568.2005.04165.xGoogle Scholar
- Vaithianathan T, Bukiya A, Liu J, Liu P, Asuncion-Chin M, Fan Z, Dopico A: Direct regulation of BK channels by phosphatidylinositol 4,5-bisphosphate as a novel signaling pathway. J Gen Physiol 2008, 132: 13-28. 10.1085/jgp.200709913Google Scholar
- Bezprozvanny I, Watras J, Ehrlich B: Bell-shaped calcium-response curves of Ins(1,4,5)P3- and calcium-gated channels from endoplasmic reticulum of cerebellum. Nature 1991, 351: 751-754. 10.1038/351751a0Google Scholar
- Guida S, Trettel F, Pagnutti S, Mantuano E, Tottene A, Veneziano L, Fellin T, Spadaro M, Stauderman K, Williams M, et al: Complete loss of P/Q calcium channel activity caused by a CACNA1A missense mutation carried by patients with episodic ataxia type 2. Am J Hum Genet 2001, 68: 759-764. 10.1086/318804Google Scholar
- Zhuchenko O, Bailey J, Bonnen P, Ashizawa T, Stockton D, Amos C, Dobyns W, Subramony S, Zoghbi H, Lee C: Autosomal dominant cerebellar ataxia (SCA6) associated with small polyglutamine expansions in the alpha 1A-voltage-dependent calcium channel. Nat Genet 1997, 15: 62-69. 10.1038/ng0197-62Google Scholar
- Yue Q, Jen J, Nelson S, Baloh R: Progressive ataxia due to a missense mutation in a calcium-channel gene. Am J Hum Genet 1997, 61: 1078-1087. 10.1086/301613Google Scholar
- Imbrici P, Cusimano A, D'Adamo M, De Curtis A, Pessia M: Functional characterization of an episodic ataxia type-1 mutation occurring in the S1 segment of hKv1.1 channels. Pflugers Arch 2003, 446: 373-379.Google Scholar
- Adachi N, Arima K, Asada T, Kato M, Minami N, Goto Yi, Onuma T, Ikeuchi T, Tsuji S, Hayashi M, Fukutani Y: Dentatorubral-pallidoluysian atrophy (DRPLA) presenting with psychosis. J Neuropsychiatry Clin Neurosci 2001, 13: 258-260. 10.1176/appi.neuropsych.13.2.258Google Scholar
- Martins S, Matamá T, Guimarães L, Vale J, Guimarães J, Ramos L, Coutinho P, Sequeiros J, Silveira I: Portuguese families with dentatorubropallidoluysian atrophy (DRPLA) share a common haplotype of Asian origin. Eur J Hum Genet 2003, 11: 808-811. 10.1038/sj.ejhg.5201054Google Scholar
- Yabe I, Sasaki H, Kikuchi S, Nonaka M, Moriwaka F, Tashiro K: Late onset ataxia phenotype in dentatorubro-pallidoluysian atrophy (DRPLA). J Neurol 2002, 249: 432-436. 10.1007/s004150200034Google Scholar
- Naito H, Oyanagi S: Familial myoclonus epilepsy and choreoathetosis: hereditary dentatorubral-pallidoluysian atrophy. Neurology 1982, 32: 798-807. 10.1212/WNL.32.8.798Google Scholar
- Bezprozvanny I: Calcium signaling and neurodegenerative diseases. Trends Mol Med 2009, 15: 89-100. 10.1016/j.molmed.2009.01.001Google Scholar
- Bezprozvanny I, Hayden M: Deranged neuronal calcium signaling and Huntington disease. Biochem Biophys Res Commun 2004, 322: 1310-1317. 10.1016/j.bbrc.2004.08.035Google Scholar
- Zhang H, Li Q, Graham R, Slow E, Hayden M, Bezprozvanny I: Full length mutant huntingtin is required for altered Ca2+ signaling and apoptosis of striatal neurons in the YAC mouse model of Huntington's disease. Neurobiol Dis 2008, 31: 80-88. 10.1016/j.nbd.2008.03.010Google Scholar
- Bezprozvanny I: Inositol 1,4,5-tripshosphate receptor, calcium signalling and Huntington's disease. Subcell Biochem 2007, 45: 323-335. 10.1007/978-1-4020-6191-2_11Google Scholar
- Stevanin G, Dürr A, Brice A: Clinical and molecular advances in autosomal dominant cerebellar ataxias: from genotype to phenotype and physiopathology. Eur J Hum Genet 2000, 8: 4-18. 10.1038/sj.ejhg.5200403Google Scholar
- Schmidt T, Lindenberg K, Krebs A, Schöls L, Laccone F, Herms J, Rechsteiner M, Riess O, Landwehrmeyer G: Protein surveillance machinery in brains with spinocerebellar ataxia type 3: redistribution and differential recruitment of 26 S proteasome subunits and chaperones to neuronal intranuclear inclusions. Ann Neurol 2002, 51: 302-310. 10.1002/ana.10101Google Scholar
- Klement I, Skinner P, Kaytor M, Yi H, Hersch S, Clark H, Zoghbi H, Orr H: Ataxin-1 nuclear localization and aggregation: role in polyglutamine-induced disease in SCA1 transgenic mice. Cell 1998, 95: 41-53. 10.1016/S0092-8674(00)81781-XGoogle Scholar
- Makarewicz D, Ziemińska E, Łazarewicz J: Dantrolene inhibits NMDA-induced 45Ca uptake in cultured cerebellar granule neurons. Neurochem Int 2003, 43: 273-278. 10.1016/S0197-0186(03)00012-3Google Scholar
- Zhao F, Li P, Chen S, Louis C, Fruen B: Dantrolene inhibition of ryanodine receptor Ca2+ release channels. Molecular mechanism and isoform selectivity. J Biol Chem 2001, 276: 13810-13816.Google Scholar
- Fruen B, Mickelson J, Louis C: Dantrolene inhibition of sarcoplasmic reticulum Ca2+ release by direct and specific action at skeletal muscle ryanodine receptors. J Biol Chem 1997, 272: 26965-26971. 10.1074/jbc.272.43.26965Google Scholar
- Krause T, Gerbershagen M, Fiege M, Weisshorn R, Wappler F: Dantrolene–a review of its pharmacology, therapeutic use and new developments. Anaesthesia 2004, 59: 364-373. 10.1111/j.1365-2044.2004.03658.xGoogle Scholar
- Gerbershagen M, Fiege M, Krause T, Agarwal K, Wappler F: Dantrolene. Pharmacological and therapeutic aspects. Anaesthesist 2003, 52: 238-245. 10.1007/s00101-003-0461-7Google Scholar
- Kobayashi S, Yano M, Suetomi T, Ono M, Tateishi H, Mochizuki M, Xu X, Uchinoumi H, Okuda S, Yamamoto T, et al: Dantrolene, a therapeutic agent for malignant hyperthermia, markedly improves the function of failing cardiomyocytes by stabilizing interdomain interactions within the ryanodine receptor. J Am Coll Cardiol 2009, 53: 1993-2005. 10.1016/j.jacc.2009.01.065Google Scholar
- Sharp AH, McPherson PS, Dawson TM, Aoki C, Campbell KP, Snyder SH: Differential immunohistochemical localization of inositol 1,4,5-trisphosphate- and ryanodine-sensitive Ca2+ release channels in rat brain. J Neurosci 1993, 13: 3051-3063.Google Scholar
- Walton PD, Airey JA, Sutko JL, Beck CF, Mignery GA, Südhof TC, Deerinck TJ, Ellisman MH: Ryanodine and inositol trisphosphate receptors coexist in avian cerebellar Purkinje neurons. J Cell Biol 1991, 113: 1145-1157. 10.1083/jcb.113.5.1145Google Scholar
- Khodakhah K, Armstrong CM: Inositol trisphosphate and ryanodine receptors share a common functional Ca2+ pool in cerebellar Purkinje neurons. Biophys J 1997, 73: 3349-3357. 10.1016/S0006-3495(97)78359-0Google Scholar
- Ln W, Li W, Yule D: Phosphorylation of type-1 inositol 1,4,5-trisphosphate receptors by cyclic nucleotide-dependent protein kinases: a mutational analysis of the functionally important sites in the S2+ and S2- splice variants. J Biol Chem 2003, 278: 45811-45817. 10.1074/jbc.M306270200Google Scholar
- Imbrici P, D'Adamo M, Cusimano A, Pessia M: Episodic ataxia type 1 mutation F184C alters Zn2 + −induced modulation of the human K + channel Kv1.4-Kv1.1/Kvbeta1.1. Am J Physiol Cell Physiol 2007, 292: C778-787.Google Scholar
- Imbrici P, D'Adamo M, Kullmann D, Pessia M: Episodic ataxia type 1 mutations in the KCNA1 gene impair the fast inactivation properties of the human potassium channels Kv1.4–1.1/Kvbeta1.1 and Kv1.4–1.1/Kvbeta1.2. Eur J Neurosci 2006, 24: 3073-3083. 10.1111/j.1460-9568.2006.05186.xGoogle Scholar
- Tonelli A, D'Angelo M, Salati R, Villa L, Germinasi C, Frattini T, Meola G, Turconi A, Bresolin N, Bassi M: Early onset, non fluctuating spinocerebellar ataxia and a novel missense mutation in CACNA1A gene. J Neurol Sci 2006, 241: 13-17. 10.1016/j.jns.2005.10.007Google Scholar