Oscillation and Repetitive Firing in Squid Axons : Comparison of experiments with computations

Space-clamped squid axons treated with low calcium and computed Hodgkin-Huxley (HH) axons were stimulated by steps of superthreshold current from 101 to 400% of the rheobasic value over a temperature range of 5–27°C. The natural frequency of sustained repetitive firing of real and computed axons depended weakly upon stimulus intensity and strongly upon temperature, with a Q₁₀ of 2.7 (experimental) and 2.6 (computed). For real axons, but not the computed axon, the intervals between the first two spikes were shorter than between subsequent spikes. Constant spike frequencies from 75 Hz at low intensities and temperatures to 330 Hz at high intensities and temperatures were soon achieved. Subthreshold and superthreshold responses were sometimes intermixed in a train of responses from a real axon responding to a constant step of current, but not predicted by HH. The time interval following a spike was always longer than that following a subthreshold oscillation in slightly decalcified real axons, as Huxley and FitzHugh also found for computed axons. There was a bias toward spikes at the beginning of the train and toward subthreshold responses later on. Some repeated patterns were found, every second, third, or fourth response being a spike. Neither the HH equations nor the computed or experimental threshold behaviors show a critical temperature to support a membrane phase transition.     

On the role of subthreshold dynamics in neuronal signaling

The role of subthreshold dynamics in neuronal signaling is examined using periodic pulse train stimulation of the Fitzhugh-Nagumo (FN) model of nerve membrane excitability and results from the squid giant axon as an experimental data base. For a broad range of stimulus conditions the first pulse in a pulse train elicited an action potential, whereas all subsequent pulses elicited subthreshold responses, both in the axon and in the FN model. These results are not well described by the Hodgkin and Huxley 1952 model. Various different patterns of subthreshold responses, including chaotic dynamics, can be observed in both systems-the FN model and the axon-depending upon stimulus conditions. For some conditions action potentials are occasionally interspersed among the subthreshold events with randomly occurring interspike intervals. The randomness is directly attributable to the underlying subthreshold chaos-deterministic chaos-rather than to a stochastic noise source. We conclude that this mechanism may contribute to multimodal interspike interval histograms which have been observed from individual neurons throughout the nervous system.     

White noise measurement of squid axon membrane impedance.

The use of white noise techniques for system identification is illustrated by the following characterization of the subthreshold membrane impedance of the squid giant axon, space-clamped in a double sucrose gap. Power spectra were also computed. Depolarization increases the resonance, shifts the resonant frequently upward and decreases the membrane's inductive reactance. Reduced external Ca⁺⁺ increases the resonance, shifts the resonant frequency downward and increases the inductive reactance.     

Spike latency and jitter of neuronal membrane patches with stochastic Hodgkin-Huxley channels

Ion channel stochasticity can influence the voltage dynamics of neuronal membrane, with stronger effects for smaller patches of membrane because of the correspondingly smaller number of channels. We examine this question with respect to first spike statistics in response to a periodic input of membrane patches including stochastic Hodgkin-Huxley channels, comparing these responses to spontaneous firing. Without noise, firing threshold of the model depends on frequency—a sinusoidal stimulus is subthreshold for low and high frequencies and suprathreshold for intermediate frequencies. When channel noise is added, a stimulus in the lower range of subthreshold frequencies can influence spike output, while high subthreshold frequencies remain subthreshold. Both input frequency and channel noise strength influence spike timing. Specifically, spike latency and jitter have distinct minima as a function of input frequency, showing a resonance like behavior. With either no input, or low frequency subthreshold input, or input in the low or high suprathreshold frequency range, channel noise reduces latency and jitter, with the strongest impact for the lowest input frequencies. In contrast, for an intermediate range of suprathreshold frequencies, where an optimal input gives a minimum latency, the noise effect reverses, and spike latency and jitter increase with channel noise. Thus, a resonant minimum of the spike response as a function of frequency becomes more pronounced with less noise. Spike latency and jitter also depend on the initial phase of the input, resulting in minimal latencies at an optimal phase, and depend on the membrane time constant, with a longer time constant broadening frequency tuning for minimal latency and jitter. Taken together, these results suggest how stochasticity of ion channels may influence spike timing and thus coding for neurons with functionally localized concentrations of channels, such as in “hot spots” of dendrites, spines or axons.     

