Globus pallidus (GP) neurons recorded in brain slices show significant variability in intrinsic electrophysiological properties. To investigate how this variability arises, we manipulated the biophysical properties of GP neurons using computer simulations. Specifically, we created a GP neuron model database with 100,602 models that had varying densities of nine membrane conductances centered on a hand-tuned model that replicated typical physiological data. To test the hypothesis that the experimentally observed variability can be attributed to variations in conductance densities, we compared our model database results to a physiology database of 146 slice recordings. The electrophysiological properties of generated models and recordings were assessed with identical current injection protocols and analyzed with a uniform set of measures, allowing a systematic analysis of the effects of varying voltage-gated and calcium-gated conductance densities on the measured properties and a detailed comparison between models and recordings. Our results indicated that most of the experimental variability could be matched by varying conductance densities, which we confirmed with additional partial block experiments. Further analysis resulted in two key observations: (1) each voltage-gated conductance had effects on multiple measures such as action potential waveform and spontaneous or stimulated spike rates; and (2) the effect of each conductance was highly dependent on the background context of other conductances present. In some cases, such interactions could reverse the effect of the density of one conductance on important excitability measures. This context dependence of conductance density effects is important to understand drug and neuromodulator effects that work by affecting ion channels.

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.

Neuronal recordings and computer simulations produce ever growing amounts of data, impeding conventional analysis methods from keeping pace. Such large datasets can be automatically analyzed by taking advantage of the well-established relational database paradigm. Raw electrophysiology data can be entered into a database by extracting its interesting characteristics (e.g., firing rate). Compared to storing the raw data directly, this database representation is several orders of magnitude higher efficient in storage space and processing time. Using two large electrophysiology recording and simulation datasets, we demonstrate that the database can be queried, transformed and analyzed. This process is relatively simple and easy to learn because it takes place entirely in Matlab, using our database analysis toolbox, PANDORA. It is capable of acquiring data from common recording and simulation platforms and exchanging data with external database engines and other analysis toolboxes, which make analysis simpler and highly interoperable. PANDORA is available to be freely used and modified because it is open-source (http://software.incf.org/software/pandora/home).

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.

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 (K v1.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.

The globus pallidus (GP) predominantly contains GABAergic projection neurons that occupy a central position in the indirect pathway of the basal ganglia. They have long dendrites that can extend through one-half the diameter of the GP in rats, potentially enabling convergence and interaction between segregated basal ganglia circuits. Because of the length and fine diameter of GP dendrites, however, it is unclear how much influence distal synapses have on spiking activity. Dendritic expression of fast voltage-dependent Na ϩ channels (NaF channels) can enhance the importance of distal excitatory synapses by allowing for dendritic spike initiation and by subthreshold boosting of EPSPs. Antibody labeling has demonstrated the presence of NaF channel proteins in GP dendrites, but the quantitative expression density of the channels remains unknown. We built a series of nine GP neuron models that differed only in their dendritic NaF channel expression level to assess the functional impact of this parameter. The models were all similar in their basic electrophysiological features; however, higher expression levels of dendritic NaF channels increased the relative effectiveness of distal inputs for both excitatory and inhibitory synapses, broadening the effective extent of the dendritic tree. Higher dendritic NaF channel expression also made the neurons more resistant to tonic inhibition and highly sensitive to clustered synchronous excitation. The dendritic NaF channel expression pattern may therefore be a critical determinant of convergence for both the striatopallidal and subthalamopallidal projections, while also dictating which spatiotemporal input patterns are most effective at driving GP neuron output.

No abstract

Rhythmic behaviors vary across individuals. We investigated the sources of this output variability across a motor system, from the central pattern generator (CPG) to the motor plant. In the bilaterally symmetric leech heartbeat system, the CPG orchestrates two coordinations in the bilateral hearts with different intersegmental phase relations (Δϕ) and periodic side-to-side switches. Population variability is large. We show that the system is precise within a coordination, that differences in repetitions of a coordination contribute little to population output variability, but that differences between bilaterally homologous cells may contribute to some of this variability. Nevertheless, much output variability is likely associated with genetic and life history differences among individuals. Variability of Δϕ were coordination-specific: similar at all levels in one, but significantly lower for the motor pattern than the CPG pattern in the other. Mechanisms that transform CPG output to motor neurons may limit output variability in the motor pattern.