Null space in the Hodgkin-Huxley Equations. A critical test.

Voltage perturbation methods based upon topological concepts are used to elicit responses from the Hodgkin-Huxley (HH) nonlinear differential equations. These responses present a critical check upon the validity of the HH model for electrical activity across squid axon membrane. It is shown that when a constant current is applied such that a stable equilibrium and rhythmic firing are present, the following predictions are inherent in the HH system of equations: (a) Small instantaneous voltage perturbations to the axon given at points along its firing spike result in phase resetting curves (when new phase versus old phase is plotted) with an average slope of 1. (b) A larger voltage perturbation (from certain points along the firing spike) results in the permanent cessation of periodic firing, with membrane voltage rapidly approaching the equilibrium value. (c) A still larger perturbation yields phase resetting curves with an average slope equal to 0. These predictions, coupled with Tasaki's experimental demonstration that squid axons in excellent condition do give repetitive firing under constant current, provide a critical test of the validity of the HH model.     

Synapse‐dependent and independent mechanisms of thalamocortical axon branching are regulated by neuronal activity

Axon branching and synapse formation are critical processes for establishing precise circuit connectivity. These processes are tightly regulated by neural activity, but the relationship between them remains largely unclear. We use organotypic coculture preparations to examine the role of synapse formation in the activity‐dependent axon branching of thalamocortical (TC) projections. To visualize TC axons and their presynaptic sites, two plasmids encoding DsRed and EGFP‐tagged synaptophysin (SYP‐EGFP) were cotransfected into a small number of thalamic neurons. Time‐lapse imaging of individual TC axons showed that most branches emerged from SYP‐EGFP puncta, indicating that synapse formation precedes emergences of axonal branches. We also investigated the effects of neuronal activity on axon branching and synapse formation by manipulating spontaneous firing activity of thalamic cells. An inward rectifying potassium channel, Kir2.1, and a bacterial voltage‐gated sodium channel, NaChBac, were used to suppress and promote firing activity, respectively. We found suppressing neural activity reduced both axon branching and synapse formation. In contrast, increasing neural activity promoted only axonal branch formation. Time‐lapse imaging of NaChBac‐expressing cells further revealed that new branches frequently appeared from the locations other than SYP‐EGFP puncta, indicating that enhancing activity promotes axonal branch formation due to an increase of branch emergence at nonsynaptic sites. These results suggest that presynaptic locations are hotspots for branch emergence, and that frequent firing activity can shift branch emergence to a synapse‐independent process. © 2015 Wiley Periodicals, Inc. Develop Neurobiol 76: 323–336, 2016     

The Mechanism of Discharge Pattern Formation in Crayfish Interneurons

Excitatory and inhibitory processes which result in the generation of output impulses were analyzed in single crayfish interneurons by using intracellular recording and membrane polarizing techniques. Individual spikes which are initiated orthodromically in axon branches summate temporally and spatially to generate a main axon spike; temporally dispersed branch spikes often pace repetitive discharge of the main axon. Hyperpolarizing IPSP's sometimes suppress axonal discharge to most of these inputs, but in other cases may interact selectively with some of them. The IPSP's reverse their polarity at a hyperpolarized level of membrane potential; they sometimes exhibit two discrete time courses indicating two different input sources. Outward direct current at the main axon near branches causes repetitive discharges which may last, with optimal current intensities, for 1 to 15 seconds. The relation of discharge frequency to current intensity is linear for an early spike interval, but above 100 to 200 impulses/sec. it begins to show saturation. In one unit the current-frequency curve exhibited two linear portions, suggesting the presence of two spike-generating sites in the axon. Current threshold measurements, using test stimuli of different durations, showed that both accommodation and "early" or "residual" refractoriness contribute to the determination of discharge rate at different frequencies.     

A fully coupled transient excited state model for the sodium channel

