|
|
||||||||
Department of Biology, Emory University, Atlanta, Georgia
Submitted 11 April 2005; accepted in final form 29 July 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Because of its small number of component neurons (14 interneurons) and well-defined synaptic connectivity, the central pattern-generating network controlling leech heartbeat is useful for exploring how rhythmic networks function (Calabrese et al. 1995
; Thompson and Stent 1976b
,c
). Within this network, two bilateral pairs of oscillator interneurons in the third and fourth ganglia pace the heartbeat rhythm. Each pair is configured as a half-center oscillator where each neuron bursts in alternation with its contralateral homolog through reciprocally inhibitory synapses (Peterson 1983a
). In isolated ganglia, pairs of oscillator interneurons maintain alternating bursting and are considered elemental half-center oscillators (Masino and Calabrese 2002a
; Peterson 1983a
). When synaptically isolated, the oscillator interneuron can display endogenous bursting (Cymbalyuk et al. 2002
). Because the rhythmic activity and membrane currents of these oscillator interneurons are well characterized (for review, see Hill et al. 2001
), they are an excellent system for exploring how modulation of specific membrane currents shapes rhythmic activity.
The leech heartbeat CPG has been shown to be modulated by a variety of sensory, neurosecretory, and motor pathways (Arbas 1984
). Modulation by the endogenous leech peptide FMRFamide has been characterized (for a review, see Calabrese et al. 1995
). Bath application of FMRFamide to the third or fourth ganglion elicits an acceleration of rhythmic bursting in the oscillator heart interneurons (Simon et al. 1992
). Voltage-clamp studies found that FMRFamide causes a negative shift in the activation and inactivation of a K current (IK1) in the oscillator interneurons (Simon et al. 1992
), activation of a persistent Na current (Schmidt et al. 1995
), and a decrease in spike-mediated synaptic transmission (Simon et al. 1994
). Additionally, FMRFamide activates a slowly activating voltage-gated K current, IKF (Nadim and Calabrese 1997
). Modeling studies identified this slowly activating K current as the most important mechanism for FMRFamide-mediated acceleration of the leech heartbeat rhythm.
Another endogenous leech peptide, leech myomodulin (Wang et al. 1998
), has also been shown to accelerate the period of oscillator interneurons. When it [or its molluscan analog myomodulin A (Cropper et al. 1987b
), subsequently referred to as myomodulin] is bath applied (1 µM) to isolated ganglia, it causes a decrease in period (Masino and Calabrese 2002b
). In this study, we demonstrate that, in addition to decreasing period, myomodulin also increases intraburst spike frequency. In synaptically isolated interneurons, it produces similar changes in burst characteristics and increases the robustness of endogenous bursting. We used voltage-clamp and current-clamp techniques to explore which membrane currents are affected by myomodulin to modulate the bursting of oscillator interneurons. We found that myomodulin enhances the hyperpolarization-activated cation current (Ih) and inhibits the electrogenic Na/K pump. By separately blocking Ih and the Na/K pump and determining the effects on burst characteristics, we demonstrate that the changes in membrane properties elicited by myomodulin may contribute to the decrease in period and increase in intraburst spike frequencies.
| METHODS |
|---|
|
|
|---|
Leeches (Hirudo medicinalis) were supplied by Leeches USA (Westbury, NY) and Biopharm (distributed by Carolina Biological Supply, Burlington, NC) and maintained in artificial pond water at 15°C. Before dissection, animals were cold anesthetized (14°C). Ganglia from midbody segments 3 and 4 were removed, individually pinned in petri dishes lined with silicone elastomer (Sylgard, Dow Corning, Midland, MI), and covered with saline. Immediately before the experiment, the glial sheath was removed with microscissors or scalpel, then the preparation was superfused continuously with saline (bath volume about 0.5 ml; flow rate about 2 ml/min). All experiments were performed at room temperature.
Unless otherwise indicated, experiments were performed in normal saline, containing (in mM): 115 NaCl, 4 KCl, 1.8 CaCl2, 10 glucose, 10 HEPES buffer; pH was adjusted to 7.4 with NaOH. For 0 Ca2+ salines, an equimolar amount of MnCl2 replaced CaCl2 (Angstadt and Calabrese 1991
, 1989
). For 0 Na+ salines, an equimolar amount of N-methyl-D-glucamine replaced NaCl, and the pH was adjusted to 7.4 with HCl. Where indicated, 2 mM CsCl was added to the saline. When used, ouabain (Sigma, St. Louis, MO) was dissolved, 10 mg/ml, in 95% ethanol, and this solution was added to saline for a concentration of 100 µM ouabain, resulting in an ethanol concentration of 0.6%. For ethanol control conditions, 95% ethanol was added to saline for an ethanol concentration of 0.6%. Where indicated, 0.5 mM bicuculline methiodide (Sigma) was added to the saline. Molluscan myomodulin (SynPep, Dublin, CA) was dissolved to a concentration of 1 mM in distilled water and frozen in 50 µl aliquots. Immediately before use, each aliquot was thawed and diluted in saline to produce a final concentration of 1 µM.
Electrophysiology
For all recordings, borosilicate microelectrodes (1 mm OD, 0.75 mm ID) were used. For extracellular recordings, cell somata were guided by gentle suction into saline-filled electrodes whose tips were pulled to about 20 µm; methods are described in Masino and Calabrese (2002a)
. For intracellular recordings, electrodes were pulled to sharp tips, and filled with 4 M K-acetate, 20 mM KCl solution, resulting in electrode resistances of 15 to 30 M
. Where indicated, TEA electrodes were filled with 1.5 M tetraethyl ammonium acetate tetrahydrate, 1.5 M Cs-acetate, 1 M K-acetate, 20 mM KCl solution. To decrease capacitance, some electrodes were coated with silicone elastomer to within nearly 500 µm of the tip, whereas all intracellular electrode tips were dipped in dimethyl polysiloxane (Sigma).
Heart interneurons were identified by position in the ganglion and by their characteristic bursting pattern (Fig. 1A1). For all intracellular recordings, neurons were impaled in normal saline, and only neurons with input resistances >60 M
, as measured with a 0.1-nA pulse, were accepted for experiments. At the end of each experiment, the bath potential was recorded, and only those experiments where it was within ±5 mV were accepted.
|
Current-clamp and voltage-clamp experiments were performed using an Axoclamp-2A or Axoclamp-2B amplifier (Axon Instruments, Union City, CA) operating in discontinuous current-clamp or discontinuous single-electrode voltage-clamp mode with a sample rate of 2.52.8 kHz. The electrode potential was monitored to ensure that it settled during each sample cycle. Output bandwidth was 0.3 kHz. Voltage-clamp gain was 0.8 to 2.0 nA/mV. Extracellular signals were amplified with a differential AC amplifier (model 1700, A-M Systems, Carlsburg, WA). Clampex software (pClamp 9.0, Axon Instruments) generated current or voltage waveforms during current-clamp and voltage-clamp recordings, respectively, and digitized and stored all electrophysiology recordings. Time-constant fitting was performed with either Clampfit (pClamp 9.0, Axon Instruments) or scripts written for MATLAB (The MathWorks, Natick, MA); all other data analysis was performed using MATLAB.
Spike detection and burst discrimination were performed as described in Masino and Calabrese (2002a)
. Briefly, spikes were detected when the voltage crossed a threshold. To prevent double counts of one spike, a 20-ms refractory period was used. Bursts were defined as groups of five or more spikes with
800 ms between each group of spikes. Burst period was calculated as the time between the median spikes of consecutive bursts (see Fig. 1A). Burst duration was defined as the time from the first to last spike of each burst, and duty cycle was defined as the ratio of burst duration and period for each cycle.
Voltage-clamp protocol and analysis
For all voltage-clamp experiments, except where indicated, protocols were first run in the control saline, indicated for each experiment below, then myomodulin or control saline was applied for 34 min and the protocol was run again. The pretreatment recordings were made to ensure the preparations in the two treatment groups were not significantly different before treatment. In all voltage-clamp protocols, at least three times the stimulus duration elapsed between stimulus repetitions, during which time the holding potential was maintained.
Synaptic current. One neuron of a pair of oscillator interneurons was recorded extracellularly to detect spikes for off-line spike-triggered averaging of postsynaptic currents. The other neuron was recorded intracellularly and held at 50 mV in single-electrode voltage-clamp mode. Synaptic currents were recorded for about 1 min (seven to ten burst cycles). Synaptic currents were averaged according to their sequence within a burst by triggering from spikes in the presynaptic neuron. For those experiments where myomodulin was applied, to ensure myomodulin had taken effect, only experiments where the period decreased were analyzed for changes in synaptic potential (two experiments out of 14 in the myomodulin group were excluded). Experiments in which the clamped neuron was spiking too much to discriminate synaptic potentials were also excluded from analysis (two experiments out of 12 in the myomodulin group were excluded).
Ramp current. For all ramp protocols, 1.8 mM Mn2+ replaced Ca2+ in the saline, to synaptically isolate the cells. Where indicated, 2 mM Cs+ and 100 µM ouabain were added to the 0 Ca2+ saline. The membrane potential was held at 40 mV and then linearly ramped to 85 mV in one second, then linearly ramped back to the 40 mV holding potential in one second. Ramp conductances were calculated as the slope of the linear fit to the IV curve between 55 and 80 mV. Average holding current was calculated as the mean of 25 ms of holding current (uncontaminated by unclamped spiking) preceding the ramp. One experiment (out of nine) was eliminated from the control group of the 0 Ca2+, 1.8 mM Mn2+ ramps because the slope was >9 SDs from the mean during the control treatment indicating an unacceptable level of damage to the neuron.
H-current.
Cells were bathed for 34 min in 0 Ca2+, 1.8 mM Mn2+ saline. Membrane potential was held at 40 mV, stepped to 60 mV for 3,040 ms, and stepped back to the 40 mV holding potential. This was performed three times, with steps to 60, 80, and 100 mV. Ih amplitude was measured as the difference between the average current between 75 and 85 ms after pulse onset and the average current between 2,990 and 2,995 ms after pulse onset. This section of the current trace was fit to a single exponential to extract the time constant of activation of Ih. The conductance of Ih was calculated by dividing the amplitude of Ih by the driving force, calculated as the voltage of the corresponding voltage step minus the Ih reversal potential, Eh = 21 mV (Angstadt and Calabrese 1989
). In a separate set of experiments measuring the effect of ZD7288 (Tocris Cookson, Ellisville, MO) on Ih, the protocol was as above except membrane potential was stepped to between 50 and 100 mV in 10-mV increments. ZD7288 (0.1 mM) was applied, and the voltage-clamp protocol was run at 5 and 8 min of application.
Ca currents.
Cells were penetrated with TEA electrodes (containing 1 M K-acetate, 1.5 M Cs-acetate, 1.5 M tetraethylammonium acetate tetrahydrate, 20 mM KCl) to block potassium currents. Cells were bathed for 34 min in 0 Na+, 5 mM Ca2+ saline. The membrane potential was held at 70 mV, stepped to 50 mV for 2 s, and back to 70 mV. This procedure was performed four times, with steps to 50, 45, 40, and 35 mV. Leak current was subtracted on-line by summing four negative prepulses, each of equal duration and one fourth the amplitude of the following protocol pulse. Linearity of leak currents was visually monitored with Axoscope (pClamp 9.0, Axon Instruments). Peak ICaF amplitude was visually ascertained as the peak current within the first 200 ms of the pulse. Peak ICaS amplitude was visually ascertained as the peak current, after the ICaF peak, within the first 500 ms of the pulse. To measure ICaS inactivation time constant for the steps to 40 and 35 mV, the voltage-clamp current was fit by four exponentials using the LevenbergMarquat fitting procedure in Clampfit. To aid fit convergence to biologically relevant values, a standard seed value for each of the exponentials was used for all fits, obtained from an oscillator heart interneuron model (based on prior published data: Hill et al. 2001
). The seed was as follows: A1 = 30 nA,
1 = 24 ms; A2 = 30 nA,
2 = 26 ms; A3 = 0.6 nA,
3 = 28 ms; A4 = 0.6 nA,
4 = 1200 ms; offset = 0.2 nA, where the A values represent the coefficients of each exponential term and the
values represent the time constants of each exponential term; offset is a constant offset added to the fitting equation. Each trace was fit twice, varying the start and end time points of the fits by 510 ms. The fit yielding the lowest SE for the slowest time constant was chosen. In one experiment (myomodulin application), the ICaS inactivation time constant for the current elicited by a step to 40 mV could not be fit successfully with our methods. This data point was excluded from the data, but the peak ICaF and peak ICaS values for this trace, as well as the inactivation time constant for the subsequent step to 35 mV, were included. To avoid erroneous measurements caused by Ca-current escape arising from imperfect space clamp, peak ICaF measurements were excluded when conductance measured at 35 mV was less than that measured at 40 mV.
K currents. No pretreatment protocols were run. Cells were bathed for 34 min in 0 Ca2+, 0 Na+ saline with myomodulin or without (control). Membrane potential was held at 35 mV, stepped to 40 mV for 3 s, and stepped back to the 35-mV holding potential. This procedure was performed six times, with steps to between 40 and 100 mV in 10-mV increments. K-current amplitude was measured as the difference between the average current between 10 and 20 ms before pulse onset and the average current between 2,040 and 2,060 ms after pulse onset.
Statistics.
Paired and unpaired t-test were performed with either Microsoft Excel or SigmaStat (Systat Software). All other statistical tests, Wilcoxon signed-rank test, Pearsons correlation, and one-way and two-way repeated-measures (RM) ANOVAs, were performed with SigmaStat. Tukeys test was used for post hoc analysis of ANOVAs. Significance of P < 0.05 is indicated on graphs by asterisks (*). All error bars indicate SE of the mean [SE = SD/
]; all values are quoted as means ± SE.
| RESULTS |
|---|
|
|
|---|
Myomodulin decreases period and increases spike frequency in oscillator heart interneurons
In accordance with previous results (Masino and Calabrese 2002b
), we found that bath-applying 1 µM myomodulin to oscillator heart interneurons in isolated ganglia significantly decreased the period by 17.3 ± 0.5% and significantly increased the mean, maximum, and minimum intraburst spike frequencies (by 13.2 ± 0.7, 14.2 ± 0.6, and 22.8 ± 2.4%, respectively, n = 20; period compared using a Wilcoxon signed-rank test attributed to nonnormality, P < = 0.001. Mean, maximum, and minimum spike frequencies compared with paired t-test, P
0.001, P
0.001, P < 0.05, respectively) (Fig. 1). Myomodulin did not change the duty cycle, (paired t-test, P = 0.30). These changes in burst characteristics were the same in the oscillator interneurons of ganglia 3 and 4, so we do not distinguish between those interneurons for the subsequent experiments in this paper. During washout of myomodulin, the period returned to premyomodulin values; mean period in wash was 3.6% greater than premyomodulin values (P = 0.7, paired t-test). The spike frequency remained elevated above premyomodulin values (mean, 10.5%; maximum, 10.4%; minimum, 7.1% of premyomodulin values; P = 0.03, P = 0.01, P = 0.50, respectively, paired t-test). Although bursting can persist for many minutes with less-invasive extracellular recordings, cells deteriorate during intracellular recordings and voltage clamp, so wash was not examined for voltage-clamp studies.
Although on average myomodulin decreased the period, its effect on period varied from a 2.5% increase to a 42.3% decrease. Similarly, myomodulins effect on mean spike frequency ranged widely from a 12.0% decrease to a 40.2% increase. We found that there was a significant negative correlation between period change and initial period, (Pearsons R = 0.629; P < 0.01), and a significant positive correlation between spike frequency change and initial spike frequency (Pearsons R = 0.493; P < 0.05) (Fig. 2, A and B). The slower the initial period, the more the rhythm sped in the presence of myomodulin; and the slower the initial intraburst spike frequency, the more it increased in the presence of myomodulin. Therefore the large range of initial burst characteristics (period = [6.1 to 13.2 s], mean spike frequency = [6.8 to 15.3 Hz]) may explain the large range of myomodulin-induced changes in the burst properties.
|
A previous study using a model of the oscillator interneurons and their synaptic interactions with coordinating neurons (simulating the condition of an isolated ganglion) demonstrated that decreasing spike-mediated inhibitory synaptic conductance increased spike frequency and decreased the period of model neurons (Hill et al. 2001
). To investigate whether myomodulin affected spike-mediated inhibitory postsynaptic currents (IPSCs), we voltage-clamped one neuron of an oscillator interneuron pair at 50 mV while recording extracellularly from its contralateral homolog (Fig. 3A) to obtain spike-triggered average IPSCs (Fig. 3B). Under these conditions, the peak amplitude of IPSCs increased during the first few spikes of each burst, stabilized between the 11th and the 30th spikes, and waned slightly at the end of the burst (Fig. 3C). To assess the effect of myomodulin on IPSC amplitude, we averaged the postsynaptic current amplitudes from the 11th to the 30th spikes, when the IPSC amplitude was most stable. Myomodulin had no significant effect on this average postsynaptic current amplitude (Fig. 3D, unpaired t-test, P = 0.12). Moreover, the average change in IPSC amplitude for the control and myomodulin applications was the same, further demonstrating no effect of myomodulin (control: 52.2 ± 20.6 pA, myomodulin: 51.9 ± 6.1 pA, unpaired t-test, P = 0.99). To assess any changes in IPSC dynamics, we fit each IPSC, from the 11th to the 30th, to a double-exponential function and extracted the rise and fall time constants (
1 and
2, respectively). Myomodulin did not change either of the averaged time constants (myomodulin:
1 = 3.1 ± 0.5 ms,
2 = 7.5 ± 0.3 ms, n = 10; control:
1 = 2.1 ± 0.2 ms,
2 = 7.0 ± 0.6 ms, n = 9; unpaired t-test,
1: P = 0.07;
2: P = 0.98).
|
In isolated ganglia 3 and 4, when one oscillator heart interneuron is voltage-clamped, the unclamped contralateral interneuron continues to receive phasic inhibition from and provide inhibition to ipsilateral axons arising from coordinating heart interneurons in ganglia 1 and 2 (Peterson 1983a
). This inhibition is functionally weaker than that from the contralateral oscillator interneuron. Under these conditions, the unclamped neuron exhibits irregular bursting that is thought to arise as a result of endogenous burst properties in the oscillator interneuron and network interactions with the coordinating axons (Cymbalyuk et al. 2002
) (Fig. 4A1). During myomodulin application, however, the unclamped neuron always exhibited regular bursting (Fig. 4A2). To quantify this effect, we calculated the coefficient of variation of period and of burst duration. Both values were smaller in the presence of myomodulin (Fig. 4, B and C; control: n = 9, myomodulin, n = 10; unpaired t-test, P < 0.05). In addition to increasing burst robustness, myomodulin significantly decreased the period (control: 12.0 ± 1.2 s; myomodulin: 6.77 ± 1.0 s; unpaired t-test, P < 0.05). Although within every preparation myomodulin increased the mean spike frequency with respect to the pretreatment condition (mean increase of 21.7 ± 4.0%); there was no significant difference between the control and myomodulin mean spike frequencies (control: 9.6 ± 0.2 Hz; myomodulin: 11.5 ± 1.3 Hz; unpaired t-test, P = 0.18). The lack of a significant myomodulin-elicited increase in spike frequencies may be a result of the initial spike frequency of the control group (10.2 ± 0.4 Hz), which is higher than that of the myomodulin group (9.4 ± 1.0 Hz). Nonetheless, these experiments demonstrated that myomodulin can enhance the robustness of bursting in an oscillator interneuron whose contralateral homolog is silenced. Furthermore, myomodulin decreased period while the contralateral homolog was voltage-clamped, supporting the idea that myomodulin has effects on burst characteristics that are independent of spike-mediated synaptic transmission between the oscillator heart interneurons.
|
700 ms between groups of at least five spikes) were analyzed. Myomodulin application significantly decreased the period by 26.6 ± 6.1% (pre-myomodulin: 6.6 ± 0.6 s; myomodulin: 5.4 ± 0.4 s; wash 6.2 ± 0.6 s, n = 5, one-way RM ANOVA, P < 0.05) and significantly increased mean spike frequency by 36.1 ± 11.3% (pre-myomodulin: 7.2 ± 1.2 Hz; myomodulin: 9.5 ± 0.9 Hz; wash: 7.2 ± 0.8 Hz, one-way RM ANOVA, P < 0.05). After washout, the period did not fully return to pre-myomodulin values, but the spike frequency did recover (P < 0.05).
|
The preceding experiments indicate that myomodulin affects membrane currents other than synaptic current. To identify which currents, we synaptically isolated interneurons with 0 Ca2+ (substituted with 1.8 mM Mn2+) saline and applied hyperpolarizing voltage-clamp ramps (40 to 85 mV) (Fig. 6A). In the presence of myomodulin, the conductance, measured as the slope of the linear current (see METHODS), during the ramps was significantly larger than that in the control group (control: 25.2 ± 1.3 nS, n = 8; myomodulin: 30.7 ± 2.3 nS, n = 9, unpaired t-test, P < 0.05), whereas the holding current was the same (control: 100.2 ± 31.7 pA; myomodulin: 229.3 ± 67.3 pA, unpaired t-test, P = 0.10). For further analysis of the myomodulin-elicited conductance, the averaged current waveforms from the control group were subtracted from those of the myomodulin group to obtain the "average difference current " (Fig. 6, B and C). Linear fits of the difference current revealed that the descending portion of the difference current had a slope of 5.5 nS with a reversal potential of 69.5 mV; the ascending portion corresponded to a slope of 3.6 nS with a reversal potential of 53.2 mV. The hysteresis in the current, indicated by the inward rectification during the down ramp, but the linearity during the up ramp, suggested a hyperpolarization activated conductance with slow activation and deactivation. Thus an increase in leak conductance seemed unlikely. It seems most likely that the myomodulin-elicited conductance was an inward-rectifier K current (IKIR) and/or a hyperpolarization-activated cation current (Ih).
|
|
|
To test whether the myomodulin enhancement of Ih could account fully for the myomodulin-elicited conductance revealed during voltage ramps, we repeated the ramp protocol, blocking Ih with 2 mM CsCl (Fig. 9A). There was no difference between the conductance of the myomodulin and control groups (control: 14.7 ± 1.9 nS, n = 9; myomodulin: 16.6 ± 1.7 nS, n = 10; unpaired t-test, P = 0.44), indicating that the myomodulin-elicited conductance revealed by the ramp protocol was indeed blocked by Cs+ and could indeed be attributed to an increase in Ih. However, there was a near constant offset between the averaged ramp current of the control and myomodulin groups, although this offset was not indicated as a statistically significant difference in holding currents (control: 100.6 ± 8.3 pA; myomodulin: 42.3 ± 7.0 pA), nor in the current at the trough of the ramp (control: 512.1 ± 19.9 pA; myomodulin: 639.8 ± 27.1 pA; two-way RM ANOVA, P = 0.09). Nonetheless, this offset indicated a further potential effect of myomodulin on interneurons. A constant offset in voltage-clamp current (not accompanied by a commensurate difference in conductance) that is present throughout all voltages could be explained by a myomodulin-induced inhibition of the Na/K pump. The absence of this offset during the ramps in 0 Ca2+ saline without Cs+ may be attributable to potential enhancement of the Na/K pump by Cs+; such enhancement by Cs+ has been observed in other leech neurons (Catarsi and Brunelli 1991
; Skou 1960
, 1965
). The effect of myomodulin on the Na/K pump may be more prominent when the pump is enhanced by Cs+.
|
Myomodulin affects the period and spike frequency in Cs+
Our results so far show that myomodulin increases Ih conductance and, in saline with 2 mM Cs+ saline (Cs+ saline), likely inhibits the Na/K pump. To determine whether myomodulin might change burst characteristics when Ih is blocked by Cs+, we recorded from oscillator interneurons extracellularly in saline, in Cs+ saline, and in Cs+ saline with myomodulin (Fig. 10A). When the superfusion was switched from saline to Cs+ saline, the period significantly increased by 23.9 ± 2.5% (n = 6; paired t-test, P < 0.001), in accordance with previous studies (Masino and Calabrese 2002b
), and the spike frequency significantly decreased by 8.7 ± 2.8% (paired t-test, P < 0.05). If increasing Ih is the only mechanism by which myomodulin changes burst properties, no further change should occur when myomodulin is added to Cs+ saline. However, application of myomodulin significantly decreased the period and increased the spike frequency (Fig. 10B, paired t-test, period: P < 0.01, spike frequency: P < 0.001). These changes demonstrate that myomodulin also affects membrane properties other than Ih. Moreover, the modulation of these other membrane properties decreases the period and increases the spike frequency. The current offset revealed by the ramp protocol in 0 Ca2+, 2 mM Cs+ saline and the subsequent disappearance of this offset in 0 Ca2+, 2 mM Cs+ saline with 100 µM ouabain suggests that, in the presence of Cs+, myomodulin inhibits the Na/K pump and thus changes the period and spike frequency.
|
Ouabain stops bursting in oscillator heart interneurons
To characterize how blocking the Na/K pump affects burst characteristics in oscillator heart interneurons, we recorded extracellularly from interneurons first in saline, then in saline plus 0.6% ethanol (the solvent vehicle for ouabain) for 4 min, followed by application of 100 µM ouabain (in 0.6% ethanol). No change in burst characteristics was seen with ethanol vehicle alone, indicating it is an innocuous vehicle for ouabain application (period change: 1.8 ± 3.1%, mean spike frequency change: 10.9 ± 5.9%, n = 7, one-way RM ANOVAs, period: P = 0.95, spike frequency: P = 0.16). Application of 100 µM ouabain resulted in the cessation of bursting within 3 min of application (n = 7). We compared burst characteristics in saline with vehicle alone to those during ouabain application (up to burst cessation). Blocking the Na/K pump with ouabain decreased period (15.9 ± 5.6%) and increased mean intraburst spike frequency (39.9 ± 6.7%) (one-way RM ANOVAs, period: P < 0.05, spike frequency: P < 0.001). These effects support our hypothesis that by inhibiting the Na/K pump, myomodulin decreases period and increases spike frequency.
Myomodulin does not alter low-threshold Ca currents
Although the results reported here are all consistent with the idea that myomodulin modulates burst characteristics by increasing Ih conductance and decreasing the Na/K pump activity, we cannot exclude other possible targets of modulation. For example, modeling studies indicate that decreasing the inactivation time constant of the slow, low-threshold Ca current (ICaS) in oscillator interneurons would decrease the period and increase spike frequency (Hill et al. 2001
). We measured the effect of myomodulin on Ca currents in 0 Na+, 5 mM Ca saline and found no change in the peak amplitudes of the slow (ICaS) or fast (ICaF) low-threshold Ca currents or the inactivation time constant of ICaS (
hCaS). Figure 11, two-way RM ANOVAs, peak ICaF: P = 0.064, peak ICaS: P = 0.52,
hCaS: P = 0.47).
|
| DISCUSSION |
|---|
|
|
|---|
Comparison of the modulation of Ih to other systems
Ih has been established as a regulator of many rhythmically active neurons (Bal and McCormick 1997
; Luthi and McCormick 1998
) and, as such, it is a common target for neuropeptide modulation (Harris-Warrick et al. 1995
; Kiehn and Harris-Warrick 1992
; Lee and Cox 2003
; Luthi and McCormick 1999
; Marder and Thirumalai 2002
). For example, in the lateral pyloric (LP) neuron of the rhythmic pyloric network of the crab stomatogastric nervous system, 5-HT was found to increase the conductance of Ih by shifting the voltage dependency of activation in a more positive direction. The resulting increase in Ih enhanced postinhibitory rebound, which phase-advanced the bursts of the LP neuron with respect to the other bursting neurons in the circuit and increased the LP burst spike frequency (Kiehn and Harris-Warrick 1992
). Whereas in many rhythmic systems the depolarizing effect of Ih causes an acceleration of rhythm, in the pre-Bötzinger complex in the mammalian brain, where respiratory rhythm is thought to be generated, blocking Ih with Cs+ or ZD 7288 increases respiratory frequency; the underlying mechanisms are as of yet unknown (Thoby-Brisson et al. 2000
). Thoby-Brisson and colleagues (2000)
found that the degree to which the frequency of the respiratory network changed by blocking Ih depended on the initial frequency, similar to our findings that the degree to which myomodulin changed period and spike frequency were dependent on their initial values (Fig. 2). Such state-dependent neuromodulation has been observed in a number of neuronal networks (Marder and Thirumalai 2002
; Nusbaum et al. 2001
).
The hyperpolarization-activated cyclic-nucleotide-gated (HCN) channels that mediate Ih have cyclic-nucleotidebinding domains (Robinson and Siegelbaum 2003
), and Ih is known to be sensitive to intracellular cyclic AMP (cAMP) (DiFrancesco 1993
; Luthi and McCormick 1998
, 1999
; Robinson and Siegelbaum 2003
). Moreover, the effects of myomodulin A in the Aplysia ARC neuromuscular system (where the peptide was first characterized; Cropper et al. 1987a
,b
) are mediated by cAMP (Brezina et al. 1994
; Hooper et al. 1994
). Thus the modulation we have observed here suggests that cAMP may mediate myomodulins effects.
Ih affects period and spike frequency in oscillator heart interneurons
Since its first characterization (Angstadt and Calabrese 1989
), Ih has been implicated as an important mechanism for burst generation and, potentially, regulation in leech oscillator heart interneurons (Calabrese 1998
; Calabrese and De Schutter 1992
; Hill et al. 2001
). Multiple generations of the single-compartment model of an oscillator interneuron (Hill et al. 2001
; Nadim et al. 1995
; Olsen et al. 1995
) have all demonstrated that changing the maximum conductance of Ih is an effective mechanism for modulating both the period and the spike frequency. Increasing the maximum conductance of Ih,
h, allows the inhibited neuron to escape earlier from synaptic inhibition and begin spiking, terminating the decelerating burst of the opposite neuron. Thus the cycle period is decreased and mean spike frequency is increased (Hill et al. 2001
). The same mechanism may contribute to the increase in both mean and minimum spike frequencies measured in the presence of myomodulin in the living neurons (Fig. 1B). Second, increasing
h augments depolarization of the initial part of the burst before Ih slowly deactivates. This increased depolarization increases the mean spike frequency in the model (Hill et al. 2001
) and may contribute to the increase in mean and maximum spike frequencies elicited by myomodulin in the living neurons (Fig. 1B). Studies of a hybrid half-center oscillator, consisting of a hardware model oscillator interneuron connected to a living oscillator interneuron, found that increasing Ih in either the hardware model or the living interneuron (through dynamic clamp; Prinz et al. 2004
; Sharp et al. 1993
) decreased the period (Sorensen et al. 2004
). These results support the conclusions reached here, by showing that our measured changes in Ih could cause the observed myomodulin-elicited changes of oscillator interneuron burst characteristics.
In an oscillator interneuron model tuned to exhibit endogenous bursting, Cymbalyuk et al. (2002)
found that increasing the maximal conductance of Ih expanded the parameter space in which endogenous bursting was possible, supporting our hypothesis that myomodulins increase of Ih may enhance the robustness of endogenous bursting in interneurons (Figs. 4 and 5).
Modulation of the Na/K pump in neurons
There is little evidence of neuromodulation of the Na/K pump in rhythmically active neurons. During pharmacologically induced rhythmic bursting, the Na/K pump has been implicated in generating rhythmic activity in neonatal rat spinal networks (Rozzo et al. 2002
), and rat midbrain dopaminergic neurons (Johnson et al. 1992
). In all examined neurons of the leech CNS, replacing extracellular Ca2+ with Co2+, Ni2+, or Mn2+ induced rhythmic oscillations that were disrupted by application of Na/K pump blockers (Angstadt and Friesen 1991
). Their computer simulations indicate that Na+ influx during a burst may activate the Na/K pump, thus causing hyperpolarization to end each burst. This mechanism is consistent with the finding that partial inhibition of the Na/K pump by low concentrations of ouabain (1025 µM) increased the cycle period of Co2+-induced oscillations (Angstadt and Friesen 1991
). However, in the normal cycling of oscillator heart interneurons, the period decreased as ouabain (100 µM) began to take effect, indicating the pump has a different role in influencing this bursting.
Modulation of the Na/K pump has been implicated in several nonrhythmic neurons, such as inhibition by DA1 receptors in cultured rat striatal neurons (Aizman et al. 2000
), and by serotonin in leech T sensory neurons (responsive to light touch) (Catarsi and Brunelli 1991
) mediated by an increase in cAMP (Catarsi et al. 1993
). Modulation of the Na/K pump by cyclic nucleotides has been demonstrated in many vertebrate tissues, including the PNS and CNS (for a review, see Therien and Blostein 2000
). Thus the modulation of the Na/K pump we have observed here again suggests that cAMP may mediate the effects of myomodulin.
Modulation of the Na/K pump in heart interneurons had not been previously observed before this study, and its functional effects have not yet been fully explored. If the myomodulin-induced inhibition of the pump is constant during myomodulin application, it would provide a voltage-independent depolarization of the neuron without an associated change in membrane conductance, similar to injecting positive current. Modeling this effect would be useful to determine whether the enhancement of the Na/K pump that we measured when blocking Ih with Cs+ might contribute to the myomodulin-induced change in burst properties measured in both Cs+ saline and normal saline. In leech T sensory neurons, Cs+ (4 mM, replacing extracellular K+) enhances the Na/K pump (Catarsi and Brunelli 1991
; Skou 1960
, 1965
), which may affect myomodulins inhibition of the Na/K pump. However, to measure the effect of myomodulins inhibition of the pump in normal saline requires an Ih blocker other than Cs+. The Ih blocker ZD7288 is ineffective in leech heart interneurons (Tobin, unpublished observations), and thus we have no way to separate, in normal saline, the effects of increasing Ih from those of decreasing the Na/K pump on heart interneuron burst characteristics.
Comparison of the modulatory effects of myomodulin to those of FMRFamide
Another well-characterized neuropeptide modulator of heart interneurons is FMRFamide (Evans et al. 1991
), which, like myomodulin, accelerates leech heartbeat (Simon et al. 1992
). FMRFamide has been shown to 1) cause a negative shift of activation and inactivation of IK1 (Simon et al. 1992
), 2) activate a persistent Na current (Schmidt et al. 1995
), 3) reduce spike-mediated synaptic transmission (Simon et al. 1994
), and 4) activate a slowly activating and deactivating K current, IKF (Nadim and Calabrese 1997
). Here we demonstrate that, although FMRFamide and myomodulin accelerate heart interneuron oscillations, they appear to act on different cellular targets. For example, we show that myomodulin increases Ih (Fig. 8), whereas FMRFamide has no effect (Angstadt and Calabrese 1989
), and myomodulin does not change spike-mediated transmission (Fig. 3E), whereas FMRFamide reduces it (Simon et al. 1994
). The changes that we have observed and that have been reported are small, however, and thus failure to observe an effect may not be definitive. We have yet to identify whether myomodulin affects K currents as FMRFamide does (Nadim and Calabrese 1997
; Simon et al. 1992
). Because these modulators both accelerate the heart rate, but apparently by different mechanisms, it has yet to be determined whether combined application of FMRFamide and myomodulin would have an effect different from that of applying either modulator alone, as has been described for muscles of the feeding apparatus in Aplysia (Brezina et al. 1996
; Nusbaum et al. 2001
) and the stomatogastric nervous system of crustaceans (Nusbaum et al. 2001
).
The neuropeptide myomodulin modulates burst properties in leech oscillator heart interneurons by modulating at least two membrane targets: increasing Ih and decreasing the Na/K pump. Although Ih has long been implicated as an important membrane current in producing the oscillatory activity of heart interneurons, this is the first demonstration of neuropeptide modulation of Ih in leech heart interneurons. Furthermore, our studies reveal a role for the Na/K pump as a modulating influence of rhythmic neuronal activity.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: R. L. Calabrese, Department of Biology, Emory University, 1510 Clifton Road N.E., Atlanta, GA 30322 (E-mail: rcalabre{at}biology.emory.edu)
| REFERENCES |
|---|
|
|
|---|
Angstadt JD and Calabrese RL. A hyperpolarization-activated inward current in heart interneurons of the medicinal leech. J Neurosci 9: 28462857, 1989.[Abstract]
Angstadt JD and Calabrese RL. Calcium currents and graded synaptic transmission between heart interneurons of the leech. J Neurosci 11: 746759, 1991.[Abstract]
Angstadt JD and Friesen WO. Synchronized oscillatory activity in leech neurons induced by calcium channel blockers. J Neurophysiol 66: 18581873, 1991.
Arbas EA. Rate modification in the heartbeat central pattern generator of the medicinal leech. J Comp Physiol A Sens Neural Behav Physiol 155: 783794, 1984.[CrossRef]
Bal T and McCormick DA. Synchronized oscillations in the inferior olive are controlled by the hyperpolarization-activated cation current I(h). J Neurophysiol 77: 31453156, 1997.
Baylor DA and Nicholls JG. After-effects of nerve impulses on signalling in the central nervous system of the leech. J Physiol 203: 571589, 1969.
Brezina V, Evans CG, and Weiss KR. Enhancement of Ca current in the accessory radula closer muscle of Aplysia californica by neuromodulators that potentiate its contractions. J Neurosci 14: 43934411, 1994.[Abstract]
Brezina V, Orekhova IV, and Weiss KR. Functional uncoupling of linked neurotransmitter effects by combinatorial convergence. Science 273: 806810, 1996.[Abstract]
Calabrese RL. Cellular, synaptic, network, and modulatory mechanisms involved in rhythm generation. Curr Opin Neurobiol 8: 710717, 1998.[CrossRef][Web of Science][Medline]
Calabrese RL and De Schutter E. Motor-pattern-generating networks in invertebrates: modeling our way toward understanding. Trends Neurosci 15: 439445, 1992.[CrossRef][Web of Science][Medline]
Calabrese RL, Nadim F, and Olsen OH. Heartbeat control in the medicinal leech: a model system for understanding the origin, coordination, and modulation of rhythmic motor patterns. J Neurobiol 27: 390402, 1995.[CrossRef][Web of Science][Medline]
Catarsi S and Brunelli M. Serotonin depresses the after-hyperpolarization through the inhibition of the Na+/K+ electrogenic pump in T sensory neurones of the leech. J Exp Biol 155: 261273, 1991.
Catarsi S, Scuri R, and Brunelli M. Cyclic AMP mediates inhibition of the Na(+)-K+ electrogenic pump by serotonin in tactile sensory neurones of the leech. J Physiol 462: 229242, 1993.
Cropper EC, Lloyd PE, Reed W, Tenenbaum R, Kupfermann I, and Weiss KR. Multiple neuropeptides in cholinergic motor neurons of Aplysia: evidence for modulation intrinsic to the motor circuit. Proc Natl Acad Sci USA 84: 34863490, 1987a.
Cropper EC, Tenenbaum R, Kolks MA, Kupfermann I, and Weiss KR. Myomodulin: a bioactive neuropeptide present in an identified cholinergic buccal motor neuron of Aplysia. Proc Natl Acad Sci USA 84: 54835486, 1987b.
Cymbalyuk GS, Gaudry Q, Masino MA, and Calabrese RL. Bursting in leech heart interneurons: cell-autonomous and network-based mechanisms. J Neurosci 22: 1058010592, 2002.
Dale N and Kuenzi F. Ionic currents, transmitters and models of motor pattern generators. Curr Opin Neurobiol 7: 790796, 1997.[CrossRef][Web of Science][Medline]
DiFrancesco D. Pacemaker mechanisms in cardiac tissue. Annu Rev Physiol 55: 455472, 1993.[CrossRef][Web of Science][Medline]
Evans BD, Pohl J, Kartsonis NA, and Calabrese RL. Identification of RFamide neuropeptides in the medicinal leech. Peptides 12: 897908, 1991.[CrossRef][Web of Science][Medline]
Gasparini S and DiFrancesco D. Action of the hyperpolarization-activated current (Ih) blocker ZD 7288 in hippocampal CA1 neurons. Pflügers Arch 435: 99106, 1997.[CrossRef][Web of Science][Medline]
Harris-Warrick RM, Coniglio LM, Levini RM, Gueron S, and Guckenheimer J. Dopamine modulation of two subthreshold currents produces phase shifts in activity of an identified motoneuron. J Neurophysiol 74: 14041420, 1995.
Harris-Warrick RM and Marder E. Modulation of neural networks for behavior. Annu Rev Neurosci 14: 3957, 1991.[CrossRef][Web of Science][Medline]
Hill AAV, Lu J, Masino MA, and Calabrese RL. A model of a segmental oscillator in the leech heartbeat neuronal network. J Comput Neurosci 10: 281302, 2001.[CrossRef][Web of Science][Medline]
Hooper SL, Probst WC, Cropper EC, Kupfermann I, and Weiss KR. Myomodulin application increases cAMP and activates cAMP-dependent protein kinase in the accessory radula closer muscle of Aplysia. Neurosci Lett 179: 167170, 1994.[CrossRef][Web of Science][Medline]
Johnson SW, Seutin V, and North RA. Burst firing in dopamine neurons induced by N-methyl-D-aspartate: role of electrogenic sodium pump. Science 258: 665667, 1992.
Katz PS. Intrinsic and extrinsic neuromodulation of motor circuits. Curr Opin Neurobiol 5: 799808, 1995.[CrossRef][Web of Science][Medline]
Katz PS. Neurons, networks, and motor behavior. Neuron 16: 245253, 1996.[CrossRef][Web of Science][Medline]
Kiehn O and Harris-Warrick RM. 5-HT modulation of hyperpolarization-activated inward current and calcium-dependent outward current in a crustacean motor neuron. J Neurophysiol 68: 496508, 1992.
Kiehn O, Johnson BR, and Raastad M. Plateau properties in mammalian spinal interneurons during transmitter-induced locomotor activity. Neuroscience 75: 263273, 1996.[CrossRef][Web of Science][Medline]
Lee SH and Cox CL. Vasoactive intestinal peptide selectively depolarizes thalamic relay neurons and attenuates intrathalamic rhythmic activity. J Neurophysiol 90: 12241234, 2003.
Luthi A and McCormick DA. H-current: properties of a neuronal and network pacemaker. Neuron 21: 912, 1998.[CrossRef][Web of Science][Medline]
Luthi A and McCormick DA. Modulation of a pacemaker current through Ca(2+)-induced stimulation of cAMP production. Nat Neurosci 2: 634641, 1999.[CrossRef][Web of Science][Medline]
Marder E and Calabrese RL. Principles of rhythmic motor pattern generation. Physiol Rev 76: 687717, 1996.
Marder E and Richards KS. Development of the peptidergic modulation of a rhythmic pattern generating network. Brain Res 848: 3544, 1999.[CrossRef][Web of Science][Medline]
Marder E and Thirumalai V. Cellular, synaptic and network effects of neuromodulation. Neural Networks 15: 479493, 2002.[CrossRef][Web of Science][Medline]
Masino MA and Calabrese RL. Phase relationships between segmentally organized oscillators in the leech heartbeat pattern generating network. J Neurophysiol 87: 15721585, 2002a.
Masino MA and Calabrese RL. Period differences between segmental oscillators produce intersegmental phase differences in the leech heartbeat timing network. J Neurophysiol 87: 16031615, 2002b.
Nadim F and Calabrese RL. A slow outward current activated by FMRFamide in heart interneurons of the medicinal leech. J Neurosci 17: 44614472, 1997.
Nadim F, Olsen OH, De Schutter E, and Calabrese RL. Modeling the leech heartbeat elemental oscillator. I. Interactions of intrinsic and synaptic currents. J Comput Neurosci 2: 215235, 1995.[CrossRef][Web of Science][Medline]
Nusbaum MP, Blitz DM, Swensen AM, Wood D, and Marder E. The roles of co-transmission in neural network modulation. Trends Neurosci 24: 146154, 2001.[CrossRef][Web of Science][Medline]
Olsen OH, Nadim F, and Calabrese RL. Modeling the leech heartbeat elemental oscillator. II. Exploring the parameter space. J Comput Neurosci 2: 237257, 1995.[CrossRef][Web of Science][Medline]
Opdyke CA and Calabrese RL. A persistent sodium current contributes to oscillatory activity in heart interneurons of the medicinal leech. J Comp Physiol A Sens Neural Behav Physiol 175: 781789, 1994.[Medline]
Peterson EL. Generation and coordination of heartbeat timing oscillation in the medicinal leech. I. Oscillation in isolated ganglia. J Neurophysiol 49: 611626, 1983a.
Prinz AA, Abbott LF, and Marder E. The dynamic clamp comes of age. Trends Neurosci 27: 218224, 2004.[CrossRef][Web of Science][Medline]
Robinson RB and Siegelbaum SA. Hyperpolarization-activated cation currents: from molecules to physiological function. Annu Rev Physiol 65: 453480, 2003.[CrossRef][Web of Science][Medline]
Rozzo A, Ballerini L, Abbate G, and Nistri A. Experimental and modeling studies of novel bursts induced by blocking Na+ pump and synaptic inhibition in the rat spinal cord. J Neurophysiol 88: 676691, 2002.
Schmidt J and Calabrese RL. Evidence that acetylcholine is an inhibitory transmitter of heart interneurons in the leech. J Exp Biol 171: 329347, 1992.
Schmidt J, Gramoll S, and Calabrese RL. Segment-specific effects of FMRFamide on membrane properties of heart interneurons in the leech. J Neurophysiol 74: 14851497, 1995.
Sharp AA, ONeil MB, Abbott LF, and Marder E. The dynamic clamp: artificial conductances in biological neurons. Trends Neurosci 16: 389394, 1993.[CrossRef][Web of Science][Medline]
Simon TW, Opdyke CA, and Calabrese RL. Modulatory effects of FMRF-NH2 on outward currents and oscillatory activity in heart interneurons of the medicinal leech. J Neurosci 12: 525537, 1992.[Abstract]
Simon TW, Schmidt J, and Calabrese RL. Modulation of high-threshold transmission between heart interneurons of the medicinal leech by FMRF-NH2. J Neurophysiol 71: 454466, 1994.
Skou JC. Further investigations on a (Mg2+ + Na+)-activated adenosine-triphosphatase possible related to the active, linked transport of Na+ and K+ across the nerve membrane. Biochim Biophys Acta 42: 623, 1960.[CrossRef]
Skou JC. Enzymatic basis for active transport of Na+ and K+ across cell membrane. Physiol Rev 45: 596617, 1965.
Sorensen M, DeWeerth S, Cymbalyuk G, and Calabrese RL. Using a hybrid neural system to reveal regulation of neuronal network activity by an intrinsic current. J Neurosci 24: 54275438, 2004.
Swensen AM and Marder E. Multiple peptides converge to activate the same voltage-dependent current in a central pattern-generating circuit. J Neurosci 20: 67526759, 2000.
Therien AG and Blostein R. Mechanisms of sodium pump regulation. Am J Physiol Cell Physiol 279: C541C566, 2000.
Thoby-Brisson M, Telgkamp P, and Ramirez JM. The role of the hyperpolarization-activated current in modulating rhythmic activity in the isolated respiratory network of mice. J Neurosci 20: 29943005, 2000.
Thompson WJ and Stent GS. Neuronal control of heartbeat in the medicinal leech. II. Intersegmental coordination of heart motor neuron activity by heart interneurons. J Comp Physiol 111: 281307, 1976b.[CrossRef]
Thompson WJ and Stent GS. Neuronal control of heartbeat in the medicinal leech. III. Synaptic relations of the heart interneurons. J Comp Physiol 111: 309333, 1976c.[CrossRef]
Wang Y, Price DA, and Sahley CL. Identification and characterization of a myomodulin-like peptide in the leech. Peptides 19: 487493, 1998.[CrossRef][Web of Science][Medline]
This article has been cited by other articles:
![]() |
A. V. Olypher and R. L. Calabrese Using Constraints on Neuronal Activity to Reveal Compensatory Changes in Neuronal Parameters J Neurophysiol, December 1, 2007; 98(6): 3749 - 3758. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. J. Norris, A. L. Weaver, A. Wenning, P. S. Garcia, and R. L. Calabrese A Central Pattern Generator Producing Alternative Outputs: Pattern, Strength, and Dynamics of Premotor Synaptic Input to Leech Heart Motor Neurons J Neurophysiol, November 1, 2007; 98(5): 2992 - 3005. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Wenning and E. P. Meyer Hemodynamics in the leech: blood flow in two hearts switching between two constriction patterns J. Exp. Biol., August 1, 2007; 210(15): 2627 - 2636. [Abstract] [Full Text] [PDF] |
||||
![]() |
A.-E. Tobin and R. L. Calabrese Endogenous and Half-Center Bursting in Morphologically Inspired Models of Leech Heart Interneurons J Neurophysiol, October 1, 2006; 96(4): 2089 - 2106. [Abstract] [Full Text] [PDF] |
||||
![]() |
A.-E. Tobin, S. D. Van Hooser, and R. L. Calabrese Creation and Reduction of a Morphologically Detailed Model of a Leech Heart Interneuron J Neurophysiol, October 1, 2006; 96(4): 2107 - 2120. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |