Cellular Sensory Mechanisms for Detecting Specific Fold-Changes in Extracellular Cues

Cellular sensory systems often respond not to the absolute levels of inputs but to the fold-changes in inputs. Such a property is called fold-change detection (FCD) and is important for accurately sensing dynamic changes in environmental signals in the presence of fluctuations in their absolute levels. Previous studies defined FCD as input-scale invariance and proposed several biochemical models that achieve such a condition. Here, we prove that the previous FCD models can be approximated by a log-differentiator. Although the log-differentiator satisfies the input-scale invariance requirement, its response amplitude and response duration strongly depend on the input timescale. This creates limitations in the specificity and repeatability of detecting fold-changes in inputs. Nevertheless, FCD with specificity and repeatability by cells has been reported in the context of Drosophila wing development. Motivated by this fact and by extending previous FCD models, we here propose two possible mechanisms to achieve FCD with specificity and repeatability. One is the integrate-and-fire type: a system integrates the rate of temporal change in input and makes a response when the integrated value reaches a constant threshold, and this is followed by the reset of the integrated value. The other is the dynamic threshold type: a system response occurs when the input level reaches a threshold, whose value is multiplied by a certain constant after each response. These two mechanisms can be implemented biochemically by appropriately combining feed-forward and feedback loops. The main difference between the two models is their memory of input history; we discuss possible ways to distinguish between the two models experimentally.     

Comparing Apples and Oranges: Fold-Change Detection of Multiple Simultaneous Inputs

Sensory systems often detect multiple types of inputs. For example, a receptor in a cell-signaling system often binds multiple kinds of ligands, and sensory neurons can respond to different types of stimuli. How do sensory systems compare these different kinds of signals? Here, we consider this question in a class of sensory systems – including bacterial chemotaxis- which have a property known as fold-change detection: their output dynamics, including amplitude and response time, depends only on the relative changes in signal, rather than absolute changes, over a range of several decades of signal. We analyze how fold-change detection systems respond to multiple signals, using mathematical models. Suppose that a step of fold F₁ is made in input 1, together with a step of F₂ in input 2. What total response does the system provide? We show that when both input signals impact the same receptor with equal number of binding sites, the integrated response is multiplicative: the response dynamics depend only on the product of the two fold changes, F₁F₂. When the inputs bind the same receptor with different number of sites n₁ and n₂, the dynamics depend on a product of power laws, . Thus, two input signals which vary over time in an inverse way can lead to no response. When the two inputs affect two different receptors, other types of integration may be found and generally the system is not constrained to respond according to the product of the fold-change of each signal. These predictions can be readily tested experimentally, by providing cells with two simultaneously varying input signals. The present study suggests how cells can compare apples and oranges, namely by comparing each to its own background level, and then multiplying these two fold-changes.     

The response of a model neurone to a white noise input

We consider a simple model of a neurone in which the input voltage is integrated to form the somatic potential, and a pulse is emitted when this reaches a threshold; the somatic potential is then reset to its resting value. We subject this model to a white-noise input, and evaluate the cross correlation between input white noise and the output pulse train; this is proportional to the small-signal impulse response of the model. Some numerical estimations are presented.     

Slope-Based Stochastic Resonance: How Noise Enables Phasic Neurons to Encode Slow Signals

