Parametric resonance and amplification in excitable membranes. The Hodgkin-Huxley model

It has previously been shown by different investigators that the excitable membrane shows a resonant sensitivity to periodic external perturbations, but its Q-factor is, as a rule, low. The present paper analyses the possible ways of increasing the membrane Q, using a model of the Hodgkin-Huxley type. It is found, in particular, that it can be increased considerably by modulating periodically the membrane capacitance or the activation and inactivation rate constants of ionic channels, with a frequency of about 2 fo (fo being the fundamental frequency of damped oscillations in the membrane), the extent of modulation not exceeding the critical value 2/Q. In this case, a significant parametric amplification of the membrane current takes place. If the modulation coefficient is above 2/Q, the membrane can display a parametric resonance that causes stable self-oscillations in the potential with a frequency approximately fo. The conditions for the realization of parametric amplification and resonance in biological membranes are discussed.     

Capacitance fluctuations causing channel noise reduction in stochastic Hodgkin-Huxley systems

Voltage-dependent ion channels determine the electric properties of axonal cell membranes. They not only allow the passage of ions through the cell membrane, but also contribute to an additional charging of the cell membrane resulting in the so-called capacitance loading. The switching of the channel gates between an open and a closed configuration is intrinsically related to the movement of gating charge within the cell membrane. At the beginning of an action potential, the transient gating current is opposite to the direction of the current of sodium ions through the membrane. Therefore, the excitability is expected to become reduced due to the influence of a gating current. Our stochastic Hodgkin-Huxley-like modeling takes into account both the channel noise-i.e. the fluctuations of the number of open ion channels-and the capacitance fluctuations that result from the dynamics of the gating charge. We investigate the spiking dynamics of membrane patches of a variable size and analyze the statistics of the spontaneous spiking. As a main result, we find that the gating currents yield a drastic reduction of the spontaneous spiking rate for sufficiently large ion channel clusters. Consequently, this demonstrates a prominent mechanism for channel noise reduction.     

Spike trains in a stochastic Hodgkin-Huxley system

We consider a standard Hodgkin-Huxley model neuron with a Gaussian white noise input current with drift parameter mu and variance parameter sigma(2). Partial differential equations of second order are obtained for the first two moments of the time taken to spike from (any) initial state, as functions of the initial values. The analytical theory for a 2-component (V,m) approximation is also considered. Let mu(c) (approximately 4.15) be the critical value of mu for firing when noise is absent. Large sample simulation results are obtained for mumu(c), for many values of sigma between 0 and 25. For the time to spike, the 2-component approximation is accurate for all sigma when mu=10, for sigma>7 when mu=5 and only when sigma>15 when mu=2. When mumu(c), most paths show similar behavior and the moments exhibit smoothly changing behavior as sigma increases. Thus there are a different number of regimes depending on the magnitude of mu relative to mu(c): one when mu is small and when mu is large; but three when mu is close to and above mu(c). Both for the Hodgkin-Huxley (HH) system and the 2-component approximation, and regardless of the value of mu, the CV tends to about 1.3 at the largest value (25) of sigma considered. We also discuss in detail the problem of determining the interspike interval and give an accurate method for estimating this random variable by decomposing the interval into stochastic and almost deterministic components.     

Ionic channel density of excitable membranes can act as a bifurcation parameter

As the maximal K⁺-conductance (or K⁺-channel density) of the Hodgkin-Huxley equations is reduced, the stable resting membrane potential bifurcates at a subcritical Hopf bifurcation into small amplitude unstable oscillations. These small amplitude solutions jump to large amplitude periodic solutions that correspond to a repetitive discharge of action potentials. Thus the specific channel density can act as a bifurcation parameter, and can control the excitability and autorhythmicity of excitable membranes.     

Application of the Hodgkin-Huxley equations to an electric field model for interaction between excitable cells

We previously described a model for the electrical transfer of excitation from one cell to the next which utilized the electric potential generated in the junctional cleft between the cells. Low-resistance connections between the cells were not used in the model, and it was assumed that the junctional membranes were excitable. This model was analyzed for the static case without capacitances and for the dynamic case in which capacitances were part of the circuit elements. For simplicity, the Na⁺ resistance (R_Na), after a threshold potential was exceeded, was allowed to decrease exponentially (to 1% of its initial value) within 0·25–1·0 ms, and possible changes in the K⁺ resistance were ignored. In this paper, we have incorporated the Hodgkin-Huxley equations into the operation of the lumped membrane units for the electrical equivalent circuit of the cell membrane. The parameters varied are the membrane capacitances, resistances, maximum Na⁺ conductance ($\text{g}{}_{\mathrm{<mtext>Na</mtext>}}$), and the radial cleft resistance (R_jc). We demonstrated that our model worked very well, i.e. the successful transfer of action potentials was achieved, with the membrane units following Hodgkin-Huxley dynamics for changes in g_Na and g_K. The calculations indicate that transmission is facilitated when the junctional units have a higher $\text{g}{}_{\mathrm{<mtext>Na</mtext>}}$ and a lower capacitance and when R_jc is elevated. Lowering the resistance of the junctional membrane units several fold, relative to the surface membrane units, also facilitated transmission; however, the absolute resistance of the junctional membrane was still well above the maximum value that would allow sufficient local-circuit current to flow to effect transmission. Thus, the electric field model provides an alternative means of cell-to-cell propagation between myocardial cells which is electrical in nature but does not require the presence of low-resistance connections between cells.     