Summary The axon membrane is simulated by standard Hodgkin-Huxley leakage and potassium channels plus a coupled transient excited state kinetic scheme for the sodium channel. This scheme for the sodium channel is as proposed previously by the author. Simultations are presented showing the form of the action potential, threshold behavior, accommodation, and repetitive firing. It is seen that the form of the individual action potential, its all-or-none nature, and its refractory period are well simulated by this model, as they are by the standard Hodgkin-Huxley model. However, the model differs markedly from the Hodgkin-Huxley model with respect to repetitive firing and accommodation to stimulating currents of slowly rising intensity, in ways that are anomn to be related to those features of the sodium inactivation which are anomalous to the H-H model. The tendency for repetitive firing is highly dependent on that parameter which primarily determintes the existence of the inactivation shift in voltage clamp experiments, in such a way that the more pronounced the inactivation shift, the less the tendency for repetitive firing,. The tendency for accommodation is highly dependent on that parameter which primarily determines the “τ_c − τ_h” separation, in such a way that the greater the separation the greater the tendency for the membrane to accommodate without firing action potentials to a slowly rising current.     

Neural repetitive firing: modifications of the Hodgkin-Huxley axon suggested by experimental results from crustacean axons.

The Hodgkin-Huxley equations for space-clamped squid axon (18 degrees C) have been modified to approximate voltage clamp data from repetitive-firing crustacean walking leg axons and activity in response to constant current stimulation has been computed. The m infinity and h infinity parameters of the sodium conductance system were shifted along the voltage axis in opposite directions so that their relative overlap was increased approximately 7 mV. Time constants tau m and tau h, were moved in a similar manner. Voltage-dependent parameters of delayed potassium conductance, n infinity and tau n, were shifted 4.3 mV in the positive direction and tau n was uniformly increased by a factor of 2. Leakage conductance and capacitance were unchanged. Repetitive activity of this modified circuit was qualitatively similar to that of the standard model. A fifth branch was added to the circuit representing a transient potassium conductance system present in the repetitive walking leg axons and in other repetitive neurons. This model, with various parameter choices, fired repetitively down to approximately 2 spikes/s and up to 350/s. The frequency vs. stimulus current plot could be fit well by a straight line over a decade of the low frequency range and the general appearance of the spike trains was similar to that of other repetitive neurons. Stimulus intensities were of the same order as those which produce repetitive activity in the standard Hodgkin-Huxley axon. The repetitive firing rate and first spike latency (utilization time) were found to be most strongly influenced by the inactivation time constant of the transient potassium conductance (tau b), the delayed potassium conductance (tau n), and the value of leakage conductance (gL). The model presents a mechanism by which stable low frequency discharge can be generated by millisecond-order membrane conductance changes.     

Digital Computer Solutions for Excitation and Propagation of the Nerve Impulse

The Hodgkin-Huxley model of the nerve axon describes excitation and propagation of the nerve impulse by means of a nonlinear partial differential equation. This equation relates the conservation of the electric current along the cablelike structure of the axon to the active processes represented by a system of three rate equations for the transport of ions through the nerve membrane. These equations have been integrated numerically with respect to both distance and time for boundary conditions corresponding to a finite length of squid axon stimulated intracellularly at its midpoint. Computations were made for the threshold strength-duration curve and for the repetitive firing of propagated impulses in response to a maintained stimulus. These results are compared with previous solutions for the space-clamped axon. The effect of temperature on the threshold intensity for a short stimulus and for rheobase was determined for a series of values of temperature. Other computations show that a highly unstable subthreshold propagating wave is initiated in principle by a just threshold stimulus; that the stability of the subthreshold wave can be enhanced by reducing the excitability of the axon as with an anesthetic agent, perhaps to the point where it might be observed experimentally; but that with a somewhat greater degree of narcotization, the axon gives only decrementally propagated impulses.     

Propagation of action potentials in squid giant axons. Repetitive firing at regions of membrane inhomogeneities