Fundamental properties of phasic firing neurons are usually characterized in a noise-free condition. In the absence of noise, phasic neurons exhibit Class 3 excitability, which is a lack of repetitive firing to steady current injections. For time-varying inputs, phasic neurons are band-pass filters or slope detectors, because they do not respond to inputs containing exclusively low frequencies or shallow slopes. However, we show that in noisy conditions, response properties of phasic neuron models are distinctly altered. Noise enables a phasic model to encode low-frequency inputs that are outside of the response range of the associated deterministic model. Interestingly, this seemingly stochastic-resonance (SR) like effect differs significantly from the classical SR behavior of spiking systems in both the signal-to-noise ratio and the temporal response pattern. Instead of being most sensitive to the peak of a subthreshold signal, as is typical in a classical SR system, phasic models are most sensitive to the signal''s rising and falling phases where the slopes are steep. This finding is consistent with the fact that there is not an absolute input threshold in terms of amplitude; rather, a response threshold is more properly defined as a stimulus slope/frequency. We call the encoding of low-frequency signals with noise by phasic models a slope-based SR, because noise can lower or diminish the slope threshold for ramp stimuli. We demonstrate here similar behaviors in three mechanistic models with Class 3 excitability in the presence of slow-varying noise and we suggest that the slope-based SR is a fundamental behavior associated with general phasic properties rather than with a particular biological mechanism.     

Threshold model for clustered firing of neurons

Two probabilistic threshold models for burst activity of cortical neurons are proposed. In model I every input impulse increases the summed effect of previous input impulses by one unit. The decay of the summed effect takes place in discrete steps of one unit. A response occurs on arrival of an input impulse, when a threshold value is attained.Although after a response the summed effect is not reset to zero, it cannot exceed the threshold either. The distribution of intervals can be resolved in two components, one for long and one for short intervals. In model II intervals of the short component are terminated by a multiple response instead of one response.     

Logarithmic and Power Law Input-Output Relations in Sensory Systems with Fold-Change Detection

Two central biophysical laws describe sensory responses to input signals. One is a logarithmic relationship between input and output, and the other is a power law relationship. These laws are sometimes called the Weber-Fechner law and the Stevens power law, respectively. The two laws are found in a wide variety of human sensory systems including hearing, vision, taste, and weight perception; they also occur in the responses of cells to stimuli. However the mechanistic origin of these laws is not fully understood. To address this, we consider a class of biological circuits exhibiting a property called fold-change detection (FCD). In these circuits the response dynamics depend only on the relative change in input signal and not its absolute level, a property which applies to many physiological and cellular sensory systems. We show analytically that by changing a single parameter in the FCD circuits, both logarithmic and power-law relationships emerge; these laws are modified versions of the Weber-Fechner and Stevens laws. The parameter that determines which law is found is the steepness (effective Hill coefficient) of the effect of the internal variable on the output. This finding applies to major circuit architectures found in biological systems, including the incoherent feed-forward loop and nonlinear integral feedback loops. Therefore, if one measures the response to different fold changes in input signal and observes a logarithmic or power law, the present theory can be used to rule out certain FCD mechanisms, and to predict their cooperativity parameter. We demonstrate this approach using data from eukaryotic chemotaxis signaling.     

Neurobiological models of two-choice decision making can be reduced to a one-dimensional nonlinear diffusion equation

The response behaviors in many two-alternative choice tasks are well described by so-called sequential sampling models. In these models, the evidence for each one of the two alternatives accumulates over time until it reaches a threshold, at which point a response is made. At the neurophysiological level, single neuron data recorded while monkeys are engaged in two-alternative choice tasks are well described by winner-take-all network models in which the two choices are represented in the firing rates of separate populations of neurons. Here, we show that such nonlinear network models can generally be reduced to a one-dimensional nonlinear diffusion equation, which bears functional resemblance to standard sequential sampling models of behavior. This reduction gives the functional dependence of performance and reaction-times on external inputs in the original system, irrespective of the system details. What is more, the nonlinear diffusion equation can provide excellent fits to behavioral data from two-choice decision making tasks by varying these external inputs. This suggests that changes in behavior under various experimental conditions, e.g. changes in stimulus coherence or response deadline, are driven by internal modulation of afferent inputs to putative decision making circuits in the brain. For certain model systems one can analytically derive the nonlinear diffusion equation, thereby mapping the original system parameters onto the diffusion equation coefficients. Here, we illustrate this with three model systems including coupled rate equations and a network of spiking neurons.     

