- Research
- Open access
- Published:
Reduction of anion reversal potential subverts the inhibitory control of firing rate in spinal lamina I neurons: towards a biophysical basis for neuropathic pain
Molecular Pain volume 2, Article number: 32 (2006)
Abstract
Background
Reduction of the transmembrane chloride gradient in spinal lamina I neurons contributes to the cellular hyperexcitability producing allodynia and hyperalgesia after peripheral nerve injury. The resultant decrease in anion reversal potential (i.e. shift in E anion to less negative potentials) reduces glycine/GABAA receptor-mediated hyperpolarization, but the large increase in membrane conductance caused by inhibitory input can nonetheless shunt concurrent excitatory input. Without knowing the relative contribution of hyperpolarization and shunting to inhibition's modulation of firing rate, it is difficult to predict how much net disinhibition results from reduction of E anion. We therefore used a biophysically accurate lamina I neuron model to investigate quantitatively how changes in E anion affect firing rate modulation.
Results
Simulations reveal that even a small reduction of E anion compromises inhibitory control of firing rate because reduction of E anion not only decreases glycine/GABAA receptor-mediated hyperpolarization, but can also indirectly compromise the capacity of shunting to reduce spiking. The latter effect occurs because shunting-mediated modulation of firing rate depends on a competition between two biophysical phenomena: shunting reduces depolarization, which translates into reduced spiking, but shunting also shortens the membrane time constant, which translates into faster membrane charging and increased spiking; the latter effect predominates when average depolarization is suprathreshold. Disinhibition therefore occurs as both hyperpolarization- and shunting-mediated modulation of firing rate are subverted by reduction of E anion. Small reductions may be compensated for by increased glycine/GABAA receptor-mediated input, but the system decompensates (i.e. compensation fails) as reduction of E anion exceeds a critical value. Hyperexcitability necessarily develops once disinhibition becomes incompensable. Furthermore, compensation by increased glycine/GABAA receptor-mediated input introduces instability into the system, rendering it increasingly prone to abrupt decompensation and even paradoxical excitation.
Conclusion
Reduction of E anion dramatically compromises the inhibitory control of firing rate and, if compensation fails, is likely to contribute to the allodynia and hyperalgesia associated with neuropathic pain. These data help explain the relative intractability of neuropathic pain and illustrate how it is important to choose therapies not only based on disease mechanism, but based on quantitative understanding of that mechanism.
Background
Neuropathic pain can arise from a multitude of pathophysiologic mechanisms or combinations thereof [1–5]. Reduced inhibition, or disinhibition, of spinal neurons constitutes an important class of these mechanisms [6–8]. Indeed, many features of neuropathic pain syndromes can be reproduced by pharmacologically blocking inhibition in the spinal cord [9–17] or through genetic changes that reduce inhibition [18]. Conversely, increasing inhibition can, in some conditions, reduce neuropathic pain [17, 19–22].
Reproducing features of neuropathic pain by reducing inhibition and relieving neuropathic pain by increasing inhibition are important observations but constitute only circumstantial evidence that disinhibition is involved in the pathogenesis of neuropathic pain. More direct evidence comes from studies showing that neuropathy can occlude the effects of pharmacologically reducing inhibition [19] and that inhibition is indeed reduced in animal models of neuropathic pain [23–32]. There is controversy whether reduction of inhibitory transmitters and/or their receptors occurs following neuropathy [33–36], but disinhibition can arise from a broad array of mechanisms.
The recent study by Coull et al. [31] suggests an alternative mechanism to explain disinhibition: reduced expression of the potassium-chloride cotransporter (KCC2) causes reduction of the chloride gradient across the neuronal membrane, which in turn leads to reduction of the anion reversal potential (i.e. E anion shifts to a less negative membrane potential) [see also [37]]. The change in driving force means that both glycine and GABAA receptor-mediated inputs produce less hyperpolarization and could even paradoxically depolarize the neuron. But even if those inputs become depolarizing on their own, their large conductances mean they may still reduce the depolarization caused by concurrent excitatory input [e.g. [38, 39]] – a phenomenon known as shunting.
Without knowing the relative contribution of shunting and hyperpolarization to firing rate modulation, it is difficult to predict how reduction of E anion will impact inhibitory control of firing rate, especially given the nonlinearity inherent in spike generation [40, 41]. Recent work has revealed a good correlation between E anion and pain threshold [42], but it remains unproven whether reductions in E anion could compromise inhibition sufficiently to produce the cellular hyperexcitability that may in turn cause the perceptual/behavioral features of neuropathic pain including allodynia (perception of pain in response to normally innocuous stimulation) and hyperalgesia (exaggerated pain perception in response to noxious stimulation).
Allodynia and hyperalgesia are typically thought to arise from hyperexcitability at the cellular and network levels [43]. It is implicit in our analysis that responses in lamina I neurons correlate with pain perception, and that hyperexcitability amongst those neurons could therefore give rise to allodynia and hyperalgesia. Although there is no doubt that lamina I neurons convey information to supraspinal targets [44, 45], it has been widely assumed that wide dynamic range cells in lamina V are solely capable of encoding stimulus intensity because lamina I cells are predominantly nociceptive specific (which is mistakenly taken to imply that they lack the capacity to modulate their response magnitude depending on stimulus intensity) and that lamina I cells do not receive low threshold input [for review see [46]]. On the contrary, multiple studies have demonstrated the capacity of lamina I projection neurons to encode stimulus intensity [47–51]; moreover, more recent work shows that the A and C fiber nociceptors that innervate lamina I can encode stimulus intensity [52]. It has also been shown that lamina I neurons receive low threshold inputs via polysynaptic pathways, but transmission through those pathways is normally suppressed by inhibition [53]. The cumulated evidence therefore indicates that lamina I projection neurons can encode nociceptive information relevant for pain perception. It logically follows that hyperexcitability of lamina I neurons may contribute to allodynia and hyperalgesia.
This study was undertaken to investigate quantitatively whether reductions in E anion could compromise inhibition sufficiently to produce cellular hyperexcitability despite compensatory changes that may develop (e.g. enhanced GABAA components to synaptic events [31]). To this end, we developed a biophysically accurate neuron model for quantitative testing under a variety of conditions. Results demonstrate that even a small reduction of E anion can cause disinhibition. But whereas compensatory changes may prevent disinhibition resulting from a small reduction of E anion, reduction of E anion at magnitudes reported by Coull et al. [31] almost certainly causes incompensable disinhibition. Disinhibition becomes incompensable when compensatory mechanisms fail; the failure of compensation means that the system decompensates and necessarily becomes hyperexcitable. Demonstration that decompensation can occur abruptly, especially in a highly compensated system, may help explain the sometimes unpredictable efficacy of therapeutic interventions and has implications for how to optimize treatment of neuropathic pain.
Results
Impact of E anionon the inhibitory control of firing rate
A multicompartment model based on available data on lamina I neurons was developed as outlined in the Methods (Fig. 1A). The model neuron was bombarded by synaptic input simulated in a biophysically realistic manner: discrete synaptic events generated inward or outward currents by opening channels in the cell membrane (rather than by simply injecting current) so that synaptic inputs influence the total membrane conductance of the postsynaptic neuron. Simulating synaptic input in this way is crucial for uncovering the role of shunting in firing rate modulation. Three varieties of inhibitory input were tested (Fig. 1B): proportional (to excitation), constant, and feedback. We also tested model neurons with different intrinsic properties (Fig. 1C): basic model containing only fast Na+ and delayed rectifier K+ channels (i.e. Hodgkin-Huxley or HH channels), tonic-spiking, and single-spiking. Tonic- and single-spiking neurons represent the opposite extremes of physiological cell types in lamina I [54]; their responses to somatic injection of constant current are illustrated in Fig. 1C. Using these models, a systematic investigation of the effects of reducing E anion on firing rate modulation was carried out.
Starting with the basic model neuron, a series of simulations was performed in which the frequency of excitatory synaptic input (f exc), the frequency of inhibitory synaptic input (f inh), and E anion were systematically varied while measuring the frequency of output spiking (f out). We first posited that both f exc and f inh increase with increasingly strong stimulation, and that the increase is proportional though not necessarily equal between excitation and inhibition; the ratio of inhibition to excitation (f inh/f exc) is reported as α and was tested at multiple values. Figure 2A shows that reduction of E anion compromises inhibition's capacity to reduce firing rate. The degree to which firing rate modulation is compromised is most readily understood by comparing the f out-f exc curves corresponding to different values of E anion (colored curves, Fig. 2A) against the f out0-f exc curve with no inhibition (i.e. α = 0; black curve, Fig. 2A). A reduction of E anion to -55 mV completely incapacitated inhibition, as evidenced by the corresponding f out-f exc curve (light green) lying very close to the no inhibition curve. Reduction of E anion to -50 or -45 mV caused paradoxical excitation, as evidenced by the corresponding f out-f exc curves (yellow and orange) lying above the no inhibition curve. Even modest reduction of E anion to -65 or -60 mV caused some disinhibition, as evidenced by the corresponding f out-f exc curves (blue and dark green) lying above the curve for E anion = -70 mV (purple) but below the no inhibition curve. The transition from modest disinhibition to complete disinhibition to paradoxical excitation as E anion was shifted from -70 mV to -45 mV occurred regardless of the ratio of inhibition to excitation (compare panels showing different values of α in Fig. 2A). The divergence the f out-f exc curves was, however, exaggerated as α was increased.
Adding feedback inhibition to proportional inhibition with α = 0.5 had an effect (Fig. 2B) similar to increasing the strength of proportional inhibition (i.e. increasing α in Fig. 2A). In other words, feedback inhibition exaggerated the effects of reducing E anion, which was manifested on the graphs as increased divergence of the f out-f exc curves. Interestingly, the divergence was more exaggerated for large reductions in E anion (i.e. -45 and -50 mV) compared with small reductions. This is explained by the fact that, when the glycine/GABAA receptor-mediated input becomes paradoxically excitatory, a positive feedback loop is born, whereas the feedback loop is negative under normal conditions. Positive feedback (i.e. feedback excitation) translates into extreme hyperexcitability.
We then tested the effect of constant inhibition. In contrast to the divergent f out-f exc curves seen for proportional inhibition and for feedback inhibition (which is also "proportional" insofar as the feedback neuron's activation is proportional to the output neuron's activity), constant inhibition caused a more parallel shift in the f out-f exc curves (Fig. 2C). The vertical separation of those curves was enhanced by increasing f inh.
Using proportional inhibition, we repeated simulations in a tonic-spiking model neuron (Fig. 2D) and in a single-spiking model neuron (Fig. 2E). While the slopes of the f out-f exc curves for these two models were different from those of the basic model (Fig. 2A), the effects of reducing E anion were very similar: in the tonic-spiking model, reduction of E anion to -55 mV caused complete disinhibition, the same as in the basic model; the single-spiking model was more resistant, requiring that E anion be reduced to around -50 mV before disinhibition becomes complete. These data demonstrate that despite the extreme differences in the intrinsic membrane properties of these cell types, firing rate modulation is altered in qualitatively the same way as E anion becomes reduced.
Data in Figure 2 thus indicate that the degree of reduction in E anion correlates with the degree to which inhibitory control of spiking is compromised: small reductions (to -65 or -60 mV) cause modest disinhibition, intermediate reductions (to -55 mV) can cause complete disinhibition (i.e. equivalent to completely removing inhibition), and large reductions (to -50 or -45 mV) cause paradoxical excitation. This correlation remains quantitatively true regardless of the circuit connectivity (feedback vs. proportional inhibition), strength of inhibition (magnitude of α or of f inh), or intrinsic neuronal properties (tonic- vs. single-spiking). In the more detailed investigation that follows, we focus on proportional inhibition in the basic model but, based on the demonstrations in Figure 2, the results can be reasonably extrapolated to other conditions.
Relative importance of shunting and hyperpolarization for firing rate modulation
The results above demonstrate that reduction of E anion significantly compromises inhibitory control of firing rate. But although a reduction in E anion should logically reduce glycine/GABAA receptor-mediated hyperpolarization, a change in E anion should not reduce glycine/GABAA receptor-mediated shunting. The results, therefore, suggest that shunting plays a relatively minor role in the modulation of firing rate, compared with the dominant role of hyperpolarization. This result was unexpected given previous reports on the importance of shunting [e.g. [55, 56]] and therefore required further investigation.
The lack of effect of shunting on firing rate contrasts its known effect on depolarization. To investigate in isolation how shunting modulates depolarization, we modified the basic model by removing HH channels, thereby preventing the model neuron from spiking and allowing us to quantify the underlying depolarization. Figure 3 shows that whereas the depol-f exc curves for E anion = -45 and -50 mV lay above the curve for no inhibition at low values of f exc, they bent downwards at higher values of f exc so that they eventually lay below the curve for no inhibition (arrow in Fig. 3). This sublinearity in the depol-f exc curves (i.e. downward bend with increasing input) is precisely what one would expect from shunting, where increasingly more excitatory input is shunted through open glycine/GABAA channels as excitation increases. These results thereby demonstrate that glycine/GABAA receptor-mediated shunting remains intact despite reduction of E anion. As expected, the sublinearity attributable to shunting was absent when there was no inhibition (black depol-f exc curve in Fig. 3 is linear).
But whereas the effects of shunting are clearly evident in the depol-f exc curves of Figure 3, they are absent from the f out-f exc curves of Figure 2 (i.e. the f out-f exc curves with inhibition are only slightly sublinear). If excitatory input drives depolarization, and depolarization drives spiking, then one should look at the relationship between depolarization and firing rate in order to identify why shunting's capacity to modulate firing rate is lost. To do this, we plotted firing rate (f out) against depolarization (Fig. 4A) rather than against f exc (as had been down in Fig. 2). Specifically, we plotted f out generated by a certain value of f exc (based on simulations in the model neuron with HH channels) against depolarization generated by the same value of f exc (based on simulations without HH channels).
Shunting can paradoxically increase firing rate
Figure 4A reveals a rather counterintuitive observation that, for the same amount of depolarization, significantly faster spiking was generated in the presence of synaptic inhibition than in its absence. This was true regardless of the value of E anion and became more evident as depolarization increased (see below). The most likely explanation for this is modulation of the membrane time constant (Ï„membrane), since Ï„membrane becomes shorter as membrane conductance increases (Fig. 4B). To test this, we increased the passive leak conductance in the model neuron to reduce Ï„membrane by half (open circle on Fig 4B); like inhibitory synaptic input, increasing the passive leak conductance caused f out to increase (dashed curve on Fig. 4A) compared with the original model neuron (solid black curve). Why does shortening Ï„membrane lead to faster spiking? When the model neuron was shunted and therefore had a short Ï„membrane, its power spectrum exhibited higher power at all but the lowest frequencies (blue curve in Fig. 4C) compared with the power spectrum for the unshunted neuron with a longer Ï„membrane (black curve). The most important consequence of this reduced filtering of high frequencies is illustrated in the inset of Figure 4C: reduced filtering allows the membrane to recharge more quickly between spikes so that higher firing rates can be achieved when the neuron is shunted than when it is not shunted.
Therefore, shunting can reduce the depolarization caused by excitatory input (Fig. 3), which reduces firing rate (Fig. 4A), but it also shortens Ï„membrane (Fig. 4B), which can increase firing rate (Fig. 4C inset). Both effects occur regardless of the value of E anion, but their relative importance for firing rate modulation depends on other factors. It is important to recall here that the depol-f exc curves with inhibition diverge from the depol-f exc curve without inhibition only when depolarization becomes suprathreshold or, at least, nearly suprathreshold (Fig. 4A). This indicates that that shortening Ï„membrane is only consequential for firing rate modulation when depolarization is suprathreshold, whereas we know from Figure 3 that shunting reduces both sub- and suprathreshold depolarization.
We have previously demonstrated that, when average depolarization is suprathreshold, spikes are generated deterministically whereas, when average depolarization is subthreshold, spikes are elicited by noisy, suprathreshold voltage fluctuations and are therefore generated probabilistic ally [41]. The rate of probabilistic spiking is reduced by shunting because, by reducing depolarization, shunting increases the difference between average depolarization and voltage threshold, which in turn reduces the likelihood of voltage fluctuations crossing threshold. If shunting can reduce suprathreshold depolarization so that it becomes subthreshold, deterministic spiking will become probabilistic and its rate will be reduced by shunting according to the above mechanism. If, on the other hand, excitation is sufficiently strong to cause net depolarization that remains suprathreshold despite shunting, spiking will be deterministic. Under those conditions, rather than being limited by the probability of threshold crossing, the interspike interval is limited by the rate of interspike membrane charging which, as demonstrated above, is sensitive to both depolarization and Ï„membrane. Shunting reduces depolarization but it also shortens Ï„membrane, where each effect has the opposite impact on the rate of membrane charging. Thus, in regard to firing rate modulation, the dual effects of increased membrane conductance counterbalance each other when spiking is deterministic, whereas the inhibitory effect via reduction of depolarization becomes dominant if spiking is probabilistic. This explanation reconciles the observations in Figures 2 and 3: shunting can reduce depolarization while at the same time having little effect on firing rate.
Since shunting becomes ineffective at modulating firing rate when average depolarization is suprathreshold, one can appreciate that by reducing the hyperpolarization caused by glycine/GABAA receptor-mediated input or, worse yet, causing that input to become depolarizing, a reduction of E anion can indirectly compromise shunting's capacity to reduce firing rate (Fig. 4D). Thus, both mechanisms through which inhibitory input normally modulates firing rate (i.e. hyperpolarization and shunting) are both susceptible (directly or indirectly) to changes in E anion. Furthermore, not only are both inhibitory mechanisms compromised by reduction of E anion, both mechanisms can contribute to paradoxical excitation if the reduction of E anion is sufficiently large.
Effects of constant vs. intermittent inhibition
We also tested whether the irregularity of inhibitory synaptic input could compromise control of spiking inasmuch as spikes can occur during inhibitory gaps (i.e. between inhibitory synaptic events). This irregularity may be expected to compromise modulation of firing rate by shunting more than it compromises modulation by hyperpolarization since inhibitory gaps are larger in the former case. The difference in size of inhibitory gaps is a direct consequence of the differential time course of inhibitory postsynaptic currents (IPSCs) and inhibitory postsynaptic potentials (IPSPs): IPSPs have a slower time course because of the low pass filtering caused by the membrane time constant, so that whereas slow IPSPs tend to overlap and produce sustained potentials, faster IPSCs (which directly parallel the time course of the synaptic conductance) overlap less and, therefore, shunting is likely to be intermittent (Fig. 5A). Quantitative testing confirmed that gaps between IPSCs close more slowly than gaps between IPSPs as f inh increases (data not shown). To test whether inhibitory gaps are functionally significant for firing rate modulation, we replaced the intermittent inhibition associated with irregular synaptic input with constant inhibition equal to the time-averaged synaptic input.
Figure 5B demonstrates that switching from intermittent to constant inhibition had virtually no effect on f out-f exc curves except at the most negative values of E anion. On initial examination, this suggests that inhibitory gaps are relatively insignificant for firing rate modulation. On closer examination however, the fact that irregular inhibition produces larger gaps in shunting than in hyperpolarization (see above and Fig. 5A) suggests that the irrelevance of switching from intermittent inhibition to constant inhibition may simply reflect the relative impotency of shunting to reduce spiking when spiking is generated in a deterministic manner (see above and Fig. 4). Following from this point, the observation that regularizing inhibition improved the effectiveness of inhibition for E anion values of -65 and -70 mV (stars in Fig. 5B) is significant, as those E anion values correspond to the same conditions under which depolarization remained subthreshold at high f exc (Fig. 3). If average depolarization remains subthreshold, then spiking is generated in a probabilistic manner and shunting should retain it capacity to reduce firing rate (see above). In short, Figure 5B shows that regularizing inhibition enhances shunting's capacity to reduce spiking when spiking is probabilistic, but not when spiking is deterministic, consistent with conclusions drawn from Figure 4 regarding the conditions under which shunting can or can not modulate spiking.
Disinhibition in the context of compensatory changes and modulation outside lamina I
Reduction of E anion does not necessarily occur in isolation. The possibility of compensatory changes complicates the conclusion that disinhibition necessarily follows from reduction of E anion. Two separate studies suggest that GABAA transmission may increase following neuropathy [31, 57]. We therefore investigated whether disinhibition develops if reduction of E anion is coupled with compensatory increases in GABA/glycine transmission.
Under normal conditions with E anion = -70 mV, even modest amounts of inhibition (e.g. α = 0.5) can reduce firing rate (Fig. 6A, dotted curve). The decrease in firing rate can be quantified by taking the ratio of the response with inhibition to the response without inhibition (f out/f out0). The f out/f out0 ratio is ~0.6 for α = 0.5, but larger values of α result in greater reduction of that ratio (Fig. 6A, dashed and solid curves). For purposes of illustration, we estimate inhibition conservatively (α = 0.5) and, by extension, consider f out/f out0 > 0.6 to represent disinhibition. If E anion were to decrease to -60 mV, inhibition with α = 0.5 would be incapable of reducing f out/f out0 to 0.6 (Fig. 6B, dotted curve), meaning disinhibition would occur under those conditions. But if α was increased to 1.2, the f out/f out0 ratio would be returned to 0.6 (Fig. 6B, blue curve). Assuming the network has the capacity to at least double α, this demonstrates that moderate reductions in E anion cause compensable disinhibition; in other words, the disinhibition that could potentially result from reduction of E anion can be compensated for by shifting the balance between excitatory and inhibitory input. However, if E anion were to decrease to -55 mV, even increasing α to 2 could not reduce f out/f out0 to 0.6 (Fig. 6C), demonstrating that disinhibition becomes incompensable if reduction of E anion becomes large. In fact, an increase in α may not only fail to reduce f out/f out0 sufficiently (as in Fig. 6C), but may even paradoxically increase f out/f out0 as in the case with E anion = -50 mV (Fig. 6D), resulting in paradoxical excitation. It is interesting to note that paradoxical excitation can occur for values of E anion below spike threshold, which is approximately -49 mV in the model neuron tested here.
Effects of E anion on the f out/f out0 ratio are summarized in Figure 6E. f out/f out0 > 0.6 represents disinhibition based on the conservative estimate of α = 0.5 prior to any compensation, while f out/f out0 > 1 represents paradoxical excitation. The relationship between α and the critical value of E anion beyond which net disinhibition necessarily develops (pink curve on Fig. 6E) shows a ceiling effect wherein increases in α, no matter how large, can not compensate for reductions in E anion beyond a certain level. For example, assuming α could increase as high as 4, disinhibition would still develop at E anion ≈ -54 mV (arrowhead in Fig. 6E). Stricter limitation on the increase in α would result in disinhibition occurring for even smaller reductions in E anion.
In a system capable of compensation, reduction of E anion does not produce disinhibition until the system decompensates; decompensation occurs when reductions in E anion outstrip compensatory changes. But although compensatory changes may prevent decompensation from occurring until a large reduction of E anion has developed, Figure 6F shows that in a system relying on strong compensation (e.g. α = 4 to maintain an f out/f out0 ratio of 0.6), small changes in E anion may cause large changes in the f out/f out0 ratio; in other words, decompensation occurs abruptly, which reflects instability within the system. This is evident from the steepness of the f out/f out0-E anion curve where it passes through f out/f out0 = 0.6 (long-dashed curve in Fig. 6F). A system relying on less compensation (e.g. α = 1) may decompensate at a lower value of E anion, but will do so gradually (dot-dashed curve in Fig. 6F), indicating that the system is more stable. Understanding the abruptness of decompensation may help guide therapeutic intervention insofar as it is preferable to reestablish a robust balance between excitation and inhibition rather than an unstable one (see below).
Pathophysiologic changes associated with neuropathic pain also occur outside lamina I. How do those changes influence activity in lamina I neurons? We have heretofore assumed a constant relationship between the strength of peripheral stimulation and the strength of synaptic input to lamina I. As a postsynaptic mechanism, disinhibition resulting from reduced E anion within lamina I neurons does not affect that relationship; on the other hand, presynaptic mechanisms including peripheral sensitization and presynaptic inhibition do (Fig. 7A): peripheral sensitization steepens that relationship by increasing the frequency of excitatory input elicited by a given stimulus, while presynaptic inhibition reduces the amplitude of individual synaptic events. The distinction between modulation of event frequency and event amplitude is irrelevant for this discussion, insofar as cumulative synaptic excitation (in Siemens per second) equals amplitude per input (in Siemens per event) multiplied by event frequency (in events per second); in short, peripheral sensitization increases cumulative synaptic excitation whereas presynaptic inhibition decreases it. Following the simple logic outlined in Figure 7B, the relationship between f out and synaptic excitation (which for all intents and purposes is equivalent to the f out-f exc relationship) can be combined with the relationship between synaptic excitation and strength of peripheral stimulation to determine the relationship between f out and strength of peripheral stimulation. Figure 7C illustrates how peripheral sensitization and presynaptic inhibition influence the relationship between f out and peripheral stimulation; note that the relationship between f out and synaptic excitation remains unchanged. If the f out/f out0 ratio is calculated using f out from the test condition and f out0 from the control condition, then the effects of peripheral sensitization and presynaptic inhibition can be visualized as leftward or rightward shifts, respectively, in the relationship between f out/f out0 and E anion (Fig. 7D). Thus, whereas peripheral sensitization will exacerbate the effects of postsynaptic disinhibition, presynaptic inhibition will mitigate the effects. Beyond that, reduced pre- or postsynaptic inhibition within polysynaptic pathways may uncover low threshold input to lamina I neurons [53], which would also alter the relationship between the strength of peripheral stimulation and synaptic input.
Implications for therapeutic interventions
It seems logical that if a reduction in E anion has compromised the strength of inhibition, interventions aimed at reestablishing inhibition's strength would be beneficial for correcting disinhibition. Ideally one would target the primary pathophysiologic cause (e.g. KCC2 expression or its direct effects on chloride extrusion from inside the cell, or upstream events including microglial activation and BDNF release [31, 42]) but with no clinically available drugs to do this, other targets must be considered. For example, the strength and kinetics of GABAA receptor-mediated input can be modulated by benzodiazepines. We therefore investigated the effects of doubling the strength (w) and τdecay of GABAergic inputs; parameters of glycinergic input were left unchanged. Figure 8A shows that augmenting individual GABAergic inputs mimics the effects of increasing α (compare with Fig. 6F): it delayed decompensation from occurring until higher values of E anion, but it did so by steepening the curve relating f out/f out0 and E anion. As explained above, this steepening means that although the intervention may rebalance the system at a normal f out/f out0 ratio, that balance is liable to be disturbed by even small changes in E anion; in other words, the system becomes increasingly unstable. Increasing GABAergic transmission also risks exacerbating paradoxical excitation if reduction of E anion is particularly large. Furthermore, although endogenous compensation may be expected to occur specifically in neurons affected by changes in E anion, a drug would not show the same specificity, suggesting the f out/f out0 ratio would be inadvertently reduced below 0.6 in normal cells (e.g. with E anion = -70 mV) exposed to the drug at the same time that the f out/f out0 ratio is potentially normalized in affected cells (e.g. with E anion = -55 mV). Therefore, although a disinhibitory mechanism intuitively suggests that inhibition should be augmented in order to offset the disinhibition, results here suggest that augmenting GABAA receptor-mediated input (or, similarly, glycine receptor-mediated input) may be unwise when disinhibition occurs through reduction of E anion. In contrast, disinhibition occurring through a different mechanism (e.g. reduced expression or function of GABAA or glycine receptors) could be successfully corrected by augmenting inhibition (see Discussion).
Although increasing inhibition ultimately fails to prevent decompensation if the reduction in E anion grows too large, and at best reestablishes a normal f out/f out0 ratio that is easily disrupted, interventions that do not target inhibition are likely to be more robust. NMDA receptors are a common pharmacological target whose blockade can reduce neuropathic pain [58]. Figure 8B illustrates that blocking NMDA receptors did not alter the slope of the curve (and thus the stability of the system) and instead caused a uniform reduction in the f out/f out0 ratio, which is similar to the effect of presynaptic changes described in Figure 7. Furthermore, reducing excitatory input did not risk exacerbating paradoxical excitation. The magnitude of effects of NMDA antagonism depends on the contribution of NMDA receptor-mediated excitation relative to AMPA receptor-mediated excitation, which is relatively small for conditions tested here (see Methods) but may increase under neuropathic conditions [59]. In any case, reducing AMPA receptor-mediated excitation had a similar effect (see below). Reducing sodium channel density, which functionally mimics the postsynaptic effects of local anesthetics and many anti-epileptic medications [60], also resulted in modulation very similar to that described above for NMDA antagonism (data not shown). Insertion of a potassium conductance such as might be activated by opioids [61] also had a similar effect (data not shown). All in all, reducing excitatory input and/or reducing intrinsic neuronal excitability reduces the f out/f out0 ratio uniformly across a broad range of E anion. Combining any of these manipulations with augmentation of inhibition had purely additive effects (Fig. 8C) rather than acting synergistically.
Based on the observation that it may be counterproductive to try to replace inhibition when disinhibition occurs through reduction of E anion, we explored whether it would be preferable to block inhibitory input altogether and balance the resulting disinhibition with reduction of excitatory input (Fig. 8D). The rationale is that although, on its own, blocking glycine/GABAA receptor-mediated input would increase the f out/f out0 ratio to 1, the appropriate reduction in excitatory input and/or intrinsic excitability could return the ratio to around 0.6 (Fig. 8D, right panel). The benefit of deliberately switching the disinhibitory mechanism (i.e. from reduction of E anion to blockade of glycine/GABAA receptors) is that blocking glycine/GABAA receptor-mediated input would uniformly adjust the f out/f out0 ratio to 1, thereby trivializing the variability of E anion and α between affected and unaffected neurons. Moreover, blocking inhibition would prevent inhibitory input from causing paradoxical excitation in the event that reduction of E anion was particularly large. One significant requirement is that the distribution and kinetics of the drugs affecting inhibitory and excitatory transmission are similar in order to avoid spatial or temporal mismatches between the two effects. Even if this approach may not be feasible in practice, it illustrates that increasing glycine/GABAA receptor-mediated input may be counterproductive in certain neuropathic conditions and that alternative approaches, possibly involving non-intuitive drug combinations, deserve consideration.
Discussion
This modeling study demonstrates that reduction of E anion in spinal lamina I neurons can result in disinhibition and hyperexcitability. Compensatory changes may successfully prevent disinhibition and maintain a normal input-output relationship, but they ultimately fail (i.e. the system decompensates) if reduction of E anion exceeds a critical value. Although incompensable disinhibition requires a relatively large reduction in E anion, the magnitude of that reduction is physiologically plausible and need not be so large as to cause glycine/GABAA receptor-mediated inputs to become paradoxically excitatory. Compensatory mechanisms eventually fail because, although glycine/GABAA receptor-mediated inputs continue to shunt excitatory inputs despite no longer causing hyperpolarization, the increase in membrane conductance that underlies shunting also shortens the membrane time constant which, under certain circumstances, allows for faster spiking. Deciphering this complex interplay between biophysical mechanisms is, ultimately, important for being able to ascribe perceptual changes (e.g. allodynia and hyperalgesia) to pathophysiologic changes at the cellular and molecular level.
Relative importance of shunting and hyperpolarization for firing rate modulation
Glycine and GABAA receptor-mediated inputs are often thought to act via shunting rather than through hyperpolarization given their large conductance and relatively small driving force. Consequently, even depolarizing glycine/GABAA receptor-mediated input can have a net inhibitory effect on spike generation by shunting concurrent excitatory input [38, 39]. In the current study, we found that shunting was far less significant than hyperpolarization when it comes to modulation of repetitive spiking. This results from increased membrane conductance having opposing effects: it reduces depolarization by shunting excitatory input (thereby decreasing firing rate) but simultaneously shortens the membrane time constant (thereby increasing firing rate). In terms of firing rate modulation, the first effect predominates when average depolarization is subthreshold and spiking is driven by suprathreshold voltage fluctuations (i.e. probabilistic spiking) but the two effects offset each other when average depolarization is suprathreshold (i.e. deterministic spiking). Shunting therefore becomes ineffective at reducing firing rate when depolarization is suprathreshold, which is especially likely to occur when reduction of E anion renders glycine/GABAA receptor-mediated input less hyperpolarizing or, worse yet, depolarizing. Kuhn et al. [62] have previously described how firing rate modulation is complicated by the dual effects of membrane conductance, noting that average depolarization was reduced by shunting while voltage fluctuations became larger because of the reduced filtering associated with a shortened membrane time constant. Their result applies to probabilistic spiking whereas the effect that we have described applies to deterministic spiking and, in that sense, is distinct.
A recent modeling study that investigated the effects of changes in E anion on firing rate modulation in a neocortical pyramidal neuron [63] concluded that, because of shunting, GABA remained inhibitory despite changes in E anion. That would appear to contradict our conclusion that shunting fails in the face of large, but physiologically plausible (as reported in Coull et al. [31]), reductions in E anion. However, Morita et al. tested only two values of E anion. According to our evaluation of the phase-planes shown in Figure 4C of their paper, a Hopf bifurcation would not have prevented repetitive spiking if a slightly lower value of E anion had been tested. There is therefore no discrepancy between our studies but, simply, a difference in the range of E anion tested.
Compensable vs. incompensable disinhibition
Given this new understanding of the biophysical mechanisms, we can predict whether the disinhibition caused by reduction of E anion can account for the hyperexcitability that is presumably responsible for the allodynia and hyperalgesia associated with neuropathic pain. We estimate that decompensation would start occurring at E anion ≈ -58 mV (see Fig. 6) assuming that, as compensation, the ratio of inhibitory to excitatory input quadruples relative to the ratio under normal conditions (i.e. α increases from 0.5 to 2). Weaker compensation (α increases to only 1) would result in decompensation starting at E anion ≈ -61 mV, whereas stronger compensation (α increases as high as 4) would prevent decompensation until E anion ≈ -54 mV. All of those estimates conservatively assume α = 0.5 under normal conditions (see Results).
Variations in compensation may explain why, with acute manipulations including intrathecal application of BDNF and activated microglia, reduction in E anion to around -62 mV caused almost an equivalent decrease in withdrawal threshold as that associated with reduction of E anion to -49 mV following peripheral nerve injury [42]: compensation that could have developed in the latter case, may not have had time to develop in the former case; it is also possible, however, that BDNF has other effects [64] that are unaccounted for in this argument. Additionally, we have not taken into account activity-dependent reduction of the chloride gradient [65–68], meaning a much smaller long-term reduction of E anion (i.e. caused by reduced KCC2 expression) may cause incompensable disinhibition once dynamic, short-term reductions of E anion (i.e. caused by activity-dependent reduction of the chloride gradient) are taken into account. Activity-dependent potassium accumulation is another fast mechanism that may exacerbate reduction of E anion leading to incompensable disinhibition [69, 70]. In any case, given that E anion reduces on average to -49 mV following peripheral nerve injury [31], the disinhibition that results is most likely incompensable; indeed, the observation that ~20% of lamina I neurons were paradoxically excited by GABA under those conditions [31] suggests that an even larger fraction were incompensably disinhibited (given that disinhibition requires less reduction in E anion than paradoxical excitation).
Regardless of the precise value of E anion at which it occurs, incompensable disinhibition approximates the conditions of pharmacologically reduced inhibition, which previous work has shown to be sufficient to produce allodynia and hyperalgesia [9–17]. It logically follows that incompensable disinhibition resulting from reduced E anion is sufficient to explain the exaggerated pain perception associated with neuropathic pain. We could not reach that conclusion if disinhibition were compensable (i.e. if compensatory mechanisms could successfully prevent disinhibition) because, in that case, even if the reduction in E anion could cause allodynia and hyperalgesia, whether or not it did would depend on the success or failure of compensation. Notably, the argument that incompensable disinhibition is sufficient to explain allodynia and hyperalgesia does not exclude other mechanisms from contributing to the aberrant perception; for example, Figure 7 illustrates how peripheral sensitization can exacerbate the hyperexcitability caused by postsynaptic disinhibition. Furthermore, under conditions in which a modest reduction of E anion occurs, such that some inhibitory capacity remains (e.g. with E anion <-55 mV), a decrease in GABAergic or glycinergic transmission (either transmitter release or receptor function; see below) would contribute to the disinhibition caused by reduction of E anion.
Reduction of Eanion vs. other mechanisms of disinhibition
Although there is general consensus that disinhibition occurs following neuropathy, the mechanism underlying disinhibition has been controversial. Studies have reported that the number of inhibitory neurons in the spinal dorsal horn decreases following peripheral nerve injury [24, 26–28], but more recent work has argued against this [34–36]. Other studies have reported a reduction of presynaptic GABA, implicating the GABA transporter GAT1 [32] and the GABA synthesizing enzyme GAD65 [30], but again this has been contradicted by the demonstration of no change in synaptosomal glycine or GABA [33]. In fact, Kontinen et al. [57] reported that GABAergic transmission was increased following neuropathy, presumably as a compensatory change, which would be consistent with the increased GABA contribution to background input to lamina I neurons reported by Coull et al. [31]. Disinhibition through reduction of E anion does not involve reduced glycinergic or GABAergic transmission but instead works at a downstream locus and controls the potency of inhibitory input. Ironically, although the system has built-in redundancy, inasmuch as it uses both glycine and GABA as fast inhibitory neurotransmitters [71–73], the reduction in E anion subverts both glycine and GABAA receptor-mediated inhibition because of the receptors' common reliance on the transmembrane chloride gradient. This contrasts the mechanism responsible for disinhibition in inflammatory pain, where prostaglandin E2-induced phosphorylation of the glycine receptor decreases glycinergic transmission [74, 75] without affecting GABAergic transmission. Under those conditions, increased GABAergic transmission can compensate for decreased glycinergic transmission [76].
In addition to reducing E anion via decreased KCC2 expression, neuropathy leads to other pathophysiologic changes in lamina I and elsewhere in the pain pathway that undoubtedly contribute to the aberrant perception associated with neuropathic pain [1, 3, 77–79]. But whereas most other changes perturb a neuron's input-output relationship directly, disinhibition acts indirectly by perturbing a modulatory process. Glycine and GABAA receptor-mediated inhibition are crucial mechanisms for the endogenous modulation of pain [80, 81], controlling the inflow of Aδ and C fiber-mediated inputs at the segmental level [32, 53, 82–85] as well as participating in descending modulation originating from the periaqueductal gray and nucleus raphe magnus [86–89]. These control mechanisms are seriously compromised by reduction of E anion. Disinhibition therefore equates with modulation of a modulatory process, or meta-modulation, and represents a higher order change that not only perturbs the system, but simultaneously compromises the control mechanisms that would otherwise correct that perturbation. This may contribute to the relative intractability of neuropathic pain compared with other types of pain. The ideal therapy would target the primary pathological change and return E anion to its normal value; this would not only reestablish a normal input-output relationship within the system, but would also restore the system's full capacity to control the input-output relationship. No such therapy currently exists.
Implications for therapeutics
The mechanism underlying disinhibition does, nonetheless, have significant implications for reestablishing inhibition through therapeutic interventions. Specifically, disinhibition through reduction of E anion means increasing glycinergic and/or GABAergic transmission may be inconsequential, and potentially even detrimental, depending on the degree of change in E anion. However, several studies have reported that allodynia was reduced by applying GABAA receptor agonists [17, 20, 22] or by transplantation of GABA-producing cells into the lumbar subarachnoid space [19]. This could be explained by an increase in presynaptic inhibition, since primary afferent terminals do not express KCC2 and are therefore not prone to disinhibition by reduced KCC2 expression [31]. However, the efficacy of increasing GABAergic transmission is controversial: Stubley et al. [21] reported that transplanted-GABA producing cell could prevent allodynia if transplanted early after nerve injury but could not reverse allodynia once it was established, other studies have reported failure of GABAA agonists to relieve neuropathic pain after ischemic spinal cord injury [90, 91]. These inconsistencies may be attributable to variation in the model of neuropathic pain, but that implicitly assumes that the underlying mechanisms are qualitatively different depending on the model. The present study suggests that variation can also be explained by quantitative differences in the degree of E anion reduction, i.e. whether the system was fully decompensated or whether residual inhibitory capacity remained.
Several authors have advocated the benefits of optimizing the treatment of neuropathic pain by choosing treatments targeted towards mechanisms implicated in the pathogenesis of that pain [92–96]. This study has approached the topic of mechanism-based therapies using quantitative modeling (see Fig. 8). By being quantitative, our results highlight how choosing the optimal therapy might not depend solely on what pathogenic mechanism is involved, but on the degree of change that has occurred, e.g. whether the system would benefit from augmenting inhibition depends on how much E anion is reduced. This may explain the variable efficacy of treatments amongst patients who are suffering from the same neuropathic pain syndrome, where the underlying mechanism is presumably the same but in whom reduction of E anion may vary.
Conclusion
Reduction of E anion can dramatically reduce the inhibitory control of spiking in spinal lamina I neurons, but whether this gives rise to the exaggerated sensitivity characterizing neuropathic pain depends quantitatively on the success or failure of compensatory mechanisms. Understanding the efficacy of therapeutic interventions also depends on this quantitative understanding. These results speak to the importance of using quantitative, biophysically accurate models to bring together the vast panoply of experimental data so that we not only identify mechanisms involved in neuropathic and other types of pain, but so that we understand those mechanisms quantitatively and fully exploit them for clinical benefit.
Methods
All simulations were performed with NEURON simulation software [97] using a compartmental model of a generic spinal lamina I neuron with resting membrane potential (V rest) = -63 mV, input resistance (R in) = 470 MΩ, and membrane time constant (τmembrane) = 31 ms, based on average values in Prescott and De Koninck [54] and Coull et al. [31]. Dendrites bifurcated up to fourth order for a total of 60 dendritic compartments (see Fig. 1A). Axial resistivity was 150 Ω·cm. An axon similar to that described by Mainen et al. [98] was attached to the soma. Fast Na+ and delayed rectifier K+ conductances based on Traub and Miles [99] were inserted at 0.1 and 0.01 S/cm2, respectively, in the soma, axon initial segment, and axon nodes. Voltage threshold for spiking was approximately -49 mV. A passive leak conductance was distributed evenly throughout the neuron and was adjusted to produce the passive membrane properties described above. Confirmatory testing (see Fig. 2D and 2E) was also performed on two other model neurons with distinct intrinsic properties (Fig. 1C). The tonic-spiking model is the same as model 2 in Figure 9 of our earlier paper [100]. The single-spiking model was derived from the tonic-spiking model by removing I Na, P, I Ca, P, I Ca, T, and I K, S, and then inserting, at 0.2 mS/cm2, a low-threshold K+ current, I K, LT, modified from the M-type K+ current described in [101]; voltage at half-maximal activation was shifted to -45 mV and kinetics were sped up 100×. For the single-spiking model, all synaptic weights (see below) were tripled in order to overcome the low intrinsic excitability of this cell type.
Synaptic conductances were modeled as a rapid exponential rise in conductance combined with a slower exponential decay in conductance described by τrise and τdecay, respectively. Synaptic current is therefore written as I syn(t) = w [1-exp(-t/τrise)]exp(-t/τdecay) (V m-E rev) where t = 0 at the onset of a synaptic event, w is synaptic weight, V m is membrane potential, E rev is reversal potential, and τrise = 0.5 ms for all synapses. Excitation was mediated by four sets of excitatory synapses distributed throughout the dendrites. Each set consisted of 5 synapses, with 2–3 AMPA synapses and the remainder being NMDA synapses. For AMPA synapses, τdecay = 5 ms and w was set so that AMPA receptor-mediated events had a peak conductance of 333 pS/synapse [100]. For NMDA synapses, τdecay = 25 ms and w was left equal to that for AMPA synapses, such that NMDA-receptor mediated events contributed 14% of the total excitatory current in voltage clamp simulations at -60 mV [102]. The voltage-dependent magnesium block of the NMDA current was modeled after Jahr and Stevens [103, 104] [see also [105]]. Reversal potential (E rev) was 0 mV for both AMPA and NMDA receptor-mediated inputs.
Eight sets of 2–3 inhibitory synapses were distributed randomly in the perisomatic region; distributing the inhibitory synapses throughout the dendrites and soma, rather than perisomatically, had no significant effect on the efficacy of inhibition, as confirmed with a separate series of simulations. Half of the inhibitory neuron sets were modeled after glycine synapses with τdecay = 12 ms and while the other half were modeled after GABA synapses with τdecay = 60 ms [31]. Synaptic weight was adjusted to produce glycine receptor-mediated events with a peak conductance of 450 pS/synapse [72]. Synaptic weight was fivefold less for GABA synapses so that, given the fivefold increase in τdecay, GABA receptor-mediated events contributed roughly the same total conductance as glycine receptor-mediated events, as reported by Coull et al. [31]. Mixed glycine/GABAA receptor-mediated inhibition simulates the conditions following neuropathy [31], but switching to full glycinergic inhibition did not significantly affect firing rate modulation; results from mixed inhibition are therefore reported throughout. Reversal potential for inhibitory events (which equates with E anion) was tested at 5 mV increments between -70 mV and -45 mV, which encompasses the range expected for normal and neuropathic conditions [31].
Simulations were intended to mimic the bombardment of a lamina I neuron with input arriving mono- and polysynaptically from multiple primary afferent fibers, where the spike train in each afferent is assumed to be semi-random (i.e. the interspike interval distribution has a coefficient of variation (CV) between 0 and 1). The summation of those spike trains gives a cumulative spike train whose distribution has a CV approaching 1 (i.e. a Poisson distribution); each set of synapses was therefore driven by an independent Poisson process. For excitatory input, the rate of each Poisson process was multiplied by the number of synaptic sets (i.e. 4) so that f exc specifies the total rate of EPSPs received by the neuron from all presynaptic cells; note that simultaneous activation of multiple excitatory synapses constitutes a single EPSP, which explains why we multiplied by the number of synaptic sets rather than by the total number of synapses in order to calculate the EPSP rate. The EPSP rate reported by Furue et al. [106] using in vivo patch clamp falls within the range we tested. Notably, each set of synapses had a slightly different spatial distribution and a different constitution of AMPA and NMDA synapses so that not all primary afferent activity elicited identical EPSPs in the postsynaptic cell. The rate of IPSPs, reported as f inh, was also calculated by multiplying the Poisson process rate by 4 so that f exc and f inh are readily comparable. In keeping with the classical view that excitation is driven mainly by input from small diameter fibers while inhibition is driven mainly by input from large diameter fibers [80], where both types of fibers can be driven by the same peripheral stimulus, we posited that f exc and f inh are proportional, although not necessarily equivalent; the constant of proportionality is reported as α, where α = f inh/f exc. Since activity in large and small diameter fibers is independent, a temporal relation between EPSPs and IPSPs is not likely to exist, except at the stimulus onset because of the differential conduction velocity in differently sized fibers and the minimum number of intervening synapses (monosynaptic excitation vs. disynaptic inhibition). We have not explicitly modeled the stimulus onset and assume, for proportional inhibition, that excitation and inhibition are temporally independent. The variable delays introduced by signal transmission through polysynaptic pathways encourage this temporal independence, and also support our decision to model inputs as Poisson processes.
In one set of simulations, the time-averaged conductance at each value of f inh was calculated and subsequently applied as a constant conductance to the soma, thereby replacing the intermittent inhibition produced by irregular synaptic input with constant inhibition. For certain simulations, feedback inhibition was introduced by setting up a simple network (see Fig. 1C) in which the neuron of interest excited another neuron (with the same intrinsic properties as the first) which, in turn, inhibited the neuron of interest. Strength of excitatory synapses was adjusted so that a spike in the presynaptic neuron typically elicited a spike in the postsynaptic neuron, meaning the firing rates of both neurons were roughly equivalent. Based on the delay between spike generation in the output neuron and spike generation in the feedback neuron, and the subsequent synaptic delay, the output neuron experiences an IPSP ~8.5 ms after each spike. Feedback inhibitory synapses were identical to inhibitory synapses described above.
Simulated temperature was 23°C, which is consistent with the kinetics we used here for voltage- and ligand-gated currents. These results can, however, be safely extrapolated to in vivo conditions (i.e. 37°C). The main concern is that kinetics will be faster at warmer temperatures, which, in the case of inhibitory synaptic input, could result in larger inhibitory gaps. Data in Figure 5 argue that such gaps are relatively unimportant, which is consistent with the absence of any appreciable effect when switching between mixed GABA/glycinergic inhibition and full glycinergic input (see above). All Simulations were 20 s long. With such long simulations, standard deviation for firing rate was only ~0.5 Hz across multiple trials; most conditions were therefore tested with a single trial.
Abbreviations
- KCC2:
-
potassium-chloride cotransporter 2
- E anion :
-
anion reversal potential
- f exc :
-
rate of excitatory synaptic input
- f inh :
-
rate of inhibitory synaptic input
- f out :
-
rate of output spiking
- GABA:
-
γ-aminobutyric acid
- HH:
-
Hodgkin-Huxley
- IPSC:
-
inhibitory postsynaptic current
- IPSP:
-
inhibitory postsynaptic potential
- Ï„membrane :
-
membrane time constant.
References
Woolf CJ, Salter MW: Neuronal plasticity: increasing the gain in pain. Science 2000, 288: 1765–1769. 10.1126/science.288.5472.1765
Woolf CJ, Mannion RJ: Neuropathic pain: aetiology, symptoms, mechanisms, and management. Lancet 1999, 353: 1959–1964. 10.1016/S0140-6736(99)01307-0
Woolf CJ: Dissecting out mechanisms responsible for peripheral neuropathic pain: implications for diagnosis and therapy. Life Sci 2004, 74: 2605–2610. 10.1016/j.lfs.2004.01.003
Koltzenburg M, Scadding J: Neuropathic pain. Curr Opin Neurol 2001, 14: 641–647. 10.1097/00019052-200110000-00014
Besson JM: The neurobiology of pain. Lancet 1999, 353: 1610–1615. 10.1016/S0140-6736(99)01313-6
Dickenson AH: Balances between excitatory and inhibitory events in the spinal cord and chronic pain. Prog Brain Res 1996, 110: 225–231.
Wiesenfeld-Hallin Z, Aldskogius H, Grant G, Hao JX, Hokfelt T, Xu XJ: Central inhibitory dysfunctions: mechanisms and clinical implications. Behav Brain Sci 1997, 20: 420–425. 10.1017/S0140525X97261480
Sandkühler J: Neurobiology of spinal nociception: new concepts. Prog Brain Res 1996, 110: 207–224.
Yaksh TL: Behavioral and autonomic correlates of the tactile evoked allodynia produced by spinal glycine inhibition: effects of modulatory receptor systems and excitatory amino acid antagonists. Pain 1989, 37: 111–123. 10.1016/0304-3959(89)90160-7
Sherman SE, Loomis CW: Morphine insensitive allodynia is produced by intrathecal strychnine in the lightly anesthetized rat. Pain 1994, 56: 17–29. 10.1016/0304-3959(94)90146-5
Sivilotti L, Woolf CJ: The contribution of GABA(A) and glycine receptors to central sensitization: disinhibition and touch-evoked allodynia in the spinal cord. J Neurophysiol 1994, 72: 169–179.
Sherman SE, Loomis CW: Strychnine-dependent allodynia in the urethane-anesthetized rat is segmentally distributed and prevented by intrathecal glycine and betaine. Can J Physiol Pharmacol 1995, 73: 1698–1705.
Sherman SE, Loomis CW: Strychnine-sensitive modulation is selective for non-noxious somatosensory input in the spinal cord of the rat. Pain 1996, 66: 321–330. 10.1016/0304-3959(96)03063-1
Sorkin LS, Puig S: Neuronal model of tactile allodynia produced by spinal strychnine: effects of excitatory amino acid receptor antagonists and a mu-opiate receptor agonist. Pain 1996, 68: 283–292. 10.1016/S0304-3959(96)03130-2
Sorkin LS, Puig S, Jones DL: Spinal bicuculline produces hypersensitivity of dorsal horn neurons: effects of excitatory amino acid antagonists. Pain 1998, 77: 181–190. 10.1016/S0304-3959(98)00094-3
Loomis CW, Khandwala H, Osmond G, Hefferan MP: Coadministration of intrathecal strychnine and bicuculline effects synergistic allodynia in the rat: an isobolographic analysis. J Pharmacol Exp Ther 2001, 296: 756–761.
Malan TP, Mata HP, Porreca F: Spinal GABA(A) and GABA(B) receptor pharmacology in a rat model of neuropathic pain. Anesthesiology 2002, 96: 1161–1167. 10.1097/00000542-200205000-00020
Ugarte SD, Homanics GE, Firestone LL, Hammond DL: Sensory thresholds and the antinociceptive effects of GABA receptor agonists in mice lacking the beta3 subunit of the GABA(A) receptor. Neuroscience 2000, 95: 795–806. 10.1016/S0306-4522(99)00481-9
Hwang JH, Yaksh TL: The effect of spinal GABA receptor agonists on tactile allodynia in a surgically-induced neuropathic pain model in the rat. Pain 1997, 70: 15–22. 10.1016/S0304-3959(96)03249-6
Eaton MJ, Plunkett JA, Martinez MA, Lopez T, Karmally S, Cejas P, Whittemore SR: Transplants of neuronal cells bioengineered to synthesize GABA alleviate chronic neuropathic pain. Cell Transplant 1999, 8: 87–101.
Stubley LA, Martinez MA, Karmally S, Lopez T, Cejas P, Eaton MJ: Only early intervention with gamma-aminobutyric acid cell therapy is able to reverse neuropathic pain after partial nerve injury. J Neurotrauma 2001, 18: 471–477. 10.1089/089771501750171092
Rode F, Jensen DG, Blackburn-Munro G, Bjerrum OJ: Centrally-mediated antinociceptive actions of GABA(A) receptor agonists in the rat spared nerve injury model of neuropathic pain. Eur J Pharmacol 2005, 516: 131–138. 10.1016/j.ejphar.2005.04.034
Woolf CJ, Wall PD: Chronic peripheral nerve section diminishes the primary afferent A-fibre mediated inhibition of rat dorsal horn neurones. Brain Res 1982, 242: 77–85. 10.1016/0006-8993(82)90497-8
Castro-Lopes JM, Tavares I, Coimbra A: GABA decreases in the spinal cord dorsal horn after peripheral neurectomy. Brain Res 1993, 620: 287–291. 10.1016/0006-8993(93)90167-L
Stiller CO, Cui JG, O'Connor WT, Brodin E, Meyerson BA, Linderoth B: Release of gamma-aminobutyric acid in the dorsal horn and suppression of tactile allodynia by spinal cord stimulation in mononeuropathic rats. Neurosurgery 1996, 39: 367–374. 10.1097/00006123-199608000-00026
Ibuki T, Hama AT, Wang XT, Pappas GD, Sagen J: Loss of GABA-immunoreactivity in the spinal dorsal horn of rats with peripheral nerve injury and promotion of recovery by adrenal medullary grafts. Neuroscience 1997, 76: 845–858. 10.1016/S0306-4522(96)00341-7
Ralston DD, Behbehani M, Sehlhorst SC, Meng XW, Ralston HJ: Decreased GABA immunoreactivity in rat dorsal horn is correlated with pain behavior: a light and electron microscope study. Proceedings of the Eighth World Congress on Pain. Edited by: Jensen TS, Turner JA and Wiesenfeld-Hallin Z. Seattle, WA, IASP Press; 1997:547–560.
Eaton MJ, Plunkett JA, Karmally S, Martinez MA, Montanez K: Changes in GAD- and GABA- immunoreactivity in the spinal dorsal horn after peripheral nerve injury and promotion of recovery by lumbar transplant of immortalized serotonergic precursors. J Chem Neuroanat 1998, 16: 57–72. 10.1016/S0891-0618(98)00062-3
Fukuoka T, Tokunaga A, Kondo E, Miki K, Tachibana T, Noguchi K: Change in mRNAs for neuropeptides and the GABA(A) receptor in dorsal root ganglion neurons in a rat experimental neuropathic pain model. Pain 1998, 78: 13–26. 10.1016/S0304-3959(98)00111-0
Moore KA, Kohno T, Karchewski LA, Scholz J, Baba H, Woolf CJ: Partial peripheral nerve injury promotes a selective loss of GABAergic inhibition in the superficial dorsal horn of the spinal cord. J Neurosci 2002, 22: 6724–6731.
Coull JA, Boudreau D, Bachand K, Prescott SA, Nault F, Sik A, De Koninck P, De Koninck Y: Trans-synaptic shift in anion gradient in spinal lamina I neurons as a mechanism of neuropathic pain. Nature 2003, 424: 938–942. 10.1038/nature01868
Miletic G, Draganic P, Pankratz MT, Miletic V: Muscimol prevents long-lasting potentiation of dorsal horn field potentials in rats with chronic constriction injury exhibiting decreased levels of the GABA transporter GAT-1. Pain 2003, 105: 347–353. 10.1016/S0304-3959(03)00250-1
Somers DL, Clemente FR: Dorsal horn synaptosomal content of aspartate, glutamate, glycine and GABA are differentially altered following chronic constriction injury to the rat sciatic nerve. Neurosci Lett 2002, 323: 171–174. 10.1016/S0304-3940(02)00157-X
Polgar E, Hughes DI, Riddell JS, Maxwell DJ, Puskar Z, Todd AJ: Selective loss of spinal GABAergic or glycinergic neurons is not necessary for development of thermal hyperalgesia in the chronic constriction injury model of neuropathic pain. Pain 2003, 104: 229–239. 10.1016/S0304-3959(03)00011-3
Polgar E, Gray S, Riddell JS, Todd AJ: Lack of evidence for significant neuronal loss in laminae I-III of the spinal dorsal horn of the rat in the chronic constriction injury model. Pain 2004, 111: 144–150. 10.1016/j.pain.2004.06.011
Polgar E, Hughes DI, Arham AZ, Todd AJ: Loss of neurons from laminas I-III of the spinal dorsal horn is not required for development of tactile allodynia in the spared nerve injury model of neuropathic pain. J Neurosci 2005, 25: 6658–6666. 10.1523/JNEUROSCI.1490-05.2005
Price TJ, Cervero F, De Koninck Y: Role of cation-chloride-cotransporters (CCC) in pain and hyperalgesia. Curr Top Med Chem 2005, 5: 547–555. 10.2174/1568026054367629
Staley KJ, Mody I: Shunting of excitatory input to dentate gyrus granule cells by a depolarizing GABA(A) receptor-mediated postsynaptic conductance. J Neurophysiol 1992, 68: 197–212.
Gulledge AT, Stuart GJ: Excitatory actions of GABA in the cortex. Neuron 2003, 37: 299–309. 10.1016/S0896-6273(02)01146-7
Chance FS, Abbott LF, Reyes AD: Gain modulation from background synaptic input. Neuron 2002, 35: 773–782. 10.1016/S0896-6273(02)00820-6
Prescott SA, De Koninck Y: Gain control of firing rate by shunting inhibition: roles of synaptic noise and dendritic saturation. Proc Natl Acad Sci USA 2003, 100: 2076–2081. 10.1073/pnas.0337591100
Coull JAM, Beggs S, Boudreau D, Boivin D, Tsuda M, Inoue K, Gravel C, Salter MW, De Koninck Y: BDNF from microglia causes the shift in neuronal anion gradient underlying neuropathic pain. Nature 2005, 438: 1017–1021. 10.1038/nature04223
Jensen TS, Baron R: Translation of symptoms and signs into mechanisms in neuropathic pain. Pain 2003, 102: 1–8. 10.1016/s0304-3959(03)00006-x
Trevino DL: The origin and projections of a spinal nociceptive and thermoreceptive pathway. In Sensory Functions of the Skin. Edited by: Zotterman Y. Oxford, Pergamon Press; 1976:367–376.
Perl ER: Pain and nociception. In Sensory processes. Edited by: Darian-Smith I. Bethesda, American Physiological Society; 1984:915–975.
Craig AD: Pain mechanisms: labeled lines versus convergence in central processing. Annu Rev Neurosci 2003, 26: 1–30. 10.1146/annurev.neuro.26.041002.131022
Craig AD, Kniffki KD: Spinothalamic lumbosacral lamina I cells responsive to skin and muscle stimulation in the cat. J Physiol 1985, 365: 197–221.
Bester H, Chapman V, Besson JM, Bernard JF: Physiological properties of the lamina I spinoparabrachial neurons in the rat. J Neurophysiol 2000, 83: 2239–2259.
Craig AD, Krout K, Andrew D: Quantitative response characteristics of thermoreceptive and nociceptive lamina I spinothalamic neurons in the cat. J Neurophysiol 2001, 86: 1459–1480.
Andrew D, Craig AD: Responses of spinothalamic lamina I neurons to maintained noxious mechanical stimulation in the cat. J Neurophysiol 2002, 87: 1889–1901.
Craig AD, Andrew D: Responses of spinothalamic lamina I neurons to repeated brief contact heat stimulation in the cat. J Neurophysiol 2002, 87: 1902–1914.
Slugg RM, Campbell JN, Meyer RA: The population response of A- and C-fiber nociceptors in monkey encodes high-intensity mechanical stimuli. J Neurosci 2004, 24: 4649–4656. 10.1523/JNEUROSCI.0701-04.2004
Torsney C, MacDermott AB: Disinhibition opens the gate to pathological pain signaling in superficial neurokinin 1 receptor-expressing neurons in rat spinal cord. J Neurosci 2006, 26: 1833–1843. 10.1523/JNEUROSCI.4584-05.2006
Prescott SA, De Koninck Y: Four cell types with distinctive membrane properties and morphologies in lamina I of the spinal dorsal horn of the adult rat. J Physiol 2002, 539: 817–836. 10.1113/jphysiol.2001.013437
Eccles JC: The Physiology of Synapses. Berlin, Springer-Verlag; 1964.
Borg-Graham LJ, Monier C, Fregnac Y: Visual input evokes transient and strong shunting inhibition in visual cortical neurons. Nature 1998, 393: 369–373. 10.1038/30735
Kontinen VK, Stanfa LC, Basu A, Dickenson AH: Electrophysiologic evidence for increased endogenous GABAergic but not glycinergic inhibitory tone in the rat spinal nerve ligation model of neuropathy. Anesthesiology 2001, 94: 333–339. 10.1097/00000542-200102000-00024
Chizh BA, Headley PM: NMDA antagonists and neuropathic pain - multiple drug targets and multiple uses. Curr Pharm Des 2005, 11: 2977–2994. 10.2174/1381612054865082
Isaev D, Gerber G, Park SK, Chung JM, Randik M: Facilitation of NMDA-induced currents and Ca2+ transients in the rat substantia gelatinosa neurons after ligation of L5-L6 spinal nerves. Neuroreport 2000, 11: 4055–4061.
Kalso E: Sodium channel blockers in neuropathic pain. Curr Pharm Des 2005, 11: 3005–3011. 10.2174/1381612054865028
Yoshimura M, North RA: Substantia gelatinosa neurones hyperpolarized in vitro by enkephalin. Nature 1983, 305: 529–530. 10.1038/305529a0
Kuhn A, Aertsen A, Rotter S: Neuronal integration of synaptic input in the fluctuation-driven regime. J Neurosci 2004, 24: 2345–2356. 10.1523/JNEUROSCI.3349-03.2004
Morita K, Tsumoto K, Aihara K: Possible effects of depolarizing GABA(A) conductance on the neuronal input-output relationship: a modeling study. J Neurophysiol 2005, 93: 3504–3523. 10.1152/jn.00988.2004
Pezet S, McMahon SB: Neurotrophins: mediators and modulators of pain. Annu Rev Neurosci 2006, 29: 507–538. 10.1146/annurev.neuro.29.051605.112929
Qian N, Sejnowski TJ: When is an inhibitory synapse effective? Proc Natl Acad Sci USA 1990, 87: 8145–8149. 10.1073/pnas.87.20.8145
Staley KJ, Soldo BL, Proctor WR: Ionic mechanisms of neuronal excitation by inhibitory GABA(A) receptors. Science 1995, 269: 977–981.
Staley KJ, Proctor WR: Modulation of mammalian dendritic GABA(A) receptor function by the kinetics of Cl- and HCO3- transport. J Physiol 1999, 519: 693–712. 10.1111/j.1469-7793.1999.0693n.x
Cordero-Erausquin M, Coull JA, Boudreau D, Rolland M, De Koninck Y: Differential maturation of GABA action and anion reversal potential in spinal lamina I neurons: impact of chloride extrusion capacity. J Neurosci 2005, 25: 9613–9623. 10.1523/JNEUROSCI.1488-05.2005
Kaila K, Lamsa K, Smirnov S, Taira T, Voipio J: Long-lasting GABA-mediated depolarization evoked by high-frequency stimulation in pyramidal neurons of rat hippocampal slice is attributable to a network-driven, bicarbonate-dependent K+ transient. J Neurosci 1997, 17: 7662–7672.
Jarolimek W, Lewen A, Misgeld U: A furosemide-sensitive K+-Cl- cotransporter counteracts intracellular Cl- accumulation and depletion in cultured rat midbrain neurons. J Neurosci 1999, 19: 4695–4704.
Todd AJ: GABA and glycine in synaptic glomeruli of the rat spinal dorsal horn. Eur J Neurosci 1996, 8: 2492–2498. 10.1111/j.1460-9568.1996.tb01543.x
Chéry N, De Koninck Y: Junctional versus extrajunctional glycine and GABA(A) receptor-mediated IPSCs in identified lamina I neurons of the adult rat spinal cord. J Neurosci 1999, 19: 7342–7355.
Keller AF, Coull JA, Chery N, Poisbeau P, De Koninck Y: Region-specific developmental specialization of GABA-glycine cosynapses in laminas I-II of the rat spinal dorsal horn. J Neurosci 2001, 21: 7871–7880.
Harvey RJ, Depner UB, Wassle H, Ahmadi S, Heindl C, Reinold H, Smart TG, Harvey K, Schutz B, Abo-Salem OM, Zimmer A, Poisbeau P, Welzl H, Wolfer DP, Betz H, Zeilhofer HU, Muller U: GlyR alpha3: an essential target for spinal PGE2-mediated inflammatory pain sensitization. Science 2004, 304: 884–887. 10.1126/science.1094925
Zeilhofer HU: The glycinergic control of spinal pain processing. Cell Mol Life Sci 2005, 62: 2027–2035. 10.1007/s00018-005-5107-2
Poisbeau P, Patte-Mensah C, Keller AF, Barrot M, Breton JD, Luis-Delgado OE, Freund-Mercier MJ, Mensah-Nyagan AG, Schlichter R: Inflammatory pain upregulates spinal inhibition via endogenous neurosteroid production. J Neurosci 2005, 25: 11768–11776. 10.1523/JNEUROSCI.3841-05.2005
Cummins TR, Dib-Hajj SD, Black JA, Waxman SG: Sodium channels and the molecular pathophysiology of pain. Prog Brain Res 2000, 129: 3–19.
Julius D, Basbaum AI: Molecular mechanisms of nociception. Nature 2001, 413: 203–210. 10.1038/35093019
Lewin GR, Lu Y, Park TJ: A plethora of painful molecules. Curr Opin Neurobiol 2004, 14: 443–449. 10.1016/j.conb.2004.07.009
Melzack R, Wall PD: Pain mechanisms: a new theory. Science 1965, 150: 971–979. 10.1126/science.150.3699.971
Game CJ, Lodge D: The pharmacology of the inhibition of dorsal horn neurones by impulses in myelinated cutaneous afferents in the cat. Exp Brain Res 1975, 23: 75–84. 10.1007/BF00238730
Baba H, Ji RR, Kohno T, Moore KA, Ataka T, Wakai A, Okamoto M, Woolf CJ: Removal of GABAergic inhibition facilitates polysynaptic A fiber-mediated excitatory transmission to the superficial spinal dorsal horn. Mol Cell Neurosci 2003, 24: 818–830. 10.1016/S1044-7431(03)00236-7
Furue H, Katafuchi T, Yoshimura M: Sensory processing and functional reorganization of sensory transmission under pathological conditions in the spinal dorsal horn. Neurosci Res 2004, 48: 361–368. 10.1016/j.neures.2003.12.005
Seagrove LC, Suzuki R, Dickenson AH: Electrophysiological characterisations of rat lamina I dorsal horn neurones and the involvement of excitatory amino acid receptors. Pain 2004, 108: 76–87. 10.1016/j.pain.2003.12.004
Buesa I, Ortiz V, Aguilera L, Torre F, Zimmermann M, Azkue JJ: Disinhibition of spinal responses to primary afferent input by antagonism at GABA receptors in urethane-anaesthetised rats is dependent on NMDA and metabotropic glutamate receptors. Neuropharmacology 2006.
McGowan MK, Hammond DL: Intrathecal GABA(B) antagonists attenuate the antinociception produced by microinjection of L-glutamate into the ventromedial medulla of the rat. Brain Res 1993, 607: 39–46. 10.1016/0006-8993(93)91487-D
Antal M, Petko M, Polgar E, Heizmann CW, Storm-Mathisen J: Direct evidence of an extensive GABAergic innervation of the spinal dorsal horn by fibres descending from the rostral ventromedial medulla. Neuroscience 1996, 73: 509–518. 10.1016/0306-4522(96)00063-2
Lin Q, Peng YB, Willis WD: Inhibition of primate spinothalamic tract neurons by spinal glycine and GABA is reduced during central sensitization. J Neurophysiol 1996, 76: 1005–1014.
Peng YB, Lin Q, Willis WD: Effects of GABA and glycine receptor antagonists on the activity and PAG-induced inhibition of rat dorsal horn neurons. Brain Res 1996, 736: 189–201. 10.1016/0006-8993(96)00668-3
Hao JX, Xu XJ, Aldskogius H, Seiger A, Wiesenfeld-Hallin Z: Allodynia-like effects in rat after ischaemic spinal cord injury photochemically induced by laser irradiation. Pain 1991, 45: 175–185. 10.1016/0304-3959(91)90186-2
Xu XJ, Hao JX, Aldskogius H, Seiger A, Wiesenfeld-Hallin Z: Chronic pain-related syndrome in rats after ischemic spinal cord lesion: a possible animal model for pain in patients with spinal cord injury. Pain 1992, 48: 279–290. 10.1016/0304-3959(92)90070-R
Woolf CJ, Bennett GJ, Doherty M, Dubner R, Kidd B, Koltzenburg M, Lipton R, Loeser JD, Payne R, Torebjork E: Towards a mechanism-based classification of pain? Pain 1998, 77: 227–229. 10.1016/S0304-3959(98)00099-2
Block BM, Hurley RW, Raja SN: Mechanism-based therapies for pain. Drug News Perspect 2004, 17: 172–186. 10.1358/dnp.2004.17.3.829015
Smith PA: Neuropathic pain: drug targets for current and future interventions. Drug News Perspect 2004, 17: 5–17. 10.1358/dnp.2004.17.1.829021
Woolf CJ: Pain: moving from symptom control toward mechanism-specific pharmacologic management. Ann Intern Med 2004, 140: 441–451.
Harden RN: Chronic neuropathic pain. Mechanisms, diagnosis, and treatment. Neurologist 2005, 11: 111–122. 10.1097/01.nrl.0000155180.60057.8e
Hines ML, Carnevale NT: The NEURON simulation environment. Neural Comput 1997, 9: 1179–1209. 10.1162/neco.1997.9.6.1179
Mainen ZF, Joerges J, Huguenard JR, Sejnowski TJ: A model of spike initiation in neocortical pyramidal neurons. Neuron 1995, 15: 1427–1439. 10.1016/0896-6273(95)90020-9
Traub RD, Miles R: Neuronal Networks of the Hippocampus. Cambridge,U.K., Cambridge; 1991.
Prescott SA, De Koninck Y: Integration time in a subset of spinal lamina I neurons is lengthened by sodium and calcium currents acting synergistically to prolong subthreshold depolarization. J Neurosci 2005, 25: 4743–4754. 10.1523/JNEUROSCI.0356-05.2005
Mainen ZF, Sejnowski TJ: Influence of dendritic structure on firing pattern in model neocortical neurons. Nature 1996, 382: 363–366. 10.1038/382363a0
Dahlhaus A, Ruscheweyh R, Sandkuhler J: Synaptic input of rat spinal lamina I projection and unidentified neurones in vitro. J Physiol 2005, 566: 355–368. 10.1113/jphysiol.2005.088567
Jahr CE, Stevens CF: A quantitative description of NMDA receptor-channel kinetic behavior. J Neurosci 1990, 10: 1830–1837.
Jahr CE, Stevens CF: Voltage dependence of NMDA-activated macroscopic conductances predicted by single-channel kinetics. J Neurosci 1990, 10: 3178–3182.
Destexhe A, Mainen ZF, Sejnowski TJ: Kinetic models of synaptic transmission. In Methods in Neuronal Modeling. Volume 1. 2nd edition. Edited by: Koch C and Segev I. Cambridge, MA, MIT; 1998:1–25.
Furue H, Narikawa K, Kumamoto E, Yoshimura M: Responsiveness of rat substantia gelatinosa neurones to mechanical but not thermal stimuli revealed by in vivo patch-clamp recording. J Physiol 1999, 521 Pt 2: 529–535. 10.1111/j.1469-7793.1999.00529.x
Acknowledgements
This research was supported by the Natural Sciences and Engineering Research Council of Canada (YDK) and by the Howard Hughes Medical Institute (TJS). YDK is a senior scholar of the Fonds de la recherche en santé du Québec. SAP was supported by postdoctoral fellowships from the Canadian Institutes of Health Research and from the Human Frontier Science Program.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The author(s) declare that they have no competing interests.
Authors' contributions
SAP designed the study and performed all simulations. YDK helped conceive the study. All authors contributed to preparation of the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Prescott, S.A., Sejnowski, T.J. & De Koninck, Y. Reduction of anion reversal potential subverts the inhibitory control of firing rate in spinal lamina I neurons: towards a biophysical basis for neuropathic pain. Mol Pain 2, 32 (2006). https://doi.org/10.1186/1744-8069-2-32
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1744-8069-2-32