Noise-Induced Coherence and Network Oscillations in a Reduced Bursting Model

The dynamics of the Hindmarsh-Rose (HR) model of bursting thalamic neurons is reduced to a system of two linear differential equations that retains the subthreshold resonance properties of the HR model. Introducing a reset mechanism after a threshold crossing, we turn this system into a resonant integrate-and-fire (RIF) model. Using Monte-Carlo simulations and mathematical analysis, we examine the effects of noise and the subthreshold dynamic properties of the RIF model on the occurrence of coherence resonance (CR). Synchronized burst firing occurs in a network of such model neurons with excitatory pulse-coupling. The coherence level of the network oscillations shows a stochastic resonance-like dependence on the noise level. Stochastic analysis of the equations shows that the slow recovery from the spike-induced inhibition is crucial in determining the frequencies of the CR and the subthreshold resonance in the original HR model. In this particular type of CR, the oscillation frequency strongly depends on the intrinsic time scales but changes little with the noise intensity. We give analytical quantities to describe this CR mechanism and illustrate its influence on the emerging network oscillations. We discuss the profound physiological roles this kind of CR may have in information processing in neurons possessing a subthreshold resonant frequency and in generating synchronized network oscillations with a frequency that is determined by intrinsic properties of the neurons. PACS 05.45.-a, 05.40.Ca, 87.18.Sn, 87.19     

Temperature Characteristics of Excitation in Space-Clamped Squid Axons

Temperature characteristics of excitability in the squid giant axon were measured for the space-clamped axon with the double sucrose gap technique. Threshold strength-duration curves were obtained for square wave current pulses from 10 µsec to 10 msec and at temperatures from 5°C to 35°C. The threshold change of potential, at which an action potential separated from a subthreshold response, averaged 17 mv at 20°C with a Q¹⁰ of 1.15. The average threshold current density at rheobase was 12 µa/cm² at 20°C with a Q¹⁰ of 2.35 compared to 2.3 obtained previously. At short times the threshold charge was 1.5·10^-8 coul/cm². This was relatively independent of temperature and occasionally showed a minimum in the temperature range. At intermediate times and all temperatures the threshold currents were less than for both the single time constant model and the two factor excitation process as developed by Hill. FitzHugh has made computer investigations of the effect of temperature on the excitation of the squid axon membrane as represented by the Hodgkin-Huxley equations. These are in general in good agreement with our experimental results.     

RECTIFICATION AND INDUCTANCE IN THE SQUID GIANT AXON

Previous measurements have shown that the electrical properties of the squid axon membrane are approximately equivalent to those of a circuit containing a capacity shunted by an inductance and a rectifier in series. Selective ion permeability of a membrane separating two electrolytes may be expected to give rise to the rectification. A quasi-crystalline piezoelectric structure of the membrane is a plausible explanation of the inductance. Some approximate calculations of behavior of an axon with these membrane characteristics have been made. Fair agreement is obtained with the observed constant current subthreshold potential and impedance during the foot of the action potential. In a simple case a formal analogy is found between the calculated membrane potential and the excitability defined by the two factor formulations of excitation. Several excitation phenomena may then be explained semi-quantitatively by further assuming the excitability proportional to the membrane potential. Some previous measurements and subthreshold potential and excitability observations are not consistent with the circuit considered and indicate that this circuit is only approximately equivalent to the membrane.     

Information transmission by isolated frog muscle spindle