The Frequency Response, Coherence, and Information Capacity of Two Neuronal Models

Two neuronal models are analyzed in which subthreshold inputs are integrated either without loss (perfect integrator) or with a decay which follows an exponential time course (leaky integrator). Linear frequency response functions for these models are compared using sinusoids, Poisson-distributed impulses, or gaussian white noise as inputs. The responses of both models show the nonlinear behavior characteristic of a rectifier for sinusoidal inputs of sufficient amplitude. The leaky integrator shows another nonlinearity in which responses become phase locked to cyclic stimuli. Addition of white noise reduces the distortions due to phase locking. Both models also show selective attenuation of high-frequency components with white noise inputs. Input, output, and cross-spectra are computed using inputs having a broad frequency spectrum. Measures of the coherence and information transmission between the input and output of the models are also derived. Steady inputs, which produce a constant “carrier” rate, and intrinsic sources, which produce variability in the discharge of neurons, may either increase or decrease coherence; however, information transmission using inputs with a broad spectrum is generally increased by steady inputs and reduced by intrinsic variability.     

Firing-rate models for neurons with a broad repertoire of spiking behaviors

Capturing the response behavior of spiking neuron models with rate-based models facilitates the investigation of neuronal networks using powerful methods for rate-based network dynamics. To this end, we investigate the responses of two widely used neuron model types, the Izhikevich and augmented multi-adapative threshold (AMAT) models, to a range of spiking inputs ranging from step responses to natural spike data. We find (i) that linear-nonlinear firing rate models fitted to test data can be used to describe the firing-rate responses of AMAT and Izhikevich spiking neuron models in many cases; (ii) that firing-rate responses are generally too complex to be captured by first-order low-pass filters but require bandpass filters instead; (iii) that linear-nonlinear models capture the response of AMAT models better than of Izhikevich models; (iv) that the wide range of response types evoked by current-injection experiments collapses to few response types when neurons are driven by stationary or sinusoidally modulated Poisson input; and (v) that AMAT and Izhikevich models show different responses to spike input despite identical responses to current injections. Together, these findings suggest that rate-based models of network dynamics may capture a wider range of neuronal response properties by incorporating second-order bandpass filters fitted to responses of spiking model neurons. These models may contribute to bringing rate-based network modeling closer to the reality of biological neuronal networks.     

The estimation of the frequency response function of a mechanoreceptor

Widespread use has been made of linear systems theory to describe the input-output relations of receptors. The frequency response function of an insect mechanoreceptor, the tactile spine of the cockroach, has been estimated by using deterministic inputs (sines and step functions), deterministic inputs added to a stochastic, auxiliary signal (band-limited white noise), and a stochastic input alone. When a stochastic input is used, spectral analysis provides methods for estimating the coherence function as well as the frequency response function. The coherence function of the tactile spine is low, suggesting that the linear frequency response function is not a good characterization of the input-output relation of the receptor. Two non-linearities, rectification and phase-locking are described. Rectification can reduce the absolute value of the frequency response measured using sine waves of all frequencies without changing its form. Phase-locking changes the form of the frequency response function at high frequencies. Use of a stochastic auxiliary signal linearizes the input-output relations of the receptor in the sense that the cycle histograms obtained with sinusoidal inputs are more sinusoidal and the form of the frequency response function agrees with that predicted from the step response over a wider range of frequencies.     

Spatiotemporal asymmetry of associative synaptic plasticity in fear conditioning pathways

Input-specific long-term potentiation (LTP) in afferent inputs to the amygdala serves an essential function in the acquisition of fear memory. Factors underlying input specificity of synaptic modifications implicated in information transfer in fear conditioning pathways remain unclear. Here we show that the strength of naive synapses in two auditory inputs converging on a single neuron in the lateral nucleus of the amygdala (LA) is only modified when a postsynaptic action potential closely follows a synaptic response. The stronger inhibitory drive in thalamic pathway, as compared with cortical input, hampers the induction of LTP at thalamo-amygdala synapses, contributing to the spatial specificity of LTP in convergent inputs. These results indicate that spike timing-dependent synaptic plasticity in afferent projections to the LA is both temporarily and spatially asymmetric, thus providing a mechanism for the conditioned stimulus discrimination during fear behavior.     

