JN Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 95: 311-322, 2006. First published September 21, 2005; doi:10.1152/jn.00904.2005
0022-3077/06 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
95/1/311    most recent
00904.2005v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wicher, D.
Right arrow Articles by Borst, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wicher, D.
Right arrow Articles by Borst, A.

Peptidergic Counter-Regulation of Ca2+- and Na+-Dependent K+ Currents Modulates the Shape of Action Potentials in Neurosecretory Insect Neurons

Dieter Wicher1, Jens Berlau1, Christian Walther2 and Alexander Borst3

1Saxon Academy of Sciences, Department of Neurohormones, Jena; 2Institute of Physiology, Department of Neuroendocrinology, Philipps University Marburg, Marburg; and 3Max-Planck Institute of Neurobiology, Department of Systems and Computational Neurobiology, Martinsried, Germany

Submitted 30 August 2005; accepted in final form 18 September 2005


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
Influx of Ca2+ and Na+ ions during an action potential can strongly affect the repolarization and the fast afterhyperpolarization (fAHP) if a neuron expresses Ca2+- and Na+-dependent K+ currents (KCa and KNa). This applies to cockroach abdominal dorsal unpaired median neurons (DUMs). Here the rapid activation of KCa depends mainly on the P/Q-type Ca2+ current. Adipokinetic hormones (AKHs)—insect counterparts to mammalian glucagon—mobilize energy reserves but also modulate neuronal activity and lead to enhanced locomotor activity. Cockroach AKH I accelerates spiking and enhances the fAHP of octopaminergic DUM neurons, and it is generally held that enhanced release of the biogenic amine from these and other neurons may lead to general arousal. AKH I modulates the voltage-gated Na+ and P/Q-type Ca2+ current and the background Ca2+ current. Upregulation of P/Q-type Ca2+ current increases the KCa current, whereas enhanced inactivation of Na+ current decreases the KNa current. We quantified the hormone-induced changes in ion currents in terms of Hodgkin-Huxley models and simulated the resulting activity of DUM neurons. Upregulation of P/Q-type Ca2+ and KCa current enhanced the hyperpolarization but had a weak effect on spiking. Downregulation of Na+ and KNa current decreased hyperpolarization and slightly accelerated spiking. Superposition of these modulations produced an increase in fAHP while the spike frequency remained unchanged. Only when the upregulation of the pacemaking Ca2+ background current was included in the simulated modulation the model reproduced the experimentally observed AKH-I-induced changes. The possible physiological relevance of this dual effect is discussed in respect to transmitter release and synaptic integration.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
Insect adipokinetic hormones (AKHs) mobilize energy reserves from storage organs such as the fat body similar to glucagon in mammals (Gäde and Auerswald 2003Go; Van der Horst 2003Go). As found recently, they also enhance locomotor activity in response to starvation (Lee and Park 2004Go). In cockroach (Periplaneta americana), there are two AKH peptides (Scarborough et al. 1984Go; Witten et al. 1984Go), one of which, AKH I (pQVNFSPNWamide, also called Neurohormone D) (Baumann et al. 1984Go; Baumann and Penelin 1984Go) accelerates spiking and alters the shape of action potentials in abdominal dorsal unpaired median (DUM) neurons. These DUM neurons are of general interest inasmuch as they release the biogenic amine octopamine. The octopaminergic system of insects, comparable to the noradrenergic system of mammals, plays an important role in arousal and affects general activity levels (e.g., Roeder 2005Go).

Previous investigations revealed that AKH I modulates a whole set of ion currents in the DUM neurons (Wicher et al. 2001Go). AKH I leads to potentiation of a P/Q-type Ca2+ current due to channel phosphorylation via protein kinase A (PKA) (Wicher 2001bGo); reduction of the fast Na+ current, again via PKA (Wicher 2001aGo); potentiation of a BK-type Ca2+–activated K+ (KCa) current (Derst et al. 2003Go; Wicher et al. 1994Go); potentiation of a Ca2+ background current via Gq-mediated stimulation of phospholipase C (PLC) (Wicher et al. 2004Go). This latter current contributes to pacemaker depolarization and spontaneous firing of the abdominal DUM neurons. Other currents investigated, such as voltage-dependent K+ currents, were not affected by the peptide (Wicher et al. 1994Go).

The present study analyzes how the complex modulation of several currents leads to the increase in spike frequency and to the altered shape of the action potential. To simulate spiking of DUM neurons, we developed a one-compartment model that includes a set of 10 ion currents described in terms of the Hodgkin-Huxley formalism. Our studies indicate that accelerated spiking is solely due to the change in the Ca2+ background current, whereas the increase in fast afterhyperpolarization (fAHP) of action potentials (Wicher et al. 1994Go) requires the change in all other modulated currents.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
Cell isolation

Isolation of dorsal unpaired (DUM) neurons from the fifth abdominal ganglion of the cockroach P. americana was performed as described previously (Wicher et al. 1994Go). Briefly, the ganglia were excised, desheathed, and incubated for 10 min at room temperature in saline (for composition see bath solution used for spike recordings) containing 0.5 mg/ml trypsin (type II, Sigma, Deisenhofen, Germany) and 0.5 mg/ml collagenase (type I, Sigma). After washing off the enzyme, the ganglia were stored in saline for ≥1 h. Prior to the measurements, the somata of DUM neurons were separated using thin metal needles.

Electrophysiology

Ion currents in HEK293 cells and in DUM neurons were measured at room temperature using whole cell patch-clamp with appropriate compensation of series resistance and of capacitive and leakage currents. Pipettes having resistances of 2–4 M{Omega} (HEK cells) or 0.5–0.8 M{Omega} (DUM cells) were pulled from borosilicate capillaries. Current measurements and data acquisition were performed using an EPC9 patch-clamp amplifier controlled by PULSE software (HEKA Elektronik).

HEK CELLS. For experiments with Periplaneta Slowpoke (pSlo) channels, the bath solution contained (in mM) 135 NaCl, 5 KCl, 1 MgCl2, 1 CaCl2, 0.33 NaH2PO4, 2 Na-pyruvate, 10 glucose, and 10 HEPES. The pipette solution contained 140 KCl, 4 NaCl, 2 Mg-ATP, and 10 HEPES. The free Ca2+ concentration was adjusted to 150 µM by adding CaCl2 and measuring [Ca2+] with calcium-sensitive electrodes (KWIK tips; WPI, Berlin). The pH of bath and pipette solution were adjusted to 7.4 and 7.3, respectively.

