Neuronal coupling by endogenous electric fields: cable theory and applications to coincidence detector neurons in the auditory brain stem

Thermal effects and ephaptic entrainment in Hodgkin-Huxley model

The brain is understood as an intricate biological system composed of numerous elements. It is susceptible to various physical and chemical influences, including temperature. The literature extensively explores the conditions that influence synapses in the context of cellular communication. However, the understanding of how the brain's global physical conditions can modulate ephaptic communication remains limited due to the poorly understood nature of ephapticity. This study proposes an adaptation of the Hodgkin and Huxley (HH) model to investigate the effects of ephaptic entrainment in response to thermal changes (HH-E). The analysis focuses on two distinct neuronal regimes: subthreshold and suprathreshold. In the subthreshold regime, circular statistics are used to demonstrate the dependence of phase differences with temperature. In the suprathreshold regime, the Inter-Spike Interval are employed to estimate phase preferences and changes in the spiking pattern. Temperature influences the model's ephaptic interactions and can modify its preferences for spiking frequency, with the direction of this change depending on specific model conditions and the temperature range under consideration. Furthermore, temperature enhance the anti-phase differences relationship between spikes and the external ephaptic signal. In the suprathreshold regime, ephaptic entrainment is also influenced by temperature, especially at low frequencies. This study reveals the susceptibility of ephaptic entrainment to temperature variations in both subthreshold and suprathreshold regimes and discusses the importance of ephaptic communication in the contexts where temperature may plays a significant role in neural physiology, such as inflammatory processes, fever, and epileptic seizures.

Impact of Extracellular Current Flow on Action Potential Propagation in Myelinated Axons

Myelinated axons conduct action potentials, or spikes, in a saltatory manner. Inward current caused by a spike occurring at one node of Ranvier spreads axially to the next node, which regenerates the spike when depolarized enough for voltage-gated sodium channels to activate, and so on. The rate at which this process progresses dictates the velocity at which the spike is conducted and depends on several factors including axial resistivity and axon diameter that directly affect axial current. Here we show through computational simulations in modified double-cable axon models that conduction velocity also depends on extracellular factors whose effects can be explained by their indirect influence on axial current. Specifically, we show that a conventional double-cable model, with its outside layer connected to ground, transmits less axial current than a model whose outside layer is less absorptive. A more resistive barrier exists when an axon is packed tightly between other myelinated fibers, for example. We show that realistically resistive boundary conditions can significantly increase the velocity and energy efficiency of spike propagation, while also protecting against propagation failure. Certain factors like myelin thickness may be less important than typically thought if extracellular conditions are more resistive than normally considered. We also show how realistically resistive boundary conditions affect ephaptic interactions. Overall, these results highlight the unappreciated importance of extracellular conditions for axon function.

Electrophysiological damage to neuronal membrane alters ephaptic entrainment

The brain is commonly understood as a complex network system with a particular organization and topology that can result in specific electrophysiological patterns. Among all the dynamic elements resulting from the circuits of the brain's network, ephapticity is a cellular communication mechanism that has received little attention. To understand the network's properties of ephaptic entrainment, we start investigating the ephaptic effect on a single neuron. In this study, we used numerical simulations to examine the relationship between alterations in ephaptic neuronal entrainment and impaired electrophysiological properties of the neuronal membrane, which can occur via spike field coherence (SFC). This change in frequency band amplitude is observed in some neurodegenerative diseases, such as Parkinson's or Alzheimer's. To further investigate these phenomena, we proposed a damaged model based on the impairment of both the resistance of the ion channels and the capacitance of the lipid membrane. Therefore, we simulated ephaptic entrainment with the hybrid neural model quadratic integrate-and-fire ephaptic (QIF-E), which mimics an ephaptic entrainment generated by an LFP (simulate a neuronal group). Our results indicate a link between peak entrainment (ephapticity) preference and a shift in frequency band when damage occurs mainly in ion channels. Finally, we discuss possible relationships between ephaptic entrainment and neurodegenerative diseases associated with aging factors.