Transmission of sensory information was calculated for the isolated frog muscle spindle receptor, using Shannon's information measure. Sinusoidal movements, random noise stretches, and sinusoids with superimposed auxiliary noise were applied as stimuli. In addition, the static prestretch level of the intrafusal muscle bundle was adjusted between resting length (L0) and L0 + 600 micron, so that the analysis of the information transmission properties covered the entire dynamic range of the sensory receptor organ. Sinusoidal stretches below 2 Hz evoked smoothly modulated cycle histograms, which were approximately linearly related to the stimulating sinewave. The transinformation rates under these conditions were generally low (5-17 bit X s-1), regardless of the amplitude of the applied movement. Increasing prestretch enhanced the modulation depth of the cycle histograms considerably, but increased the transinformation rates by less than 10 bit X s-1. By contrast, sinusoids above 2 Hz evoked clearly nonlinear cycle histograms, because each action potential was firmly phase-locked to a small segment of the stretch cycle. Under these conditions the transinformation rates grew larger with increasing stimulus frequency and approached 130 bit X s-1 at 60 Hz. Small amplitude sinusoidal stretches, however, evoked considerable transinformation rates in the high frequency region only then, when the spindle receptor was extended to higher prestretch levels. Random stretches evoked transinformation rates between 5 and 30 bit X s-1 depending on both the prestretch level and the intensity of the noise stimulus. The linear response components carried only about 25% of the transinformation rates transmitted by both the linear and nonlinear response components. Auxiliary noise stimuli greatly improved the information transmission of sinusoidal stretches. For example, a pure sinusoid evoked 5 bit X s-1. Adding a noise signal with equal energy to the sinusoidal movement elicited 20 bit X s-1. This facilitation effect of auxiliary noise was restricted to low frequency sinusoidal stimuli. The present results are discussed with respect to the information transmission properties of various sensory systems evaluated by either the same or different information processing procedure as that used in the present study. The functional significance of high transinformation rates sent by the muscle spindle to the central nervous system is discussed with respect to motor control.     

Assessment of TTX-s and TTX-r Action Potential Conduction along Neurites of NGF and GDNF Cultured Porcine DRG Somata

Nine isoforms of voltage-gated sodium channels (NaV) have been characterized and in excitable tissues they are responsible for the initiation and conduction of action potentials. For primary afferent neurons residing in dorsal root ganglia (DRG), individual neurons may express multiple NaV isoforms extending the neuron’s functional capabilities. Since expression of NaV isoforms can be differentially regulated by neurotrophic factors we have examined the functional consequences of exposure to either nerve growth factor (NGF) or glial cell line-derived neurotrophic factor (GDNF) on action potential conduction in outgrowing cultured porcine neurites of DRG neurons. Calcium signals were recorded using the exogenous intensity based calcium indicator Fluo-8®, AM. In 94 neurons, calcium signals were conducted along neurites in response to electrical stimulation of the soma. At an image acquisition rate of 25 Hz it was possible to discern calcium transients in response to individual electrical stimuli. The peak amplitude of electrically-evoked calcium signals was limited by the ability of the neuron to follow the stimulus frequency. The stimulus frequency required to evoke a half-maximal calcium response was approximately 3 Hz at room temperature. In 13 of 14 (93%) NGF-responsive neurites, TTX-r NaV isoforms alone were sufficient to support propagated signals. In contrast, calcium signals mediated by TTX-r NaVs were evident in only 4 of 11 (36%) neurites from somata cultured in GDNF. This establishes a basis for assessing action potential signaling using calcium imaging techniques in individual cultured neurites and suggests that, in the pig, afferent nociceptor classes relying on the functional properties of TTX-r NaV isoforms, such as cold-nociceptors, most probably derive from NGF-responsive DRG neurons.     

Spike-frequency adaptation generates intensity invariance in a primary auditory interneuron

Adaptation of the spike-frequency response to constant stimulation, as observed on various timescales in many neurons, reflects high-pass filter properties of a neuron's transfer function. Adaptation in general, however, is not sufficient to make a neuron's response independent of the mean intensity of a sensory stimulus, since low frequency components of the stimulus are still transmitted, although with reduced gain. We here show, based on an analytically tractable model, that the response of a neuron is intensity invariant, if the fully adapted steady-state spike-frequency response to constant stimuli is independent of stimulus intensity. Electrophysiological recordings from the AN1, a primary auditory interneuron of crickets, show that for intensities above 60 dB SPL (sound pressure level) the AN1 adapted with a time-constant of approximately 40 ms to a steady-state firing rate of approximately 100 Hz. Using identical random amplitude-modulation stimuli we verified that the AN1's spike-frequency response is indeed invariant to the stimulus' mean intensity above 60 dB SPL. The transfer function of the AN1 is a band pass, resulting from a high-pass filter (cutoff frequency at 4 Hz) due to adaptation and a low-pass filter (100 Hz) determined by the steady-state spike frequency. Thus, fast spike-frequency adaptation can generate intensity invariance already at the first level of neural processing.     