DUM NEURONS. Spiking of neurons was recorded under current-clamp conditions without current injection. The bath solution contained (in mM) 190 NaCl, 5 KCl, 5 CaCl2, 2 MgCl2, and 10 HEPES (pH = 7.4), and patch pipettes (resistance >1.5 M{Omega}) were filled with solution composed of (in mM) 190 K-gluconate, 5 NaCl, 2 Mg-ATP, 1 CaCl2, 3 EGTA, and 10 HEPES (pH = 7.25). Between recordings (duration: 1 s) the cells were held at –70 mV under voltage clamp.

Na+ currents were separated as described (Wicher 2001aGo) with pipette solution composed of (in mM) 5 NaCl, 100 choline methylsulfate, 30 TEA-Br, 3 CsCl, 60 CsOH, 2 Mg-ATP, 1 CaCl2, 5 EGTA, and 10 HEPES and bath solution containing (in mM) 60 Na isethionate, 90 choline methylsulfate, 40 TEA-Br, 7 MgCl2, 1 CdCl2, and 10 HEPES. For Ca2+ current measurements (Wicher and Penzlin 1997Go), the pipette solution contained (in mM) 100 choline methylsulfate, 30 TEA-Br, 8 CsCl, 60 CsOH, 2 Mg-ATP, 1 CaCl2, 5 EGTA, and 10 HEPES, and the bath solution (in mM) had 5 CaCl2, 190 choline methylsulfate, and 10 HEPES and 5 µM tetrodotoxin.

K+ currents in DUM neurons were measured with pipette solution containing (in mM) 180 K-gluconate, 10 NaCl, 2 Mg-ATP, 1 CaCl2, 5 mannitol, 3 EGTA, and 10 HEPES. The bath solution to measure KCa currents contained (in mM) 190 NaCl, 5 KCl, 5 CaCl2, 2 MgCl2, and 10 HEPES and 0.5 µM tetrodotoxin (TTX). KCa currents were separated from other currents by applying the BK-channel blocker iberiotoxin (100 nM). The purely voltage gated delayed rectifier (DR) and the A-type current were measured with the above bath solution containing in addition 2 CdCl2. To isolate the DR, the A-type current was depressed by 3 mM 4-aminopyridine. The bath solution to measure KNa currents contained (in mM) 190 NaCl, 5 KCl, 1 CaCl2, 2 MgCl2, 2 CdCl2, and 10 HEPES. KNa currents, together with Na currents, were separated by comparing recordings performed in the absence and presence of the Na+ channel blocker TTX (0.5 µM). The pH of bath solutions was adjusted to 7.4 and that of pipette solutions to 7.25. Liquid junction potentials were compensated.

Application or washout of blocking agents was performed by transferring the cell (attached to the pipette tip) within a glass tube into the various solutions. A complete and fast solution change within a few milliseconds was achieved by sucking a small amount of solution into the tube.

Data analysis

Results are given as means ± SE (n = number of cells). Fitting procedures were performed with the software Prism2 (Graph Pad, San Diego, CA) and with IGOR (WaveMetrics, Lake Oswego, OR).

Chemicals

Salts, tetraethylammonium, tetrodotoxin, and verapamil were purchased from Sigma (Deisenhofen, Germany). {omega}-Agatoxin IVA and {omega}-conotoxin GVIA were obtained from Alomone Labs (Jerusalem, Israel) and iberiotoxin from Calbiochem (Schwalbach, Germany).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
Relationship between voltage-gated Ca currents and KCa currents in DUM neurons

Periplaneta DUM neurons have previously been shown to express pSlo, the {alpha} subunit of large conductance KCa (BK) channels (Derst et al. 2003Go). While pSlo channels heterologously expressed in HEK293 cells give rise to a purely noninactivating KCa current, the KCa current in DUM neurons is composed of a transient (KCa,t) and a sustained (KCa,s) component (Derst et al. 2003Go; Grolleau and Lapied 1995Go; Wicher et al. 1994Go). pSlo channels are sensitive to the BK channel blocker iberiotoxin (IbTx; IC50 = 45 nM). The KCa channels in DUM neurons seem to be more sensitive to IbTx because both KCa current components are blocked by 100 nM IbTx (Derst et al. 2003Go). This is illustrated in Fig. 1A, which shows sample traces of K+ currents, measured prior and after application of IbTx and the KCa current obtained by subtraction of these currents. There is an initial rapid KCa current rise (the time to peak ranges from 7 ms at –20 mV to 3 ms at ≥20 mV) and a fast decay (decay to 50% takes ~2 ms) of the first part of the KCa current which is dominated by KCa,t. It then reaches a sustained level, i.e., KCa,s, some 5 ms after stepping the voltage. Activation of both components of the KCa current starts around –50 mV, i.e., at about the same voltage as that of the start of the Ca2+ current (Fig. 1B, top). The time-to-peak of the KCa current is comparable to that of the voltage-gated Ca2+ current (Fig. 1D).



View larger version (30K):
[in this window]
[in a new window]
 
FIG. 1. KCa currents in dorsal unpaired median (DUM) neurons and their relation to Ca2+ currents. A: separation of the KCa current by iberiotoxin (IbTx). Currents activated by the voltage jump indicated were recorded before (control) and 2 min after application of 100 nM IbTx (top) and were digitally subtracted yielding the IbTx-sensitive KCa current (bottom) that is composed of a transient and a sustained component. B: I-V relationships for the transient KCa current (IK,Ca,t), the sustained KCa current (IK,Ca,s, measured 20 ms after start of depolarization) and for voltage-gated Ca2+ peak current in DUM neurons. KCa currents were separated as described in the preceding text (n = 7). All currents were normalized to the maximal Ca2+ peak current (V = 0 mV). The total voltage-gated Ca2+ current comprises 3 components shown in the lower I-V relation. The currents shown are those that were blocked by the L-type channel blocker verapamil (10 µM), the N-type channel blocker {omega}-conotoxin GVIA ({omega}-CgTx, 1 µM), or the P/Q-type channel blocker {omega}-agatoxin IVA ({omega}-AgaTx, 50 nM). Currents were normalized as in the preceding text. C: separation of Ca2+ current components by {omega}-AgaTx and verapamil. Top: currents activated by the voltage jump indicated (control), 1 min after application of {omega}-AgaTx, and 1 min after application of verapamil (in the continuous presence of {omega}-AgaTx). Subtraction yields the {omega}-AgaTx-sensitive P/Q-type current and the verapamil-sensitive L-type current (bottom). D: time to peak for KCa,t and total Ca current (n = 7).

 
It should be pointed out that there is no pharmacological tool to separate KCa,t and KCa,s. However, these components differ in their Ca2+ sensitivity: When voltage-gated currents are blocked and the intracellular [Ca2+] is raised by perfusing the cell with Ca2+–rich pipette solution, the first current evolving—as seen on repeated application of depolarizing voltage protocols—is KCa,t. KCa,s appears ~1 min later, indicating that KCa,s requires a higher [Ca2+]i than KCa,t to be activated by depolarization (Derst et al. 2003Go). The fact that KCa,t inactivated in spite of a constant [Ca2+]i indicates that inactivation is an intrinsic channel property perhaps produced by accessory {beta} subunits. It cannot be excluded that in our present measurements there is some overlap of currents, i.e., that there is a small contribution of KCa,s to the KCa peak current at least at higher depolarizations. Because, however, such error is expected to be small and will not have consequences for the results obtained in the following text, we will refer the KCa peak current as KCa,t.