On some properties of stochastic information processes in neurons and neuron populations

The information in the nervous spike trains and its processing by neural units are discussed. In these problems, our attention is focused on the stochastic properties of neurons and neuron populations. There are three subjects in this paper, which are the spontaneous type neuron, the forced type neuron and the reciprocal inhibitory pairs.

1. The spontaneous type neuron produces spikes without excitatory inputs. The mathematical model has the following assumptions. The neuron potential (NP) has the fluctuation and obeys the Ornstein-Uhlenbeck process, because the N P is not so perfectly random as that of the Wiener process but has an attraction to the rest value. The threshold varies exponentially and the NP has the constant lower limit. When the NP reaches the threshold, the neuron fires and the NP is reset to a certain position. After a firing, an absolute refractory period exists. In discussing the stochastic properties of neurons, the transition probability density function and the first passage time density function are the important quantities, which are governed by the Kolmogorov's equations. Although they can be set up easily, we can rarely obtain the analytical solutions in time domain. Moreover, they cover only simple properties. Hence the numerical analysis is performed and a good deal of fair results are obtained and discussed.
2. The forced type neuron has input pulse trains which are assumed to be based on the Poisson process. Other assumptions and methods are almost the same as above except the diffusion approximation of the stochastic process. In this case, we encounter the inhomogeneous process due to the pulse-frequency-modulation, whose first passage time density reveals the multimodal distribution. The numerical analysis is also tried, and the output spike interval density is further discussed in the case of the periodic modulation.
3. Two types of reciprocal inhibitory pairs are discussed. The first type has two excitatory driving inputs which are mutually independent. The second type has one common excitatory input but it advances in two ways, one of which has a time lag. The neuron dynamics is the same as that of the forced type neuron and each neuron has an identical structure. The inputs are assumed to be based on the Poisson process and the inhibition occurs when the companion neuron fires. In this case, the equations of the probability density functions are not obtained. Hence the computer simulation is tried and it is observed that the stochastic rhythm emerges in spite of the temporally homogeneous inputs. Furthermore, the case of inhomogeneous inputs is discussed.

     

A Single Compartment Quantal Response Model with a Renewal Type Input Process and IID Releases

We deal with a single compartment quantal response model, where unlike the previous models, which do not have any input after the administration of a single dose Z(0)=z at time t=0, we allow inputs of doses after time t=0. More precisely, the system uses the (s, S) input policy as in inventory models, and has IID releases. Also when the amount of dose in the subject reaches 0, there is a probability p to stop having input thereafter. Among other results, the probability that the subject never responds and the expressions for some quantities of interest are obtained.     

Direction selectivity of excitation and inhibition in simple cells of the cat primary visual cortex

Direction selectivity in simple cells of primary visual cortex, defined from their spike responses, cannot be predicted using linear models. It has been suggested that the shunting inhibition evoked by visual stimulation is responsible for the nonlinear component of direction selectivity. Cortical inhibition would suppress a neuron's firing when stimuli move in the nonpreferred direction, but would allow responses to stimuli in the preferred direction. Models of direction selectivity based solely on input from the lateral geniculate nucleus, however, propose that the nonlinear response is caused by spike threshold. By extracting excitatory and inhibitory components of synaptic inputs from intracellular records obtained in vivo, we demonstrate that excitation and inhibition are tuned for the same direction, but differ in relative timing. Further, membrane potential responses combine in a linear fashion. Spike threshold, however, quantitatively accounts for the nonlinear component of direction selectivity, amplifying the direction selectivity of spike output relative to that of synaptic inputs.     