Effects of reduction in potassium conductance on impulse conduction were studied in squid giant axons. Internal perfusion of axons with tetraethylammonium (TEA) ions reduces G K and causes the duration of action potential to be increased up to 300 ms. This prolongation of action potentials does not change their conduction velocity. The shape of these propagating action potentials is similar to membrane action potentials in TEA. Axons with regions of differing membrane potassium conductances are obtained by perfusing the axon trunk and one of its two main branches with TEA after the second branch has been filled with normal perfusing solution. Although the latter is initially free of TEA, this ion diffuses in slowly. Up until a large amount of TEA has diffused into the second branch, action potentials in the two branches have very different durations. During this period, membrane regions with prolonged action potentials are a source of depolarizing current for the other, and repetitive activity may be initiated at transitional regions. After a single stimulus in either axon region, interactions between action potentials of different durations usually led to rebound, or a short burst, of action potentials. Complex interactions between two axon regions whose action potentials have different durations resembles electric activity recorded during some cardiac arrhythmias.     

Factors affecting slow regular firing in the suprachiasmatic nucleus in vitro

In isolated slices of hypothalamus, suprachiasmatic nucleus (SCN) neurons were recorded intracellularly. Blockade of Ca++ channels increased spike duration, eliminating an early component of the afterhyperpolarization (AHP) that followed evoked spikes. The duration and reversal potential of AHPs were, however, unaffected, suggesting that only an early, fast component of the AHP was Ca(++)-dependent. Unlike other central neurons that exhibit pacemaker activity, therefore, SCN neurons do not display a pronounced, long-lasting Ca(++)-dependent AHP. Extracellular Ba++ and intracellular Cs+ both revealed slow depolarizing potentials evoked either by depolarizing current injection, or by repolarization following large hyperpolarizations. They had different effects on the shape of spikes and the AHPs that followed them, however. Cs+, which blocks almost all K+ channels, dramatically reduced resting potential, greatly increased spike duration (to tens of milliseconds), and blocked AHPs completely. In contrast, Ba++ had little effect on resting potential and produced only a small increase in spike duration, depressing an early Ca(++)-dependent component and a later Ca(++)-independent component of the AHP. The relatively weak pacemaker activity of SCN neurons appears to involve voltage-dependent activation of at least one slowly inactivating inward current, which brings the cells to firing threshold and maintains tonic firing; both Ca(++)-dependent and Ca(++)-independent K+ channels, which repolarize cells after spikes and maintain interspike intervals; and Ca++ channels, which contribute to activation of Ca(++)-activated K+ currents and may also contribute to slow depolarizing potentials. In the absence of powerful synaptic inputs, SCN neurons express a pacemaker activity that is sufficient to maintain an impressively regular firing pattern. Slow, repetitive activation of optic input, however, increases local circuit activity to such an extent that the normal pacemaker potentials are overridden and firing patterns are altered. Since SCN neurons are very small and have large input resistances, they are particularly susceptible to synaptic input.     

Excitation properties of the squid axon membrane and model systems with current stimulation. Statistical evaluation and comparison.

The space-clamped squid axon membrane and two versions of the Hodgkin-Huxley model (the original, and a strongly adapting version) are subjected to a first order dynamic analysis. Stable, repetitive firing is induced by phase-locking nerve impulses to sinusoidal currents. The entrained impulses are then pulse position modulated by additional, small amplitude perturbation sinusoidal currents with respect to which the frequencies response of impulse density functions are measured. (Impulse density is defined as the number of impulses per unit time of an ensemble of membranes with each membrane subject to the same stimulus). Two categories of dynamic response are observed: one shows clear indications of a corner frequency, the other has the corner frequency obscured by dynamics associated with first order conductance perturbations in the interspike interval. The axon membrane responds with first order perturbations whereas the unmodified Hodgkin-Huxley model does not. Quantitative dynamic signatures suggest that the relaxation times of axonal recovery excitation variables are twice as long as those of the corresponding model variables. A number of other quantitative differences between axon and models, including the values of threshold stimuli are also observed.     

Phase resetting and bifurcation in the ventricular myocardium.