The Ca2+ current activating positive to –50 mV is composed of three subtypes (Wicher and Penzlin 1997Go): P/Q-type current sensitive to {omega}-agatoxin IVA ({omega}-AgaTx), N-type current sensitive to {omega}-conotoxin GVIA ({omega}-CgTx), and L-type current sensitive to verapamil (Fig. 1B, bottom). The effect of {omega}-AgaTx (50 nM) and verapamil (10 µM) on Ca2+ current is illustrated in Fig. 1C; the concentrations used have been previously shown to block the respective currents within 1 min (Wicher and Penzlin 1997Go).

We applied the Ca2+ channel blockers to learn to what extent L-, N-, and P/Q-type Ca2+ currents might supply the Ca2+ for activation of the KCa current. The three Ca2+ currents cause a Ca2+ influx of different size and activation kinetics (Fig. 1, B, bottom, and C), and their suppression affected the KCa currents differently: {omega}-AgaTx (50 nM) led to reduction of both the transient and the sustained component of KCa current (Fig. 2, A and C), {omega}-CgTx (1 µM) had no effect on KCa current (n = 5, not shown), whereas verapamil (10 µM) had a clearly weaker effect on the transient component than {omega}-AgaTx and suppressed the sustained KCa current to a similar degree (Fig. 2, B and D). The lack of effect of N-type channel block might indicate a spatial separation of N-type Ca2+ channels and KCa channels so that Ca2+ entering through N-type channels hardly diffuses to KCa channels. On the other hand, the weak inhibition of the transient KCa current component by the block of the L-type channel may be largely due to the rather slow activation kinetics of L-type channel. As illustrated in Fig. 2B, the peak of L-type current is attained after the peak of the verapamil-sensitive KCa current. The rise in intracellular [Ca2+] due to the L-type current might thus overlap with the intrinsic inactivation of the channel conducting KCa,t. By contrast, the P/Q-type current activates more rapidly and reaches the peak earlier than the KCa current sensitive to {omega}-AgaTx (Fig. 2A). Thus the fast activating P/Q-type current seems to supply most of the Ca2+ required for the activation of the transient KCa current component. On the other hand, the sustained KCa current may rely on Ca2+ entering via both types of Ca2+ channels because it is reduced by about the same extent by blocking P/Q- or L-type channels (Fig. 2, C and D).



View larger version (22K):
[in this window]
[in a new window]
 
FIG. 2. Block of P/Q-type (A) and L-type (B) Ca2+ currents reduces KCa currents differently. A: example of a KCa current and a Ca2+ current blocked by {omega}-AgaTx (50 nM). Currents activated by the voltage jump indicated were recorded before and 1 min after application of {omega}-AgaTx. The displayed currents are the difference currents, i.e., control–{omega}-AgaTx. B: recordings of KCa current and Ca2+ current blocked by verapamil (10 µM). Currents activated by the voltage jump indicated were recorded before and 1 min after application of verapamil; displayed are the difference currents. The measurements of Ca2+ and KCa currents for A and B had to be carried out in separate cells and representative traces are co-displayed. C: I-V relationships for the transient KCa current (IK,Ca,t) and the sustained KCa current (IK,Ca,s) blocked by {omega}-AgaTx. Currents were separated as described in A, n = 5. D: I-V relationships for K,Ca,t and K,Ca,s current blocked by verapamil. Currents were separated as described in A and B, n = 6.

 
Peptidergic upregulation of KCa current

The Periplaneta adipokinetic hormone AKH I, a peptide released from the corpora cardiaca, is known to enhance KCa currents (Wicher et al. 1994Go). AKH I also potentiates the P/Q-type Ca2+ but not the N- and L-type Ca2+ currents in DUM neurons (Wicher 2001bGo). Figure 3 demonstrates the AKH-induced effect on KCa currents: the total outward current produced by a voltage step to 0 mV before and after application of 10 nM AKH I, and the difference current representing the peptide-induced KCa current are shown in Fig. 3A, 1 and 2, respectively. AKH I thus upregulates both the transient and the sustained KCa component. This may be plausibly attributed to the upregulation of the P/Q-type current as this current was seen in the preceding text to provide Ca2+ for the activation of both KCa components. The threshold concentration for the AKH I-induced potentiation of the P/Q-type Ca2+ current was 1 pM (Wicher 2001bGo). At this concentration, there was already some increase in the KCa current (Fig. 3A3), which, however, was mainly restricted to the transient component (Fig. 3B1). Increasing [AKH I] from 1 pM to 10 nM induced progressively larger KCa currents (Fig. 3C). Compared with the KCa current under control conditions (isolated with IbTx), the current produced by 10 nM AKH I differed in that the transient component activated at more negative potentials. The nonlinearity in the I-V relation for this component around –10 mV (Fig. 3B2) is especially remarkable. It seems to reflect the increase in P/Q-type current by AKH I that is most pronounced at voltages ranging from –20 to 0 mV (Wicher 2001bGo). The possibility that AKH I affects purely voltage-gated K+ currents such as the delayed rectifier (IK DR) and the A-type current (IK A) could be ruled out because in the presence of 1 mM Cd2+, which completely blocks the Ca2+ currents (Wicher and Penzlin 1997Go), AKH I failed to affect K+ currents (n = 5 cells, not shown). Taken together, these results are compatible with the assumption that the upregulation of the KCa current is caused by the increase in P/Q-type current.



View larger version (23K):
[in this window]
[in a new window]
 
