On the Role of Theory and Modeling in Neuroscience.

In recent years, the field of neuroscience has gone through rapid experimental advances and a significant increase in the use of quantitative and computational methods. This growth has created a need for clearer analyses of the theory and modeling approaches used in the field. This issue is particularly complex in neuroscience because the field studies phenomena that cross a wide range of scales and often require consideration at varying degrees of abstraction, from precise biophysical interactions to the computations they implement. We argue that a pragmatic perspective of science, in which descriptive, mechanistic, and normative models and theories each play a distinct role in defining and bridging levels of abstraction, will facilitate neuroscientific practice. This analysis leads to methodological suggestions, including selecting a level of abstraction that is appropriate for a given problem, identifying transfer functions to connect models and data, and the use of models themselves as a form of experiment.

Could simplified stimuli change how the brain performs visual search tasks? A deep neural network study.

Visual search is a complex behavior influenced by many factors. To control for these factors, many studies use highly simplified stimuli. However, the statistics of these stimuli are very different from the statistics of the natural images that the human visual system is optimized by evolution and experience to perceive. Could this difference change search behavior? If so, simplified stimuli may contribute to effects typically attributed to cognitive processes, such as selective attention. Here we use deep neural networks to test how optimizing models for the statistics of one distribution of images constrains performance on a task using images from a different distribution. We train four deep neural network architectures on one of three source datasets-natural images, faces, and x-ray images-and then adapt them to a visual search task using simplified stimuli. This adaptation produces models that exhibit performance limitations similar to humans, whereas models trained on the search task alone exhibit no such limitations. However, we also find that deep neural networks trained to classify natural images exhibit similar limitations when adapted to a search task that uses a different set of natural images. Therefore, the distribution of data alone cannot explain this effect. We discuss how future work might integrate an optimization-based approach into existing models of visual search behavior.

Distal spike initiation zone location estimation by morphological simulation of ionic current filtering demonstrated in a novel model of an identified Drosophila motoneuron.

Studying ion channel currents generated distally from the recording site is difficult because of artifacts caused by poor space clamp and membrane filtering. A computational model can quantify artifact parameters for correction by simulating the currents only if their exact anatomical location is known. We propose that the same artifacts that confound current recordings can help pinpoint the source of those currents by providing a signature of the neuron's morphology. This method can improve the recording quality of currents initiated at the spike initiation zone (SIZ) that are often distal to the soma in invertebrate neurons. Drosophila being a valuable tool for characterizing ion currents, we estimated the SIZ location and quantified artifacts in an identified motoneuron, aCC/MN1-Ib, by constructing a novel multicompartmental model. Initial simulation of the measured biophysical channel properties in an isopotential Hodgkin-Huxley type neuron model partially replicated firing characteristics. Adding a second distal compartment, which contained spike-generating Na+ and K+ currents, was sufficient to simulate aCC's in vivo activity signature. Matching this signature using a reconstructed morphology predicted that the SIZ is on aCC's primary axon, 70 µm after the most distal dendritic branching point. From SIZ to soma, we observed and quantified selective morphological filtering of fast activating currents. Non-inactivating K+ currents are filtered ∼3 times less and despite their large magnitude at the soma they could be as distal as Na+ currents. The peak of transient component (NaT) of the voltage-activated Na+ current is also filtered more than the magnitude of slower persistent component (NaP), which can contribute to seizures. The corrected NaP/NaT ratio explains the previously observed discrepancy when the same channel is expressed in different cells. In summary, we used an in vivo signature to estimate ion channel location and recording artifacts, which can be applied to other neurons.

The transcription factors islet and Lim3 combinatorially regulate ion channel gene expression.

Expression of appropriate ion channels is essential to allow developing neurons to form functional networks. Our previous studies have identified LIM-homeodomain (HD) transcription factors (TFs), expressed by developing neurons, that are specifically able to regulate ion channel gene expression. In this study, we use the technique of DNA adenine methyltransferase identification (DamID) to identify putative gene targets of four such TFs that are differentially expressed in Drosophila motoneurons. Analysis of targets for Islet (Isl), Lim3, Hb9, and Even-skipped (Eve) identifies both ion channel genes and genes predicted to regulate aspects of dendritic and axonal morphology. Significantly, some ion channel genes are bound by more than one TF, consistent with the possibility of combinatorial regulation. One such gene is Shaker (Sh), which encodes a voltage-dependent fast K(+) channel (Kv1.1). DamID reveals that Sh is bound by both Isl and Lim3. We used body wall muscle as a test tissue because in conditions of low Ca(2+), the fast K(+) current is carried solely by Sh channels (unlike neurons in which a second fast K(+) current, Shal, also contributes). Ectopic expression of isl, but not Lim3, is sufficient to reduce both Sh transcript and Sh current level. By contrast, coexpression of both TFs is additive, resulting in a significantly greater reduction in both Sh transcript and current compared with isl expression alone. These observations provide evidence for combinatorial activity of Isl and Lim3 in regulating ion channel gene expression.