With the dynamic differential equations of Beeler, G. W., and H. Reuter (1977, J. Physiol. [Lond.]. 268:177-210), we have studied the oscillatory behavior of the ventricular muscle fiber stimulated by a depolarizing applied current I app. The dynamic solutions of BR equations revealed that as I app increases, a periodic repetitive spiking mode appears above the subthreshold I app, which transforms to a periodic spiking-bursting mode of oscillations, and finally to chaos near the suprathreshold I app (i.e., near the termination of the periodic state). Phase resetting and annihilation of repetitive firing in the ventricular myocardium were demonstrated by a brief current pulse of the proper magnitude applied at the proper phase. These phenomena were further examined by a bifurcation analysis. A bifurcation diagram constructed as a function of I app revealed the existence of a stable periodic solution for a certain range of current values. Two Hopf bifurcation points exist in the solution, one just above the lower periodic limit point and the other substantially below the upper periodic limit point. Between each periodic limit point and the Hopf bifurcation, the cell exhibited the coexistence of two different stable modes of operation; the oscillatory repetitive firing state and the time-independent steady state. As in the Hodgkin-Huxley case, there was a low amplitude unstable periodic state, which separates the domain of the stable periodic state from the stable steady state. Thus, in support of the dynamic perturbation methods, the bifurcation diagram of the BR equation predicts the region where instantaneous perturbations, such as brief current pulses, can send the stable repetitive rhythmic state into the stable steady state.     

Poisson process stimulation of an excitable membrane cable model.

The convergence of multiple inputs within a single-neuronal substrate is a common design feature of both peripheral and central nervous systems. Typically, the result of such convergence impinges upon an intracellularly contiguous axon, where it is encoded into a train of action potentials. The simplest representation of the result of convergence of multiple inputs is a Poisson process; a general representation of axonal excitability is the Hodgkin-Huxley/cable theory formalism. The present work addressed multiple input convergence upon an axon by applying Poisson process stimulation to the Hodgkin-Huxley axonal cable. The results showed that both absolute and relative refractory periods yielded in the axonal output a random but non-Poisson process. While smaller amplitude stimuli elicited a type of short-interval conditioning, larger amplitude stimuli elicited impulse trains approaching Poisson criteria except for the effects of refractoriness. These results were obtained for stimulus trains consisting of pulses of constant amplitude and constant or variable durations. By contrast, with or without stimulus pulse shape variability, the post-impulse conditional probability for impulse initiation in the steady-state was a Poisson-like process. For stimulus variability consisting of randomly smaller amplitudes or randomly longer durations, mean impulse frequency was attenuated or potentiated, respectively. Limitations and implications of these computations are discussed.     

Nonlinear cable equations for axons. II. Computations and experiments with external current electrodes

We have investigated the steady-state potential and current distributions resulting from current injection into a close-fitting channel into which a squid axon is placed. Hybrid computer solutions of the cable equations, using the Hodgkin-Huxley equations to give the membrane current density, were in good agreement with experimental observations. A much better fit was obtained when the Hodgkin-Huxley leakage conductance was reduced fivefold.     

Spike initiation and propagation on axons with slow inward currents

We investigate spike initiation and propagation in a model axon that has a slow regenerative conductance as well as the usual Hodgkin-Huxley type sodium and potassium conductances. We study the role of slow conductance in producing repetitive firing, compute the dispersion relation for an axon with an additional slow conductance, and show that under appropriate conditions such an axon can produce a traveling zone of secondary spike initiation. This study illustrates some of the complex dynamics shown by excitable membranes with fast and slow conductances.     

Firing properties of the soma and axon of the abdominal stretch receptor neurons in the crayfish (Astacus leptodactylus)