Piezoelectric-excited millimeter-sized cantilever (PEMC) sensors detect Bacillus anthracis at 300 spores/mL

Piezoelectric-excited millimeter-sized cantilever (PEMC) sensors consisting of a piezoelectric and a borosilicate glass layer with a sensing area of 2.48 mm2 were fabricated. Antibody specific to Bacillus anthracis (BA, Sterne strain 7702) spores was immobilized on PEMC sensors, and exposed to spores (300 to 3x10(6) spores/mL). The resonant frequency decreased at a rate proportional to the spore concentration and reached a steady state frequency change of 5+/-5 Hz (n=3), 92+/-7 Hz (n=3), 500+/-10 Hz (n=3), 1030+/-10 Hz (n=2), and 2696+/-6 Hz (n=2) corresponding to 0, 3x10(2), 3x10(3), 3x10(4), and 3x10(6) spores/mL, respectively. The reduction in resonant frequency is proportional to the change in cantilever mass, and thus the observed changes are due to the attachment of spores on the sensor surface. Selectivity of the antibody-functionalized sensor was determined with samples of BA (3x10(6)/mL) mixed with Bacillus thuringiensis (BT; 1.5x10(9)/mL) in various volume ratios that yielded BA:BT ratios of 1:0, 1:125, 1:250, 1:500 and 0:1. The corresponding resonance frequency decreases were, respectively, 2345, 1980, 1310, 704 and 10 Hz. Sample containing 100% BT spores (1.5x10(9)/mL and no BA) gave a steady state frequency decrease of 10 Hz, which is within noise level of the sensor, indicating excellent selectivity. The observed binding rate constant for the pure BA and BT-containing samples ranged from 0.105 to 0.043 min-1 in the spore concentration range 300 to 3x10(6)/mL. These results show that detection of B. anthracis spore at a very low concentration (300 spores/mL) and with high selectivity in presence of another Bacillus spore (BT) can be accomplished using piezoelectric-excited millimeter-sized cantilever sensors.     

Current noise parameters derived from voltage noise and impedance in embryonic heart cell aggregates.

We have recorded membrane impedance and voltage noise in the pacemaker range of potentials (-70 to -59 mV) from spheroidal aggregates of 7-d embryonic chick ventricle cells made quiescent by exposure to tetrodotoxin in medium containing 4.5 mM K+. The input capacitance is proportional to aggregate volume and therefore to total membrane area. The specific membrane capacitance is 1.24 microF/cm2. The input resistance at constant potential is inversely proportional to aggregate volume and therefore to total membrane area. The specific membrane resistance in 18 k omega . cm2 at -70 mV and increases to 81 k omega . cm2 at -59 mV. The RC time constant is 22 ms at -70 mV and increases to 146 ms at -59 mV. The aggregate transmembrane small-signal impedance can be represented by a parallel RC circuit itself in parallel with an inductive branch consisting of a resistor (rL) and an inductor (L) in series. The time constant of the inductive branch (L/rL) is 340 ms, and is only weakly dependent on potential. Correlation functions of aggregate voltage noise and the impedance data were modeled by a population of channels with simple open-close kinetics. The time constant of a channel (tau s) derived from the noise analysis is 300 ms. The low frequency limit of the pacemaker current noise (SI[0]), derived from the voltage noise and impedance, increases from 10(-20) A2/Hz . cm2 at -67 mV to 10(-19) A2/Hz . cm2 at -61 mV.     

Membrane potential synchrony in primary visual cortex during sensory stimulation

