Regulation of spike timing in visual cortical circuits

A train of action potentials (a spike train) can carry information in both the average firing rate and the pattern of spikes in the train. But can such a spike-pattern code be supported by cortical circuits? Neurons in vitro produce a spike pattern in response to the injection of a fluctuating current. However, cortical neurons in vivo are modulated by local oscillatory neuronal activity and by top-down inputs. In a cortical circuit, precise spike patterns thus reflect the interaction between internally generated activity and sensory information encoded by input spike trains. We review the evidence for precise and reliable spike timing in the cortex and discuss its computational role.     

Biophysical mechanism of spike threshold dependence on the rate of rise of the membrane potential by sodium channel inactivation or subthreshold axonal potassium current

Spike threshold filters incoming inputs and thus gates activity flow through neuronal networks. Threshold is variable, and in many types of neurons there is a relationship between the threshold voltage and the rate of rise of the membrane potential (dVm/dt) leading to the spike. In primary sensory cortex this relationship enhances the sensitivity of neurons to a particular stimulus feature. While Na⁺ channel inactivation may contribute to this relationship, recent evidence indicates that K⁺ currents located in the spike initiation zone are crucial. Here we used a simple Hodgkin-Huxley biophysical model to systematically investigate the role of K⁺ and Na⁺ current parameters (activation voltages and kinetics) in regulating spike threshold as a function of dVm/dt. Threshold was determined empirically and not estimated from the shape of the Vm prior to a spike. This allowed us to investigate intrinsic currents and values of gating variables at the precise voltage threshold. We found that Na⁺ inactivation is sufficient to produce the relationship provided it occurs at hyperpolarized voltages combined with slow kinetics. Alternatively, hyperpolarization of the K⁺ current activation voltage, even in the absence of Na⁺ inactivation, is also sufficient to produce the relationship. This hyperpolarized shift of K⁺ activation allows an outward current prior to spike initiation to antagonize the Na⁺ inward current such that it becomes self-sustaining at a more depolarized voltage. Our simulations demonstrate parameter constraints on Na⁺ inactivation and the biophysical mechanism by which an outward current regulates spike threshold as a function of dVm/dt.     

EGF converts transit-amplifying neurogenic precursors in the adult brain into multipotent stem cells

We studied the transformation of sensory input as it progresses from vibrissa primary sensor (S1) to motor (M1) cortex. Single-unit activity was obtained from alert adult rats that did not to whisk upon application of punctate, rhythmic stimulation of individual vibrissae. The spike response of units in S1 cortex largely reproduced the shape of the stimulus. In contrast, the spiking output of units in M1 cortex were modulated solely as a sinusoid at the repetition rate of the stimulus for frequencies between 5 and 15 Hz; this range corresponds to that of natural whisking. Thus, the S1 to M1 transformation extracts the fundamental frequency from a spectrally rich stimulus. We discuss our results in terms of a band-pass filter with a center frequency that adapts to the change in stimulation rate.     

Learning of Precise Spike Times with Homeostatic Membrane Potential Dependent Synaptic Plasticity

Precise spatio-temporal patterns of neuronal action potentials underly e.g. sensory representations and control of muscle activities. However, it is not known how the synaptic efficacies in the neuronal networks of the brain adapt such that they can reliably generate spikes at specific points in time. Existing activity-dependent plasticity rules like Spike-Timing-Dependent Plasticity are agnostic to the goal of learning spike times. On the other hand, the existing formal and supervised learning algorithms perform a temporally precise comparison of projected activity with the target, but there is no known biologically plausible implementation of this comparison. Here, we propose a simple and local unsupervised synaptic plasticity mechanism that is derived from the requirement of a balanced membrane potential. Since the relevant signal for synaptic change is the postsynaptic voltage rather than spike times, we call the plasticity rule Membrane Potential Dependent Plasticity (MPDP). Combining our plasticity mechanism with spike after-hyperpolarization causes a sensitivity of synaptic change to pre- and postsynaptic spike times which can reproduce Hebbian spike timing dependent plasticity for inhibitory synapses as was found in experiments. In addition, the sensitivity of MPDP to the time course of the voltage when generating a spike allows MPDP to distinguish between weak (spurious) and strong (teacher) spikes, which therefore provides a neuronal basis for the comparison of actual and target activity. For spatio-temporal input spike patterns our conceptually simple plasticity rule achieves a surprisingly high storage capacity for spike associations. The sensitivity of the MPDP to the subthreshold membrane potential during training allows robust memory retrieval after learning even in the presence of activity corrupted by noise. We propose that MPDP represents a biophysically plausible mechanism to learn temporal target activity patterns.     

