# Influence of extracellular oscillations on neural communication: a computational perspective

Zoran Tiganj, Sylvain Chevallier, Eric Monacelli

► To cite this version:

Zoran Tiganj, Sylvain Chevallier, Eric Monacelli. Influence of extracellular oscillations on neural communication: a computational perspective. Frontiers in Computational Neuroscience, 2014, 8 (9), 10.3389/fncom.2014.00009. hal-01214044

HAL Id: hal-01214044

https://hal.uvsq.fr/hal-01214044

Submitted on 9 Oct 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Influence of extracellular oscillations on neural communication: a computational perspective

Zoran Tiganj*, Sylvain Chevallier and Eric Monacelli

LISV, Université de Versailles Saint-Quentin-en-Yvelines, Vélizy-Villacoublay, France

**Edited by:**
Terrence J. Sejnowski, The Salk Institute for Biological Studies, USA

**Reviewed by:**
G. B. Ermentrout, University of Pittsburgh, USA
Benjamin Torben-Nielsen, Hebrew University of Jerusalem, Israel

**\*Correspondence:**
Zoran Tiganj, Theoretical Cognitive Neuroscience Lab, Department of Psychological and Brain Sciences, Center for Memory and Brain, Boston University, 2 Cummington Street, Boston, MA 02215, USA
e-mail: zorant@bu.edu

Neural communication generates oscillations of electric potential in the extracellular medium. In feedback, these oscillations affect the electrochemical processes within the neurons, influencing the timing and the number of action potentials. It is unclear whether this influence should be considered only as noise or it has some functional role in neural communication. Through computer simulations we investigated the effect of various sinusoidal extracellular oscillations on the timing and number of action potentials. Each simulation is based on a multicompartment model of a single neuron, which is stimulated through spatially distributed synaptic activations. A thorough analysis is conducted on a large number of simulations with different models of CA3 and CA1 pyramidal neurons which are modeled using realistic morphologies and active ion conductances. We demonstrated that the influence of the weak extracellular oscillations, which are commonly present in the brain, is rather stochastic and modest. We found that the stronger fields, which are spontaneously present in the brain only in some particular cases (e.g., during seizures) or that can be induced externally, could significantly modulate spike timings.

**Keywords:** extracellular oscillations, local field potentials, ephaptic coupling, nonsynaptic communication, multi-compartment model, NEURON simulation environment, pyramidal neurons

## 1. INTRODUCTION

Ion exchange between neurons and the surrounding extracellular medium is essential for neural communication. In the extracellular medium, this exchange causes oscillations of electric potential, commonly referred to as Local Field Potentials (LFPs) (Ebersole and Milton, 2003; Buzsáki et al., 2012). The LFPs often exhibit significant levels of temporal and spatial coherence and are linked with a number of cognitive phenomena (see Ward, 2003; Buzsáki and Draguhn, 2004; Uhlhaas et al., 2010; Wang, 2010 for reviews). For instance, theta and gamma oscillations, particularly in the hippocampus, are often linked with memory formation and retrieval (Schack et al., 2002; Sederberg et al., 2003; Osipova et al., 2006; Nyhus and Curran, 2010), while attention is associated with a reduced alpha rhythm and increased gamma rhythm (Thut et al., 2006; Lisman and Buzsáki, 2008; Capotosto et al., 2009; Deco and Thiele, 2009; Schroeder and Lakatos, 2009). Different rhythms are related with processing in early sensory systems (Koepsell et al., 2010) and in speech generation and processing (Peelle and Davis, 2012). Furthermore, it has recently been proposed that some neurons code information by firing with respect to the phase of the ongoing oscillations (Fries et al., 2007; Montemurro et al., 2008; Kayser et al., 2009).

The extracellular oscillations are commonly considered to be an epiphenomena, a sporadic consequence of neural dynamics. Contrary to this hypothesis, it has been speculated that the periodic oscillations in the extracellular medium in some cases have a functional role, that they support the neural communication. This alternative hypothesis, based on *ephaptic transmission*, assumes that the neurons might not communicate only

by exchanging information through their synaptic connections, but additionally through the extracellular medium (Fröhlich and McCormick, 2010; Weiss and Faber, 2010; Anastassiou et al., 2011). Even though the LFPs are generally too weak to trigger an action potential, they can still affect the spike number and timing (Parra and Bikson, 2004; Radman et al., 2007; Anastassiou et al., 2010).

If the hypotheses about the functional role of the extracellular oscillations are proven right, it will provide new insights for understanding the neural code, indicating mechanisms that can contribute to the widely observed synchrony of the neural activity. This finding could have impact on theoretical and clinical research, including neuropharmacology and neuroprosthetics. Additionally, it would be beneficial for developing new approaches in computational models of the brain, and more generally, models of artificial neural networks, since the existing models (e.g., Izhikevich, 2006; Markram, 2006) do not address the ephaptic transmission.