Action potentials (APs) and impulse responses in the soma and axon of the rapidly and slowly adapting (SA) abdominal stretch receptor neurons of the crayfish (Astacus leptodactylus) were recorded with single microelectrode current-clamp technique. Impulse frequency response to constant current injection was almost constant in the SA neuron while the response decayed completely in the rapidly adapting (RA) neuron. Mean impulse frequency responses to current stimulations were similar in the receptor neuron pairs. In the RA neuron additional current steps evoked additional impulses while a sudden drop in the current amplitude caused adaptation. Impulse duration was dependent on the rate of rise when current ramps were used. Adaptation was facilitated when calculated receptor current was used. Exposing the neuron to 3 mmol/l TEA or scorpion venom resulted in partly elongated impulse responses. SA neuron could continuously convert the current input into impulse frequency irrespective of previous stimulation conditions. Exposing the SA neuron to 3 mmol/l TEA or 1 mmol/l Lidocaine reduced impulse duration to large current stimulations. The SA neuron fired spontaneously if it was exposed to 5-10 mmol/l Lidocaine or 10(-2) mg/ml Leiurus quinquestriatus venom. The action potential (AP) amplitudes in the RA soma, RA axon, SA soma, and SA axon were significantly different between components of all pairs. Duration of the AP in the axon of the RA neuron was significantly shorter than those in the RA soma, SA soma, and SA axon. Diameter of the RA axon was larger than that of the SA axon. Non-adapting impulse responses were promptly observed only in the SA axons. The results indicate that the RA neuron is a sort of rate receptor transducing the rapid length changes in the receptor muscle while the SA neuron is capable of transducing the maintained length changes in the receptor muscle. The differences in firing properties mainly originate from the differences in the active and passive properties of the receptor neurons.     

Phase resetting of the rhythmic activity of embryonic heart cell aggregates. Experiment and theory.

Injection of a current pulse of brief duration into an aggregate of spontaneously beating chick embryonic heart cells resets the phase of the activity by either advancing or delaying the time of occurrence of the spontaneous beat subsequent to current injection. This effect depends upon the polarity, amplitude, and duration of the current pulse, as well as on the time of injection of the pulse. The transition from prolongation to shortening of the interbeat interval appears experimentally to be discontinuous for some stimulus conditions. These observations are analyzed by numerical investigation of a model of the ionic currents that underlie spontaneous activity in these preparations. The model consists of: Ix, which underlies the repolarization phase of the action potential, IK2, a time-dependent potassium ion pacemaker current, Ibg, a background or time-independent current, and INa, an inward sodium ion current that underlies the upstroke of the action potential. The steady state amplitude of the sum of these currents is an N-shaped function of potential. Slight shifts in the position of this current-voltage relation along the current axis can produce either one, two, or three intersections with the voltage axis. The number of these equilibrium points and the voltage dependence of INa contribute to apparent discontinuities of phase resetting. A current-voltage relation with three equilibrium points has a saddle point in the pacemaker voltage range. Certain combinations of current-pulse parameters and timing of injection can shift the state point near this saddle point and lead to an interbeat interval that is unbounded . Activation of INa is steeply voltage dependent. This results in apparently discontinuous phase resetting behavior for sufficiently large pulse amplitudes regardless of the number of equilibrium points. However, phase resetting is fundamentally a continuous function of the time of pulse injection for these conditions. These results demonstrate the ionic basis of phase resetting and provide a framework for topological analysis of this phenomenon in chick embryonic heart cell aggregates.     

Modelling the Effects of Electrical Coupling between Unmyelinated Axons of Brainstem Neurons Controlling Rhythmic Activity

Gap junctions between fine unmyelinated axons can electrically couple groups of brain neurons to synchronise firing and contribute to rhythmic activity. To explore the distribution and significance of electrical coupling, we modelled a well analysed, small population of brainstem neurons which drive swimming in young frog tadpoles. A passive network of 30 multicompartmental neurons with unmyelinated axons was used to infer that: axon-axon gap junctions close to the soma gave the best match to experimentally measured coupling coefficients; axon diameter had a strong influence on coupling; most neurons were coupled indirectly via the axons of other neurons. When active channels were added, gap junctions could make action potential propagation along the thin axons unreliable. Increased sodium and decreased potassium channel densities in the initial axon segment improved action potential propagation. Modelling suggested that the single spike firing to step current injection observed in whole-cell recordings is not a cellular property but a dynamic consequence of shunting resulting from electrical coupling. Without electrical coupling, firing of the population during depolarising current was unsynchronised; with coupling, the population showed synchronous recruitment and rhythmic firing. When activated instead by increasing levels of modelled sensory pathway input, the population without electrical coupling was recruited incrementally to unpatterned activity. However, when coupled, the population was recruited all-or-none at threshold into a rhythmic swimming pattern: the tadpole “decided” to swim. Modelling emphasises uncertainties about fine unmyelinated axon physiology but, when informed by biological data, makes general predictions about gap junctions: locations close to the soma; relatively small numbers; many indirect connections between neurons; cause of action potential propagation failure in fine axons; misleading alteration of intrinsic firing properties. Modelling also indicates that electrical coupling within a population can synchronize recruitment of neurons and their pacemaker firing during rhythmic activity.     