Knösche T. Modelling the effect of ephaptic coupling on spike propagation in peripheral nerve fibres

Experimental and theoretical studies have shown that ephaptic coupling leads to the synchronisation and slowing down of spikes propagating along the axons within peripheral nerve bundles. However, the main focus thus far has been on a small number of identical axons, whereas realistic peripheral nerve bundles contain numerous axons with different diameters. Here, we present a computationally efficient spike propagation model, which captures the essential features of propagating spikes and their ephaptic interaction, and facilitates the theoretical investigation of spike volleys in large, heterogeneous fibre bundles. We first lay out the theoretical basis to describe how the spike in an active axon changes the membrane potential of a passive axon. These insights are then incorporated into the spike propagation model, which is calibrated with a biophysically realistic model based on Hodgkin-Huxley dynamics. The fully calibrated model is then applied to fibre bundles with a large number of axons and different types of axon diameter distributions. One key insight of this study is that the heterogeneity of the axonal diameters has a dispersive effect, and that a higher level of heterogeneity requires stronger ephaptic coupling to achieve full synchronisation between spikes.

Ephaptic entrainment in hybrid neuronal model

In recent decades, there has been a growing interest in the impact of electric fields generated in the brain. Transmembrane ionic currents originate electric fields in the extracellular space and are capable of affecting nearby neurons, a phenomenon called ephaptic neuronal communication. In the present work, the Quadratic Integrated-and-Fire model (QIF-E) underwent an adjustment/improvement to include the ephaptic entrainment behavior between neurons and electric fields. Indeed, the aim of our study is to validate the QIF-E model, which is a model to estimate the influence of electric fields on neurons. For this purpose, we evaluated whether the main properties observed in an experiment by Anastassiou et al. (Nat Neurosci 14:217-223, 2011), which analyzed the effect of an electric field on cortical pyramidal neurons, are reproduced with the QIF-E model. In this way, the analysis tools are employed according to the neuronal activity regime: (i) for the subthreshold regime, the circular statistic is used to describe the phase differences between the input stimulus signal (electrode) and the modeled membrane response; (ii) in the suprathreshold regime, the Population Vector and the Spike Field Coherence are used to estimate phase preferences and the entrainment intensity between the input stimulus and Action Potentials. The results observed are (i) in the subthreshold regime the values of the phase differences change with distinct frequencies of the input stimulus; (ii) in the supra-threshold regime the preferential phase of Action Potentials changes for different frequencies. In addition, we explore other parameters of the model, such as noise and membrane characteristic-time, in order to understand different types of neurons and extracellular environment related to ephaptic communication. Such results are consistent with results observed in empirical experiments based on ephaptic phenomenon. In addition, the QIF-E model allows further studies on the physiological importance of ephaptic communication in the brain, and its simplicity may open a door to simulate the ephaptic response in neuronal networks and assess the impact of ephaptic communication in such scenarios.

An electrodiffusive neuron-extracellular-glia model for exploring the genesis of slow potentials in the brain

Within the computational neuroscience community, there has been a focus on simulating the electrical activity of neurons, while other components of brain tissue, such as glia cells and the extracellular space, are often neglected. Standard models of extracellular potentials are based on a combination of multicompartmental models describing neural electrodynamics and volume conductor theory. Such models cannot be used to simulate the slow components of extracellular potentials, which depend on ion concentration dynamics, and the effect that this has on extracellular diffusion potentials and glial buffering currents. We here present the electrodiffusive neuron-extracellular-glia (edNEG) model, which we believe is the first model to combine compartmental neuron modeling with an electrodiffusive framework for intra- and extracellular ion concentration dynamics in a local piece of neuro-glial brain tissue. The edNEG model (i) keeps track of all intraneuronal, intraglial, and extracellular ion concentrations and electrical potentials, (ii) accounts for action potentials and dendritic calcium spikes in neurons, (iii) contains a neuronal and glial homeostatic machinery that gives physiologically realistic ion concentration dynamics, (iv) accounts for electrodiffusive transmembrane, intracellular, and extracellular ionic movements, and (v) accounts for glial and neuronal swelling caused by osmotic transmembrane pressure gradients. The edNEG model accounts for the concentration-dependent effects on ECS potentials that the standard models neglect. Using the edNEG model, we analyze these effects by splitting the extracellular potential into three components: one due to neural sink/source configurations, one due to glial sink/source configurations, and one due to extracellular diffusive currents. Through a series of simulations, we analyze the roles played by the various components and how they interact in generating the total slow potential. We conclude that the three components are of comparable magnitude and that the stimulus conditions determine which of the components that dominate.