Phase-of-firing coding of natural visual stimuli in primary visual cortex

We investigated the hypothesis that neurons encode rich naturalistic stimuli in terms of their spike times relative to the phase of ongoing network fluctuations rather than only in terms of their spike count. We recorded local field potentials (LFPs) and multiunit spikes from the primary visual cortex of anaesthetized macaques while binocularly presenting a color movie. We found that both the spike counts and the low-frequency LFP phase were reliably modulated by the movie and thus conveyed information about it. Moreover, movie periods eliciting higher firing rates also elicited a higher reliability of LFP phase across trials. To establish whether the LFP phase at which spikes were emitted conveyed visual information that could not be extracted by spike rates alone, we compared the Shannon information about the movie carried by spike counts to that carried by the phase of firing. We found that at low LFP frequencies, the phase of firing conveyed 54% additional information beyond that conveyed by spike counts. The extra information available in the phase of firing was crucial for the disambiguation between stimuli eliciting high spike rates of similar magnitude. Thus, phase coding may allow primary cortical neurons to represent several effective stimuli in an easily decodable format.     

Millisecond-Scale Motor Encoding in a Cortical Vocal Area

Studies of motor control have almost universally examined firing rates to investigate how the brain shapes behavior. In principle, however, neurons could encode information through the precise temporal patterning of their spike trains as well as (or instead of) through their firing rates. Although the importance of spike timing has been demonstrated in sensory systems, it is largely unknown whether timing differences in motor areas could affect behavior. We tested the hypothesis that significant information about trial-by-trial variations in behavior is represented by spike timing in the songbird vocal motor system. We found that neurons in motor cortex convey information via spike timing far more often than via spike rate and that the amount of information conveyed at the millisecond timescale greatly exceeds the information available from spike counts. These results demonstrate that information can be represented by spike timing in motor circuits and suggest that timing variations evoke differences in behavior.     

Predicting spike timing of neocortical pyramidal neurons by simple threshold models

Neurons generate spikes reliably with millisecond precision if driven by a fluctuating current—is it then possible to predict the spike timing knowing the input? We determined parameters of an adapting threshold model using data recorded in vitro from 24 layer 5 pyramidal neurons from rat somatosensory cortex, stimulated intracellularly by a fluctuating current simulating synaptic bombardment in vivo. The model generates output spikes whenever the membrane voltage (a filtered version of the input current) reaches a dynamic threshold. We find that for input currents with large fluctuation amplitude, up to 75% of the spike times can be predicted with a precision of ±2 ms. Some of the intrinsic neuronal unreliability can be accounted for by a noisy threshold mechanism. Our results suggest that, under random current injection into the soma, (i) neuronal behavior in the subthreshold regime can be well approximated by a simple linear filter; and (ii) most of the nonlinearities are captured by a simple threshold process.     

Rhythm-induced spike-timing patterns characterized by 1D firing maps

We explore patterns in the spike timing of neurons receiving periodic inputs, with an emphasis on stable characteristics which are realized in both models and in-vitro whole-cell recordings. We report on whole-cell recordings of pyramidal CA1 cells from rat hippocampus and entorhinal cortex and compare this data to model simulations. Cells were injected with a constant current to induce a steady firing rate and then a modest rhythm was added which altered the spike times and their corresponding phases relative to the rhythm. For both experiment and theory the relationship between consecutive spike phases is characterized by a probability distribution with peaks concentrated near a one-dimensional firing map. As is well-known, stable fixed points of this map correspond to the neuron phase-locking to the rhythm. We show that the interaction between noise and sufficiently steep maps can also cause a new kind of spike-time organization, in which consecutive spike time pairs organize into discrete clusters, with transitions between these clusters proceeding in a fixed sequence. This structure is not just a vestige of the noise-free dynamics. This slow dynamics and temporal organization in the relationship between consecutive spike phases is not evident in either the neuron’s voltage traces or single phase or interspike interval histograms. Furthermore, the consecutive spike relationship is also evident in consecutive ISIs, and hence this ordering can be observed without detailed knowledge of the rhythm (e.g. without concurrent LFP recordings).     

The role of individual spikes and spike patterns in population coding of stimulus location in rat somatosensory cortex

This report addresses the nature of population coding in sensory cortex by applying information theoretic analysis to data recorded simultaneously from neuron pairs located in primary somatosensory cortex of anaesthetised rats. We studied how cortical spike trains code for the location of a whisker stimulus on the rat's snout. We found that substantially more information was conveyed by 10 ms precision spike timing compared with that conveyed by the number of spikes counted over a 40 ms response interval. Most of this information was accounted for by the timing of individual spikes. In particular, it was the first post-stimulus spikes that were crucial. Spike patterns within individual cells played a smaller role; spike patterns across cells were negligible. This pattern of results was robust both to the exact nature of the stimulus set and to the precision at which spikes were binned.     