FIG. 3. Modulation of KCa currents by the peptide AKH I. A1: current traces recorded on jumps to 0 mV before (control) and 2 min after application of 10 nM AKH I. A2: KCa current induced by 10 nM AKH I (difference current of traces shown in A1). A3: KCa current induced by 1 pM AKH I (difference current: 2 min AKH I–control); note the onset of the Ca2+ inward current on which the KCa current is superimposed. B1: I-V relationships for transient and sustained KCa current induced by 1 pM AKH I. Currents were separated as described in A1, n = 4. B2: I-V relationships for transient and sustained KCa current induced by 10 nM AKH I. Currents were separated as described in A1, n = 5. C: concentration dependence of AKH I effect on transient KCa current activated by depolarization to –20 and +20 mV. Currents were separated as described in A1, n = 4.

 
The increase in KCa current is expected to increase the fast afterhyperpolarization (fAHP) of the action potential as has been previously observed in DUM neurons in the presence of AKH I (Wicher et al. 1994Go). Conversely, a reduction of KCa current was seen to prolong action potentials and to reduce the hyperpolarization (Derst et al. 2003Go).

Activation of PKA does not affect IKCa

Some Slo channels are known to be modulated by phosphorylation via protein kinases such as PKA (Schubert and Nelson 2001Go). AKH I in fact activates PKA, which is a necessary step in the upregulation of the P/Q-type Ca2+ current (Wicher 2001bGo). In principle the regulation of the KCa currents by AKH I might be dual, i.e., be partly due to KCa channel phosphorylation. To evaluate whether PKA can affect the pSlo channel, we coexpressed pSlo together with the Periplaneta AKH receptor (pAKHR) in HEK293 cells. Stimulation of heterologously expressed pAKHR with AKH I was previously shown to increase [cAMP] via activating Gs proteins in HEK293 cells (Wicher, unpublished observation). Application of AKH I (10 nM), after loading the cell with a Ca2+-rich pipette solution, had, however, no effect on pSlo currents activated by depolarization (Fig. 4A, left). There was neither a change in the size of currents nor in the kinetics. In another series of experiments, we induced a rise in [cAMP] by bath application of the membrane-permeable cAMP analogue 8-bromo-cAMP. Again, there was no change in the properties of pSlo currents (Fig. 4A, right). To exclude the possibility that PKA might be fully activated under control conditions, we tested whether inhibition of PKA would affect pSlo currents. Application of 10 µM KT5720, which was previously shown to completely abolish the AKH I effect on P/Q-type Ca2+ current (Wicher 2001bGo), did not change pSlo currents. These results, which are summarized in Fig. 4B, do not support the possibility that PKA regulates pSlo. Recently direct binding of a PKA and a Src tyrosin kinase to domains in the C-terminus of the Drosophila Slo channel (dSlo) and phosphorylation of the dSlo channel protein by both kinases has been demonstrated (Wang et al. 1999Go). However, no difference in peak current amplitude nor in voltage dependence of dSlo gating was observed after coexpression of dSlo with either protein kinase (Wang et al. 1999Go). On the other hand, the free catalytic subunit of PKA binds to dSlo and leads to downregulation of channel activity (Zhou et al. 2002Go). This modulation, however, does not involve phosphorylation of the only consensus PKA-substrate site in the C-terminal domain of dSlo.



View larger version (15K):
[in this window]
[in a new window]
 
FIG. 4. Test for a direct effect of AKH I on KCa channel. pSlo, the {alpha} subunit of Periplaneta KCa channel, was coexpressed with the pAKH receptor in HEK293 cells. A: sample traces of pSlo currents activated by jumps from –90 to +80 mV before (control) and 2 min after application of AKH I (10 nM) or 8-bromo-cAMP (8-br-cAMP, 2 µM). Note that neither AKH I nor 8-br-cAMP affected pSlo currents. B: effect of AKH I (10 nM), 8-br-cAMP (2 µM), and KT5720 (10 µM) on pSlo currents. Currents measured in the presence of these agents were normalized to control currents. n = 5.

 
Relationship between voltage-gated Na currents and KNa currents in DUM neurons

DUM neurons express Na+-dependent K+ channels (Wicher et al. 1994Go; Grolleau and Lapied 1994Go). Therefore blocking the Na+ current with tetrodotoxin (TTX) leads to disappearance of a transient outward current component (Fig. 5A), and the total TTX-sensitive current is the sum of the Na+ current and the Na+-dependent K+ (KNa) current (Fig. 5B). Due to the lack of a specific blocker, it is impossible to separate the KNa current from the Na+ current (Grolleau and Lapied 1994Go). Although the TTX-sensitive outward current is contaminated with the Na+ current, we will refer to it as the KNa current. It activates very rapidly, its peak following that of the Na+ current within <1 ms (Fig. 5C). Because KNa channels are not voltage-gated, the I-V relation of the KNa peak current mirrors the I-V relation of the Na+ current (Fig. 5D).



View larger version (20K):
[in this window]
[in a new window]
 
FIG. 5. KNa currents and Na+ currents in DUM neurons. A: sample traces of currents activated by depolarizations as indicated in B before (control) and 2 min after application of TTX (0.5 µM). B: biphasic current, i.e., the difference of currents shown in A results from superposition of KNa current and Na+ current. For comparison, a Na+ current is shown that was recorded when K+ currents are suppressed. C: voltage dependence of the time to peak for KNa and Na+ currents (n = 5). D: I-V relationships for voltage-gated Na+ peak current and for KNa peak current in DUM neurons. KNa and Na+ currents were separated by means of TTX as demonstrated in A and B (n = 7). The size of Na+ (and KNa) currents is thus underestimated, yet this procedure was adequate for comparing the data with the currents suppressed by AKH I. E: KNa plus Na+ current suppressed by 10 nM AKH I (difference current: control–2 min AKH I). F: I-V relationships for KNa peak current and for Na+ peak current suppressed by 10 nM AKH I (control–2 min AKH I, n = 7).

 
Peptidergic upregulation of KNa current

AKH I accelerates the inactivation of the Na+ current in DUM neurons, thereby reducing both peak size and duration and thus also the net Na+ influx (Wicher 2001aGo). Because activation of the KNa current requires a high intracellular Na+ concentration (Dryer 1994Go), any reduction of Na+ influx should entail a reduction of the KNa current. An example of a current suppressed by AKH I (i.e., the difference of the current measured under control conditions and in presence of AKH I) is shown in Fig. 5E. The relation between the AKH-I-sensitive Na+ current and the AKH-I-sensitive KNa current (Fig. 5F) resembles that between the TTX-sensitive Na+ current and KNa current under control conditions (Fig. 5D). Thus the reduction of KNa peak current was proportional to the reduction of Na+ peak current, and the downregulation of KNa current by AKH I, therefore is probably solely due to the reduction of the Na+ influx. A test of a possible effect of PKA on the KNa channel (Slo2 or slack) (Yuan et al. 2003Go) could not be performed because this channel has not yet been cloned in Periplaneta. Furthermore, an experimental analysis of the role of the peptidergic modulation of Na+ current and KNa current on action potential shape is impaired by the lack of tools to separate the currents: Li+ permeates Na+ channels but fails to activate KNa channels in some preparations. By contrast, in DUM neurons Li+ also activates the KNa current, and it can thus not be used as tool to separate the Na+ current from KNa currents (Grolleau and Lapied 1994Go).