Ephaptic coupling in white matter fibre bundles modulates axonal transmission delays

Axonal connections are widely regarded as faithful transmitters of neuronal signals with fixed delays. The reasoning behind this is that extracellular potentials caused by spikes travelling along axons are too small to have an effect on other axons. Here we devise a computational framework that allows us to study the effect of extracellular potentials generated by spike volleys in axonal fibre bundles on axonal transmission delays. We demonstrate that, although the extracellular potentials generated by single spikes are of the order of microvolts, the collective extracellular potential generated by spike volleys can reach several millivolts. As a consequence, the resulting depolarisation of the axonal membranes increases the velocity of spikes, and therefore reduces axonal delays between brain areas. Driving a neural mass model with such spike volleys, we further demonstrate that only ephaptic coupling can explain the reduction of stimulus latencies with increased stimulus intensities, as observed in many psychological experiments.

Axonal fibre bundles that connect distant cortical areas contain millions of densely packed axons. When synchronous spike volleys travel through such fibre bundles, the extracellular potential within the bundles is perturbed. We use computer simulations to examine the magnitude and shape of this perturbation, and demonstrate that it is sufficiently strong to affect axonal transmission speeds. Since most spikes within a spike volley are positioned in an area where the extracellular potential is negative (relative to a distant reference), the resulting depolarisation of the axonal membranes accelerates the spike volley on average. This finding is in contrast to previous studies of ephaptic coupling effects between axons, where ephaptic coupling was found to slow down spike propagation. Our finding has consequences for information transmission and synchronisation between cortical areas.

A mathematical model of ephaptic interactions in neuronal fiber pathways: Could there be more than transmission along the tracts?

While numerous studies of ephaptic interactions have focused on either axons of peripheral nerves or on cortical structures, no attention has been given to the possibility of ephaptic interactions in white matter tracts. Inspired by the highly organized, tightly packed geometry of axons in fiber pathways, we aim to investigate the potential effects of ephaptic interactions along these structures that are resilient to experimental probing. We use axonal cable theory to derive a minimal model of a sheet of N ephaptically coupled axons. Numerical solutions of the proposed model are explored as ephaptic coupling is varied. We demonstrate that ephaptic interactions can lead to local phase locking between adjacent traveling impulses and that, as coupling is increased, traveling impulses trigger new impulses along adjacent axons, resulting in finite size traveling fronts. For strong enough coupling, impulses propagate laterally and backwards, resulting in complex spatiotemporal patterns. While common large-scale brain network models often model fiber pathways as simple relays of signals between different brain regions, our work calls for a closer reexamination of the validity of such a view. The results suggest that in the presence of significant ephaptic interactions, the brain fiber tracts can act as a dynamic active medium.

Starting from the FitzHugh-Nagumo cable model, we derive a system of nonlinear coupled partial differential equations (PDEs) to model a sheet of N ephaptically coupled axons. We also present a continuous limit approximation transforming the model into a two-dimensional field equation. We numerically solve the equations exploring the dynamics as coupling strength is varied. We observe phase locking of adjacent impulses and coordination of subthreshold dynamics. Strong enough coupling generates complex spatiotemporal patterns as new impulses form traveling fronts propagating laterally and backwards. The transition between different dynamic regimes happens abruptly at critical values of parameter. The results put into question the validity of assuming the role of fiber pathways to be that of mere interneuronal transmission and call for further investigation of the matter.