Feedback control of variability in the cycle period of a central pattern generator.

We address how feedback to a bursting biological pacemaker with intrinsic variability in cycle length can affect that variability. Specifically, we examine a hybrid circuit constructed of an isolated crab anterior burster (AB)/pyloric dilator (PD) pyloric pacemaker receiving virtual feedback via dynamic clamp. This virtual feedback generates artificial synaptic input to PD with timing determined by adjustable phase response dynamics that mimic average burst intervals generated by the lateral pyloric neuron (LP) in the intact pyloric network. Using this system, we measure network period variability dependence on the feedback element's phase response dynamics and find that a constant response interval confers minimum variability. We further find that these optimal dynamics are characteristic of the biological pyloric network. Building upon our previous theoretical work mapping the firing intervals in one cycle onto the firing intervals in the next cycle, we create a theoretical map of the distribution of all firing intervals in one cycle to the distribution of firing intervals in the next cycle. We then obtain an integral equation for a stationary self-consistent distribution of the network periods of the hybrid circuit, which can be solved numerically given the uncoupled pacemaker's distribution of intrinsic periods, the nature of the network's feedback, and the phase resetting characteristics of the pacemaker. The stationary distributions obtained in this manner are strongly predictive of the experimentally observed distributions of hybrid network period. This theoretical framework can provide insight into optimal feedback schemes for minimizing variability to increase reliability or maximizing variability to increase flexibility in central pattern generators driven by pacemakers with feedback.

Degradation of extracellular chondroitin sulfate delays recovery of network activity after perturbation.

Chondroitin sulfate proteoglycans (CSPGs) are widely studied in vertebrate systems and are known to play a key role in development, plasticity, and regulation of cortical circuitry. The mechanistic details of this role are still elusive, but increasingly central to the investigation is the homeostatic balance between network excitation and inhibition. Studying a simpler neuronal circuit may prove advantageous for discovering the mechanistic details of the cellular effects of CSPGs. In this study we used a well-established model of homeostatic change after injury in the crab Cancer borealis to show first evidence that CSPGs are necessary for network activity homeostasis. We degraded CSPGs in the pyloric circuit of the stomatogastric ganglion with the enzyme chondroitinase ABC (chABC) and found that removal of CSPGs does not influence the ongoing rhythm of the pyloric circuit but does limit its capacity for recovery after a networkwide perturbation. Without CSPGs, the postperturbation rhythm is slower than in controls and rhythm recovery is delayed. In addition to providing a new model system for the study of CSPGs, this study suggests a wider role for CSPGs, and perhaps the extracellular matrix in general, beyond simply plastic reorganization (as observed in mammals) and into a foundational regulatory role of neural circuitry.

Differential effects of conductances on the phase resetting curve of a bursting neuronal oscillator.

The intrinsically oscillating neurons in the crustacean pyloric circuit have membrane conductances that influence their spontaneous activity patterns and responses to synaptic activity. The relationship between the magnitudes of these membrane conductances and the response of the oscillating neurons to synaptic input has not yet been fully or systematically explored. We examined this relationship using the phase resetting curve (PRC), which summarizes the change in the cycle period of a neuronal oscillator as a function of the input's timing within the oscillation. We first utilized a large database of single-compartment model neurons to determine the effect of individual membrane conductances on PRC shape; we found that the effects vary across conductance space, but on average, the hyperpolarization-activated and leak conductances advance the PRC. We next investigated how membrane conductances affect PRCs of the isolated pacemaker kernel in the pyloric circuit of Cancer borealis by: (1) tabulating PRCs while using dynamic clamp to artificially add varying levels of specific conductances, and (2) tabulating PRCs before and after blocking the endogenous hyperpolarization-activated current. We additionally used a previously described four-compartment model to determine how the location of the hyperpolarization-activated conductance influences that current's effect on the PRC. We report that while dynamic-clamp-injected leak current has much smaller effects on the PRC than suggested by the single-compartment model, an increase in the hyperpolarization-activated conductance both advances and reduces the noisiness of the PRC in the pacemaker kernel of the pyloric circuit in both modeling and experimental studies.