DNA vaccination induced protective immunity against SARS CoV-2 infection in hamsterss

The development of efficient vaccines against COVID-19 is an emergent need for global public health. The spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a major target for the COVID-19 vaccine. To quickly respond to the outbreak of the SARS-CoV-2 pandemic, a nucleic acid-based vaccine is a novel option, beyond the traditional inactivated virus vaccine or recombinant protein vaccine. Here, we report a DNA vaccine containing the spike gene for delivery via electroporation. The spike genes of SARS-CoV and SARS-CoV-2 were codon optimized for mammalian cell expression and then cloned into mammalian cell expression vectors, called pSARS-S and pSARS2-S, respectively. Spike protein expression was confirmed by immunoblotting after transient expression in HEK293T cells. After immunization, sera were collected for antigen-specific antibody and neutralizing antibody titer analyses. We found that both pSARS-S and pSARS2-S immunization induced similar levels of antibodies against S2 of SARS-CoV-2. In contrast, only pSARS2-S immunization induced antibodies against the receptor-binding domain of SARS-CoV-2. We further found that pSARS2-S immunization, but not pSARS-S immunization, could induce very high titers of neutralizing antibodies against SARS-CoV-2. We further analyzed SARS-CoV-2 S protein-specific T cell responses and found that the immune responses were biased toward Th1. Importantly, pSARS2-S immunization in hamsters could induce protective immunity against SARS-CoV-2 challenge in vivo. These data suggest that DNA vaccination could be a promising approach for protecting against COVID-19.     

The Possible Role of Spike Patterns in Cortical Information Processing

When the same visual stimulus is presented across many trials, neurons in the visual cortex receive stimulus-related synaptic inputs that are reproducible across trials (S) and inputs that are not (N). The variability of spike trains recorded in the visual cortex and their apparent lack of spike-to-spike correlations beyond that implied by firing rate fluctuations, has been taken as evidence for a low S/N ratio. A recent re-analysis of in vivo cortical data revealed evidence for spike-to-spike correlations in the form of spike patterns. We examine neural dynamics at a higher S/N in order to determine what possible role spike patterns could play in cortical information processing. In vivo-like spike patterns were obtained in model simulations. Superpositions of multiple sinusoidal driving currents were especially effective in producing stable long-lasting patterns. By applying current pulses that were either short and strong or long and weak, neurons could be made to switch from one pattern to another. Cortical neurons with similar stimulus preferences are located near each other, have similar biophysical properties and receive a large number of common synaptic inputs. Hence, recordings of a single neuron across multiple trials are usually interpreted as the response of an ensemble of these neurons during one trial. In the presence of distinct spike patterns across trials there is ambiguity in what would be the corresponding ensemble, it could consist of the same spike pattern for each neuron or a set of patterns across neurons. We found that the spiking response of a neuron receiving these ensemble inputs was determined by the spike-pattern composition, which, in turn, could be modulated dynamically as a means for cortical information processing.     

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.     

Akt1 interacts with epidermal growth factor receptors and hedgehog signaling to increase stem/transit amplifying cells in the embryonic mouse cortex

A subset of precursors in the embryonic mouse cortex and in neurospheres expresses a higher level of the serine/threonine kinase Akt1 than neighboring precursors. We reported previously that the functional significance of high Akt1 expression was enhanced Akt1 activity, resulting in an increase in survival, proliferation, and self-renewal of multipotent stem/transit amplifying cells. Akt1 can interact with a number of signaling pathways, but the extrinsic factors that are required for specific effects of elevated Akt1 expression have not been identified. In this study we addressed the contributions of signaling via epidermal growth factor (EGF) and hedgehog (Hh) receptors. In EGF receptor-null precursors or following transient inhibition of EGF receptor tyrosine kinase activity, elevating Akt1 by retroviral transduction could still increase survival and proliferation but could not increase self-renewal. We also found that elevated Akt1 expression induced the expression of EGF receptors (EGFRs) in wild-type precursors. Several extrinsic factors, including Shh, can induce EGFR expression by cortical precursors, and we found that elevating Akt1 allowed them to respond to a subthreshold concentration of Shh to induce EGFRs. In precursors that lack the Hh receptor smoothened, however, elevating Akt1 did not increase EGFR expression or self-renewal, though it could still stimulate proliferation. These findings suggest that a subset of precursors in the embryonic cortex that express an elevated level of Akt1 can respond to lower concentrations of Shh than neighboring precursors, resulting in an increase in their expression of EGFRs. Signaling via EGFRs is required for their self-renewal.     