When the primary visual cortex (V1) is activated by sensory stimulation, what is the temporal correlation between the synaptic inputs to nearby neurons? This question underlies the origin of correlated activity, the mechanism of how visually evoked activity emerges and propagates in cortical circuits, and the relationship between spontaneous and evoked activity. Here, we have recorded membrane potential from pairs of V1 neurons in anesthetized cats and found that visual stimulation suppressed low-frequency membrane potential synchrony (0-10 Hz), and often increased synchrony at high frequencies (20-80 Hz). The increase in high-frequency synchrony occurred for neurons with similar orientation preferences and for neurons with different orientation preferences and occurred for a wide range of stimulus orientations. Thus, while only a subset of neurons spike in response to visual stimulation, a far larger proportion of the circuit is correlated with spiking activity through subthreshold, high-frequency synchronous activity that crosses functional domains.     

‘Communication’ in weakly electric fish, Gnathonemus niger (Mormyridae) I. Variation of electric organ discharge (EOD) frequency elicited by controlled electric stimuli

The weakly electric fish, Gathonemus niger, discharged with a frequency of 4 to 8 Hz during the day and 10 to 16 Hz during the night. The frequency of superimposed burst discharges (32 to 56 Hz) was independent of diurnal factors. The variation of the electric organ discharge frequency during the day was investigated in response to controlled electric stimulus patterns: (a) A free running stimulus frequency of 4 Hz, simulating the resting frequency of another fish, and different stimulus intensities, simulating different distances between two fish. (b) Free running frequencies of 4, 8, 16, …, 128 Hz and two particular stimulus intensities. (c) Discharge coupled stimuli (each discharge triggered an electric stimulus with a fixed delay) and different stimulus intensities.All three kinds of stimuli elicited defined and predictable response discharge patterns supporting the assumption that an electric fish would respond to a particular discharge pattern of another fish also in a similar and predictable manner. Low stimulus intensities (0·04 to 0·2 mV per cm) caused cessation of the discharge activity, a ‘hiding’ or ‘listening’ response. The discharge rate increased linearly with the logarithm of the stimulus intensity. The fish was particularly sensitive to stimulus frequencies which simulated its burst activity (32 to 56 Hz). Discharge coupled stimuli showed that the fish responded to about eight times lower stimulus intensities if the stimulus occurred between two discharges (15 to 30 m-s after the fish's discharge) than if the stimulus occurred within or immediately after the discharge. All suprathreshold stimuli elicited a typical discharge pattern: The irregular resting discharge activity became significantly regular. The degree of regularity was even improved during maintained stimulation. The regularisation of the discharge activity is thought to be involved in the fish's electrolocating system whereas frequency variations are considered as being involved in both the locating system and as communication signals among weakly electric fish.     

A behavioral study of vibrational sensitivity in the pigeon (Columba livia)

Summary The pigeon (Columba livia) has a well-developed ability to detect weak vibrations. Using the method of heart-rate conditioning the vibrational sensitivity was determined for four pigeons at an error probability of P<0.025. The threshold-frequency relationships indicate that the greatest sensitivity to vibrational stimuli is found in the frequency range from 300 to 1,000 Hz with thresholds of about 0.1 m; lowest threshold is 0.04 m at 500 Hz (Fig. 4). Pigeons can respond not only to the frequency of a stimulus, but also to its intensity. The interval decrement (in %) of ECG is a positive correlative function of the stimulus intensity, the calculated values being approximately 4–5% per order of magnitude of the stimulus amplitude (in m) at best frequencies (Fig. 5). The value of vibration detection for birds is discussed.Abbreviation ECG electrocardiogram     

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.     

Electrolocation in the presence of jamming signals: electroreceptor physiology

Responses of ampullary and tuberous electroreceptor afferents were studied using moving electrolocation targets and electrical modulations of the animal's electric organ discharge as stimuli. The ability of the electroreceptors to encode these stimuli was measured with and without various forms of electrical jamming signals. The goal of this study was to measure the deterioration in electroreceptor responses due to the jamming signals, and to compare these results with the behavioral measures of electrolocation under the same conditions of jamming as described in the preceding report (Bastian 1987). 1. Three types of jamming stimuli were used to interfere with the tuberous electroreceptor afferents' ability to respond to the test stimuli mentioned above: Broad-band noise, high frequency stimuli consisting of a sinusoidal waveform having a frequency maintained at a chosen difference frequency (DF) from the EOD frequency of the fish being studied, and 5 or 50 Hz sinusoidal stimuli. 2. The tuberous receptor afferents' spontaneous frequency was sensitive to continuous presentation of all but the 5 Hz jamming signals. The 4 Hz DF signal caused the largest increase in spontaneous activity, the 50 Hz stimulus was intermediate in effectiveness, and the noise stimulus caused the smallest increase. Estimates of the variability of the ongoing receptor activity were also made, and both the 4 Hz DF and the 50 Hz stimuli reduced the coefficient of variation of the receptor activity, but noise had no significant effect on this parameter. Noise, 4 Hz DF, and 50 Hz jamming signals also reduced the tuberous receptors' responses to a 100 ms EOD amplitude modulation, and the 5 Hz stimulus was again ineffective. 3. Noise and 4 Hz DF jamming were also effective in reducing tuberous receptor afferents' responses to a moving metal electrolocation target. The 4 Hz DF stimulus was most effective in reducing the receptor's ability to encode information about the target. Receptor responses showed about a three-fold larger decrease per 10 dB increase in DF jamming amplitude as compared to similar sized increases in noise amplitude. Threshold target distances were also determined with and without noise and DF jamming, and again, the noise stimulus was less effective in reducing the distance at which electrolocation targets were just detectable. 4. Recordings from ampullary receptor afferents confirmed that the galvanic potentials produced by metal electrolocation targets stimulate these receptors while EOD distortions caused by such objects probably do not.(ABSTRACT TRUNCATED AT 400 WORDS)     

Electrical Resonance in the θ Frequency Range in Olfactory Amygdala Neurons

The cortical amygdala receives direct olfactory inputs and is thought to participate in processing and learning of biologically relevant olfactory cues. As for other brain structures implicated in learning, the principal neurons of the anterior cortical nucleus (ACo) exhibit intrinsic subthreshold membrane potential oscillations in the θ-frequency range. Here we show that nearly 50% of ACo layer II neurons also display electrical resonance, consisting of selective responsiveness to stimuli of a preferential frequency (2–6 Hz). Their impedance profile resembles an electrical band-pass filter with a peak at the preferred frequency, in contrast to the low-pass filter properties of other neurons. Most ACo resonant neurons displayed frequency preference along the whole subthreshold voltage range. We used pharmacological tools to identify the voltage-dependent conductances implicated in resonance. A hyperpolarization-activated cationic current depending on HCN channels underlies resonance at resting and hyperpolarized potentials; notably, this current also participates in resonance at depolarized subthreshold voltages. KV7/KCNQ K⁺ channels also contribute to resonant behavior at depolarized potentials, but not in all resonant cells. Moreover, resonance was strongly attenuated after blockade of voltage-dependent persistent Na⁺ channels, suggesting an amplifying role. Remarkably, resonant neurons presented a higher firing probability for stimuli of the preferred frequency. To fully understand the mechanisms underlying resonance in these neurons, we developed a comprehensive conductance-based model including the aforementioned and leak conductances, as well as Hodgkin and Huxley-type channels. The model reproduces the resonant impedance profile and our pharmacological results, allowing a quantitative evaluation of the contribution of each conductance to resonance. It also replicates selective spiking at the resonant frequency and allows a prediction of the temperature-dependent shift in resonance frequency. Our results provide a complete characterization of the resonant behavior of olfactory amygdala neurons and shed light on a putative mechanism for network activity coordination in the intact brain.     

Neonatal heart rate variability and its relation to respiration

The heart rate and respiration signals from nine healthy full term neonates were studied using autoregressive spectral analysis and cross-correlation techniques. The heart rate spectra could be divided into three regions of activity: a very low frequency (VLF) region from 0-0.04 Hz; a low frequency (LF) band from 0.04-0.20 Hz; and a high frequency (HF) region above 0.20 Hz. The newborns exhibited very little respiratory sinus arrhythmia in their heart rate variability in contrast to the situation for adults and older infants. However, variations in heart rate correlated strongly with changes in the breath amplitude, leading to what may be termed a breath amplitude sinus arrhythmia. The neonatal heart rate behaviour under stable conditions of oscillation could be simulated with a nonlinear control model provided the delay time in the baroreceptor loop of the model was taken to be approximately 2 seconds longer than in adults. This is consistent with the immature neurological status of neonates.