Phase maintenance in a rhythmic motor pattern during temperature changes in vivo.

Central-pattern-generating neural circuits function reliably throughout an animal's life, despite constant molecular turnover and environmental perturbations. Fluctuations in temperature pose a problem to the nervous systems of poikilotherms because their body temperature follows the ambient temperature, thus affecting the temperature-dependent dynamics of various subcellular components that constitute neuronal circuits. In the crustacean stomatogastric nervous system, the pyloric circuit produces a triphasic rhythm comprising the output of the pyloric dilator, lateral pyloric, and pyloric constrictor neurons. In vitro, the phase relationships of these neurons are maintained over a fourfold change in pyloric frequency as temperature increases from 7°C to 23°C. To determine whether these temperature effects are also found in intact crabs, in the presence of sensory feedback and neuromodulator-rich environments, we measured the temperature dependence of the pyloric frequency and phases in vivo by implanting extracellular electrodes into Cancer borealis and Cancer pagurus and shifting tank water temperature from 11°C to 26°C. Pyloric frequency in the intact crab increased significantly with temperature (Q10 = 2-2.5), while pyloric phases were generally conserved. For a subset of the C. borealis experiments, animals were subsequently dissected and the stomatogastric ganglion subjected to a similar temperature ramp in vitro. We found that the maximal frequency attained at high temperatures in vivo is lower than it is under in vitro conditions. Our results demonstrate that, over a wide temperature range, the phases of the pyloric rhythm in vivo are generally preserved, but that the frequency range is more restricted than it is in vitro.

Activity-dependent alternative splicing increases persistent sodium current and promotes seizure.

Activity of voltage-gated Na channels (Na(v)) is modified by alternative splicing. However, whether altered splicing of human Na(v)s contributes to epilepsy remains to be conclusively shown. We show here that altered splicing of the Drosophila Na(v) (paralytic, DmNa(v)) contributes to seizure-like behavior in identified seizure mutants. We focus attention on a pair of mutually exclusive alternate exons (termed K and L), which form part of the voltage sensor (S4) in domain III of the expressed channel. The presence of exon L results in a large, non-inactivating, persistent I(Nap). Many forms of human epilepsy are associated with an increase in this current. In wild-type (WT) Drosophila larvae, ∼70-80% of DmNa(v) transcripts contain exon L, and the remainder contain exon K. Splicing of DmNa(v) to include exon L is increased to ∼100% in both the slamdance and easily-shocked seizure mutants. This change to splicing is prevented by reducing synaptic activity levels through exposure to the antiepileptic phenytoin or the inhibitory transmitter GABA. Conversely, enhancing synaptic activity in WT, by feeding of picrotoxin is sufficient to increase I(Nap) and promote seizure through increased inclusion of exon L to 100%. We also show that the underlying activity-dependent mechanism requires the presence of Pasilla, an RNA-binding protein. Finally, we use computational modeling to show that increasing I(Nap) is sufficient to potentiate membrane excitability consistent with a seizure phenotype. Thus, increased synaptic excitation favors inclusion of exon L, which, in turn, further increases neuronal excitability. Thus, at least in Drosophila, this self-reinforcing cycle may promote the incidence of seizure.

Modifying spiking precision in conductance-based neuronal models.

The temporal precision of a neuron's spiking can be characterized by calculating its "jitter," defined as the standard deviation of the timing of individual spikes in response to repeated presentations of a stimulus. Sub-millisecond jitters have been measured for neurons in a variety of experimental systems and appear to be functionally important in some instances. We have investigated how modifying a neuron's maximal conductances affects jitter using the leaky integrate-and-fire (LIF) model and an eight-conductance Hodgkin-Huxley type (HH8) model. We observed that jitter can be largely understood in the LIF model in terms of the neuron's filtering properties. In the HH8 model we found the role of individual conductances in determining jitter to be complicated and dependent on the model's spiking properties. Distinct behaviors were observed for populations with slow (<11.5 Hz) and fast (>11.5 Hz) spike rates and appear to be related to differences in a particular channel's activity at times just before spiking occurs.

Temporal dynamics of Na/K pump mediated memory traces: insights from conductance-based models of Drosophila neurons.