A review of the integrate-and-fire neuron model: II. Inhomogeneous synaptic input and network properties

The integrate-and-fire neuron model describes the state of a neuron in terms of its membrane potential, which is determined by the synaptic inputs and the injected current that the neuron receives. When the membrane potential reaches a threshold, an action potential (spike) is generated. This review considers the model in which the synaptic input varies periodically and is described by an inhomogeneous Poisson process, with both current and conductance synapses. The focus is on the mathematical methods that allow the output spike distribution to be analyzed, including first passage time methods and the Fokker–Planck equation. Recent interest in the response of neurons to periodic input has in part arisen from the study of stochastic resonance, which is the noise-induced enhancement of the signal-to-noise ratio. Networks of integrate-and-fire neurons behave in a wide variety of ways and have been used to model a variety of neural, physiological, and psychological phenomena. The properties of the integrate-and-fire neuron model with synaptic input described as a temporally homogeneous Poisson process are reviewed in an accompanying paper (Burkitt in Biol Cybern, 2006).     

Enhanced Sensitivity to Rapid Input Fluctuations by Nonlinear Threshold Dynamics in Neocortical Pyramidal Neurons

The way in which single neurons transform input into output spike trains has fundamental consequences for network coding. Theories and modeling studies based on standard Integrate-and-Fire models implicitly assume that, in response to increasingly strong inputs, neurons modify their coding strategy by progressively reducing their selective sensitivity to rapid input fluctuations. Combining mathematical modeling with in vitro experiments, we demonstrate that, in L5 pyramidal neurons, the firing threshold dynamics adaptively adjust the effective timescale of somatic integration in order to preserve sensitivity to rapid signals over a broad range of input statistics. For that, a new Generalized Integrate-and-Fire model featuring nonlinear firing threshold dynamics and conductance-based adaptation is introduced that outperforms state-of-the-art neuron models in predicting the spiking activity of neurons responding to a variety of in vivo-like fluctuating currents. Our model allows for efficient parameter extraction and can be analytically mapped to a Generalized Linear Model in which both the input filter—describing somatic integration—and the spike-history filter—accounting for spike-frequency adaptation—dynamically adapt to the input statistics, as experimentally observed. Overall, our results provide new insights on the computational role of different biophysical processes known to underlie adaptive coding in single neurons and support previous theoretical findings indicating that the nonlinear dynamics of the firing threshold due to Na⁺-channel inactivation regulate the sensitivity to rapid input fluctuations.     

Threshold-dominated regulation hides genetic variation in gene expression networks

Background  

In dynamical models with feedback and sigmoidal response functions, some or all variables have thresholds around which they regulate themselves or other variables. A mathematical analysis has shown that when the dose-response functions approach binary or on/off responses, any variable with an equilibrium value close to one of its thresholds is very robust to parameter perturbations of a homeostatic state. We denote this threshold robustness. To check the empirical relevance of this phenomenon with response function steepnesses ranging from a near on/off response down to Michaelis-Menten conditions, we have performed a simulation study to investigate the degree of threshold robustness in models for a three-gene system with one downstream gene, using several logical input gates, but excluding models with positive feedback to avoid multistationarity. Varying parameter values representing functional genetic variation, we have analysed the coefficient of variation (CV) of the gene product concentrations in the stable state for the regulating genes in absolute terms and compared to the CV for the unregulating downstream gene. The sigmoidal or binary dose-response functions in these models can be considered as phenomenological models of the aggregated effects on protein or mRNA expression rates of all cellular reactions involved in gene expression.     

Lateral input of particulate organic matter from bank slopes surpasses direct litter fall in the uppermost reaches of a headwater stream in Hokkaido, Japan