Bifurcation of periodic solutions of Hodgkin-Huxley model for the squid giant axon

The Hodgkin-Huxley model of the space-clamped squid giant axon is shown to admit unstable periodic solutions for current stimuli less than the stimulus at which the rest state becomes linearly unstable. The periodic solutions are demonstrated both by bifurcation theory and by numerical integration. The presence of subcritical unstable oscillations explains the discontinuous behaviour of the amplitude of the repetitive response as a function of current stimulus     

Nonlinear cable equations for axons. I. Computations and experiments with internal current injection

Steady-state potential and current distributions resulting from internal injection of current in the squid giant axon have been measured experimentally and also computed from nonlinear membrane cable equation models by numerical methods, using the Hodgkin-Huxley equations to give the membrane current density. The solutions obtained by this method satisfactorily reproduce experimental measurements of the steady-state distribution of membrane potential. Computations of the input current-voltage characteristic for a nonlinear cable were in excellent agreement with measurements on axons. Our results demonstrate the power of Cole's equation to extract the nonlinear membrane characteristics simply from measurement of the input resistance.     

Spontaneous Activity in Crustacean Neurons

Single units which discharged with regular spontaneous rhythms without intentional stimulation were observed in the ventral nerve cord by intracellular recording close to the sixth abdominal ganglion. These units were divided into two groups: group A units in which interspike intervals varied less than 10 msec.; group B units in which interspike intervals varied within a range of 10 to 30 msec. Group A units maintained "constant" interspike intervals and could not be discharged by sensory inputs, while the majority of group B units could be discharged by appropriate sensory nerve stimulation. Both group A and B units discharged to direct stimulation when the stimulating and recording electrodes were placed in the same ganglionic intersegment, and directly evoked single spikes reset the spontaneous rhythm. In group B units, presynaptic volleys reset the spontaneous rhythm of some units; but in others, synaptically evoked spikes were interpolated within the spontaneous rhythm without resetting. The phenomenon of enhancement could also be demonstrated in spontaneously active units as a result of repetitive stimulation. It is concluded that endogenous pacemaker activity is responsible for much of the regular spontaneous firing observed in crayfish central neurons, and that interaction of evoked responses with such pacemaker sites can produce a variety of effects dependent upon the anatomical relationships between pacemaker and synaptic regions.     

Suppression of pacemaker activity by rapid repetitive phase delay

Spontaneous activity of pacemaker cells or structures may be suppressed by rapid repetitive stimulation. Conditions are that the oscillator's phase reset curve, characterizing the phase resetting effect of single stimuli, has a phase delay part and that the interval between the stimuli falls within a range of values, determined by the form oo the phase reset curve. Under these conditions, which appeared the same as those for stable underdrive pacing, the pacemaker becomes stably entrained to the stimuli without firing, i.e. it is kept within a certian part of its limit cycle because the pulses repeatedly delay the next coming action potential. This rapid stimulation suppression of pacemaker activity is demonstrated experimentally on a simple electronic pacemaker cell model for two types of phase reset curves, a biphasic one for depolarizing and a monophasic one for hyperpolarizing pulses. Computer simulations of coupled pacemaker cells, interacting by phase reset curves, illustrate how this type of pacemaker suppression may protect a population of pacemaker cells like the sinus node in the heart against arrhythmias.Supported by the Netherlands Foundation for Medical Research FUNGOSupported by the Netherlands Organization for the Advancement of Pure Research ZWO