Sodium potassium ATPases (Na/K pumps) mediate long-lasting, dynamic cellular memories that can last tens of seconds. The mechanisms controlling the dynamics of this type of cellular memory are not well understood and can be counterintuitive. Here, we use computational modeling to examine how Na/K pumps and the ion concentration dynamics they influence shape cellular excitability. In a Drosophila larval motor neuron model, we incorporate a Na/K pump, a dynamic intracellular Na+ concentration, and a dynamic Na+ reversal potential. We probe neuronal excitability with a variety of stimuli, including step currents, ramp currents, and zap currents, then monitor the sub- and suprathreshold voltage responses on a range of time scales. We find that the interactions of a Na+-dependent pump current with a dynamic Na+ concentration and reversal potential endow the neuron with rich response properties that are absent when the role of the pump is reduced to the maintenance of constant ion concentration gradients. In particular, these dynamic pump-Na+ interactions contribute to spike rate adaptation and result in long-lasting excitability changes after spiking and even after sub-threshold voltage fluctuations on multiple time scales. We further show that modulation of pump properties can profoundly alter a neuron's spontaneous activity and response to stimuli by providing a mechanism for bursting oscillations. Our work has implications for experimental studies and computational modeling of the role of Na/K pumps in neuronal activity, information processing in neural circuits, and the neural control of animal behavior.

Co-variation of ionic conductances supports phase maintenance in stomatogastric neurons.

Neuronal networks produce reliable functional output throughout the lifespan of an animal despite ceaseless molecular turnover and a constantly changing environment. Central pattern generators, such as those of the crustacean stomatogastric ganglion (STG), are able to robustly maintain their functionality over a wide range of burst periods. Previous experimental work involving extracellular recordings of the pyloric pattern of the STG has demonstrated that as the burst period varies, the inter-neuronal delays are altered proportionally, resulting in burst phases that are roughly invariant. The question whether spike delays within bursts are also proportional to pyloric period has not been explored in detail. The mechanism by which the pyloric neurons accomplish phase maintenance is currently not obvious. Previous studies suggest that the co-regulation of certain ion channel properties may play a role in governing neuronal activity. Here, we observed in long-term recordings of the pyloric rhythm that spike delays can vary proportionally with burst period, so that spike phase is maintained. We then used a conductance-based model neuron to determine whether co-varying ionic membrane conductances results in neural output that emulates the experimentally observed phenomenon of spike phase maintenance. Next, we utilized a model neuron database to determine whether conductance correlations exist in model neuron populations with highly maintained spike phases. We found that co-varying certain conductances, including the sodium and transient calcium conductance pair, causes the model neuron to maintain a specific spike phase pattern. Results indicate a possible relationship between conductance co-regulation and phase maintenance in STG neurons.

Model calcium sensors for network homeostasis: sensor and readout parameter analysis from a database of model neuronal networks.

In activity-dependent homeostatic regulation (ADHR) of neuronal and network properties, the intracellular Ca(2+) concentration is a good candidate for sensing activity levels because it is correlated with the electrical activity of the cell. Previous ADHR models, developed with abstract activity sensors for model pyloric neurons and networks of the crustacean stomatogastric ganglion, showed that functional activity can be maintained by a regulation mechanism that senses activity levels solely from Ca(2+). At the same time, several intracellular pathways have been discovered for Ca(2+)-dependent regulation of ion channels. To generate testable predictions for dynamics of these signaling pathways, we undertook a parameter study of model Ca(2+) sensors across thousands of model pyloric networks. We found that an optimal regulation signal can be generated for 86% of model networks with a sensing mechanism that activates with a time constant of 1 ms and that inactivates within 1 s. The sensor performed robustly around this optimal point and did not need to be specific to the role of the cell. When multiple sensors with different time constants were used, coverage extended to 88% of the networks. Without changing the sensors, it extended to 95% of the networks by letting the sensors affect the readout nonlinearly. Specific to this pyloric network model, the sensor of the follower pyloric constrictor cell was more informative than the pacemaker anterior burster cell for producing a regulatory signal. Conversely, a global signal indicating network activity that was generated by summing the sensors in individual cells was less informative for regulation.

Identifiable cells in the crustacean stomatogastric ganglion.

Neural circuits rely on slight physiological differences between the component cells for proper function. When any circuit is analyzed, it is important to characterize the features that distinguish one cell type from another. This review describes the methods used to identify the neurons of the crustacean stomatogastric ganglion.

Responses of a bursting pacemaker to excitation reveal spatial segregation between bursting and spiking mechanisms.