In forested streams, surrounding riparian forests provide essential supplies of organic matter to aquatic ecosystems. We focused on two pathways of particulate organic matter inputs: direct input from upper riparian forests and indirect lateral input from bank slopes, for which there are limited quantitative data. We investigated the inputs of coarse particulate organic matter (CPOM) and carbon and nitrogen in the CPOM into the uppermost reaches of a headwater stream with steep bank slopes in Hokkaido, Japan. CPOM collected by litter traps was divided into categories (e.g., leaves, twigs) and weighed. Monthly nitrogen and carbon inputs were also estimated. The annual direct input of CPOM (ash-free dry mass) was 472 g m⁻², a common value for temperate riparian forests. The annual lateral CPOM input was 353 g m⁻¹ and 941 g m⁻² when they were converted to area base. This value surpassed the direct input. Organic matter that we could not separate from inorganic sediments contributed to the total lateral input from the bank slopes (124 g m⁻¹); this organic matter contained relatively high amounts of nitrogen and carbon. At uppermost stream reaches, the bank slope would be a key factor to understanding the carbon and nitrogen pathways from the surrounding terrestrial ecosystem to the aquatic ecosystem.     

Effects of active versus passive dendritic membranes on the transfer properties of a simulated neuron

In a simulated neuron with a dendritic tree, the relative effects of active and passive dendritic membranes on transfer properties were studied. The simulations were performed by means of a digital computer. The computations calculated the changes in transmembrane voltages of many compartments over time as a function of other biophysical variables. These variables were synaptic input intensity, critical firing threshold, rate of leakage of current across the membrane, and rate of longitudinal current spread between compartments. For both passive and active dendrites, the transfer properties of the soma studied for different rates of longitudinal current spread. With low rates of current spread, graded changes in firing threshold produced correspondingly graded changes in output discharge. With high rates of current spread, the neuron became a bistable operator where spiking was enhanced if the threshold was below a certain level and suppressed if the threshold was above that level. Since alterations in firing threshold were shown to have the same effect on firing rate as alterations in synaptic input intensity, the neuron can be said to change from graded to contrast-enhancing in its response to stimuli of different intensities. The presence or absence of dendritic spiking was found to have a significant effect on the integrative properties of the simulated neuron. In particular, contrast enhancement was considerably more pronounced in neurons with passive than with active dendrites in that somatic spike rates reached a higher maximum when dendrites were passive. With active dendrites, a less intense input was needed to initiate somatic spiking than with passive dendrites because a distal dendritic spike could easily propagate by means of longitudinal current spread to the soma. Once somatic spiking was initiated, though, spike rates tended to be lower with active than with passive dendrites because the soma recovered more slowly from its post-spike refractory period if it was also influenced by refractory periods in the dendrites. The experiment of comparing neurons with active and passive dendrites was repeated at a different, higher value of synaptic input. The same differences in transfer properties between the active and passive cases emerged as before. Spiking patterns in neurons with active dendrites were also affected by the time distribution of synaptic inputs. In a previous study, inputs had been random over both space and time, varying about a predetermined mean, whereas in the present study, inputs were random over space but uniform over time. When inputs were made uniform over time, spiking became more difficult to initiate and the transition from graded to bistable response became less sharp.     

The response of a nerve cylinder to spatially distributed white noise inputs

The response of a passive nerve cylinder (or dendritic tree in the equivalent cylinder representation) to random white noise input currents is determined. Results for the mean, variance and covariance of the depolarization are obtained for an arbitrary number of independent spatially distributed inputs. The case of a cylinder with sealed ends is considered in detail. The differences that arise when the input currents are distributed over a small but finite region of space instead of concentrated at a point are investigated. In the case of distributed inputs, the expectation is smoother near the stimulus and the variance becomes finite over the entire cable length including the region of the applied stimulus. Away from the stimulus, there are no appreciable differences between the responses for the two cases. The interaction between an excitatory input and an inhibitory input at various locations is examined and one case of more than two inputs is also analysed to study effects which could not have been discerned from point models for a neuron with random inputs.