Hyperpolarization-Activated Current Induces Period-Doubling Cascades and Chaos in a Cold Thermoreceptor Model

Dynamics of antiphase bursting modulated by the inhibitory synaptic and hyperpolarization-activated cation currents

Antiphase bursting related to the rhythmic motor behavior exhibits complex dynamics modulated by the inhibitory synaptic current (Isyn), especially in the presence of the hyperpolarization-activated cation current (Ih). In the present paper, the dynamics of antiphase bursting modulated by the Ih and Isyn is studied in three aspects with a theoretical model. Firstly, the Isyn and the slow Ih with strong strength are the identified to be the necessary conditions for the antiphase bursting. The dependence of the antiphase bursting on the two currents is different for low (escape mode) and high (release mode) threshold voltages (Vth) of the inhibitory synapse. Secondly, more detailed co-regulations of the two currents to induce opposite changes of the bursting period are obtained. For the escape mode, increase of the Ih induces elevated membrane potential of the silence inhibited by a strong Isyn and shortened silence duration to go beyond Vth, resulting in reduced bursting period. For the release mode, increase of the Ih induces elevated tough value of the former part of the burst modulated by a nearly zero Isyn and lengthen burst duration to fall below Vth, resulting in prolonged bursting period. Finally, the fast-slow dynamics of the antiphase bursting are acquired. Using one-and two-parameter bifurcations of the fast subsystem of a single neuron, the burst of the antiphase bursting is related to the stable limit cycle, and the silence modulated by a strong Isyn to the stable equilibrium to a certain extent. The Ih mainly modulates the dynamics within the burst and quiescent state. Furthermore, with the fast subsystem of the coupled neurons, the silence is associated with the unstable equilibrium point. The results present theoretical explanations to the changes in the bursting period and fast-slow dynamics of the antiphase bursting modulated by the Isyn and Ih, which is helpful for understanding the antiphase bursting and modulating rhythmic motor patterns.

Transient and Steady-State Properties of Drosophila Sensory Neurons Coding Noxious Cold Temperature

Coding noxious cold signals, such as the magnitude and rate of temperature change, play essential roles in the survival of organisms. We combined electrophysiological and computational neuroscience methods to investigate the neural dynamics of Drosophila larva cold-sensing Class III (CIII) neurons. In response to a fast temperature change (-2 to -6°C/s) from room temperature to noxious cold, the CIII neurons exhibited a pronounced peak of a spiking rate with subsequent relaxation to a steady-state spiking. The magnitude of the peak was higher for a higher rate of temperature decrease, while slow temperature decrease (-0.1°C/s) evoked no distinct peak of the spiking rate. The rate of the steady-state spiking depended on the magnitude of the final temperature and was higher at lower temperatures. For each neuron, we characterized this dependence by estimating the temperature of the half activation of the spiking rate by curve fitting neuron's spiking rate responses to a Boltzmann function. We found that neurons had a temperature of the half activation distributed over a wide temperature range. We also found that CIII neurons responded to decrease rather than increase in temperature. There was a significant difference in spiking activity between fast and slow returns from noxious cold to room temperature: The CIII neurons usually stopped activity abruptly in the case of the fast return and continued spiking for some time in the case of the slow return. We developed a biophysical model of CIII neurons using a generalized description of transient receptor potential (TRP) current kinetics with temperature-dependent activation and Ca²⁺-dependent inactivation. This model recapitulated the key features of the spiking rate responses found in experiments and suggested mechanisms explaining the transient and steady-state activity of the CIII neurons at different cold temperatures and rates of their decrease and increase. We conclude that CIII neurons encode at least three types of cold sensory information: the rate of temperature decrease by a peak of the firing rate, the magnitude of cold temperature by the rate of steady spiking activity, and direction of temperature change by spiking activity augmentation or suppression corresponding to temperature decrease and increase, respectively.

Visualization of currents in neural models with similar behavior and different conductance densities

Conductance-based models of neural activity produce large amounts of data that can be hard to visualize and interpret. We introduce visualization methods to display the dynamics of the ionic currents and to display the models’ response to perturbations. To visualize the currents’ dynamics, we compute the percent contribution of each current and display them over time using stacked-area plots. The waveform of the membrane potential and the contribution of each current change as the models are perturbed. To represent these changes over a range of the perturbation control parameter, we compute and display the distributions of these waveforms. We illustrate these procedures in six examples of bursting model neurons with similar activity but that differ as much as threefold in their conductance densities. These visualization methods provide heuristic insight into why individual neurons or networks with similar behavior can respond widely differently to perturbations.

The nervous system contains networks of neurons that generate electrical signals to communicate with each other and the rest of the body. Such electrical signals are due to the flow of ions into or out of the neurons via proteins known as ion channels. Neurons have many different kinds of ion channels that only allow specific ions to pass. Therefore, for a neuron to produce an electrical signal, the activities of several different ion channels need to be coordinated so that they all open and close at certain times.

Researchers have previously used data collected from various experiments to develop detailed models of electrical signals in neurons. These models incorporate information about how and when the ion channels may open and close, and can produce numerical simulations of the different ionic currents. However, it can be difficult to display the currents and observe how they change when several different ion channels are involved.