Modelling the effects of ephaptic coupling on selectivity and response patterns during artificial stimulation of peripheral nerves

Artificial electrical stimulation of peripheral nerves for sensory feedback restoration can greatly benefit from computational models for simulation-based neural implant design in order to reduce the trial-and-error approach usually taken, thus potentially significantly reducing research and development costs and time. To this end, we built a computational model of a peripheral nerve trunk in which the interstitial space between the fibers and the tissues was modelled using a resistor network, thus enabling distance-dependent ephaptic coupling between myelinated axons and between fascicles as well. We used the model to simulate a) the stimulation of a nerve trunk model with a cuff electrode, and b) the propagation of action potentials along the axons. Results were used to investigate the effect of ephaptic interactions on recruitment and selectivity stemming from artificial (i.e., neural implant) stimulation and on the relative timing between action potentials during propagation. Ephaptic coupling was found to increase the number of fibers that are activated by artificial stimulation, thus reducing the artificial currents required for axonal recruitment, and it was found to reduce and shift the range of optimal stimulation amplitudes for maximum inter-fascicular selectivity. During propagation, while fibers of similar diameters tended to lock their action potentials and reduce their conduction velocities, as expected from previous knowledge on bundles of identical axons, the presence of many other fibers of different diameters was found to make their interactions weaker and unstable.

Finite Element Simulation of Ionic Electrodiffusion in Cellular Geometries

Mathematical models for excitable cells are commonly based on cable theory, which considers a homogenized domain and spatially constant ionic concentrations. Although such models provide valuable insight, the effect of altered ion concentrations or detailed cell morphology on the electrical potentials cannot be captured. In this paper, we discuss an alternative approach to detailed modeling of electrodiffusion in neural tissue. The mathematical model describes the distribution and evolution of ion concentrations in a geometrically-explicit representation of the intra- and extracellular domains. As a combination of the electroneutral Kirchhoff-Nernst-Planck (KNP) model and the Extracellular-Membrane-Intracellular (EMI) framework, we refer to this model as the KNP-EMI model. Here, we introduce and numerically evaluate a new, finite element-based numerical scheme for the KNP-EMI model, capable of efficiently and flexibly handling geometries of arbitrary dimension and arbitrary polynomial degree. Moreover, we compare the electrical potentials predicted by the KNP-EMI and EMI models. Finally, we study ephaptic coupling induced in an unmyelinated axon bundle and demonstrate how the KNP-EMI framework can give new insights in this setting.

Action potential propagation and synchronisation in myelinated axons

With the advent of advanced MRI techniques it has become possible to study axonal white matter non-invasively and in great detail. Measuring the various parameters of the long-range connections of the brain opens up the possibility to build and refine detailed models of large-scale neuronal activity. One particular challenge is to find a mathematical description of action potential propagation that is sufficiently simple, yet still biologically plausible to model signal transmission across entire axonal fibre bundles. We develop a mathematical framework in which we replace the Hodgkin-Huxley dynamics by a spike-diffuse-spike model with passive sub-threshold dynamics and explicit, threshold-activated ion channel currents. This allows us to study in detail the influence of the various model parameters on the action potential velocity and on the entrainment of action potentials between ephaptically coupled fibres without having to recur to numerical simulations. Specifically, we recover known results regarding the influence of axon diameter, node of Ranvier length and internode length on the velocity of action potentials. Additionally, we find that the velocity depends more strongly on the thickness of the myelin sheath than was suggested by previous theoretical studies. We further explain the slowing down and synchronisation of action potentials in ephaptically coupled fibres by their dynamic interaction. In summary, this study presents a solution to incorporate detailed axonal parameters into a whole-brain modelling framework.

