Presented here is a method for investigation of the roles of the Na+/K+ pump and persistent Na+ current in leech heart interneurons using dynamic clamp.
The Na+/K+ pump, often thought of as a background function in neuronal activity, contributes an outward current (Ipump) that responds to the internal concentration of Na+ ([Na+]i). In bursting neurons, such as those found in central pattern generator (CPG) neuronal networks that produce rhythmic movements, the [Na+]i and therefore the Ipump, can be expected to vary throughout the burst cycle. This responsiveness to electrical activity, combined with independence from membrane potential, endow Ipump with dynamical properties not common to channel-based currents (e.g., voltage- or transmitter-gated or leak channels). Moreover, in many neurons, the pump's activity is modulated by a variety of modulators, further expanding the potential role of Ipump in rhythmic bursting activity. This paper shows how to use a combination of modeling and dynamic clamp methods to determine how Ipump and its interaction with persistent Na+ current influence rhythmic activity in a CPG. Specifically, this paper will focus on a dynamic clamp protocol and computational modeling methods in heart interneurons of medicinal leeches.
Heartbeat in leeches is driven by a CPG consisting of 9 bilateral pairs of heart interneurons (HNs) distributed across as many mid-body segmental ganglia. At the core of the CPG are mutually inhibitory pairs of interneurons located in the 3rd and 4th segmental ganglia that form half-center oscillators (HCOs) (Figure 1A). These neurons continue to burst when synaptically isolated pharmacologically using bicuculline1. Others, such as the pair in the 7th segmental ganglia (the focus of this protocol), are also bursters, capable of producing bursting activity when synaptically isolated. They are not mutually connected and receive only descending input, and thus are easily isolated by severing the ganglion from the rest of the nerve cord. This independent bursting activity is sensitive to introduced leak current caused by penetration with sharp microelectrodes for recording but vigorously burst when recorded with loose patch methods1.
Both individual HN neurons and HN HCOs have been modeled (Hodgkin-Huxley-based single isopotential compartment models of HN neurons containing all experimentally identified voltage-gated and synaptic currents), and all the burst characteristics of the living system have been successfully captured2. Myomodulin, an endogenous neuropeptide in leeches, markedly decreases the period (T) of the burst rhythm of isolated HN neurons and HN HCOs. This modulator acts to increase h-current (hyperpolarization-activated inward current, Ih) and to decrease Ipump3. This observation led to the exploration of how Ipump interacts with Ih, and how their co-modulation contributes to the rhythmic activity of HN neurons. Activation of the pump by increasing [Na+]i (using the ionophore monensin) speeds the HN burst rhythm in both HN HCOs and isolated HN neurons4. This speed-up was dependent on Ih. When Ih was blocked (2 mM Cs+), the burst period was not altered by this method of pump activation; however, the burst duration (BD) was curtailed, and the interburst interval (IBI) increased in both HN HCOs and isolated HN neurons4.
For this protocol, all the currents of a living HN(7) neuron, including the pump current, Ipump, are incorporated in the HN model as follows:
(1)
where C is the membrane capacitance (in nF), V is the membrane potential (in V), t is time (in s). Detailed ionic current descriptions and equations have been described elsewhere2,4. The complete HN model neuron runs in real time (Figure 2). The software will be made available on GitHub upon publication and will be suitable to run on the digital signal processing board described in the Table of Materials. Here, the focus of enquiry is the Na+/K+ pump current (Ipump) and the voltage-gated currents contributing significant Na+ flux: a fast Na+ current (INa) and a persistent Na+ current (IP). The maximal conductances of these currents are respectively. The Na+/K+ pump exchanges three intracellular Na+ ions for two extracellular K+ ions, thus producing a net outward current. Importantly, it pumps 3 times as much Na+ out of the neuron as this current indicates, which is important for calculating the intracellular Na+ concentration.
The Na+/K+ pump current depends on intracellular Na+ concentrations and is expressed by the following sigmoidal function:
(2)
where [Na]i is the intracellular Na+ concentration, is the maximal Na+/K+ pump current, [Na]ih is the intracellular Na+ concentration for the half-activation of the Na+/K+ pump, and [Na]is the sensitivity of the Na+/K+ pump to [Na]i. [Na]i builds as a result of the Na+ influxes carried by IP and INa and is diminished by the Na+ efflux of the Na+/K+ pump. The contribution of Ih and ILeak to the total Na+ flux is small and is not considered in the real-time model.
(3)
where, v is the volume (~6.7 pL) of the intracellular Na+ reservoir, F is Faraday's constant, and the extracellular Na+ concentration is kept constant.
Voltage-gated and leak conductances have been differentiated-these respond to membrane potential – from the pump current, which is regulated by the calculated intracellular Na+ concentration ([Na+]i). [Na+]i is built up through Na+ entry via the fast Na+ current (INa) that produces action potentials (spikes) and the persistent Na+ current (IP) that provides the depolarization to support spiking. [Na+]i is, in turn, reduced by the action of the pump through the extrusion of Na+. Baseline living HN values of (5nS) and (150 nS) have been assumed, and we take account of any added dynamic clamp .
The goal of the protocol described here is to manipulate Ipump precisely and reversibly in real time to discover how it interacts with voltage-gated currents (persistent Na+ current in the current protocol) to control rhythmic bursting in single HNs. To accomplish this goal, dynamic clamp was used, which artificially introduces, upon command, a precise amount of any current that can be calculated as the model is running. This method has advantages over pharmacological manipulation of the pump, which affects the entire tissue, can have off-target effects that are often hard to reverse, and cannot be precisely manipulated. Dynamic clamp5,6 reads the voltage of a recorded neuron in real time (Figure 1B) and calculates and injects, in real time, the amount of any current based on model equations and the set values of any or . Similar methods can easily be applied to any neuron that can be recorded intracellularly. However, parameters will have to be rescaled to the neuron chosen, and the neuron should be isolated from synaptic inputs, e.g., pharmacologically.
NOTE: Invertebrate animal experimental subjects are not regulated by the NIH or Emory and Georgia State Universities. All measures were nevertheless taken to minimize the suffering of the leeches used in this work.
1. Prepare isolated ganglion 7 from the leech nerve cord
2. Identify and record leech heart interneurons with sharp microelectrodes
3. Build a real-time HN or another model neuron
4. Implement and vary dynamic clamp conductances/currents
Modeling with the addition of Ipump4 brought the experimental findings presented in the introduction section into sharper focus and began to explain the pump-assisted mechanism of bursting. The real-time model demonstrated here has been tuned ( and parameters chosen) so that it produces regular rhythmic activity falling within the bounds of normal activity as observed in experiments – f, IBI, BD, T - and continues to produce such activity when the myomodulin-modulated parameters (the maximal pump current) and (maximal conductance of h-current) are varied or co-varied in the model. The parameter values determined can be used as a benchmark or canonical set for modeling experiments. In these model instances, Ipump oscillates throughout the burst cycle as [Na+]i around a baseline level. Ipump contributes to burst termination during the burst phase, and the hyperpolarization it produces activates Ih during the IBI; notice the maximal level of Ih near burst initiation (Figure 2).
Although the real-time HN model has all implemented currents2,4 available for dynamic clamping, the focus here was on and , which are available for changes while the model is running in the dynamic clamp GUI (Figure 3). Dynamic clamp allows the experimenter to add (or subtract with a negative or ) any conductance or current into a neuron artificially that mimics the voltage and ionic dependence of a real conductance or current. Thus, it is possible to fully explore how a particular conductance/current interacts with the endogenous conductances/currents inside cells (Figure 1). The real-time HN model indicates that the persistent Na+ current (IP) in HN neurons contributes much of the Na+ entry strongly affecting [Na+]i (Figure 2) and thus, Ipump. Because IP is active at relatively negative membrane potentials, it opposes Ipump even during the IBI.
These observations indicate that it is instructive to explore interactions between and in isolated HN neurons with dynamic clamp as discussed previously8,9,10. These experiments (ongoing) are performed with sharp microelectrode recordings in single, synaptically isolated HN(7) neurons (seventh ganglion severed from the nerve cord). To date, these experiments show that robust bursting is restored in tonically active HN neurons (due to microelectrode penetration introduced leak) by co-addition of IP and Ipump with dynamic clamp (Figure 4). This is an important observation indicating that a bursting mechanism is available in these neurons (even when leak is compromised) that results from the interaction of Ipump and IP. Preliminary results indicate their strong complicated interaction, which can be explored in the model and experiments (Figure 5).
In conclusion, Ipump in response to periodic increases in [Na+]i during bursting activity contributes to the burst rhythm through burst termination (decreasing BD). The interaction of IP and Ipump constitutes a mechanism that is sufficient to support endogenous bursting activity; this mechanism can reinstate robust bursting in HN interneurons recorded intracellularly in ganglion 7. The interaction between IP and Ipump through [Na+]i affects the HN burst period non-monotonically and ensures robustness of autonomous bursting. These conclusions are in line with experiments and modeling in vertebrate systems11,12.
Figure 1: Leech heart interneuron electrical activity and implementation of Ipump and IP with dynamic clamp. (A) Normal bursting activity simultaneously recorded, extracellularly (top) and intracellularly (bottom), in a leech heartbeat HCO from a third ganglion, a schematic of the recorded neurons and their mutually inhibitory synaptic connections at right. (B) Dynamic clamp schematic when recording a HN(7) interneuron in an isolated ganglion 7; note there is no synaptic interaction between the two HN(7) interneurons. (Ci) Bursting in a leak-compromised HN(7) interneuron. (Cii) More robust bursting can be produced by adding dynamic clamp Ipump ( = 0.1 nA), which makes up for the microelectrode induced leak, but depresses excitability, and (1 nS), which increases excitability. Black dashed lines indicate baseline values. Abbreviations: HN = heart interneuron; HCO = half-center oscillator; Ipump = outward current; IP = persistent Na+ current; = maximal Na+/K+ pump current; = maximal conductance of the persistent Na+ current; Vm = membrane potential; [Na+]i = internal concentration of Na+. Please click here to view a larger version of this figure.
Figure 2: Single HN interneuron model showing traces for membrane potential (Vm), Ih, Ipump, [Na+]i, and IP. Outward hyperpolarizing currents are negative, and inward depolarizing currents are positive. Black dashed lines indicate baseline values. Abbreviations: HN = heart interneuron; Ipump = outward current; IP = persistent Na+ current; = maximal Na+/K+ pump current; Ih = hyperpolarization-activated inward current; = maximal conductance of the persistent Na+ current; = maximal conductance of the hyperpolarization-activated inward current; Vm = membrane potential; [Na+]i = internal concentration of Na+. Please click here to view a larger version of this figure.
Figure 3: Graphical user interface of real-time heart interneuron (HN) model and dynamic clamp implemented on a digital signal processing board. Upper Left: Red Math boxes are user-determined parameter boxes for the real-time model, whereas Blue Live boxes are user-determined parameter boxes used in the dynamic clamp. El = the reversal potential of the leak current; Gl = leak conductance; Gh = h-current maximal conductance; Gp = P current maximal conductance; GCaS = slow calcium current maximal conductance; PumpMax = pump maximal current; [GSyn2 maximal synaptic conductance to the respective neuron; ThreshSyn2 spike crossing threshold for mediating a synaptic potential – these used to make a hybrid (living/model) half-center oscillator not illustrated here.]. Lower Left for Dynamic Clamp. At the very left are 5 computed values of dynamic clamp variables: Ipump = pump current injected; Ih = h-current injected (not used here); IP = P current injected; NaI = calculated internal Na+ concentration; ENa = calculated sodium reversal potential. Lower Left for Dynamic Clamp. To the right of the computed variables are 6 user determined parameter boxes: GNa = assumed endogenous fast sodium maximal conductance use to calculate Na+ flux associated with action potentials; PumpMaxL = maximal pump current to be injected by dynamic clamp; Naih see equation (2); Gh = maximal conductance to determine h-current to be injected by dynamic clamp; Gp = assumed endogenous P current maximal conductance use to calculate Na+ flux associated with endogenous P current; GpinHNLive = maximal conductance to determine P current to be injected by dynamic clamp. Please click here to view a larger version of this figure.
Figure 4: Dynamic clamp analysis of independent HN(7) bursting. Upregulation of from (A) 4.0 nS to (B) 9.0 nS slows down the independent HN burst rhythm. Experimental traces show rhythmic bursting in isolated HN(7) neuron with dynamic clamp. The ranges of oscillation of [Na+]i and Vm increase with upregulated . Traces top to bottom: recorded Vm, injected Ipump, calculated [Na+]i, and injected IP. Black dashed lines indicate baseline values. Abbreviations: HN = heart interneuron; Ipump = outward current; IP = persistent Na+ current; = maximal Na+/K+ pump current; = maximal conductance of the persistent Na+ current; Vm = membrane potential; [Na+]i = internal concentration of Na+. Please click here to view a larger version of this figure.
Figure 5: Dynamic clamp analysis of independent HN(7) bursting. Upregulation of tends to decrease, followed by an increased HN burst period. In individual experiments (points connected by lines) using dynamic clamp, values were swept while was held constant. Colors represent different constant levels of added used in different experiments. Abbreviations: HN = heart interneuron; = maximal Na+/K+ pump current; = maximal conductance of the persistent Na+ current. Please click here to view a larger version of this figure.
Modeling, dynamic clamp, and the resulting analyses that they enable are useful techniques for exploring how individual and groups of ionic conductance/currents contribute to the electrical activity of neurons (Figure 1, Figure 2,Figure 4, and Figure 5). The use of these techniques shows how the Na+/K+ pump current (Ipump) interacts with voltage-gated currents, particularly the persistent Na+ current (IP), to promote robust bursting in the leech heartbeat pattern generator's core HNs. By combining dynamic clamp experiments and modeling, it is possible to test models more directly than is possible with ordinary voltage recording and current clamp techniques. The results gathered from the dynamic clamp experiments (Figure 5) will be used to further refine the HN model. The basic method of dynamic clamping demonstrated here can be customized to reflect the properties of any neuron under study if a mathematical model of neuronal currents can be determined with voltage clamp experiments.
Successful completion of the experiments of the type shown here requires careful impalement of an HN or other neuron when using a sharp microelectrode, because strong bursting is curtailed by electrode penetration1. (Whole-cell patch recording techniques, which minimize introduced leak, are also applicable to other neurons, but do not work well on leech neurons.) It is critical that the impalement of the HN neuron causes minimal damage to the neuron (added leak), and input resistance should be monitored and must be in the range of 60-100 MOhms for successful experiments4.
Dynamic clamp is a powerful technique, but it has limitations imposed by neuronal geometry because the artificial conductances are implemented at the site of the recording electrode-usually the cell body-not at the site where rhythm-generating currents are usually localized5,6,10. In leech HN neurons, the cell body is electrically close to the integration zone (main neurite) of the neuron where most active currents are localized, and spikes are initiated.
The authors have nothing to disclose.
We thank Christian Erxleben for preliminary dynamic clamp experiments on HN(7) neurons that demonstrated their bursting capabilities. Angela Wenning aided the experiments with expert advice. We acknowledge the NIH for funding this work through Grant 1 R21 NS111355 to GSC and RLC.
ANIMALS | |||
Hirudo verbana | Leech.com, https://www.leech.com/collections/live-leeches | live leeches 2-3 grams | |
CHEMICALS | |||
ARTIFICIAL POND WATER | |||
CaCl2 | Sigma Aldrich | C5670-100G | 1.8 mM add last after adjusting pH |
glucose | Sigma Aldrich | G7021-100G | 10 mM |
HEPES | Sigma Aldrich | H4034-100G | 10 mM |
Instant Ocean (sea salt ) | Spectrum Brands Inc., Madison, WI | 0.05% (w/v) diluted in deionized water | |
KCl | Sigma Aldrich | P9333-500G | 4 mM |
NaCl | Sigma Aldrich | S7653-250G | 115 mM |
NaOH 0.1 N Solution | Sigma Aldrich | 2105-50ML | Adjust to pH 7.4 with NaOH |
MICROELECTRODES | |||
K Acetate | Sigma Aldrich | P1190-100G | 2 M |
KCl | Sigma Aldrich | P9333-500G | 20 mM |
SALINE | |||
EQUIPMENT | |||
#5 Forceps | Fine Science Tools Dumont | 11251-30 OR 11251-20 | For general leech dissection |
AxoClamp 2A/2B DCC electrometer | Axon Instruments Molecular Devices | 2A/2B | For recording of neuronal membrane potential and discontinuous current clamp |
Black resin | Dow Sylguard | 170 | Lines general dissect dish |
Capilary glass 1 mm outer diameter, 0.75 mm inner diameter | A-M Systems | 615000 | For fabricating sharp microelectrodes |
Clear resin | Dow Sylguard | 184 | Lines Petri dish used to mount ganglion for electrophysilogy |
Dark field condenser | Nikon | Dry 0.95-0.80 MBL 1210 | For illuminating the ganglion preparation during cell impalement |
Digidata 1440A | Axon CNS Molecular Devices | 1440A | Performs A to D and D to A for data acquisition and stimulation during electrophysiology |
Digital signal processing board | dSpace | CLP1104 | Our software implements all the conductances/currents in our model HN neuron on a DS1103 dSPACE PPC Controller Board in real-time at a rate of 20 kHz with a ControlDesk GUI (dSPACE, Paderborn, Germany)9. |
Falming/Brown Microelectrode Puller | Sutter Instruments | P-97 | For fabricating sharp microelectrodes |
Fiber-Lite high intensity illuminator | Dolan Jenner Industries | 170D | For illuminating the general dissection and for illuminating the ganglion preparation during cell impalement |
Headstage amplifier for AxoClamp 2A | Axon Instruments | HS-2A Gain:0.1LU | Now part of Molecular Devices for recording of neuronal membrane potential and discontinuous current clamp |
Light guide | Dolan Jenner Industries | Rev R 38 08 3729107 | For illuminating the general dissection and for illuminating the ganglion preparation during cell impalement |
Micromanipulator | Sutter Instruments | MPC-385 | Micromanipulator for cell impalement with microelectrodes |
Micromanipulator controller | Sutter Instruments | MPC-200 | Controls micromanipulators for cell impalement with microelectrodes |
Minuten pins | BioQuip | 0.15 mm diameter 1208SA | Should be shortened by curtting to ~5 mm |
Optical Breadboard 3' x 5' x 8" | Newport | Obsolete | With the 4 pneumatic Isolators below used to construct a vibration free workspace for electrophysiology |
Oscilloscope | HAMEG Instruments | HM303-6 | To monitor electrode setteling during DCC |
Pascheff-Wolff spring scissors | Moria | Supplied by Fine Science Tools (Foster City, CA) catalog # 15371-92 | |
pClamp 9 Software | Axon Instruments | 9 | Now part of Moleculear Devices uses the Digidata 1440 for data acquisition and stimulation during electrophysiology |
Pneumatic Isolators 28" | Newport | Obsolete | With optical breadboard used to construct a vibration free workspace for electrophysiology |
Simulink / MATLAB software | MathWorks | 2006 (Obsolete) | Implements dynamic clamp on the digital signal processing board |
Stereomicroscope | Wild | M5A | 10x Eye Pieces used for dissecting the leech and removingand desheathing ganglia |
Steromicroscope | Wild | M5 | 20x Eye Pieces used in electrophysiologcal station to visualize neuron for microelectrode penetration |
Student Vannas Spring Scissors | Fine Science Tools | 91500-09 | For general leech dissection |