Modeling the effect of peptidergic counter-regulation of Ca2+- and Na+-dependent K+ currents on the action potential shape

We performed a modeling study of the differential modulation of the K+ currents with the aim of answering two main questions: Does modulation of these currents affect the pacemaker depolarization and does superposition of the modulatory effects explain quantitatively the observed increase in fast AHP (fAHP) of action potentials?

To simulate the endogeneous spiking of DUM neurons, we included 10 currents in a one-compartment model of a DUM neuron: one Na+ current, four Ca2+ currents (the background current, a low-voltage activated current, and 2 high-voltage activated currents, i.e., P/Q-type and non-P/Q-type current), and five K+ currents (the delayed rectifier, KDR, the A-type current, KA, the KNa current and the KCa currents, KCa,t and KCa,s). The currents were described in terms of the Hodgkin-Huxley formalism (cf. APPENDIX, Fig. A1 and Table 1). Figure 6 compares recorded spikes from a DUM neuron with the spike pattern generated by the model, using the parameters given in Table 1. The model reproduced the experimentally observed characteristics of the action potential such as threshold, overshoot, and AHP as well as the resting firing frequency.



View larger version (30K):
[in this window]
[in a new window]
 
FIG. A1. Description of DUM cell ion currents in terms of Hodgkin-Huxley models. Current separation was performed as described in METHODS. Symbols represent means of 5–9 cells, bars denote SE. The curves are fits of Hodgkin-Huxley models to the data described by parameters given in Table 1. Columns show steady-state activation (minf) and inactivation (minf), activation time constant {tau}m, and inactivation time constant {tau}h for Na+ current (A) (cf. Wicher 2001aGo), P/Q- and non-P/Q-type Ca2+ currents (B), delayed rectifier (DR) K+ current (C), A-type K+ current (D) (cf. Grolleau and Lapied 1995Go), KNa current (E) (cf. Grolleau and Lapied 1994Go), and transient and sustained KCa current (F; IKCa,t and IKCa,s).

 

View this table:
[in this window]
[in a new window]
 
TABLE 1. Parameters for voltage- and ion-dependent membrane currents

 


View larger version (17K):
[in this window]
[in a new window]
 
FIG. 6. Measured (A) and modeled (B) spiking of a DUM neuron. A: current-clamp recording (no current injected). B: calculated activity produced by the model described in the APPENDIX using the parameters given in Table. 1.

 
To get insight into the role of the KNa current and the transient KCa current in the regulation of neuronal activity, we investigated the effects of varying the size of these currents. The KNa current was found to be essential for the stability of repetitive activity (Fig. 7A). Both the KNa and the KCa,t current affected the spike frequency (Fig. 7A) and the fAHP (Fig. 7C) although with different efficiency. On reduction of the KNa current, there was a robust increase in spike frequency, whereas an increase produced only a weak reduction (Fig. 7A). Changing the KCa,t current led to inversely proportional changes in spike frequency. When KCa,t was increased by >30%, spiking stopped because the action potential threshold was no longer reached (Fig. 7A). On the other hand, spiking did not become irregular even when the KCa,t conductance was reduced by 90%. Reduction of the KNa conductance by >40% caused irregular spiking. Reduction of KCa,t caused a slight attenuation of the fAHP. Complete block of KCa,t was estimated to reduce the AHP by 4 mV, which was somewhat less than found experimentally, i.e., after blocking KCa with iberiotoxin (9 mV) (Wicher 2001bGo). The dependence of fAHP on IKNa was similar to that on IKCa,t but considerably steeper (Fig. 7B).



View larger version (27K):
[in this window]
[in a new window]
 
FIG. 7. Changes of spike frequency (A and B) and fast afterhyperpolarization, fAHP (C and D), on varying the strengths of the currents INa, I K,Na, ICa, and IK,Ca,t. Changes produced by varying either IK,Na or IK,Ca,t are shown in A and C. B and D: demonstration of what happens when the Na+ + KNa current and the P/Q-type Ca2+ and KCa currents change concomitantly in the manner observed on application of 10 nM AKH I. Grel = 1 corresponds to the Gmax values given in Table. 1. The change in Grel of INa was produced by changing the inactivation according to the modulation by AKH I, i.e., by increasing {tau} h slope up and decreasing {tau} h slope dn (see RESULTS) until the amplitude of INa corresponded to a given figure of Grel. A direct change of Gmax (with no change in inactivation kinetics)—as performed with the Ca2+ current—would lead to drastically different effects (cf. text). When 1 set of parameters was varied the others were fixed at Grel = 1. *, instability of spiking.

 
We further investigated how the combined variation of the inward currents, the Na current and the P/Q-type Ca current and the coupled K+ currents, i.e., KNa and KCa, affect spiking. Increasing the Na+ conductance unexpectedly produced a rise in frequency (by 100% on increasing G by 30%) and attenuated the AHP (by 4 mV on increase by 30%) while reducing the Na+ conductance terminated spiking because the pacemaker activity became too weak (not shown). We then adjusted the Na+ conductance solely by changing the inactivation parameters thus imitating the Na+ current-modulation caused by AKH I (Fig. 7, B and D). Under these conditions, the variation of Na+ current produced changes more similar to those found for varying the KNa conductance alone (compare Fig. 7, B with A). Reduction of Na+ conductance caused a steep increase in frequency and a decrease in hyperpolarization until spiking became irregular at >40% reduction (Fig. 7, B and D). A rise in Na+ conductance did not affect both the spike frequency and the AHP (Fig. 7, B and D). Variation of the P/Q Ca2+ current had no effect on spike frequency but on AHP (Fig. 7, B and D). With G reduced by >50% the AHP is only –50 mV, and it saturates at –62 mV when G is increased by >50%.

The main conclusions about the role of the KNa current and the transient KCa current in the regulation of spiking are there is a critical size of KNa current that is required for stable spiking (Fig. 7, A and B) and the AHP is determined by up- and downregulation of the KCa current by the P/Q-type Ca2+ current and by downregulation of the KNa current by the Na+ current.

Application of 10 nM AKH I produces faster spiking of DUM neurons and affects the action potential shape. The spike frequency increase ranges from 15 to 34%, the mean is 25.2 ± 4.0% (n = 7, Fig. 8A). The hyperpolarization increases by 2–6 mV, on average by 4.0 ± 1.0 mV (n = 7, Fig. 8B). The question, then, is whether our model reproduces these changes and whether this requires that in fact all observed AKH I-induced changes in ion currents are implemented, i.e., upregulation of P/Q-type Ca2+ current and KCa currents, downregulation of Na+ current, and KNa current as well as upregulation of Ca2+ background current. In terms of our model (cf. Table 1), this means the following changes: G (Ca P/Q) from 209 to 287 nS, {tau} h Slope up (Na) from +11 to +13 mV, {tau} h Slope dn (Na) from –22 to –18 mV, and G (Ca, back) from 0.023 to 0.07 nS. These changes produced indeed an accelerated spiking. In the example shown in Fig. 8C, the increase in spike frequency amounts to 25%; this is in accordance with the mean AKH I effect. Furthermore, the hyperpolarization is increased by 3 mV (Fig. 8D). The model also predicts a slight reduction in overshoot (by 3 mV, Fig. 8C), which was, however, not observed with DUM neurons (Fig. 8A).



View larger version (27K):
[in this window]
[in a new window]
 
FIG. 8. Measured (A and B) and calculated (C and D) effect of AKH I (10 nM) on the activity (A and C) and the shape of an action potential (black traces, control; red traces, 2 min after application of 10 nM AKH I). AKH I accelerates spiking both in the neuron (A) and in the model (C). Also, AKH I leads to stronger hyperpolarization of measured (B) and calculated (D) action potentials. The data used for simulation of the AKH effect are given in the text.

 
To assess whether the initial size of KNa and KCa currents may affect the simulation of the AKH effect, we varied their conductances by ±10% and repeated the simulation. Variations of the KNa current did not affect the spike frequency increase. When the initial G was increased by 10% the only effect was a slight reduction of hyperpolarization-increase by 1 mV. A similar effect on hyperpolarization was seen when the initial KCa currents were reduced by 10%. While this variation had no effect on the increase in spike frequency, the variation of initial KCa currents by +10% led to an AKH-induced acceleration by 34% which matches the upper limit of the AKH I effect observed in DUM neurons. Thus the model sufficiently and stably reproduces the changes in neuronal activity by AKH I.

We finally simulated a situation where AKH I would only lead to upregulation of the P/Q-type Ca2+ current and both KCa currents (Fig. 9A). This caused an increase in hyperpolarization (by 5 mV), which is in line with the measured AKH I effect (4 mV), whereas the spike frequency was hardly affected (reduction by 2%). Downregulation of Na+ current and KNa current accelerated spiking (by 7%) and slightly reduced the hyperpolarization (by 2 mV; Fig. 9B). Another effect was a slight reduction of the overshoot (by 3 mV). When the modulations of P/Q-type Ca2+ current and KCa currents were combined with those of Na+ and KNa currents, the change of the hyperpolarization was in accord with that observed experimentally (increase by 3 mV; Fig. 9C). The spike frequency, however, remained virtually unaffected (increase by 2%). When in addition the upregulation of Ca2+ background current was implemented, the result was an increase in spike frequency by 25% with no further effects on the action potentials i.e., the observed AKH I effects were now correctly reproduced as already mentioned in the preceding text (Fig. 9D).



View larger version (21K):
[in this window]
[in a new window]
 
FIG. 9. Simulation of the effects on spiking obtained by AKH-I-induced modulations of ion currents. Black traces, activity under control conditions; red traces, activity changed due to the indicated modulations. A: upregulation of P/Q-type Ca2+ and KCa current. B: downregulation of Na+ and KNa current. C: combined modulations of P/Q-type Ca2+ and KCa current and downregulation of Na+ and KNa current. D: like C, but the upregulation of Ca2+ background current is included. In this situation, i.e., when all AKH-I-induced modulations are implemented, the model reproduces the experimentally observed effects (for further explanations, cf. text).

 
In conclusion, the modulation of two voltage-gated currents (P/Q-type Ca2+ current and Na+ current) together is involved in shaping the action potential, but it does not affect spike frequency. In line with previous findings (Wicher et al. 2004Go), the acceleration of spiking by AKH I solely relies on the upregulation of the Ca2+ background current.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
Octopamine, the chemical messenger released by the DUM neurons investigated here, is known to upregulate general activity levels (Roeder 2005Go; Saraswati et al. 2004Go). In Drosophila, it has recently been shown that, on starvation, AKH-release produces enhanced locomotor activity (Lee and Park 2004Go). Conceivably, the physiological relevance of the AKH I-induced modulation of spiking lies in the enhanced release of octopamine both to local targets (e.g., respiratory muscles) and to the hemolymph thereby contributing to a rise of systemic octopamine level.

Neuropeptides, in general, affect not a single but several ionic currents in the excitable cells they act on, e.g., neurons and skeletal muscles of invertebrates (Brezina et al. 2005Go; Nusbaum et al. 2001Go; Wicher et al. 2001Go) or neurons of vertebrates (e.g., Pena and Ramirez 2004Go). In the present study on cockroach DUM neurons, we have completed knowledge of the neuro-modulatory actions on ion currents of the hormone AKH I, and we have tested, in a model of this neuron, how the interplay of these effects leads to the observed changes in firing and time course of action potentials.

To understand the role of different ion currents in the discharge patterns of neurons, various computational models have been developed for well studied examples from invertebrates and vertebrates (e.g., Dale 1995Go; Golowasch et al. 1992Go). This approach has also been used in analyzing the modulation of currents by chemical messengers (e.g., Harris-Warrick et al. 1995Go). Neuronal models designed to simulate spontaneous activity either implement KCa currents (Buchholtz et al. 1992Go), KNa currents (Dale 1995Go), or KNa and slowly activating KCa currents (Dale and Kuenzi 1997Go). In our model, the KNa current and the KCa,t current, two rapidly activating currents contributing to action potential repolarization, act in parallel and are both subject to modulation. In contrast to various other models, ours lacks the hyperpolarization-activated cation current, Ih. This type of current, which has not been found in cockroach DUM neurons (Grolleau and Lapied 2000Go) but in Drosophila CNS (Marx et al. 1999Go) and in bee olfactory receptor neurons (Gisselmann et al. 2003Go), is present in some but not all spontaneously active insect neurons (Wicher et al. 2001Go). Generally, for the neurons to spontaneously spike in the tonic mode, only some persistent inward current plus a low-threshold Ca2+ current seem to be required. In the case of cockroach DUMs, a Ca2+ background current, a low-voltage-activated Ca2+ current and a persistent Na+ current contribute to pacemaking (Grolleau and Lapied 2000Go).