Alonso and Marder used simple mathematical concepts to develop new methods to display ionic currents in computational models of neurons. These tools use color to capture changes in ionic currents and provide insights into how the opening and closing of ion channels shape electrical signals.

The methods developed by Alonso and Marder could be adapted to display the behavior of biochemical reactions or other topics in biology and may, therefore, be useful to analyze data generated by computational models of many different types of cells. Additionally, these methods may potentially be used as educational tools to illustrate the coordinated opening and closing of ion channels in neurons and other fundamental principles of neuroscience that are otherwise hard to demonstrate.

Chaos versus noise as drivers of multistability in neural networks

The multistable behavior of neural networks is actively being studied as a landmark of ongoing cerebral activity, reported in both functional Magnetic Resonance Imaging (fMRI) and electro- or magnetoencephalography recordings. This consists of a continuous jumping between different partially synchronized states in the absence of external stimuli. It is thought to be an important mechanism for dealing with sensory novelty and to allow for efficient coding of information in an ever-changing surrounding environment. Many advances have been made to understand how network topology, connection delays, and noise can contribute to building this dynamic. Little or no attention, however, has been paid to the difference between local chaotic and stochastic influences on the switching between different network states. Using a conductance-based neural model that can have chaotic dynamics, we showed that a network can show multistable dynamics in a certain range of global connectivity strength and under deterministic conditions. In the present work, we characterize the multistable dynamics when the networks are, in addition to chaotic, subject to ion channel stochasticity in the form of multiplicative (channel) or additive (current) noise. We calculate the Functional Connectivity Dynamics matrix by comparing the Functional Connectivity (FC) matrices that describe the pair-wise phase synchronization in a moving window fashion and performing clustering of FCs. Moderate noise can enhance the multistable behavior that is evoked by chaos, resulting in more heterogeneous synchronization patterns, while more intense noise abolishes multistability. In networks composed of nonchaotic nodes, some noise can induce multistability in an otherwise synchronized, nonchaotic network. Finally, we found the same results regardless of the multiplicative or additive nature of noise.

Synchronization transition in neuronal networks composed of chaotic or non-chaotic oscillators

Chaotic dynamics has been shown in the dynamics of neurons and neural networks, in experimental data and numerical simulations. Theoretical studies have proposed an underlying role of chaos in neural systems. Nevertheless, whether chaotic neural oscillators make a significant contribution to network behaviour and whether the dynamical richness of neural networks is sensitive to the dynamics of isolated neurons, still remain open questions. We investigated synchronization transitions in heterogeneous neural networks of neurons connected by electrical coupling in a small world topology. The nodes in our model are oscillatory neurons that – when isolated – can exhibit either chaotic or non-chaotic behaviour, depending on conductance parameters. We found that the heterogeneity of firing rates and firing patterns make a greater contribution than chaos to the steepness of the synchronization transition curve. We also show that chaotic dynamics of the isolated neurons do not always make a visible difference in the transition to full synchrony. Moreover, macroscopic chaos is observed regardless of the dynamics nature of the neurons. However, performing a Functional Connectivity Dynamics analysis, we show that chaotic nodes can promote what is known as multi-stable behaviour, where the network dynamically switches between a number of different semi-synchronized, metastable states.

Dynamical Mechanism of Hyperpolarization-Activated Non-specific Cation Current Induced Resonance and Spike-Timing Precision in a Neuronal Model

Hyperpolarization-activated cyclic nucleotide-gated cation current (I_h) plays important roles in the achievement of many physiological/pathological functions in the nervous system by modulating the electrophysiological activities, such as the rebound (spike) to hyperpolarization stimulations, subthreshold membrane resonance to sinusoidal currents, and spike-timing precision to stochastic factors. In the present paper, with increasing g_h (conductance of I_h), the rebound (spike) and subthreshold resonance appear and become stronger, and the variability of the interspike intervals (ISIs) becomes lower, i.e., the enhancement of spike-timing precision, which are simulated in a conductance-based theoretical model and well explained by the nonlinear concept of bifurcation. With increasing g_h, the stable node to stable focus, to coexistence behavior, and to firing via the codimension-1 bifurcations (Hopf bifurcation, saddle-node bifurcation, saddle-node bifurcations on an invariant circle, and saddle homoclinic orbit) and codimension-2 bifurcations such as Bogdanov-Takens (BT) point related to the transition between saddle-node and Hopf bifurcations, are acquired with 1- and 2-parameter bifurcation analysis. The decrease of variability of ISIs with increasing g_h is induced by the fast decrease of the standard deviation of ISIs, which is related to the increase of the capacity of resisting noisy disturbance due to the firing becomes far away from the bifurcation point. The enhancement of the rebound (spike) with increasing g_h builds up a relationship to the decrease of the capacity of resisting disturbance like the hyperpolarization stimulus as the resting state approaches the bifurcation point. The “typical”-resonance and non-resonance appear in the parameter region of the stable focus and node far away from the bifurcation points, respectively. The complex or “strange” dynamics, such as the “weak”-resonance for the stable node near the transition point between the stable node and focus and the non-resonance for the stable focus close to the codimension-1 and −2 bifurcation points, are discussed.