With more and more data becoming available on white-matter tracts, the need arises to develop modelling frameworks that incorporate these data at the whole-brain level. This requires the development of efficient mathematical schemes to study parameter dependencies that can then be matched with data, in particular the speed of action potentials that cause delays between brain regions. Here, we develop a method that describes the formation of action potentials by threshold activated currents, often referred to as spike-diffuse-spike modelling. A particular focus of our study is the dependence of the speed of action potentials on structural parameters. We find that the diameter of axons and the thickness of the myelin sheath have a strong influence on the speed, whereas the length of myelinated segments and node of Ranvier length have a lesser effect. In addition to examining single axons, we demonstrate that action potentials between nearby axons can synchronise and slow down their propagation speed.

ELFENN: A Generalized Platform for Modeling Ephaptic Coupling in Spiking Neuron Models

The transmembrane ionic currents that underlie changes in a cell's membrane potential give rise to electric fields in the extracellular space. In the context of brain activity, these electric fields form the basis for extracellularly recorded signals, such as multiunit activity, local field potentials and electroencephalograms. Understanding the underlying neuronal dynamics and localizing current sources using these signals is often challenging, and therefore effective computational modeling approaches are critical. Typically, the electric fields from neural activity are modeled in a post-hoc form, i.e., a traditional neuronal model is used to first generate the membrane currents, which in turn are then used to calculate the electric fields. When the conductivity of the extracellular space is high, the electric fields are weak, and therefore treating membrane currents and electric fields separately is justified. However, in brain regions of lower conductivity, extracellular fields can feed back and significantly influence the underlying transmembrane currents and dynamics of nearby neurons—this is often referred to as ephaptic coupling. The closed-loop nature of ephaptic coupling cannot be modeled using the post-hoc approaches implemented by existing software tools; instead, electric fields and neuronal dynamics must be solved simultaneously. To this end, we have developed a generalized modeling toolbox for studying ephaptic coupling in compartmental neuron models: E_LFENN (E_Lectric Field Effects in Neural Networks). In open loop conditions, we validate the separate components of E_LFENN for modeling membrane dynamics and associated field potentials against standard approaches (NEURON and LFPy). Unlike standard approaches however, E_LFENN enables the closed-loop condition to be modeled as well, in that the field potentials can feed back and influence membrane dynamics. As an example closed-loop case, we use E_LFENN to study phase-locking of action potentials generated by a population of axons running parallel in a bundle. Being able to efficiently explore ephaptic coupling from a computational perspective using tools, such as E_LFENN will allow us to better understand the physical basis of electric fields in the brain, as well as the conditions in which these fields may influence neuronal dynamics in general.

How does the presence of neural probes affect extracellular potentials?

OBJECTIVE

Mechanistic modeling of neurons is an essential component of computational neuroscience that enables scientists to simulate, explain, and explore neural activity. The conventional approach to simulation of extracellular neural recordings first computes transmembrane currents using the cable equation and then sums their contribution to model the extracellular potential. This two-step approach relies on the assumption that the extracellular space is an infinite and homogeneous conductive medium, while measurements are performed using neural probes. The main purpose of this paper is to assess to what extent the presence of the neural probes of varying shape and size impacts the extracellular field and how to correct for them.

APPROACH

We apply a detailed modeling framework allowing explicit representation of the neuron and the probe to study the effect of the probes and thereby estimate the effect of ignoring it. We use meshes with simplified neurons and different types of probe and compare the extracellular action potentials with and without the probe in the extracellular space. We then compare various solutions to account for the probes' presence and introduce an efficient probe correction method to include the probe effect in modeling of extracellular potentials.

MAIN RESULTS

Our computations show that microwires hardly influence the extracellular electric field and their effect can therefore be ignored. In contrast, multi-electrode arrays (MEAs) significantly affect the extracellular field by magnifying the recorded potential. While MEAs behave similarly to infinite insulated planes, we find that their effect strongly depends on the neuron-probe alignment and probe orientation.

SIGNIFICANCE

Ignoring the probe effect might be deleterious in some applications, such as neural localization and parameterization of neural models from extracellular recordings. Moreover, the presence of the probe can improve the interpretation of extracellular recordings, by providing a more accurate estimation of the extracellular potential generated by neuronal models.