Electroreceptor mechanisms of the ampullae of lorenzini in skates

Responses of the receptor epithelium of single electrically isolated ampullae of Lorenzini and spike responses of nerve fibers connected to them to electrical stimulation under voltage clamping conditions were studied in experiments on the Black Sea skateRaja clavata. The preparations had an input resistance of 200–800 KΩ, a transepithelial resting potential of between 0 and ?2 mV, and the usual spontaneous spike activity in their afferent fibers. Thresholds of the electroreceptors were 2–10 µV (current about 10^?11 A). Within the working range of the electroreceptors (up to ±500 µV, current up to 1 nA) the current-voltage characteristic curve of the epithelium was linear without any evidence of spike or wave activity of the receptor cells. With negative currents of over 1–10 nA a regenerative spike appeared on the epithelium and was accompanied by an uncharacteristic pattern of spike discharge of the nerve fiber. It is concluded that, contrary to Bennett's hypothesis, spike or oscillatory activity bears no relationship to normal working of electroreceptors. It is postulated that "secondarily sensitive" receptor cells share a common functional organization, which is based on a chemical synapse with high electrical sensitivity.     

Neurotoxic action of some antiarrhythmic agents: comparative effects of propafenone,lidocaine and amiodarone on leech Retzius nerve cell

A series of antiarrhythmic drugs was studied on spontaneous spike activity and depolarizing outward potassium current in leech Retzius nerve cells. Propafenone (0.7 μM/ml) produced a cardiac-like action potential with a rapid depolarization followed by a sustained depolarization or plateau, which is terminated after 250 msec by a rapid repolarization. The effect of lidocaine (0.7 μM/ml) on spontaneous spike activity was less pronounced, and early afterdepolarization has been recorded. Amiodarone at the same and much higher concentrations (3 μM/ml) did not generate either a cardiac-like action potential or an early afterdepolarization. In the voltage clamp experiments, fast and slow calcium-activated outward potassium currents were suppressed with propafenone and lidocaine but not with amiodarone. These results suggest that the antiarrhythmic drugs, propafenone and lidocaine modulate calcium-activated potassium channels in leech Retzius nerve cells.     

Patterns of spontaneous unit preoptic neurosecretory cell discharges in the rainbow trout, Salmo gairdneri

Extracellular antidromic potentials recorded from the neurosecretory cell body were characterized by the following criteria: constant latency, the ability to follow a high frequency rate of stimulation and the collision test. The latency of the antidromic potentials ranged from 12 to 24 ms (17.46 +/- 3.10 SD) which gave a mean conduction velocity of 0.19 m/s, typical of unmyelinated nerve fibers. Two components could be clearly distinguished in the antidromic potential. A small "A" spike which showed constant latency and a large "B" spike with a variable latency and amplitude. A delay of 6.5 ms between the two spikes could occur and sometimes the "B" spike was blocked leaving only the "A" spike. Four patterns of spontaneous activity seem to emerge: Type I (26% of units, M +/- SD = 0.77 +/- 0.32 sp/s) corresponds to a slow and irregular pattern of activity; Type II (28% of units, M = 1.58 +/- 0.47 sp/s) is hard to classify and may be related to an irregular bursting pattern of activity; Type III (28% of units, M = 2.59 +/- 1.19 sp/s) corresponds to a continuous pattern of activity; Type IV (18% of units) represents a rhythmic pattern of activity with an active phase of about 3 min (M = 2.42 +/- 0.90 min), a silent phase of about 4 min (M = 3.89 +/- 3.02 min) and a maximal frequency of unit discharge in the range 2-18 sp/s. No statistical differences exist for the mean dorsal aortic pressure (DAP) between the four types of neurosecretory cell activity.     

Neuronal Spike Timing Adaptation Described with a Fractional Leaky Integrate-and-Fire Model