How AKH I leads to changes of Ca2+- and Na+-dependent K+ currents

Because AKH I produces changes in Ca2+ and Na+ currents, the KCa and the KNa current must also change. For IKCa, we experimentally ruled out the possibility of an additional direct modulation of the respective channels. For IKNa, this is also rather unlikely because our model calculations returned realistic results without implementing any direct modulation of this current. The observed slight discrepancy with respect to the amplitude of the action potential (i.e., reduction in overshoot by 3 mV) indicates that in the real neurons the effect of Na+ current reduction on overshoot is compensated for by the reduction of KNa current. There might be a more sophisticated relationship between Na+ and KNa current than included in our model.

Modulation of the KCa current

Although three types of voltage-gated Ca2+ channels are expressed in DUM neurons (Fig. 1B), only P/Q-type channels appear to provide Ca2+ for activation of the transient KCa current component. The sustained KCa component can additionally be activated by Ca2+ influx through L-type channels (Fig. 2). Coupling of BK channels to L- and Q-type channels has been reported for rat chromaffin cells (Prakriya and Lingle 1999Go). In mouse sympathetic neurons, however, BK channels are linked to N-type channels, whereas L- and P-type channels activate Ca2+-activated Cl currents (Martinez-Pinna et al. 2000Go). In hippocampal neurons, N-type channels activate BK channels, L-type channels activate small conductance (SK) KCa channels and P/Q-type channels fail to activate either type of channel (Marrion and Tavalin 1998Go).

Blocking the BK currents in DUM neurons under current-clamp conditions demonstrated that only the rapidly activating, transient component KCa,t affects the shape of the action potential. KCa,t controls its duration (for further K+ currents involved in this process, cf. Wicher et al. 2001Go). In addition KCa,t determines the size of fAHP (Derst et al. 2003Go). These roles of the BK current are comparable to those in various mammalian neurons (Faber and Sah 2002Go). The sustained KCa component, according to the predictions of our model, controls the interspike interval similar to the purely voltage-dependent K currents KA and KDR although, of course, in Ca2+-dependent fashion (not shown).

Modulation of the KNa current

KNa channels require a relatively high intracellular Na+ concentration (~50 mM) to become activated. In some preparations, prolonged discharges or long-lasting depolarizations are required to activate IK,Na (Dryer 1994Go). Co-localization of Na+ and KNa channels may, however, allow a sufficient rise in [Na+] to activate the KNa current already during a single action potential (Koh et al. 1994Go). In hippocampal CA1 neurons, for example, the KNa current is responsible for the fAHP after a single action potential (Liu and Stan Leung 2004Go).

In bursting neocortical neurons, the KNa current is the main cause of the postexcitatory hyperpolarization (Franceschetti et al. 2003Go). In these cells, the KCa current, too, makes an albeit small contribution to the hyperpolarization. In DUM neurons, the fAHP is about equally dependent on KNa and KCa current, the former activating slightly more rapidly than the latter (which accords to the somewhat different activation kinetics of the Na+ and the Ca2+ currents). The modulations of the P/Q-type Ca2+ current and the Na+ current by AKH I change the initial proportion of KCa to KNa current in favor of the KCa current.

Why counterregulation of KCa and KNa currents?

Would one considers acceleration of spiking as the main effect of AKH I on the DUM neurons, then it remains enigmatic why a whole combination of currents (P/Q-type Ca2+ current, Na+ current, KCa and KNa current) is modulated by the peptide because this does not produce a change in spike frequency. In this light, it would be sufficient to upregulate the Ca2+ background current to get the observed acceleration of spiking.

Spike frequency, however, is not the sole parameter involved in controlling secretion. It is known from vertebrate neurons that presynaptic voltage-gated Ca2+ currents (N and P/Q types) play a central role in triggering transmitter release (Reid et al. 2003Go). Although we have no evidence that the upregulation of P/Q-type current in the soma of DUM neurons observed with AKH I also occurs in the release zones of its terminals, this would be a possible means to enhance octopamine release. That a modulator may produce the same changes in the soma and the terminals of a neuron is shown by the example of neurons in rat superior cervical ganglion where application of norepinephrine reduces both neurotransmission and somatic Ca2+ currents (Stephens and Mochida 2005Go). Another important principle involved in the control of transmitter release is the regulation of the BK channel-mediated fAHP. On the one hand, the fAHP can limit transmitter release (Raffaelli et al. 2004Go). On the other hand, it is necessary for an effective repolarization of action potentials, thus supporting the terminal's capacity to produce full-size action potentials on repetitive firing (e.g., Sausbier et al. 2004Go). In the case of the DUM neurons, it seems possible that the downregulation of Na+ current and KNa current is necessary to limit the effect of Ca2+ and KCa current modulation on hyperpolarization.

Conclusion

We have demonstrated a considerable complexity behind the "simple" increase in a neuron's firing rate induced by a hormone. Such complexity may be of functional relevance in the terminals of the neuron. Alternatively—or additionally—it could enable the neuron to respond differently to a (more or less) constant synaptic input. Such kind of mechanism was, for example, found for dopamine modulation in lobster stomatogastric neurons (Harris-Warrick et al. 1995Go). It might also be involved in the adaptation of insect neurons to different behavioral situations as observed, for example, in locust thoracic DUM neurons - ranging from patterned activity to complete inhibitory silencing (Pflüger 1999Go).


    APPENDIX
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
Modeling spike activity of DUM neurons

We simulated a one-compartment model with each current being described by a Hodgkin-Huxley formalism (Yamada et al. 1998Go). Beside the leak current, there are 10 currents implemented: One Na+ current (INa), four Ca2+ currents, and 5 K+ currents. The Ca2+ currents are the voltage-independent background current ICa,back, the low-voltage-activated (LVS) current ICa LVA, and the two high-voltage-activated currents ICa P/Q (P/Q-type current) and ICa nP/Q (non-P/Q-type current = total HVA current –P/Q-type current). The K+ currents are the voltage-activated delayed rectifier IKDR and the A-type current IKA, the voltage-independent Na+-activated current IKNa, and the Ca2+-activated current composed of the transient IKCa,t and the sustained IKCa,s. In general, each current I(t) is modeled according to the following equation [G(t) is the conductance, Veq the equilibrium potential, Vm(t) the membrane potential]