Multimodal Modeling of Neural Network Activity: Computing LFP, ECoG, EEG, and MEG Signals With LFPy 2.0

Recordings of extracellular electrical, and later also magnetic, brain signals have been the dominant technique for measuring brain activity for decades. The interpretation of such signals is however nontrivial, as the measured signals result from both local and distant neuronal activity. In volume-conductor theory the extracellular potentials can be calculated from a distance-weighted sum of contributions from transmembrane currents of neurons. Given the same transmembrane currents, the contributions to the magnetic field recorded both inside and outside the brain can also be computed. This allows for the development of computational tools implementing forward models grounded in the biophysics underlying electrical and magnetic measurement modalities. LFPy (LFPy.readthedocs.io) incorporated a well-established scheme for predicting extracellular potentials of individual neurons with arbitrary levels of biological detail. It relies on NEURON (neuron.yale.edu) to compute transmembrane currents of multicompartment neurons which is then used in combination with an electrostatic forward model. Its functionality is now extended to allow for modeling of networks of multicompartment neurons with concurrent calculations of extracellular potentials and current dipole moments. The current dipole moments are then, in combination with suitable volume-conductor head models, used to compute non-invasive measures of neuronal activity, like scalp potentials (electroencephalographic recordings; EEG) and magnetic fields outside the head (magnetoencephalographic recordings; MEG). One such built-in head model is the four-sphere head model incorporating the different electric conductivities of brain, cerebrospinal fluid, skull and scalp. We demonstrate the new functionality of the software by constructing a network of biophysically detailed multicompartment neuron models from the Neocortical Microcircuit Collaboration (NMC) Portal (bbp.epfl.ch/nmc-portal) with corresponding statistics of connections and synapses, and compute in vivo-like extracellular potentials (local field potentials, LFP; electrocorticographical signals, ECoG) and corresponding current dipole moments. From the current dipole moments we estimate corresponding EEG and MEG signals using the four-sphere head model. We also show strong scaling performance of LFPy with different numbers of message-passing interface (MPI) processes, and for different network sizes with different density of connections. The open-source software LFPy is equally suitable for execution on laptops and in parallel on high-performance computing (HPC) facilities and is publicly available on GitHub.com.

Ephaptic Coupling Promotes Synchronous Firing of Cerebellar Purkinje Cells

Correlated neuronal activity at various timescales plays an important role in information transfer and processing. We find that in awake-behaving mice, an unexpectedly large fraction of neighboring Purkinje cells (PCs) exhibit sub-millisecond synchrony. Correlated firing usually arises from chemical or electrical synapses, but, surprisingly, neither is required to generate PC synchrony. We therefore assessed ephaptic coupling, a mechanism in which neurons communicate via extracellular electrical signals. In the neocortex, ephaptic signals from many neurons summate to entrain spiking on slow timescales, but extracellular signals from individual cells are thought to be too small to synchronize firing. Here we find that a single PC generates sufficiently large extracellular potentials to open sodium channels in nearby PC axons. Rapid synchronization is made possible because ephaptic signals generated by PCs peak during the rising phase of action potentials. These findings show that ephaptic coupling contributes to the prevalent synchronization of nearby PCs.

Senomic view of the cell: Senome versus Genome

In the legacy of Thomas Henry Huxley, and his 'epigenetic' philosophy of biology, cells are proposed to represent a trinity of three memory-storing media: Senome, Epigenome, and Genome that together comprise a cell-wide informational architecture. Our current preferential focus on the Genome needs to be complemented by a similar focus on the Epigenome and a here proposed Senome, representing the sum of all the sensory experiences of the cognitive cell and its sensing apparatus. Only then will biology be in a position to embrace the whole complexity of the eukaryotic cell, understanding its true nature which allows the communicative assembly of cells in the form of sentient multicellular organisms.

Alternating chimeras in networks of ephaptically coupled bursting neurons