A possible model for receptors in excitable membranes

A model is proposed for receptors in excitable membranes based on the following assumptions. The receptor site and the process it excites in the membrane are located close to each other. The change of the electrostatic potential in the neighbourhood of the receptor site on the adsorption of a molecule (or ion) influences a potential dependent process in the membrane, such as ion permeability, rate of enzymatic reactions, ion binding etc. A comment is also made about the connection between measured physiological activity of a molecule and its ?real” physical activity.     

Effect of pyromecaine on ionic currents in electrically excitable membranes

The effects of the local anesthetic pyromecaine on the action potential parameters of the guinea-pig heart ventricular cardiomyocytes as well as on Na and K ionic currents of the neurons in the rat spinal ganglions have been studied. Pyromecaine was shown to reduce preferentially the first derivative of the ascending phase of cardiomyocyte action potential suggesting the blocking action of the anesthetic on the fast sodium inward current. The experiments on the isolated neurons in the rat spinal ganglions have shown that interaction of pyromecaine with inactivated Na+ channels makes a considerable contribution to the blocking effect of anesthesia.     

On the mechanism of spiking and bursting in excitable cells

A mathematical model previously developed to explain beta-cell membrane potential oscillations has been modified to accommodate the external variation of K+, Na+ and Ca2+ concentrations. Our model, which is applicable to excitable cells, incorporates the barrier kinetics. Hodgkin-Huxley-type gating mechanism, and an electrogenic Na+-K+ pump. Numerical solutions of our model are in agreement with many of the experimental results reported in the literature on excitable cells.     

Electrically induced phase transitions via the dipole model in excitable membranes

Calculations based upon the membrane dipole model have been carried out and indicate the possibility of cooperative behavior during the membrane excitation process. An explicit expression is obtained for the critical electric field,E, which must be present to initiate the cooperative structural transition assumed to occur during membrane excitation. An hypothesis concerning the occurrence of two distinct phase transitions in the membrane resulting in the rapid influx of sodium ions and sodium ion inactivation, respectively, is presented.     

Noise effects on spike propagation in the stochastic Hodgkin-Huxley models

Effects of membrane current noise on spike propagation along a nerve fiber are studied. Additive current noise and channel noise are considered by using stochastic versions of the Hodgkin-Huxley model. The results of computer simulation show that the membrane noise causes considerable variation of the propagation time of a spike (thus changes in interspike intervals) for a small unmyelinated fiber of radius 0.1 approximately 1 micron.     

Mechanisms of use-dependent block of sodium channels in excitable membranes by local anesthetics

Many local anesthetics promote reduction in sodium current during repetitive stimulation of excitable membranes. Use-, frequency-, and voltage-dependent responses describe patterns of peak INa when pulse width, pulse frequency, and pulse amplitude are varied. Such responses can be viewed as reflecting voltage-sensitive shifts in equilibrium between conducting, unblocked channels and nonconducting, blocked channels. The modulated-receptor hypothesis postulates shifts in equilibrium as the result of a variable-affinity receptor and modified inactivation gate kinetics in drug-complexed channels. An alternative view considers drug blocking in the absence of these two features. We propose that drug binds to a constant-affinity channel receptor where receptor access is regulated by the channel gates. Specifically, we view channel binding sites as guarded by the channel gate conformation, so that unlike receptors where ligands have continuous access, blocking agent access is variable during the course of an action potential. During the course of an action potential, the m and h gates change conformation in response to transmembrane potential. Conducting channels with both gates open leave the binding site unguarded and thus accessible to drug, whereas nonconducting channels, with gates in the closed conformation, act to restrict drug access to unbound receptors and possibly to trap drug in drug-complexed channels. We develop analytical expressions characterizing guarded receptors as "apparently" variable-affinity binding sites and predicting shifts in "apparent" channel inactivation in the hyperpolarizing direction. These results were confirmed with computer simulations. Furthermore, these results are in quantitative agreement with recent investigations of lidocaine binding in cardiac sodium channels.     

One of the possible regimes of bursting activity in the Hodgkin-Huxley model

A mathematical model of the Hodgkin--Huxley neuron membrane that comprises a "fast" second-order system and two "slow" equations is considered. The criterion for a self-oscillatory solution similar to the bursting activity in pacemaker neurons is found. The results obtained agree well with the computations carried out for the initial model.     

Intermittency in the Hodgkin-Huxley model

We show that action potentials in the Hodgkin-Huxley neuron model result from a type I intermittency phenomenon that occurs in the proximity of a saddle-node bifurcation of limit cycles. For the Hodgkin-Huxley spatially extended model, describing propagation of action potential along axons, we show the existence of type I intermittency and a new type of chaotic intermittency, as well as space propagating regular and chaotic diffusion waves. Chaotic intermittency occurs in the transition from a turbulent regime to the resting regime of the transmembrane potential and is characterised by the existence of a sequence of action potential spikes occurring at irregular time intervals.