The voltage trace of neuronal activities can follow multiple timescale dynamics that arise from correlated membrane conductances. Such processes can result in power-law behavior in which the membrane voltage cannot be characterized with a single time constant. The emergent effect of these membrane correlations is a non-Markovian process that can be modeled with a fractional derivative. A fractional derivative is a non-local process in which the value of the variable is determined by integrating a temporal weighted voltage trace, also called the memory trace. Here we developed and analyzed a fractional leaky integrate-and-fire model in which the exponent of the fractional derivative can vary from 0 to 1, with 1 representing the normal derivative. As the exponent of the fractional derivative decreases, the weights of the voltage trace increase. Thus, the value of the voltage is increasingly correlated with the trajectory of the voltage in the past. By varying only the fractional exponent, our model can reproduce upward and downward spike adaptations found experimentally in neocortical pyramidal cells and tectal neurons in vitro. The model also produces spikes with longer first-spike latency and high inter-spike variability with power-law distribution. We further analyze spike adaptation and the responses to noisy and oscillatory input. The fractional model generates reliable spike patterns in response to noisy input. Overall, the spiking activity of the fractional leaky integrate-and-fire model deviates from the spiking activity of the Markovian model and reflects the temporal accumulated intrinsic membrane dynamics that affect the response of the neuron to external stimulation.     

Regional Difference in Responsiveness of Norepinephrine-Sensitive Cyclic AMP-Generating Systems of Rat Cerebral Cortex with Iron-Induced Epileptic Activity

Responsiveness of norepinephrine (NE)-sensitive cyclic AMP (cAMP)-generating systems was determined in slices from different areas of the rat cerebral cortex in which FeCl2 solution was injected unilaterally into the sensorimotor cortex to induce epileptic activity. In anterior cortical areas of rats in which the appearance of electrographic isolated spikes was dominant either ipsilaterally or contralaterally to the injection site 8-10 days after the injection, the cAMP accumulations elicited by NE and an NE-phentolamine combination were greater on the side of dominant spike activity than on the other. In anterior cortical areas of rats showing dominant spike activity on either side of the cortex 31-60 days after the injection, the cAMP accumulation elicited by NE was smaller on the dominant side than on the other. In anterior cortical areas of rats showing nearly equal spike activity on the two sides 31-60 days after the injection, the cAMP accumulations elicited by NE and an NE-phentolamine combination were greater on the side ipsilateral to the injection site than on the other. In anterior and posterior cortical areas of rats in which the appearance of spike and wave complexes, as well as isolated spikes, was detected 31-60 days after the injection, the cAMP accumulations elicited by NE and combinations of NE and phentolamine or propranolol were greater on the side ipsilateral to the injection site than on the other. The elicitation by an NE-propranolol combination, but not by an NE-phentolamine combination, of cAMP accumulation was almost completely inhibited by 8-phenyltheophylline.(ABSTRACT TRUNCATED AT 250 WORDS)     

Dynamics of action potential initiation in the GABAergic thalamic reticular nucleus in vivo

Understanding the neural mechanisms of action potential generation is critical to establish the way neural circuits generate and coordinate activity. Accordingly, we investigated the dynamics of action potential initiation in the GABAergic thalamic reticular nucleus (TRN) using in vivo intracellular recordings in cats in order to preserve anatomically-intact axo-dendritic distributions and naturally-occurring spatiotemporal patterns of synaptic activity in this structure that regulates the thalamic relay to neocortex. We found a wide operational range of voltage thresholds for action potentials, mostly due to intrinsic voltage-gated conductances and not synaptic activity driven by network oscillations. Varying levels of synchronous synaptic inputs produced fast rates of membrane potential depolarization preceding the action potential onset that were associated with lower thresholds and increased excitability, consistent with TRN neurons performing as coincidence detectors. On the other hand the presence of action potentials preceding any given spike was associated with more depolarized thresholds. The phase-plane trajectory of the action potential showed somato-dendritic propagation, but no obvious axon initial segment component, prominent in other neuronal classes and allegedly responsible for the high onset speed. Overall, our results suggest that TRN neurons could flexibly integrate synaptic inputs to discharge action potentials over wide voltage ranges, and perform as coincidence detectors and temporal integrators, supported by a dynamic action potential threshold.     

Design of recombinant protein-based SARS-CoV entry inhibitors targeting the heptad-repeat regions of the spike protein S2 domain

Entry of SARS-CoV into a target cell is initiated by binding of the S1 domain of spike protein to a receptor, followed by conformational changes of the spike protein S2 domain, resulting in the formation of a six-helix bundle by the heptad-repeat (HR1 and HR2) regions. Our previous studies have demonstrated that peptides derived from HR2 region could inhibit SARS-CoV entry. However, synthesis of these peptides is at high cost. In this study, we designed two recombinant proteins, one containing two HR1 and one HR2 peptides (denoted HR121), and the other consisting of two HR2 and one HR1 peptides (designated HR212). These two proteins could be easily purified with the low cost of production, exhibiting high stability and potent inhibitory activity on entry of the HIV/SARS pseudoviruses with IC(50) values of 4.13 and 0.95muM, respectively. These features suggest that HR121 and HR212 can serve as potent inhibitors of SARS-CoV entry.