The distinctive phenomenon of the chimera state has been explored in neuronal systems under a variety of different network topologies during the last decade. Nevertheless, in all the works, the neurons are presumed to interact with each other directly with the help of synapses only. But, the influence of ephaptic coupling, particularly magnetic flux across the membrane, is mostly unexplored and should essentially be dealt with during the emergence of collective electrical activities and propagation of signals among the neurons in a network. Through this article, we report the development of an emerging dynamical state, namely, the alternating chimera, in a network of identical neuronal systems induced by an external electromagnetic field. Owing to this interaction scenario, the nonlinear neuronal oscillators are coupled indirectly via electromagnetic induction with magnetic flux, through which neurons communicate in spite of the absence of physical connections among them. The evolution of each neuron, here, is described by the three-dimensional Hindmarsh-Rose dynamics. We demonstrate that the presence of such non-locally and globally interacting external environments induces a stationary alternating chimera pattern in the ensemble of neurons, whereas in the local coupling limit, the network exhibits a transient chimera state whenever the local dynamics of the neurons is of the chaotic square-wave bursting type. For periodic square-wave bursting of the neurons, a similar qualitative phenomenon has been witnessed with the exception of the disappearance of cluster states for non-local and global interactions. Besides these observations, we advance our work while providing confirmation of the findings for neuronal ensembles exhibiting plateau bursting dynamics and also put forward the fact that the plateau pattern actually favors the alternating chimera more than others. These results may deliver better interpretations for different aspects of synchronization appearing in a network of neurons through field coupling that also relaxes the prerequisite of synaptic connectivity for realizing the chimera state in neuronal networks.

Propagation of beta/gamma rhythms in the cortico-basal ganglia circuits of the parkinsonian rat

Much of the motor impairment associated with Parkinson’s disease is thought to arise from pathological activity in the networks formed by the basal ganglia (BG) and motor cortex. To evaluate several hypotheses proposed to explain the emergence of pathological oscillations in parkinsonism, we investigated changes to the directed connectivity in BG networks following dopamine depletion. We recorded local field potentials (LFPs) in the cortex and basal ganglia of rats rendered parkinsonian by injection of 6-hydroxydopamine (6-OHDA) and in dopamine-intact controls. We performed systematic analyses of the networks using a novel tool for estimation of directed interactions (nonparametric directionality, NPD). We used a “conditioned” version of the NPD analysis that reveals the dependence of the correlation between two signals on a third reference signal. We find evidence of the dopamine dependency of both low-beta (14–20 Hz) and high-beta/low-gamma (20–40 Hz) directed network interactions. Notably, 6-OHDA lesions were associated with enhancement of the cortical “hyperdirect” connection to the subthalamic nucleus (STN) and its feedback to the cortex and striatum. We find that pathological beta synchronization resulting from 6-OHDA lesioning is widely distributed across the network and cannot be located to any individual structure. Furthermore, we provide evidence that high-beta/gamma oscillations propagate through the striatum in a pathway that is independent of STN. Rhythms at high beta/gamma show susceptibility to conditioning that indicates a hierarchical organization compared with those at low beta. These results further inform our understanding of the substrates for pathological rhythms in salient brain networks in parkinsonism.

NEW & NOTEWORTHY We present a novel analysis of electrophysiological recordings in the cortico-basal ganglia network with the aim of evaluating several hypotheses concerning the origins of abnormal brain rhythms associated with Parkinson’s disease. We present evidence for changes in the directed connections within the network following chronic dopamine depletion in rodents. These findings speak to the plausibility of a “short-circuiting” of the network that gives rise to the conditions from which pathological synchronization may arise.

Signatures of Somatic Inhibition and Dendritic Excitation in Auditory Brainstem Field Potentials