The main motivations for the proposed study are to deepen the understanding of this prospective mode of neural communication and to investigate whether the externally induced electrical fields (e.g., through Transcranial Magnetic Stimulation (TMS) (Barker et al., 1985; Hallett, 2000; Walsh and Cowey, 2000; Allen et al., 2007), Transcranial Electrical Stimulation (TEC) (Deans et al., 2007; Kirov et al., 2009; Ozen et al., 2010; Reato et al., 2010) and Deep Brain Stimulation (DBS) (Benabid et al., 2009; Joucla and Yvert, 2012) might promote the ephaptic coupling. Modifying the ephaptic coupling could alter synchrony of the neural assemblies and potentially help to correct pathological network dynamics.
Through a computational study, we investigated what changes in the timing and number of action potentials can be induced by the temporal and spatial oscillations of the extracellular potential. Our study is conducted at a single neuron level and is done through several simulation setups, hard to implement experimentally. Contrary to the common experimental protocols, where the stimulation is done by current injection into the soma (see for example Anastassiou et al., 2011), we used PostSynaptic Potentials (PSPs) which are applied at different spatial locations of the dendrites. We generated a sufficient number of PSPs for the neuron to exhibit different firing patterns. Moreover, we analyzed the influence of extracellular oscillations with respect to the spatial part of the dendritic tree receiving most of the inputs. We modeled the dendrites with active ion conductances and considered different dendritic morphologies. Since the considered neuronal model is multi-compartmental, we computed the strength of the extracellular oscillations for each compartment separately, so that it depends on the spatial location of the compartment.

Focusing on models of the CA1 and the CA3 neurons, we investigated the influence of the extracellular oscillations at a single neuron level depending on:

1. Properties of the active ion conductances, i.e., we used different sets of active ion channels.

2. Spatial distribution of the dendritic inputs, i.e., we altered the locations of the synaptic inputs between the basal and apical branches.

3. Morphological properties, i.e., we used different shapes of the dendritic tree.

For each of these three cases, we examined the following properties of the extracellular oscillations (as shown in **Figure 1**): the temporal and spatial frequencies, the amplitude, and the phase. In this study, we always addressed the influence of the oscillations on an individual neuron, not on neural assemblies. Through randomly delivered PSPs the membrane potential was set to fluctuate around the threshold level, with a certain probability of crossing it. We studied how the above properties of the extracellular oscillations influence the number and the timing of generated spikes.

We did not investigate directly the influence of the spontaneous LFPs or the artificially induced oscillations, but the influence of the extracellular oscillations on a more general level, regardless of their origin. The link between this general level and the spontaneous LFPs of a particular temporal and spatial properties can still be made and will be discussed here.

It is important to mention that the proposed study is difficult to perform experimentally, either *in vivo* or *in vitro*. It is particularly difficult to produce repetitive stimulation with the PSPs of controlled amplitude, timing and spatial distribution. Since such repetitive stimulation is essential for our investigation, we find it suitable to approach the problem through a computational study. Complex interactions of nonlinear ionic currents within branched dendritic trees can be described through a set of equations, based on the Hodgkin–Huxley model, which cannot be solved analytically. We used the simulator NEURON (Hines and

<figure><img src="image_3.png" /><figcaption>FIGURE 1 | Flowchart of the computational model.</figcaption></figure>

Carnevale, 1997) to numerically solve the equations (Mainen and Sejnowski, 1998). This study is not a replacement for *in vivo* or *in vitro* experiments, but rather a complementary approach that can help in focusing future research.

Common use of the Hodgkin–Huxley model assumes that the extracellular medium is equipotential and it is typically considered as an electrical ground. Here we assume that the electric potential in the extracellular space varies through time and space. If we consider the neural membrane to consist of the ohmic conductance (active and passive) and the capacitance, the above assumption implies that the voltage on the conductance and capacitance is not simply equal to the intracellular membrane potential, which is the case in the Hodgkin–Huxley model.
Instead, it is set to be equal to the difference between the intra-
cellular membrane potential and the instantaneous amplitude of
the extracellular oscillations.

The paper is organized as follows. Parameters of the model
are closely analyzed in the following section. In particular, we
describe the neural morphologies and give details on the proper-
ties of the ion channels in the model. We also describe the model
of PSPs and the extracellular oscillations. In section 3 we present
the results of the simulations: first, when the stimulation is done
with the current injection, and then when it is done with the PSPs.
The results illustrate the importance of the amplitude, the tempo-
ral and spatial frequency and the spatial phase of the extracellular
oscillations. Finally, the results are discussed in section 4.

## 2. MATERIALS AND METHODS

This section provides details on the computational models that
are used to obtain the results presented in the paper. The mod-
els are implemented in neural simulator NEURON (Hines and
Carnevale, 1997) that we run on a Linux-based cluster, consisting
of 75 dual-core computers. We used multi-compartment models
for two different neuron types: the CA1 and the CA3 rat hip-
 pocampal pyramidal neurons. The time step of the simulation
was 25µs. The stability of the models with respect to the time step
was tested by verifying that for smaller times steps the influence
of the extracellular oscillations does not change significantly. To
account for the artifacts caused by the initial conditions we added
an extra period of 100 ms at the beginning of each simulation (not
displayed in the results) that allowed the membrane potential to
grow from the initial –65 mV to the value controlled only by the
PSPs.

2.1. MORPHOLOGY

Each neuron in our model consists of the dendritic tree, the soma
and the short axon. The length of the axon was determined exper-
imentally (100µm), in a such way that simulating longer axon did
not quantitatively change the results. For each of the two neuron
types, we investigated five different morphologies of the dendritic
tree. The morphologies of the CA3¹ are taken from Samsonovich
and Ascoli (2003) and those of the CA1 cells<fnref n="2" /> are taken from
Pyapali et al. (1998).

All the morphology files are given in SWC format. In this for-
mat, each segment is described as a cylinder and is associated with
a type (soma, axon, basal, or apical dendrite), a spatial position
of the end point, a diameter, and the information about which
segment is connected to the given one on the path to the soma.
Thus, from such files we were able to precisely reconstruct 3D
shapes of the given morphologies. Each morphology file describes
a branched dendritic tree, with the number of compartments
spanning from 151 to 503. The size of the compartments was
computed using the dossy rule proposed by Carnevale and
Hines (2006).

## 2.2. PROPERTIES OF THE ACTIVE ION CONDUCTANCES

To study how the influence of the extracellular oscillations on the
spike number and timing depends on the active ion conductances
in the soma and the dendrites, two different sets of active ion con-
ductances were used, one corresponding to the CA1 and the other
to the CA3 pyramidal neurons. Notice that adding active conduc-
tances in the dendrites was necessary for the PSPs to propagate
in a significant amount from the apical dendrites to the soma. In
such way the PSPs could contribute to the generation of the den-
dritic as well as the somatic spikes (Kim et al., 2012). Properties of
the ion channels used in the simulations have been described in
Migliore et al. (2005) for the CA1 neurons and in Hemond et al.
(2008) for the CA3 neurons<fnref n="3" />:

* CA1: Uniform passive properties were $\tau_m$ = 28 ms (membrane time constant), $R_m$=28 kΩcm<fnref n="2" /> (radial resistance of the membrane), $R_a$ = 150 Ωcm (axial resistance of the membrane), while the active properties included sodium, DR-, and A-type potassium conductances and a hyperpolarization-activated $I_h$ current. The kinetics for all the channels and the spatial distribution of the channels were identical to those in Migliore et al. (2005).

* CA3: Passive: $\tau_m$ = 35 ms, $R_m$ = 25 Ωcm<fnref n="2" />, $R_a$ = 150 Ωcm. The set of active channel properties included sodium conductances, a repertoire of potassium conductances: delayed rectifier, M-current, fast-inactivating A-type, and slowly-inactivating D-type current, three voltage-gated $Ca^{2+}$ conductances ($CaV$ N-, L-, and T-type), two calcium-dependent potassium conductances (KC and KAHP), and $I_h$ current.

These two sets give rather different profiles of ion conductances.
While for the CA1 neurons only the main conductances involved
in the action potential generation and regulation were included,
for the CA3 neurons we included a much broader set of ion
conductances.

2.3. MODEL OF POSTSYNAPTIC POTENTIALS

The spatial location of each synapse and their corresponding
PSPs times were determined randomly for each simulation. We
repeated each simulation for several different properties of the
extracellular field, e.g., we varied the amplitude of the extracel-
lular field. To study the influence of the extracellular oscillations,
we kept the spatial locations and the PSPs times fixed for each
different property of the extracellular field.

To emphasize spatiotemporal summation of the PSPs, we bal-
anced the contributions from different synapses by pacing them
only onto dendritic compartments thinner than 0.35µm. The
probability of a synapse being placed onto any compartment thin-
ner than 0.35µm was set to be equal, unless otherwise specified.
The exact spatial location within the compartment was also deter-
mined randomly, with equal probability that the synapse will be
placed onto any of the sections within the compartment. One

<footnote n="1">¹We used morphologies named: cell_a, cell_c, cell_m,
cell_t, and cell_v; they are available for download on
http://krasnowl.gmu.edu/cn3/L-Neuron/database/.</footnote>

<footnote n="2">²We used morphologies named: n121CNG, n129CNG, n141CNG,
n142CNG, and n145CNG; they are available for download on
http://www.neuromorpho.org.</footnote>

<footnote n="3">³Both models are available in ModelDB on the following addresses:
http://senselab.med.yale.edu/modeldb/showmodel.asp?model=101629 and
http://senselab.med.yale.edu/modeldb/showmodel.asp?model=55035.</footnote>
example of synaptic locations generated in such way is given in **Figure 2A**.

Temporal profile of the PSPs was modeled as a Poisson process. The number of the PSPs per second was adjusted so that the neuron exhibits a firing rate *F* within a desired range. We define three ranges of the firing rate: low for *F* ≤ 5 Hz, medium for 5 < *F* < 15 Hz and high for *F* ≥ 15 Hz. To achieve these firing rates we altered the of the PSPs to be between 200 and 7000 Hz.

The synaptic conductance is modeled for $t \ge t_{onset}$ as an alpha function:

$$g = g_{\max} \frac{t - t_{\text{onset}}}{\tau} e^{-\frac{1}{\tau}(t - t_{\text{onset}} - \tau)}, \quad (1)$$

while for $t < t_{onset}$, $g = 0$. The parameters within the alpha function were the same for all the synapses: $\tau = 0.1$ ms, $g_{max} = 1\mu S$.

## 2.4. MODEL OF THE EXTRACELLULAR OSCILLATIONS

The voltage oscillations were simulated in the extracellular medium by adding a voltage source to each compartment<fnref n="4" />, as illustrated on **Figure 2B**. We modeled the oscillations as a temporal and spatial fluctuation. The change across space was only occurring along the somatodendritic axis, which we denote as the *x* axis. The other two axes of the Cartesian coordinate system are denoted as *y* and *z*. Consequently, at any time, the amplitude of the extracellular oscillations was different for different spatial

points along the somatodendritic axis, but the same for any spatial point along the other two axis (*y* and *z*). This is depicted on a 2D plane in **Figure 3**. Shape of the oscillations is described by a simple sinusoidal function $V_e(x, y, z, t)$. Then, the extracellular potential at the spatial point (*x*, *y*, *z*) and at the time *t* writes as follows:

$$V_e(x, y, z, t) = V_o \sin(2\pi f_t t) \sin(2\pi f_s x + \phi_s) \quad (2)$$

where $V_o$ is the amplitude of the oscillations, $f_t$ is the temporal frequency, $f_s$ is the spatial frequency and $\phi_s$ is the spatial phase.

In the context of the Hodgkin–Huxley model, the equation for the current $I_k$ going through the *k*th ion channel is:

$$I_k = g_k(V_i - V_e - V_k) \quad (3)$$

where $V_i$ is the intracellular potential, $V_k$ is the reversal potential of the *k*th ion channel and $g_k$ is the channel conductance.

To have better understanding of the influence of the extracellular oscillations on neural activity, the simulations were repeated for different values of $V_o$, $f_t$, $f_s$ and $\phi_s$.

The amplitude of the extracellular oscillations depends on many factors, from the animal species, to the regions of the brain and the cognitive states (Buzsáki et al., 2012). In the band below 100 Hz, the amplitude rarely exceeds 1 mV (Anastassiou et al., 2010; Brette and Destexhe, 2011; Buzsáki et al., 2012). We performed the simulations for the following set of amplitudes $V_o$: 0.25, 0.5, 1, 2 and 4 mV. This wide range of values allowed us to investigate the significance of the amplitude in a rather general

<figure><img src="image_4.png" /></figure>

<footnote n="4">⁴Within the simulator NEURON, we used mechanism extracellular to drive the voltage oscillations into each compartment.</footnote>
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | The amplitude of the extracellular oscillations across one full temporal cycle of 100 ms, i.e., $f_t = 10$ Hz. In this example $V_o = 4$ mV, $f_s = 0.8$ mm$^{-1}$ and $\phi_s = \frac{3\pi}{10}$ rad. The parameter $\phi_s$ is expressed at the bottom of the plots as the reference.</figcaption></figure>

context, including the spontaneous and the artificially induced oscillations.

The extracellular oscillations that exhibit a significant spatial extent commonly have temporal frequencies below 100 Hz, although there are notable exceptions, such as the hippocampal sharp waves (Buzsáki, 1986). Even though the high amplitude oscillations are known to occur in several particular bands, e.g., theta waves (Kahana et al., 2001), we performed the simulations for a relatively broad set of frequency values: 1, 5, 10, 20, 50, and 100 Hz. These values have been chosen to investigate whether the values observed in reality are in some sense optimal for high or low influence on neural activity.

Regarding the spatial frequency, sources and sinks were located at the opposite sides of the dendritic tree, e.g., if the sources were located around the apical dendrites, then the sinks were located around the basal dendrites and vice versa. Therefore, we set $f_s = 0.8$ mm$^{-1}$ as illustrated in **Figure 3**. The influence of the extracellular oscillations when sinks and sources were more distant from each other was also investigated by setting the spatial frequency $f_s = 0.4$ mm$^{-1}$. In this case, to ensure the consistency of the results, the strength of the electrical field was kept constant: when the distance between the sink and the source was doubled, the amplitude of the extracellular oscillations was also doubled. The spatial distributions of the oscillations for both cases are shown in **Figure 4A**.

The spatial phase of the oscillations was set so that the peak amplitude was reached around the middle of the basal dendrites and around the upper end of the apical dendrites, as these areas often contain the largest number of branches. The exact value of the phase was $\phi_s = \frac{3\pi}{10}$ rad, as shown on the left hand side plot in **Figure 4B**. Additionally, we investigated a case when the peak amplitude was reached around the soma. In that case $\phi_s = \frac{\pi}{2}$ rad was used (see the right hand side plot in **Figure 4B**).

# 3. RESULTS

In this section we first present the results of a simulation conducted with current injection as the stimulus input. We then present the results from a set of simulations in which the PSPs, instead of current injection, were used to stimulate the neuron.

## 3.1. STIMULATION WITH CURRENT INJECTION

As an initial validation of the models used in this study we performed a simulation that is similar to the experimental protocol described in Anastassiou et al. (2011). The goal of the simulation was to analyze the influence of extracellular oscillations on spike timing during current injection into soma. The results are displayed in **Figure 5**.

For the CA3 neuron model, the injection of a 0.35 nA constant current into the soma led to a tonic spiking, such that each action potential was fired periodically, approximately every 75 ms. Adding the extracellular oscillations affected this periodicity, degrading the period stability. The time shift of the spikes increases with the amplitude of the oscillations. However, the increase is not linear, due to the complex nonlinear dynamics of the voltage-gated ion channels in the soma and in the dendritic tree.

The results displayed in **Figure 5** are qualitatively similar to those reported in the experiments presented in Anastassiou et al. (2011). This type of simulation demonstrates that the extracellular oscillations can affect the timing of action potentials. Nevertheless, the influence of the low amplitude oscillations, which are observable in different cortical structures, seems to be rather modest. The temporal shift of the firing instants is very small, e.g., the influence of $V_o = 0.25$ mV oscillations results in 0.5 ms change in the timing of the action potential shown in the zoomed part on the top of **Figure 5**.

## 3.2. STIMULATION WITH POSTSYNAPTIC POTENTIALS

Current injection into the soma gives an illustration of the influence of extracellular oscillations on the timing of action potentials. However, this is not a very realistic setup and. Stimulation mediated by the PSPs is more biologically realistic. The PSPs properties were set as described in section 2.3 and the results of the simulation are presented for different properties of the extracellular oscillations according to the parameters presented in section 2.4.

### 3.2.1. Influence of the amplitude of the extracellular oscillations

The simulations were first performed to explore the influence of the amplitude of the extracellular oscillations on the spike
<figure><img src="image_6.png" /></figure>

**FIGURE 4** | (A) The amplitude of the extracellular oscillations for two spatial frequencies, $f_s = 0.8$ mm (the left side plot) and $0.4 mm^{-1}$ (the right side plot).
(B) The amplitude of the extracellular oscillations for two spatial phases, $\phi_s = \frac{3\pi}{10}$ rad (left side plot) and $\frac{\pi}{2}$ rad (right side plot).

<figure><img src="image_7.png" /></figure>

**FIGURE 5** | Influence of the extracellular oscillations on the timing of the action potentials during current injection into the soma.

The top boxes highlight the influence of the extracellular oscillations on the spike timing. The baseline (no extracellular oscillations) is

displayed in blue. The dashed line represents the extracellular oscillations. We display the maximum amplitude used in the simulations, $V_o = 4$ mV, scaled up by the factor of two for better visibility.

number and timing. The temporal frequency was set to 10 Hz, the spatial frequency to $0.8 \text{ mm}^{-1}$ and the spatial phase to $-\frac{3\pi}{10}$ rad. The plots on the top of **Figure 6** show the membrane potential at the soma of a CA3 neuron for different values of the amplitude of the oscillations. Hereafter we will refer to *baseline* as the evolution of the soma potential evoked by the PSPs in absence of the extracellular oscillations, i.e., it corresponds to the simulations with $V_o = 0$ mV. The baseline is shown in blue in **Figure 6**. In this example, the PSPs frequency was tuned so that the firing rate *F* stays below 5 Hz. The zooms of the four different parts of the plot illustrate the small changes of the membrane potential, as well as the small time shifts in the spike timing, both due to the influence of the extracellular oscillations.

The bottom plot in **Figure 6** further illustrates how the membrane potential depends on the amplitude of the extracellular oscillations. The plot shows the membrane potential at the

soma in the presence of 4 mV extracellular oscillations subtracted from the membrane potential at the soma with no extracellular oscillations (black line). To illustrate the phase relationship, on the same plot we also displayed the extracellular oscillations (blue line). The phase and amplitude of the difference of the two membrane potentials follow those of the extracellular oscillations, but due to a large nonlinearity of the voltage-gated channels, large fluctuations of the amplitude occur frequently.

To have a better understanding of the correlation between the extracellular oscillations and the resulting changes in the timing and the number of action potentials, we performed an ensemble of simulations: simulations on 5 CA1 and 5 CA3 models were repeated 48 times (called *iterations*), resulting in a total number of 480 simulations. As described in section 2.3, in each of the 480 simulations the time onsets of the PSPs and the spatial locations of
<figure><img src="image_8.png" /><figcaption>**FIGURE 6** | **Top:** Influence of the extracellular oscillations on the timing of the action potentials during the stimulation with PSPs. **Bottom:** The extracellular oscillations at the soma ($V_o$ = 4 mV) and the difference between the membrane potential</figcaption></figure>

of the soma without and with the extracellular oscillations of 4 mV amplitude. The vertical lines on the **bottom** plot that exceed the plot boundaries results from action potentials visible on the **top** plot.

the synapses were randomly chosen. **Figure 7A** displays the raster plots obtained for the medium firing rate, $5 < F < 15$ Hz, of the CA3 neuron models (5 cell morphologies, each used for 48 iterations, resulting in 240 simulations). Differences in the timing and the numbers of action potentials are evident, but they are difficult to quantify directly from the raster plots.

To provide a meaningful visualization of the changes caused by the extracellular oscillations, several processing steps were made. First, we created a peristimulus time histogram: the spikes were binned, that is, the time axis was divided into 10 ms bins and the spikes were gathered to obtain a spike count per bin. Then, the spike count of each bin was subtracted from the corresponding baseline spike count. The next step was to generate cumulative results, where all the bins corresponding to the same phase of the oscillations were summed. For example, when 10 Hz sinusoidal oscillations were applied during a 1 s long simulation, the spike counts at different phases of the oscillations were summed. The resulting graphs for basal and apical stimulations of the CA3 and the CA1 neural models, together with the plot of the phase of the oscillations at the soma, are shown in **Figure 7B**. Comprehensive figures of all the intermediate steps are provided in supplementary material.

The stimulations of the basal dendrites, on either the CA1 or the CA3 neurons, yielded a significant negative correlation for the high amplitude oscillations, $V_o \ge 1$ mV. On the contrary, the correlation was positive for the apical dendrites, with a smaller absolute value than the correlation observed for the basal dendrites. These trends, anti-correlation for the basal dendrites and correlation for the apical, are evident in **Figure 7C** where

the cross-correlations between the bins and the phase of the oscillations for three different firing rates *F* are plotted.

The negative correlation during the stimulation of the basal dendrites is caused by the fact that an increase of the extracellular oscillations induces a decrease in the difference between the intracellular and the extracellular potential as written in Equation (3). Thus, the positive half-wave (with respect to the phase of the oscillations at the soma) of the extracellular oscillations drives the membrane potential away from the firing threshold, resulting in the negative correlation and vice versa.

The opposite signs of the correlation coefficients for the stimulation of the basal and the apical dendrites is a direct consequence of the asymmetrical morphology of the pyramidal neurons. As shown in **Figure 3**, for the spatial frequency of $f_s = 0.8$ mm$^{-1}$, when the positive half-wave of the extracellular oscillations is around the soma, the negative half-wave will be around the apical dendrites. The negative half-wave of the oscillations facilitates the spike emission, raising the membrane potential toward the firing threshold. When the overall positive correlation is observed, it indicates that the distal oscillations have a higher influence than those around the soma. This is observed when the neuron is stimulated in the apical part of the dendritic tree, as shown on right hand side plot in **Figure 7C**.

To explain the difference between the absolute values of the correlation coefficients, which are higher for the stimulations of the basal dendrites, we speculate that the difference is influenced by the morphology of the pyramidal neurons: the apical dendrites are widely spread and the basal dendrites are gathered around the soma. Only a small part of the apical tree was affected by the peak
<figure><img src="image_9.png" /><figcaption>**FIGURE 7** | (A) The raster plots showing firing times of CA3 neurons for different amplitudes of the extracellular oscillations, $V_o = 0, 0.25, 0.5, 1, 2$ and $4$ mV, across 48 iterations for 5 morphologies, resulting in 240 spike trains. The phase of the extracellular oscillations is given with red dashed lines. (B) The average (across all 240 iterations) difference between the firing frequencies with and without the extracellular oscillations. Additionally, an observation of 1 s for each</figcaption></figure>

iteration is cumulated and averaged in a single period of the extracellular oscillations (100 ms). All the plots are shown for the medium firing rate, i.e., $5 < F < 15$ Hz. The sin wave, indicating the phase of the extracellular oscillations, is plotted with a dotted line and with a dimensionless amplitude. (C) Correlation coefficients between bins of 1 ms and the phase of the extracellular oscillations for the CA3 and the CA1 neurons for different firing rates F.
of the oscillations, while most of the basal dendrites are affected by the peak. Furthermore, the soma and a nonnegligible part of the apical tree proximal to the soma were stimulated by oscillations with the phase opposite to the oscillations taking place in the distal apical dendrites.

The firing rate plays also a role in the sensitivity of the mem-
brane potential to the influence of the extracellular oscillations.
For CA3 models, it appears that when the basal dendrites were
stimulated the largest change of the correlation coefficients is
observed for low firing rate, then for medium and eventually for
high firing rate. This result is inverted for the CA1 neurons, as
illustrated in **Figure 7B**, suggesting that the link between the fir-
ing rate and the sensitivity to the extracellular oscillations highly
depends on the properties of the ion channels. In our simula-
tions the difference might come from the calcium and potassium
conductances, since they differ in the two models.

After observing how timing of the action potentials is influ-
enced by the extracellular oscillations, we investigated the overall
change in the number of action potentials as a function of the
extracellular oscillations. This indicator is given in **Figure 8A**. The
overall number of action potentials is computed separately for
the positive and the negative period of the extracellular oscilla-
tions. The number of action potentials for the positive period is
shown with red dashed lines and for the negative period with blue
dot-dashed lines.

For the stimulation of the basal dendrites, both the CA3 and
the CA1 neurons exhibit a decrease in the number of action
potentials during the positive period and an increase during the
negative period. The change is larger for the CA3 than for the
CA1 neurons, suggesting that the dynamical properties of the ion
channels are important for this kind of analysis. The effects are
inverted for the apical dendrites and the variations are smaller
than for the basal dendrites: during the positive period there
is an increase in the number of action potentials, while during
the negative period it decreases or does not show a significant
change. This is in accordance with the previous observations

where the proximal and distal contributions from the apical den-
drites are balancing each other out. The average number of action
potentials during both the positive and the negative period is
shown with the solid black line. The average number slightly
increases with the increase of the amplitude of the extracellular
oscillations.

The last indicator we explored assesses the influence of the
extracellular oscillations on the precise timing of action poten-
tials. **Figure 8B** displays the mean timing and its standard devi-
ation. The mean timing tends to stay near the center of the half
period, which is 25 ms for the positive half-wave and 75 ms for
the negative. In cases where the mean value was initially biased,
i.e., significantly different from 25 or 75 ms without the extra-
cellular oscillations, the extracellular oscillations bring the mean
value closer to the middle of the half period. This is evident in
the results of the stimulation of the apical dendrites in **Figure 8B**.
This result is expected since 25 and 75 ms correspond to the peak
of the amplitude of the extracellular oscillations.

**3.2.2. Influence of the temporal frequency of the extracellular oscillations**

To investigate the influence of the temporal frequency of
the extracellular oscillations, we conducted simulations with
V₀ = 4 mV and the following set of the temporal frequencies
f_t = 1, 5, 10, 20, 50, and 100 Hz. The results are shown in
**Figure 9A** and displayed as histograms representing the differ-
ence in the number of action potentials caused by the oscillations.
**Figure 9A** is given for the stimulation of the CA1 basal dendrites,
and the figures that display the results of the stimulation of the
CA1 apical dendrites and the CA3 basal and apical dendrites are
given in the supplemental. The results are cumulated over a single
period of oscillations. The phase of the oscillations at the soma
is indicated with dashed red lines. The stimulation with the PSPs
was set so that the firing frequency was in the medium firing range
interval, i.e., 5 < *F* <15 Hz. Spatial phase of the oscillations was
3π/10 rad.

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

**FIGURE 8** | (A) The number of action potentials as a function of the amplitude of the extracellular oscillations, with the frequency of 10 Hz and for the medium firing rate (5 < *F* <15 Hz). The results for the CA3 neuron models are shown in the first row and for the CA1 in the second row. The curves are obtained from results of 240 simulations (48 iterations for 5 morphologies) with 1 s duration. The results during the positive period w.r.t. the potential at the soma, are

shown with red dashed lines, and during the negative period with blue
dot-dashed lines and the average number of the action potentials is
shown with solid black lines. (B) The mean timing of the action
potentials with the standard deviations. The x-axis denotes the
amplitude of the 10Hz extracellular oscillations. To facilitate the visual
comparison, we subtracted 50 ms from the negative half-wave mean
timing (blue dash-dotted lines).
<figure><img src="image_13.png" /></figure>

**FIGURE 9** | (A) Influence of the frequency of the extracellular oscillations on the number of action potentials for the stimulation of the basal dendrites of the CA3 neuron model. From top to bottom the plots correspond to the following frequencies: $f_t = 1, 5, 10, 20, 50$ and 100 Hz. The bins are created by

<figure><img src="image_14.png" /></figure>

cumulating the results over 1 s of simulation. (B) Cross-correlation between the difference of the bins (cumulated over a single period of the extracellular oscillations) and the extracellular oscillations. The cross-correlation is computed for six different frequencies of the extracellular oscillations.

To evaluate the correlation between the change in the number of action potentials and the frequency of the extracellular oscillations, we computed the cross-correlation coefficient for each frequency of the oscillations. The results for the basal and the apical stimulation of the CA1 and the CA3 neurons are given in **Figure 9B**. The stimulation of the basal dendrites led to anti-correlation with the phase of the oscillations at the soma, which tends to be larger for higher frequencies, except for 100 Hz where it gets out of the phase. When the apical dendrites were stimulated, the relative change in the number of action potentials is still of the same order as when the basal dendrites were stimulated, but the level of correlation with the oscillations at the soma significantly varies for different frequencies. This suggests that the oscillations have important influence, not only at the soma, but at the dendrites as well.

### 3.2.3. Influence of the spatial frequency and the spatial phase of the extracellular oscillations

We performed a sensitivity analysis of the neural activity to the spatial frequency and the spatial phase of the extracellular oscillations. The methodology detailed in section 3.2.1 was applied with a different spatial frequency, $f_s = 0.4 \text{ mm}^{-1}$ instead of $0.8 \text{ mm}^{-1}$ as shown in **Figure 4A**, and with a different spatial phase, $\phi_s = \frac{\pi}{2}$ rad instead of $\frac{3\pi}{10}$ rad as shown in **Figure 4B**. The cross-correlation coefficients, computed as for **Figure 7C**, are shown in **Figure 10A**. The resulting number of action potentials is shown in **Figure 10B** for $f_s = 0.8 \text{ mm}^{-1}$ and $\phi_s = \frac{\pi}{2}$ rad and in **Figure 10C** for $f_s = 0.4 \text{ mm}^{-1}$ and $\phi_s = \frac{\pi}{2}$ rad. Even though the results are different for each parameter, especially for the higher amplitudes of the extracellular oscillations, the difference often appears to be rather stochastic, reflecting the highly nonlinear processing defined by the activation and inactivation functions of the active ion channels.

When $\phi_s = \frac{\pi}{2}$ rad, the peak of the extracellular oscillations is around the soma, while for $\phi_s = \frac{3\pi}{10}$ rad the peak occurs in the basal part of the dendritic tree (see **Figure 4**). The lack of significant difference between the two cases indicates that the influence of the extracellular oscillations at the soma is not more important than the influence at the dendritic tree. This suggests that for an accurate modeling of the nonsynaptic communication, the active properties of the dendritic tree should not be neglected.

The results for the apical dendrites show a great diversity, resulting mainly from the fact that the apical dendrites occupy a larger volume than the basal dendrites, thus the locations of the PSPs inputs are more spread for the apical dendrites. Also, variations across morphologies were more significant for the apical dendrites than for the basal ones, in particular for the CA1 neurons which have highly branched apical part of their dendritic tree along the entire somatodendritic axis (Spruston, 2008). This could account for the variability of the results when the stimulation is given through the apical dendrites.

## 4. DISCUSSION

Studies that investigate how the LFPs influence timing and number of action potentials are commonly based on stimulating neurons by injection of a constant current into the soma (Parra and Bikson, 2004; Radman et al., 2007; Anastassiou et al., 2010, 2011). In this paper, we show that the neural models with active conductances and dendrites stimulated by spatially distributed alpha function-shaped excitatory PSPs (EPSPs) have highly stochastic firing dynamics, which can to some degree be influenced by the parameters of the extracellular oscillations.

The weak extracellular oscillations studied in this paper correspond to those existing spontaneously in the brain. The strength of the extracellular oscillations in the mammalian brain is
<figure><img src="image_15.png" /><figcaption>**FIGURE 10** | (A) The correlation coefficients between the bins of 1 ms width and the phase of the extracellular oscillations. Extracellular oscillations with several different spatial properties for the CA3 and the CA1 neuron models are</figcaption></figure>

used. (B) The number of action potentials as a function of the amplitude when $f_s = 0.8 \text{ mm}^{-1}$ and $\phi_s = \frac{\pi}{2}$ rad. (C) The number of action potentials as a function of the amplitude when $f_s = 0.4 \text{ mm}^{-1}$ and $\phi_s = \frac{\pi}{2}$ rad.

commonly considered to go up to 3 mV mm$^{-1}$. This corresponds to 1 mV oscillations in the present study. Our results suggest that these weak oscillations have a modest influence on the timing of the action potentials. Larger fields, above 6 mV mm$^{-1}$ (which correspond to 2 and 4 mV oscillations in our study) have a significant influence on the timing. Such strong fields are seldom observed in the brain, e.g., during the hippocampal sharp

waves (Buzsáki, 1986; Sullivan et al., 2011). Another source of strong oscillations could be external stimulations, for example techniques such as TES, TMS or DBS. These methods could cause significant and coherent changes in the neural activity. We speculate that inducing a synchronous neural activity in the brain using strong extracellular electric fields could help to restore the normal neural activity.
It is important to emphasize that it has been shown that small correlations, which could be induced by the extracellular oscillations, can be amplified through network interactions and lead to increase in synchronization among the neurons (Parra and Bikson, 2004; Radman et al., 2007). The neuronal synchrony has been suggested to have a computational purpose, in particular to help maintain neural rhythmic activity (Weiss and Faber, 2010).

Recent findings suggest that the extracellular oscillations of the magnitude spontaneously present in the brain usually change the membrane potential only by 1–2 mV (Radman et al., 2007; Anastassiou et al., 2010; Buzsáki et al., 2012). This is rather small with respect to the depolarization required to reach the firing threshold level from the resting potential, which is usually about 15 mV. Nevertheless, since neurons are constantly bombarded with synaptic inputs of various amplitudes, their membrane potentials are generally significantly different from the resting potential and stay in a balanced state often close to the threshold level (Wilson, 2008). When the membrane potential is close to the firing threshold, the extracellular oscillations might affect neural firing more strongly. Our results suggest that the changes in the overall firing rate are rather small, generally below 5% even for the oscillations of 4 mV amplitude. The weaker fields usually do not cause an overall change greater than 1%, see **Figures 8A, 10B,C**.

Changes in the firing rates show that, for a given simulation setup, the influence of the extracellular oscillations is more significant when the synaptic inputs are received through the basal than through the apical dendrites. We speculate that this is due to the spatial properties of the extracellular oscillations. The phase of the oscillations is the same along the entire basal part of the dendritic tree, as well as along the soma. On the contrary, the phase of the oscillations along the apical part of the tree is the opposite to that around the soma, thus the resulting influence of the oscillations when synaptic inputs are on the apical tree gets balanced out.

We also observed that the correlation of the firing times with the phase of the extracellular oscillations for the CA3 neuron models is stronger than for the CA1 models. This emphasizes the importance of the ion channels which have a nonlinear dependence on the voltage and alter significantly the influence of the oscillations. The presented study did not investigate all the possible scenarios through which the extracellular oscillations could have a functional role. For instance, one possible functional role could be in gating some particular (hypothetical) ion channels that could be sensitive primarily to some properties of the extracellular oscillations (e.g., voltage or slope). Such ion channels could transfer this stochastic influence of the extracellular oscillations on the membrane potential into a more structured one, which could alter the neural activity in some functional way. In our study, the CA3 model, which contained more potassium and calcium channels, was generally influenced by the extracellular oscillations in a more coherent way.

The functional role of the spontaneous LFPs is an open question, but we have shown that even if the mean firing rate is not affected, the spontaneous LFPs could influence spike timings for up to several milliseconds. This could affect the neural synchronization and thus the entire network dynamics. This is important in the context of development of biologically inspired artificial

neural networks. A way to test whether an artificial neural network behaves in the same way as a biological one could be to evaluate its robustness to small changes in the spike timing.

It has been observed that the amplitude of extracellular oscillations is larger for smaller brains, i.e., it decreases from rat to cat, and from cat to primate (Buzsáki et al., 2012). One possible explanation for this is that, due to the brain size, bodies of the pyramidal neurons are vertically the best aligned in the rat brain, then in the cat brain, and in the human brain. Since studies, including this one, suggest that the amplitude of the oscillations is crucial for their coherent influence on firing properties, a decrease in the amplitude in larger brains can be an argument against the functional role of the oscillations.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fncom.2014.00009/abstract

**Figure S1 | Spikes from Figure 9 in the main part of the paper binned into 10 ms time bins (only first 300 out of 1000 ms are shown).** Bins corresponding with different amplitudes of the extracellular oscillations are color coded. Phase of the extracellular oscillations at the soma is indicated with dashed black line.

**Figure S2 | Difference of the time bins from Figure S1.**

**Figure S3 | Influence of the frequency of the extracellular oscillations on the number of action potentials for the stimulation of the apical dendrites of CA3 neuron model.** This figure is produced in analogous way as **Figure 9A** in the main part of the paper. (A) The stimulation of the CA1 apical dendrites. (B) The stimulation of the CA3 basal dendrites. (C) The stimulation of the CA3 apical dendrites.

## REFERENCES

Allen, E. A., Pasley, B. N., Duong, T., and Freeman, R. D. (2007). Transcranial magnetic stimulation elicits coupled neural and hemodynamic consequences. *Science* 317, 1918–1921. doi: 10.1126/science.1146426

Anastassiou, C. A., Montgomery, S. M., Barahona, M., Buzsáki, G., and Koch, C. (2010). The effect of spatially inhomogeneous extracellular electric fields on neurons. *J. Neurosci.* 30, 1925–1936. doi: 10.1523/JNEUROSCI.3635-09.2010

Anastassiou, C. A., Perin, R., Markram, H., and Koch, C. (2011). Ephaptic coupling of cortical neurons. *Nat. Neurosci.* 14, 217–223. doi: 10.1038/nn.2727

Barker, A. T., Jalinous, R., and Freeston, I. L. (1985). Non-invasive magnetic stimulation of human motor cortex. *Lancet* 1, 1106–1107. doi: 10.1016/S0140-6736(85)92413-4

Benabid, A. L., Chabardes, S., Mitrofanis, J., and Pollak, P. (2009). Deep brain stimulation of the subthalamic nucleus for the treatment of Parkinson's disease. *Lancet Neurol.* 8, 67–81. doi: 10.1016/S1474-4422(08)70291-6

Brette, R., and Destexhe, A. (2011). *Handbook of Neural Activity Measurement*. Cambridge: Cambridge University Press.

Buzsáki, G. (1986). Hippocampal sharp waves: their origin and significance. *Brain Res.* 398, 242–252. doi: 10.1016/0006-8993(86)91483-6

Buzsáki, G., Anastassiou, C. A., and Koch, C. (2012). The origin of extracellular fields and currents – EEG, ECOG, LFP and spikes. *Nat. Rev. Neurosci.* 13, 407–420. doi: 10.1038/nrn3241

Buzsáki, G., and Draguhn, A. (2004). Neuronal oscillations in cortical networks. *Science* 304, 1926–1929. doi: 10.1126/science.1099745

Capotosto, P., Babiloni, C., Romani, G. L., and Corbetta, M. (2009). Frontoparietal cortex controls spatial attention through modulation of anticipatory alpha rhythms. *J. Neurosci.* 29, 5863–5872. doi: 10.1523/JNEUROSCI.0539-09.2009

Carnevale, N. T., and Hines, M. L. (2006). *The NEURON Book*. Cambridge: Cambridge University Press. doi: 10.1017/CBO9780511541612
Deans, J. K., Powell, A. D., and Jefferys, J. G. (2007). Sensitivity of coherent oscillations in rat hippocampus to AC electric fields. *J. Physiol.* 583, 555–565. doi: 10.1113/jphysiol.2007.137711

Deco, G., and Thiele, A. (2009). Attention: oscillations and neuropharmacology. *Eur. J. Neurosci.* 30, 347–354. doi: 10.1111/j.1460-9568.2009.06833.x

Ebersole, J. S., and Milton, J. (2003). “The electroencephalogram (EEG): a measure of neural synchrony,” in *Epilepsy as a Dynamic Disease* (Berlin; Heidelberg: Springer), 51–68. doi: 10.1007/978-3-662-05048-4_5

Fries, P., Nikolic, D., and Singer, W. (2007). The gamma cycle. *Trends Neurosci.* 30, 309–316. doi: 10.1016/j.tins.2007.05.005

Fröhlich, F., and McCormick, D. A. (2010). Endogenous electric fields may guide neocortical network activity. *Neuron* 67, 129–143. doi: 10.1016/j.neuron.2010.06.005

Hallett, M. (2000). Transcranial magnetic stimulation and the human brain. *Nature* 406, 147–150. doi: 10.1038/35018000

Hemond, P., Epstein, D., Boley, A., Migliore, M., Ascoli, G. A., and Jaffe, D. B. (2008). Distinct classes of pyramidal cells exhibit mutually exclusive firing patterns in hippocampal area CA3b. *Hippocampus* 18, 411–424. doi: 10.1002/hipo.20404

Hines, M. L., and Carnevale, N. T. (1997). The NEURON simulation environment. *Neural Comput.* 9, 1179–1209. doi: 10.1162/neco.1997.9.6.1179

Izhikevich, E. M. (2006). *Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting*. Cambridge: The MIT Press.

Joucla, S., and Yvert, B. (2012). Modeling extracellular electrical neural stimulation: from basic understanding to mea-based applications. *J. Physiol.* 106, 146–158. doi: 10.1016/j.jphysparis.2011.10.003

Kahana, M. J., Seelig, D., and Madsen, J. R. (2001). Theta returns. *Curr. Opin. Neurobiol.* 11, 739–744. doi: 10.1016/S0959-4388(01)00278-1

Kayser, C., Montemurro, M. A., Logothetis, N. K., and Panzeri, S. (2009). Spike-phase coding boosts and stabilizes information carried by spatial and temporal spike patterns. *Neuron* 61, 597–608. doi: 10.1016/j.neuron.2009.01.008

Kim, S., Guzman, S. J., H., H., and Jonas, P. (2012). Active dendrites support efficient initiation of dendritic spikes in hippocampal ca3 pyramidal neurons. *Nat. Neurosci.* 15, 600–606. doi: 10.1038/nn.3060

Kirov, R., Weiss, C., Siebner, H. R., Born, J., and Marshall, L. (2009). Slow oscillation electrical brain stimulation during waking promotes EEG theta activity and memory encoding. *Proc. Natl. Acad. Sci. U.S.A.* 106, 15460–15465. doi: 10.1073/pnas.0904438106

Koepsell, K., Wang, X., Hirsch, J. A., and Sommer, F. T. (2010). Exploring the function of neural oscillations in early sensory systems. *Front. Neurosci.* 4:53. doi: 10.3389/neuro.01.010.2010

Lisman, J., and Buzsáki, G. (2008). A neural coding scheme formed by the combined function of gamma and theta oscillations. *Schizophr. Bull.* 34, 974–980. doi: 10.1093/schbul/sbn060

Mainen, Z. F., and Sejnowski, T. J. (1998). “Modeling active dendritic processes in pyramidal neurons,” in *Methods in Neuronal Modeling, 2nd Edn.*, eds C. Koch and I. Segev (Cambridge: MIT Press), 170–209.

Markram, H. (2006). The blue brain project. *Nat. Rev. Neurosci.* 9, 153–160. doi: 10.1038/nrn1848

Migliore, M., Ferrante, M., and Ascoli, G. A. (2005). Signal propagation in oblique dendrites of cal pyramidal cells. *J. Neurophysiol.* 94, 4145–4155. doi: 10.1152/in.00521.2005

Montemurro, M. A., Rasch, M. J., Murayama, Y., Logothetis, N. K., and Panzeri, S. (2008). Phase-of-firing coding of natural visual stimuli in primary visual cortex. *Curr. Biol.* 18, 375–380. doi: 10.1016/j.cub.2008.02.023

Nyhus, E., and Curran, T. (2010). Functional role of gamma and theta oscillations in episodic memory. *Neurosci. Biobehav. Rev.* 34, 1023–1035. doi: 10.1016/j.neubiorev.2009.12.014

Osipova, D., Takashima, A., Oostenveld, R., Fernandez, G., Maris, E., and Jensen, O. (2006). Theta and gamma oscillations predict encoding and retrieval of declarative memory. *J. Neurosci.* 26, 7523–7531. doi: 10.1523/JNEUROSCI.1948-06.2006

Ozen, S., Sirota, A., Belluscio, M. A., Anastassiou, C. A., Stark, E., Koch, C., et al. (2010). Transcranial electric stimulation entrains cortical neuronal populations in rats. *J. Neurosci.* 30, 11476–11485. doi: 10.1523/JNEUROSCI.5252-09.2010

Parra, L. C., and Bikson, M. (2004). Model of the effect of extracellular fields on spike time coherence. *Conf. Proc. IEEE Eng. Med. Biol. Soc.* 6, 4584–4587. doi: 10.1109/IEMBS.2004.1404271

Peelle, J. E., and Davis, M. H. (2012). Neural oscillations carry speech rhythm through to comprehension. *Front. Psychol.* 3:320. doi: 10.3389/fpsyg.2012. 00320

Pyapali, G., Sik, A., Penttonen, M., Buzsáki, G., and Turner, D. (1998). Dendritic properties of hippocampal cal pyramidal neurons in the rat: intracellular staining *in vivo* and *in vitro*. *J. Comp. Neurol.* 391, 335–352. doi: 10.1002/(SICI)1096-9861(19980216)391:3<335::AID-CNE4>3.0.CO;2-2

Radman, T., Su, Y., An, J. H., Parra, L. C., and Bikson, M. (2007). Spike timing amplifies the effect of electric fields on neurons: implications for endogenous field effects. *J. Neurosci.* 27, 3030–3036. doi: 10.1523/JNEUROSCI.0095-07.2007

Reato, D., Rahman, A., Bikson, M., and Parra, L. C. (2010). Low-intensity electrical stimulation affects network dynamics by modulating population rate and spike timing. *J. Neurosci.* 30, 15067–15079. doi: 10.1523/JNEUROSCI.2059-10.2010

Samsonovich, A. V., and Ascoli, G. A. (2003). Statistical morphological analysis of hippocampal principal neurons indicates cell-specific repulsion of dendrites from their own cell. *J. Neurosci. Res.* 71, 173–187. doi: 10.1002/jnr.10475

Schack, B., Vath, N., Petsche, H., Geissler, H. G., and Möller, E. (2002). Phase-coupling of theta-gamma EEG rhythms during short-term memory processing. *Int. J. Psychophysiol.* 44, 143–163. doi: 10.1016/S0167-8760(01)00199-4

Schroeder, C. E., and Lakatos, P. (2009). Low-frequency neuronal oscillations as instruments of sensory selection. *Trends Neurosci.* 32, 9–18. doi: 10.1016/j.tins.2008.09.012

Sederberg, P. B., Kahana, M. J., Howard, M. W., Donner, E. J., and Madsen, J. R. (2003). Theta and gamma oscillations during encoding predict subsequent recall. *J. Neurosci.* 23, 10809–10814.

Spruston, N. (2008). Pyramidal neurons: dendritic structure and synaptic integration. *Nat. Rev. Neurosci.* 9, 206–221. doi: 10.1038/nrn2286

Sullivan, D., Csicsvari, J., Mizuseki, K., Montgomery, S., Diba, K., and Buzsáki, G. (2011). Relationships between hippocampal sharp waves, ripples, and fast gamma oscillation: influence of dentate and entorhinal cortical activity. *J. Neurosci.* 31, 8605–8616. doi: 10.1523/JNEUROSCI.0294-11.2011

Thut, G., Nietzel, A., Brandt, S. A., and Pascual-Leone, A. (2006). alpha-band electroencephalographic activity over occipital cortex indexes visuospatial attention bias and predicts visual target detection. *J. Neurosci.* 26, 9494–9502. doi: 10.1523/JNEUROSCI.0875-06.2006

Uhlhaas, P., Roux, F., Rodriguez, E., Rotarska-Jagiela, A., and Singer, W. (2010). Neural synchrony and the development of cortical networks. *Trends Cogn. Sci.* 14, 72–80. doi: 10.1016/j.tics.2009.12.002

Walsh, V., and Cowey, A. (2000). Transcranial magnetic stimulation and cognitive neuroscience. *Nat. Rev. Neurosci.* 1, 73–79. doi: 10.1038/35036239

Wang, X.-J. J. (2010). Neurophysiological and computational principles of cortical rhythms in cognition. *Physiol. Rev.* 90, 1195–1268. doi: 10.1152/phys-rev.00035.2008

Ward, L. M. (2003). Synchronous neural oscillations and cognitive processes. *Trends Cogn. Sci.* 7, 553–559. doi: 10.1016/j.tics.2003.10.012

Weiss, S. A., and Faber, D. S. (2010). Field effects in the CNS play functional roles. *Front. Neural Circuits* 4:15. doi: 10.3389/fncir.2010.00015

Wilson, C. (2008). Up and down states. *Scholarpedia* 3, 1410. doi: 10.4249/scholarpedia.1410

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Received: 03 July 2013; accepted: 16 January 2014; published online: 07 February 2014.

Citation: Tiganj Z, Chevallier S and Monacelli E (2014) Influence of extracellular oscillations on neural communication: a computational perspective. *Front. Comput. Neurosci.* 8:9. doi: 10.3389/fncom.2014.00009

This article was submitted to the journal *Frontiers in Computational Neuroscience*.

Copyright © 2014 Tiganj, Chevallier and Monacelli. This is an open-access article distributed under the terms of the Creative Commons Attribution License

(CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No

use, distribution or reproduction is permitted which does not comply with these terms.