(A1)
with

(A2)
The exponents for m (mexp) and h (hexp) ranged from 0 to 4 and are specified for each of the current. m(t) and h(t) are described by the following first-order kinetics

(A3)
and

(A4)
The voltage dependence of the steady-state values for m and h are calculated according to the following equations

(A5)
and

(A6)
The time constants for m and h were also modeled as voltage dependent with the following equations

(A7)
and

(A8)
This description results in a total of 18 parameters specifying each of the currents. These parameters given in Table 1 were determined by describing the various currents in terms of Hodgkin-Huxley models (Fig. A1). The currents to be fitted according to the above equations were separated as outlined in METHODS. For technical reasons, i.e., to get an adequate voltage clamp, the currents were measured using solutions with reduced concentration of permeating ions. For example, the Na+ current was measured with a bath solution containing 60 mM Na+ (instead of 190 mM in standard bath solution). In such cases, the maximum conductance Gmax was corrected, the data given in Table 1 correspond to standard conditions (solutions as used for current-clamp measurements).

The three ion-activated currents, i.e., IKNa and both the transient and the sustained IKCa, deviated from the preceding description in the following aspects.

First, to model IKNa, the actual Na+-inward current carried by INa was calculated and multiplied with an empirical determined scaling factor called sensitivity factor. This approximation appears to be justified at least for the KNa peak current, which was found to be proportional to the Na current (Fig. 5D). Moreover, using this approximation, we were able to simulate a KNa current with amplitude and kinetic properties similar to the current obtained in the cells. The Na+-inward current times the sensitivity factor was used in Eqs. A5 and A6 instead of the membrane voltage to calculate m and h, respectively. Because the time constants of the KNa current are virtually voltage independent (Fig. A1, E2 and E3), they were fixed, i.e., Eqs. A7 and A8 did not apply. However, the current had a particular delay of ~2 ms with respect to the Na-inward current. The sensitivity factor and the delay were adjusted to reproduce time course and size of the modeled current according to the measured current at different voltages (for size, cf. Fig. 5D).

Second, to model IKCa,t and IKCa,s, we had to make some simplifying assumptions. We first considered these currents purely voltage dependent. To separate the currents, we modeled KCa,t with a m2h kinetics and subtracted it from the total KCa current. Under these conditions, we determined the voltage dependence of parameters as given in Table 1 and Fig. A1, F1–F3. These parameters, of course, reflect a "mean" Ca2+ supply through voltage-gated Ca2+ channels. To link the Ca2+ to the KCa current, we determined the effect of increase and decrease of the Ca2+ current (using 10 and 2 mM Ca2+ in the bath solution, respectively) on the KCa current. With this information, we calculated a sensitivity factor linking the amplitudes of KCa,t and KCa,s to the Ca2+ current amplitude (Table 1). For modeling the KCa currents, we made a further simplifying assumption in that we calculated these currents from the actual Ca2+ current. According to the results shown in Fig. 2, we coupled the KCa currents to the P/Q-type Ca2+ current. Although this neglects the fact that also the L-type Ca2+ current provides Ca2+ for the activation of the KCa,s component (Fig. 2, B and D), it is more important for the purpose of this study to investigate the effect of P/Q-type current modulation on KCa currents and spiking. To reproduce the kinetics and voltage dependence of the KCa currents measured in neurons by the modeled currents, we had to introduce a delay (~1.5 ms) between Ca2+ current and KCa currents.

The background Ca2+ current, ICa,back, was modeled with a constant, voltage- and ion-independent conductance (cf. Wicher et al. 2004Go).

The total number of parameters including leak conductance, leak potential and membrane capacitance amounted to 156. The simulation could be run in either the current- or voltage-clamp mode at a variable time-resolution ranging from 1 to 50 µs. With a maximum number of 20,000 time steps, this corresponded to a total simulation time ranging from 20 ms to 1 s. All parameters could be saved in a parameter file in ASCII format, which allowed for easy editing.

In voltage-clamp, the total current, i.e., the sum of all voltage- and ion-activated currents, was calculated as leak-subtracted and without capacitive current

(A9)
In current clamp, the actual membrane voltage Vm at time t was updated from its previous value at time t-1 according to

(A10)
This equation can be derived from fact that the sum of all currents (active, leak, injected and capacitive) equal to 0



The simulation software was written in IDL programming language (RSI Research Systems, Boulder, CO) with a menu-driven GUI. The compiled code including the parameter list is available on request from the authors and can be run under the "IDL virtual machine" environment, which can be downloaded at no cost from the RSI Website (http://www.rsinc.com/idlvm/).


    GRANTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
This work was supported by Deutsche Forschungsgemeinschaft Grant Wi 1422/2–4, 5).


    FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Address for reprint requests and other correspondence: D. Wicher, Saxon Academy of Sciences, Dept. Neurohormones, Erbertstr. 1, 07743 Jena, Germany (E-mail: b6widi{at}pan.zoo.uni-jena.de)


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 REFERENCES
 
Baumann E and Penzlin H. Sequence analysis of neurohormone D, a neuropeptide of an insect, Periplaneta americana. Biomed Biochim Acta 43: K13–16, 1984.

Brezina V, Horn CC, Weiss KR, and Orekhova IV. Modeling neuromuscular modulation in Aplysia. III. Interaction of central motor commands and peripheral modulatory state for optimal behavior. J Neurophysiol 93: 1523–1556, 2005.[Abstract/Free Full Text]

Buchholtz F, Golowasch J, Epstein IR, and Marder E. Mathematical model of an identified stomatogastric ganglion neuron. J Neurophysiol 67: 332–340, 1992.[Abstract/Free Full Text]

Dale N. Experimentally derived model for the locomotor pattern generator in the Xenopus embryo. J Physiol 489: 489–510, 1995.[Abstract/Free Full Text]

Dale N and Kuenzi FM. Ion channels and the control of swimming in the Xenopus embryo. Prog Neurobiol 53: 729–756, 1997.[CrossRef][ISI][Medline]

Derst C, Messutat S, Walther C, Eckert M, Heinemann SH, and Wicher D. The large conductance Ca2+-activated potassium channel (pSlo) of the cockroach Periplaneta americana: structure, localisation in neurones and electrophysiology. Eur J Neurosci 17: 1197–1212, 2003.[CrossRef][ISI][Medline]

Dryer SE. Na+-activated K+