Extracellular voltage recordings (V_e ; field potentials) provide an accessible view of in vivo neural activity, but proper interpretation of field potentials is a long-standing challenge. Computational modeling can aid in identifying neural generators of field potentials. In the auditory brainstem of cats, spatial patterns of sound-evoked V_e can resemble, strikingly, V_e generated by current dipoles. Previously, we developed a biophysically-based model of a binaural brainstem nucleus, the medial superior olive (MSO), that accounts qualitatively for observed dipole-like V_e patterns in sustained responses to monaural tones with frequencies >∼1000 Hz (Goldwyn et al., 2014). We have observed, however, that V_e patterns in cats of both sexes appear more monopole-like for lower-frequency tones. Here, we enhance our theory to accurately reproduce dipole and non-dipole features of V_e responses to monaural tones with frequencies ranging from 600 to 1800 Hz. By applying our model to data, we estimate time courses of paired input currents to MSO neurons. We interpret these inputs as dendrite-targeting excitation and soma-targeting inhibition (the latter contributes non-dipole-like features to V_e responses). Aspects of inferred inputs are consistent with synaptic inputs to MSO neurons including the tendencies of inhibitory inputs to attenuate in response to high-frequency tones and to precede excitatory inputs. Importantly, our updated theory can be tested experimentally by blocking synaptic inputs. MSO neurons perform a critical role in sound localization and binaural hearing. By solving an inverse problem to uncover synaptic inputs from V_e patterns we provide a new perspective on MSO physiology.SIGNIFICANCE STATEMENT Extracellular voltages (field potentials) are a common measure of brain activity. Ideally, one could infer from these data the activity of neurons and synapses that generate field potentials, but this "inverse problem" is not easily solved. We study brainstem field potentials in the region of the medial superior olive (MSO); a critical center in the auditory pathway. These field potentials exhibit distinctive spatial and temporal patterns in response to pure tone sounds. We use mathematical modeling in combination with physiological and anatomical knowledge of MSO neurons to plausibly explain how dendrite-targeting excitation and soma-targeting inhibition generate these field potentials. Inferring putative synaptic currents from field potentials advances our ability to study neural processing of sound in the MSO.

An Evaluation of the Accuracy of Classical Models for Computing the Membrane Potential and Extracellular Potential for Neurons

Two mathematical models are part of the foundation of Computational neurophysiology; (a) the Cable equation is used to compute the membrane potential of neurons, and, (b) volume-conductor theory describes the extracellular potential around neurons. In the standard procedure for computing extracellular potentials, the transmembrane currents are computed by means of (a) and the extracellular potentials are computed using an explicit sum over analytical point-current source solutions as prescribed by volume conductor theory. Both models are extremely useful as they allow huge simplifications of the computational efforts involved in computing extracellular potentials. However, there are more accurate, though computationally very expensive, models available where the potentials inside and outside the neurons are computed simultaneously in a self-consistent scheme. In the present work we explore the accuracy of the classical models (a) and (b) by comparing them to these more accurate schemes. The main assumption of (a) is that the ephaptic current can be ignored in the derivation of the Cable equation. We find, however, for our examples with stylized neurons, that the ephaptic current is comparable in magnitude to other currents involved in the computations, suggesting that it may be significant-at least in parts of the simulation. The magnitude of the error introduced in the membrane potential is several millivolts, and this error also translates into errors in the predicted extracellular potentials. While the error becomes negligible if we assume the extracellular conductivity to be very large, this assumption is, unfortunately, not easy to justify a priori for all situations of interest.

Dipolar extracellular potentials generated by axonal projections

Extracellular field potentials (EFPs) are an important source of information in neuroscience, but their physiological basis is in many cases still a matter of debate. Axonal sources are typically discounted in modeling and data analysis because their contributions are assumed to be negligible. Here, we established experimentally and theoretically that contributions of axons to EFPs can be significant. Modeling action potentials propagating along axons, we showed that EFPs were prominent in the presence of terminal zones where axons branch and terminate in close succession, as found in many brain regions. Our models predicted a dipolar far field and a polarity reversal at the center of the terminal zone. We confirmed these predictions using EFPs from the barn owl auditory brainstem where we recorded in nucleus laminaris using a multielectrode array. These results demonstrate that axonal terminal zones can produce EFPs with considerable amplitude and spatial reach.