Central pattern generators (CPGs) frequently include bursting neurons that serve as pacemakers for rhythm generation. Phase resetting curves (PRCs) can provide insight into mechanisms underlying phase locking in such circuits. PRCs were constructed for a pacemaker bursting complex in the pyloric circuit in the stomatogastric ganglion of the lobster and crab. This complex is comprised of the Anterior Burster (AB) neuron and two Pyloric Dilator (PD) neurons that are all electrically coupled. Artificial excitatory synaptic conductance pulses of different strengths and durations were injected into one of the AB or PD somata using the Dynamic Clamp. Previously, we characterized the inhibitory PRCs by assuming a single slow process that enabled synaptic inputs to trigger switches between an up state in which spiking occurs and a down state in which it does not. Excitation produced five different PRC shapes, which could not be explained with such a simple model. A separate dendritic compartment was required to separate the mechanism that generates the up and down phases of the bursting envelope (1) from synaptic inputs applied at the soma, (2) from axonal spike generation and (3) from a slow process with a slower time scale than burst generation. This study reveals that due to the nonlinear properties and compartmentalization of ionic channels, the response to excitation is more complex than inhibition.

Conductance ratios and cellular identity.

Recent experimental evidence suggests that coordinated expression of ion channels plays a role in constraining neuronal electrical activity. In particular, each neuronal cell type of the crustacean stomatogastric ganglion exhibits a unique set of positive linear correlations between ionic membrane conductances. These data suggest a causal relationship between expressed conductance correlations and features of cellular identity, namely electrical activity type. To test this idea, we used an existing database of conductance-based model neurons. We partitioned this database based on various measures of intrinsic activity, to approximate distinctions between biological cell types. We then tested individual conductance pairs for linear dependence to identify correlations. Contrary to experimental evidence, in which all conductance correlations are positive, 32% of correlations seen in this database were negative relationships. In addition, 80% of correlations seen here involved at least one calcium conductance, which have been difficult to measure experimentally. Similar to experimental results, each activity type investigated had a unique combination of correlated conductances. Finally, we found that populations of models that conform to a specific conductance correlation have a higher likelihood of exhibiting a particular feature of electrical activity. We conclude that regulating conductance ratios can support proper electrical activity of a wide range of cell types, particularly when the identity of the cell is well-defined by one or two features of its activity. Furthermore, we predict that previously unseen negative correlations and correlations involving calcium conductances are biologically plausible.

Geometry and dynamics of activity-dependent homeostatic regulation in neurons.

To maintain activity in a functional range, neurons constantly adjust membrane excitability to changing intra- and extracellular conditions. Such activity-dependent homeostatic regulation (ADHR) is critical for normal processing of the nervous system and avoiding pathological conditions. Here, we posed a homeostatic regulation problem for the classical Morris-Lecar (ML) model. The problem was motivated by the phenomenon of the functional recovery of stomatogastric neurons in crustaceans in the absence of neuromodulation. In our study, the regulation of the ionic conductances in the ML model depended on the calcium current or the intracellular calcium concentration. We found an asymptotic solution to the problem under the assumption of slow regulation. The solution provides a full account of the regulation in the case of correlated or anticorrelated changes of the maximal conductances of the calcium and potassium currents. In particular, the solution shows how the target and parameters of the regulation determine which perturbations of the conductances can be compensated by the ADHR. In some cases, the sets of compensated initial perturbations are not convex. On the basis of our analysis we formulated specific questions for subsequent experimental and theoretical studies of ADHR.

Phase resetting curves allow for simple and accurate prediction of robust N:1 phase locking for strongly coupled neural oscillators.

Existence and stability criteria for harmonic locking modes were derived for two reciprocally pulse coupled oscillators based on their first and second order phase resetting curves. Our theoretical methods are general in the sense that no assumptions about the strength of coupling, type of synaptic coupling, and model are made. These methods were then tested using two reciprocally inhibitory Wang and Buzsáki model neurons. The existence of bands of 2:1, 3:1, 4:1, and 5:1 phase locking in the relative frequency parameter space was predicted correctly, as was the phase of the slow neuron's spike within the cycle of the fast neuron in which it occurred. For weak coupling the bands are very narrow, but strong coupling broadens the bands. The predictions of the pulse coupled method agreed with weak coupling methods in the weak coupling regime, but extended predictability into the strong coupling regime. We show that our prediction method generalizes to pairs of neural oscillators coupled through excitatory synapses, and to networks of multiple oscillatory neurons. The main limitation of the method is the central assumption that the effect of each input dies out before the next input is received.

Current compensation in neuronal homeostasis.

How do neurons maintain stable intrinsic properties over long periods of time as the channels that govern excitability turn over in the membrane? In this issue of Neuron, MacLean et al. argue that homeostatic regulation of intrinsic activity can occur by an activity-independent mechanism.