# Pharmacodynamics of the Glutamate Receptor Antagonists in the Rat Barrel Cortex

Daria Vinokurova, Andrey Zakharov, Julia Lebedeva, Gulshat F Burkhanova, Kseniya A Chernova, Nailya Lotfullina, Rustem Khazipov, Guzel Valeeva

► To cite this version:

Daria Vinokurova, Andrey Zakharov, Julia Lebedeva, Gulshat F Burkhanova, Kseniya A Chernova, et al. Pharmacodynamics of the Glutamate Receptor Antagonists in the Rat Barrel Cortex. Frontiers in Pharmacology, 2018, 9, 10.3389/fphar.2018.00698. [hal-01963826](https://hal HAL Id:)

Submitted on 21 Dec 2018

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>

Check for
updates

# Pharmacodynamics of the Glutamate Receptor Antagonists in the Rat Barrel Cortex

Daria Vinokurova<fnref n="1" />,²*, Andrey V. Zakharov<fnref n="1" />†, Julia Lebedeva<fnref n="1" />, Gulshat F. Burkhanova<fnref n="1" />, Kseniya A. Chernova<fnref n="1" />, Nailya Lotfullina<fnref n="1" />,², Rustem Khazipov<fnref n="1" />,²* and Guzel Valeeva<fnref n="1" />*

Epipial application is one of the approaches for drug delivery into the cortex. However, passive diffusion of epipially applied drugs through the cortical depth may be slow, and different drug concentrations may be achieved at different rates across the cortical depth. Here, we explored the pharmacodynamics of the inhibitory effects of epipially applied ionotropic glutamate receptor antagonists CNQX and dAPV on sensory-evoked and spontaneous activity across layers of the cortical barrel column in urethane-anesthetized rats. The inhibitory effects of CNQX and dAPV were observed at concentrations that were an order higher than in slices *in vitro*, and they slowly developed from the cortical surface to depth after epipial application. The level of the inhibitory effects also followed the surface-to-depth gradient, with full inhibition of sensory evoked potentials (SEPs) in the supragranular layers and L4 and only partial inhibition in L5 and L6. During epipial CNQX and dAPV application, spontaneous activity and the late component of multiple unit activity (MUA) during sensory-evoked responses were suppressed faster than the short-latency MUA component. Despite complete suppression of SEPs in L4, sensory-evoked short-latency multiunit responses in L4 persisted, and they were suppressed by further addition of lidocaine suggesting that spikes in thalamocortical axons contribute ~20% to early multiunit responses. Epipial CNQX and dAPV also completely suppressed sensory-evoked very fast (~500 Hz) oscillations and spontaneous slow wave activity in L2/3 and L4. However, delta oscillations persisted in L5/6. Thus, CNQX and dAPV exert inhibitory actions on cortical activity during epipial application at much higher concentrations than *in vitro*, and the pharmacodynamics of their inhibitory effects is characterized by the surface-to-depth gradients in the rate of development and the level of inhibition of sensory-evoked and spontaneous cortical activity.

**Keywords:** glutamate receptor antagonists, barrel cortex, epipial application, drug delivery, sensory-evoked potential, fast oscillations, slow wave activity

**OPEN ACCESS**

**Edited by:**
Luc Zimmer,
Université Claude Bernard Lyon 1,
France

**Reviewed by:**
Werner Kilb,
Johannes Gutenberg Universität
Mainz, Germany
Dirk Isbrandt,
Deutsches Zentrum für
Neurodegenerative Erkrankungen
(DZNE), Germany

*Correspondence:
Rustem Khazipov
roustem.khazipov@inserm.fr
Guzel Valeeva
gurvaleeva@kpfu.ru

**Specialty section:**
This article was submitted to
Experimental Pharmacology
and Drug Discovery,
a section of the journal
Frontiers in Pharmacology

Received: 30 March 2018

Accepted: 08 June 2018

Published: 03 July 2018

**Citation:**
Vinokurova D, Zakharov AV, Lebedeva J, Burkhanova GF, Chernova KA, Lotfullina N, Khazipov R and Valeeva G (2018) Pharmacodynamics of the Glutamate Receptor Antagonists in the Rat Barrel Cortex. *Front. Pharmacol.* 9:698.
doi: 10.3389/fphar.2018.00698

## INTRODUCTION

The epipial application of pharmacological agents is widely used for drug delivery into the cortex (Pockberger et al., 1984; Andreasen et al., 1989; Hablitz and Sutor, 1990; Kohling et al., 1993; Conti and Minelli, 1994; Ikeda et al., 2002). This approach appears to be particularly useful for compounds, which are not very permeable through the brain–blood barrier (BBB), and therefore

<footnote n="1">¹ Laboratory of Neurobiology, Kazan Federal University, Kazan, Russia, ² Mediterranean Institute of Neurobiology - National Institute of Health and Medical Research, Aix-Marseille University, UMR1249, Marseille, France</footnote>

<footnote>† These authors have contributed
equally to this work.</footnote>
are not effective in systemic delivery. While in a clinical setting the use of epipial drug delivery is limited, this approach may be of interest during drug development as the BBB-impermeable molecule can easily be tested in various disease models with future development of BBB-permeable precursor forms or other delivery methods if the compound shows its efficacy. However, the primary requirement for using epipial application is to know how well the drug penetrates into the brain tissue after application, what the time course of penetration at different depths is, and what drug concentrations should be used during epipial application for effective concentrations at the cortical depth.

During epipial application, drugs passively diffuse through the cortical tissue mainly via extracellular space compartments (Ohata and Marmarou, 1992; Groothuis et al., 2007). However, the diffusion of molecules within the cortical tissue is more complex than diffusion in a free medium. The diffusion within brain tissue is known to be constrained by many factors including tortuosity of extracellular space, direction of axon bundles, uptake into cells, transient binding to receptors or transporters, interaction with negative-charged elements of the extracellular matrix, trapping of molecules in dead-space microdomains, solubility in lipids and electrical charge of diffusing molecules (Sykova and Nicholson, 2008). Therefore, the effective diffusion coefficients for compounds in brain tissue are, on average, half of that in the free medium and may vary depending on brain region and animal age (ibid). In keeping with the relatively slow diffusion rates and the thickness of the cortex, drug concentration may develop at different rates and attain different levels at different depths of the cortex. This raises the question of the pharmacodynamics of drugs during epipial application. To elucidate this question, we have chosen a combination of the glutamate ionotropic AMPA/kainate and NMDA receptor antagonists CNQX and dAPV as a model drug, using the spontaneous and sensory-evoked activity in different cortical layers as a readout of their penetration into the cortex.

AMPA/kainate and NMDA receptors are two types of ionotropic glutamate receptors at glutamatergic synapses, pharmacological blockade of which with the selective antagonists CNQX (Honore et al., 1988) and dAPV (Davies et al., 1981) completely suppresses spontaneous and evoked excitatory glutamatergic synaptic transmission (Neuman et al., 1988; Andreasen et al., 1989; Hablitz and Sutor, 1990; Long et al., 1990; Llano et al., 1991; Traynelis et al., 2010). Glutamatergic synapses are critical for generation of the network driven activities involving local and large-scale intracortical excitatory connections and thalamic glutamatergic inputs. In the barrel cortex, which has highly organized excitatory connectivity within and between cortical layers, internally generated spontaneous activity and sensory evoked responses mediated by the thalamus display patterned cross-layer network dynamics (Kandel and Buzsaki, 1997; Steriade, 2001; Buzsaki, 2006; Sakata and Harris, 2009). For example, in the rodent barrel cortex, whisker driven thalamic inputs that mainly project to L4 neurons, are sequentially processed via the canonical L4 – L2/3 – L5 glutamatergic pathway [for review, (Douglas and Martin, 2004; Brecht, 2007; Petersen, 2007; Feldmeyer, 2012;

Feldmeyer et al., 2013)]. Also, direct input from the thalamus drives activity in the infragranular layers (Bureau et al., 2006; Brecht, 2007; Meyer et al., 2010; Oberlaender et al., 2012; Constantinople and Bruno, 2013; Crocker-Buque et al., 2015). Activation of the thalamic synaptic inputs and intracortical synapses, as well as voltage-gated conductances in the pre- and postsynaptic neurons underlie transmembrane currents that are at the origin of the extracellular local field potentials (LFPs) recorded at different depths of the cortical column during spontaneous and sensory-evoked activities (Buzsaki et al., 2012).

In the present study, we attempted to determine the effective concentration, exposure times, and inhibition levels produced by epipially applied CNQX/dAPV in different layers of the cortical barrel column. Also, analyzing CNQX/dAPV effects through the time course of their penetration into the cortex we made several observations that could be of interest in cortical barrel network physiology.

# MATERIALS AND METHODS

## Ethical Approval

This work has been carried out in accordance with EU Directive 2010/63/EU for animal experiments, and all animal-use protocols were approved by the French National Institute of Health and Medical Research (INSERM, protocol N007.08.01) and Kazan Federal University on the use of laboratory animals (ethical approval by the Institutional Animal Care and Use Committee of Kazan State Medical University N9-2013).

## Surgery

Wistar rats of either sex from postnatal day P19–39 (P0 = day of birth) were used. The rat cerebral cortex acquires the diffusion characteristics of the adult brain by P20–P23 (Lehmenkuhler et al., 1993; Vorisek and Sykova, 1997b; Mazel et al., 2002). In addition, we observed no correlation between animal age and the magnitude of the CNQX/dAPV effect, both on L4 sensory evoked potential (SEP) amplitude (r = 0.09, p = 0.80) and L4 evoked multiple unit activity (MUA) (r = -0.47, p = 0.14) within a given age range. Therefore, P19–39 animals were pooled into one group. Surgery was performed under isoflurane anesthesia (4% for induction, 2% for maintenance), and urethane (1 g/kg, i.p.) was injected at the end of surgery. The skull of the animal was cleaned of skin and periosteum, dried and covered with cyanacrylamide glue except for a 4–9 mm² window above the left barrel cortex. Then the skull was covered by dental cement (Grip Cement, Dentsply Sirona, Milford, DE, United States). A metal ring (15 mm inner diameter) was fixed to the rat’s head by dental cement. The area inside the ring was left cement-free. After surgery, the animals were warmed, and left for 1 h to recover. The head was attached to a ball-joint holder by the metal ring. During recordings animals were heated via a thermal pad (35–36°C). A chlorided silver wire, placed in the visual cortex, served as the ground electrode.

A cranial window ~3–4 mm in diameter was drilled above the barrel cortex (window edges: AP –0.5 to –4.5 mm; lateral 2.5–6.5 mm from bregma) (Paxinos and Watson, 2007; Khazipov
et al., 2015), and a 1–2 mm-long dura incision was made in the middle of the cranial window using a 27G needle. Dura dissection was done carefully to avoid bleeding and the formation of blood clots on the brain surface, which could obstruct drug diffusion. The cranial window was encircled with a 1 mm-high cement wall to form an epipial chamber. Artificial cerebrospinal fluid (ACSF) during control recordings and CNQX/dAPV dissolved in ACSF were applied to the epipial chamber at 20–40 µl volumes every 5–10 min.

## Extracellular Recordings

Extracellular LFPs and MUA were recorded from a single barrel column using 16-site linear silicon probes (100 µm separation distance between recording sites, NeuroNexus Technologies, Ann Arbor, MI, United States). The probe was inserted vertically to the cortical surface to a depth of 1300–1900 µm (depending on animal age). The signals from extracellular recordings were amplified and filtered (×10,000; 0.15 Hz–10 kHz) using a Digital Lynx (Neuralynx, United States) amplifier, digitized at 32 kHz and saved on a PC for *post hoc* analysis. The whiskers were trimmed to a length of 0.8–1.5 mm and were stimulated using a piezoelectric bending actuator (Noliac, Denmark) using 200 ms square pulses with 5–10 s intervals. A needle (22G) was glued to the end of piezo actuator and the tip of the whisker was inserted into the blunt tip of the needle. The principal whisker (PW) was identified by the shortest latency MUA responses in layer 4 evoked by single whisker deflection.

## Preliminary Selection of the Concentrations and the Mode of Application

To estimate the concentrations of CNQX and dAPV which efficiently inhibit cortical activity, we first applied drugs using the “superfused cortex” preparation (Khazipov and Holmes, 2003; Minlebaev et al., 2007), where the dura is completely removed within the entire cranial window and the cortical surface is continuously superfused. Superfusion of the cortex with 40 µM CNQX/100 µM dAPV, which is more than twice the concentration that completely suppress cortical activity in L4 in neonatal rats (Minlebaev et al., 2007, 2009), caused reduction of the L4 SEP slope to only 60 ± 17% (n = 4; *p* = 0.13), and suppression of spontaneous L4 MUA to only 27 ± 7% (n = 4; *p* = 0.13) of control values as estimated at 45–60 min after drug application. These moderate effects likely reflect a deceleration of drug diffusion with age (Lehmenkuhler et al., 1993; Vorisek and Sykova, 1997b). A further increase of drug concentrations to 80 µM CNQX/200 µM dAPV also failed to completely inhibit activity in L4 (not shown). As further increase of concentration would be extremely wasteful if applied via perfusion, we next tested two concentrations using static epipial drug application, 170 µM CNQX + 700 µM dAPV and 0.5 mM CNQX + 2 mM dAPV. We found that the reduction of L4 SEP slope to 18 ± 17% (n = 4) of control values produced by 170 µM CNQX/700 µM dAPV did not significantly differ (*p* = 0.83) from SEP slope reduction to 2 ± 2% (n = 9) produced by 0.5 mM CNQX/2 mM dAPV. Similarly, the inhibition of spontaneous L4

MUA to 7 ± 6% and 3 ± 1% of control values caused by these two combinations of drug concentrations, respectively, were not significantly different (*p* = 0.70). Hereafter we describe the effects of epipially applied CNQX and dAPV at these two concentrations with the results pooled together in group statistics.

## Histology

After recordings the animals were deeply anesthetized with urethane (3 g/kg, intraperitoneally) and perfused through the left cardiac ventricle with 4% paraformaldehyde and 1% glutaraldehyde. The brains were removed and left for fixation in the solution containing glutar- and paraformaldehyde at room temperature. Then the brains were rinsed in PBS and cut into 100 µm-thick slices using a Vibratome (Leica Biosystems, Wetzlar, Germany). The slices, cut in the coronal plane, were used for identification of probe position. Electrode positions were identified from the DiI tracks overlaid on microphotographs of the sections in oblique light or after cresyl-violet staining and depth of the short-latency spikes during sensory-evoked responses.

## Data Analysis and Statistics

Raw data were preprocessed using custom-developed routines in MATLAB environment (MathWorks, United States). The wide-band signal was down-sampled to 1000 Hz and used as the LFP signal. Positive polarity is shown as up in all figures. For action potential detection, the raw wide-band signal was filtered (bandpass 300–5000 Hz) and negative events exceeding five standard deviations calculated over the most silent 1 s length segment of the filtered trace were considered as spikes (MUA). SEPs were detected as the first LFP troughs of sensory evoked responses. LFPs and extracellular spikes were analyzed by custom-written, MATLAB-based programs. Current-source density (CSD) analysis across cortical depth was used to eliminate volume conduction and localize synaptic currents. CSD was computed according to a differential scheme for the second spatial derivative along recording sites. The up-states were detected by deep channels located in L6. Negative LFP deflections exceeding 100 µV and accompanied by MUA were considered as up-states. The early and late components of evoked MUA response were calculated individually for each animal because of slight variations in the temporal properties of sensory response across animals. The early evoked MUA was defined as a time interval where the density of L4 spikes persisting after the full blockade of SEP by CNQX/dAPV exceeded 10% of the level of peak MUA density. Thus, the average onset and duration of the early MUA period were of 7.0 ± 0.2 and 4.3 ± 0.5 ms, respectively (n = 11). The time point corresponding to a negative extremum of derivative from evoked MUA density in control was chosen as the onset of the late evoked MUA period, and, on average, accounted for 5.0 ± 0.3 ms after the end of the early evoked MUA period. The late evoked MUA was assessed within the following 20 ms time window.

Statistical analysis was performed using the Matlab Statistics toolbox. The two-side Wilcoxon rank sum test was performed to assess the significance of differences between groups of data with the level of significance kept at *p* < 0.05. Pearson's correlation
coefficient (r) was used to describe linear relation between variables.

## Drugs

Urethane, CNQX (6-Cyano-7-nitroquinoxaline-2,3-dione), dAPV (D(-)-2-Amino-5-phosphonopentanoic acid) and lidocaine N-ethyl bromide solution were purchased from Sigma. Isoflurane was purchased from Baxter.

## RESULTS

We addressed the pharmacodynamics of the ionotropic glutamate receptor antagonists (170–500 µM CNQX and 0.7–2 mM dAPV) after their epipial administration using linear multisite silicone probe recordings of the LFP signals and MUA at different depths of the cortical barrel column. As a readout of the drug penetration we used several parameters of sensory-evoked responses to brief mechanical stimulation of the PW and spontaneous activity in different cortical layers. These parameters included: (i) slope and amplitude of the SEP; (ii) early and late components of the sensory-evoked MUA response; (iii) characteristics of spontaneous delta-wave activity and frequency of spontaneous MUA. Each of these parameters, as well as the time for attaining half-maximal effect (T₁/₂) and the steady state levels of inhibition were calculated for different cortical layers through the time course of 1 h long drug application.

## SEP

Under control conditions, PW deflection evoked SEP with a characteristic CSD profile with two early sinks in L4 and the L5/6 border followed by sinks in L2/3 and L5 (Figures 1A, 3A; Di et al., 1990; Castro-Alamancos and Oldford, 2002; Roy et al., 2011; Reyes-Puerta et al., 2015). The largest sensory response was observed in L4, where average L4 SEP amplitude and slope were of 724 ± 145 µV and 197 ± 41 µV/ms, respectively (n = 11). After epipial combined CNQX and dAPV application (170–500 µM CNQX and 0.7–2 mM dAPV), the inhibitory effect of the blockers on SEP was first observed in superficial layers 2/3, and then it successively involved L4 and L5/6 (Figures 1, 2). The T₁/₂ of SEP slope inhibition in L2/3 and L4 were 3.9 ± 0.5 and 9.4 ± 0.5 min after epipial application of blockers, respectively, while in L5/6 it accounted for 28.0 ± 1.9 min (n = 11; Figures 2A,C). The level of SEP inhibition at 45–60 min of application was stronger at the superficial than deep layers attaining –4 ± 1%, 0 ± 1%, and 32 ± 7% of the control values in L2/3, L4, and L5/6 (L5 and L6 pooled together), respectively (n = 11; Figure 2A). It is noteworthy that the level of SEP inhibition in the infragranular layers also showed a surface to depth gradient from L5a to 6B (Figure 1B, right panel).

## Sensory-Evoked MUA

In parallel with a reduction of SEP, epipial CNQX/dAPV application induced an inhibition of sensory-evoked MUA, which also developed earlier and was in general greater at the superficial layers (Figure 2B). Thus, the T₁/₂ of sensory-evoked MUA inhibition was 2.3 ± 0.2, 5.3 ± 0.8, and 13.4 ± 2.2 min in L2/3,

L4, and L5/6, respectively (n = 11). The steady state level of inhibition of sensory-evoked MUA was also achieved earlier in the superficial layers attaining 5 ± 3, 20 ± 4, and 24 ± 6% of control values in L2/3, L4, and L5/6, respectively. Interestingly, the late component of the sensory-evoked MUA response was inhibited faster and more strongly than the early component (Figures 1B, 2C). Also, the early MUA component was almost completely suppressed in L2/3 (p = 0.014) and L5a (p = 0.002), but it was less affected in L4 (attaining 20 ± 4% of control values compared to 5 ± 3% in L2/3, p = 0.001; and to 9 ± 6% in L5a, p = 0.002; Figures 1A,B).

Paradoxically, the early MUA component in L4 persisted even after complete suppression of the L4 SEP implying that the spikes at the onset of the sensory-evoked response can occur despite complete blockade of thalamocortical synaptic transmission (Figures 1, 3). In individual experiments, early evoked L4 MUA density in the presence of CNQX/dAPV varied from 19 to 54% of the control values with the average value of 38 ± 3% (average early MUA density was 2.7 ± 0.4 spikes/ms in control conditions and 1.1 ± 0.2 spikes/ms after drug application; n = 11; p = 0.003; Figure 3C). The majority of the CNQX/dAPV insensitive L4 spikes fired mainly within a time window corresponding to the rising phase of evoked MUA in control peaking at 7.8 ± 0.3 ms after the stimulus that was 0.5 ± 0.1 ms earlier than the peak of evoked MUA response in control conditions (n = 11; p = 0.006; Figure 3D). We hypothesized that this CNQX/dAPV-resistant component of the sensory evoked early MUA response could be generated by spikes in thalamocortical axons. In keeping with this hypothesis, after addition of the voltage gated sodium channel blocker lidocaine (2%) to the epipial solution, the peak of the CNQX/dAPV-resistant early evoked L4 MUA response was suppressed from 31 ± 5% of control values to 8 ± 5% during the 15–25 min after lidocaine application with T1/2 of 8 ± 2 min (n = 4; p = 0.03; Figures 3B–E).

## Fast Oscillations

In control conditions, sensory evoked responses were characterized by fast oscillations (FOs) as described earlier (Kandel and Buzsaki, 1997; Jones and Barth, 1999; Jones et al., 2000; Barth, 2003). FOs were observed in all cortical layers, were locked to the stimulus onset and could be clearly seen on the average multiunit response and LFP oscillations after high-pass filtering (Figures 1A, 3A, 4). A population average L4 MUA power spectrum revealed a peak frequency of FOs at 390 ± 45 Hz (n = 11; Figure 4B). During CNQX/dAPV application, along with a suppression of SEP and MUA response, L4 FOs reduced in power and duration and were completely eliminated after full blockade of SEP (Figure 4B). It is noteworthy, that in individual animals, the power spectrum of L4 multiunit responses displayed two prominent peaks at 451 ± 49 and 645 ± 40 Hz (n = 11), both were also suppressed by CNQX/dAPV (data not shown).

## Spontaneous Activity

Spontaneous activity in the cortical barrel column was characterized by slow wave oscillations. Up-states of slow oscillations typically involved the entire column, and were characterized by sinks and an increase in MUA in all layers with
<figure><img src="image_3.png" /></figure>

**FIGURE 1** | Depth profile of CNQX and dAPV action on barrel cortex activity. (A) Top, Stimulus-triggered LFP averages (black traces) of sensory responses evoked by principal whisker deflection across cortical depths of corresponding barrel column overlaid on color-coded current source density plot (CSD) in control, 7 and 25 min after 500 µM CNQX and 2 mM dAPV application (25 responses were averaged for each CSD plot). Bottom, corresponding stimulus-triggered averages for MUA. Stimulus onset is indicated by the black arrowhead and vertical magenta line. The cortical layer borders are shown left of the CSD and MUA density plots, and recording electrode depths are shown on the right. Red square brackets indicate the electrodes of the recording that were used for building group data plots on panel (B). The gray bars above MUA density plots limit the time windows for early and late MUA density calculation. (B) Group data (from 11 animals and two concentrations, 170 µM CNQX/700 µM dAPV and 0.5 mM CNQX/2 mM dAPV) on the magnitude of CNQX and dAPV effect on early and late components of the sensory-evoked MUA, spontaneous MUA and SEP slope at the beginning of drug application (4–9 min, the time point corresponding to SEP slope reduction by 20% of control values) and at the end of 1 h of application (45–60 min). Note the gradual development of CNQX/dAPV effect from superficial toward deep cortical layers, and the persistence of early evoked L4 spikes after full blockade of SEP and nearly complete block of evoked MUA in L2/3 and L5a.
<figure><img src="image_4.png" /></figure>

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

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

**FIGURE 2** | The time course of CNQX and dAPV effect on sensory evoked responses through the depth of the barrel cortex. The decrease in SEP slope (A) and MUA during SEP (B) in different layers of cortical barrel column following application of CNQX and dAPV. The application time is shown with a gray bar above the plot. White circles on the curves indicate a half-maximal effect time (T¹/²). The colored areas surrounding each curve show SE bands. (C) Statistical data on T¹/² of CNQX and dAPV effect on SEP slope (**left**), early evoked MUA (**middle**) and late evoked MUA (**right**) across cortical layers. The open circles correspond to individual experiments. The medians of boxplots are shown by red lines. (A–C), pooled data from 11 animals and two concentrations, 170 µM CNQX/700 µM dAPV and 0.5 mM CNQX/2 mM dAPV.

their onset in the infragranular layers (Figure 5) as described previously (Sakata and Harris, 2009). Similarly to sensory-evoked responses the spontaneous activity was gradually suppressed by epipial CNQX/dAPV in all layers of the cortical barrel column following a surface-to-depth gradient with faster and stronger inhibition at the superficial layers (Figures 1B, 5, 6). In L2/3 and L4, where complete SEP inhibition by CNQX/dAPV was observed, LFP deflections during the up-states were first inhibited and then switched polarity from the negative to positive direction, and the up-state related sinks in these layers were converted to sources (Figures 5, 6A). Also, MUA density during up-states decreased in L2/3 from 41.8 ± 10.0 to 0.3 ± 0.2 spikes/s after drug application (p = 0.001) and in L4 from 106 ± 17 to 4 ± 2 spikes/s (p = 0.001) (n = 11; Figures 6B,E). During 1 h of CNQX/dAPV application, up-states were reduced yet persisted, maintaining negative LFP polarity, sinks, and MUA activation in the infragranular layers (Figure 5). In L6, the up-states frequency decreased from 1.1 ± 0.1 s⁻¹ in control conditions to 0.2 ± 0.1 s⁻¹ (n = 11; p = 0.001; Figures 6C,F), the up-states amplitude reduced from 376 ± 63 to 236 ± 55 µV, and the up-states duration reduced from 238 ± 19 to 133 ± 12 ms (n = 11; p = 0.001; Figures 6A,D). The MUA frequency during up-states in L5/6 was suppressed by CNQX/dAPV from 106 ± 13 to 15 ± 3 spikes/s (n = 11;

p = 0.001; Figures 6B,E). These findings are consistent with the importance of intracortical excitatory connections in the generation of slow wave activity and the initiation of up-states in the infragranular layers (Silva et al., 1991; Sanchez-Vives and McCormick, 2000; Timofeev et al., 2000; Sakata and Harris, 2009).

DISCUSSION

In the present study, we assessed the time course and the magnitude of inhibitory effects produced by epipially applied ionotropic glutamate receptor antagonists on sensory evoked and spontaneous network activity through the depth of the rat barrel cortex *in vivo*. The described pharmacodynamics of these drugs may also allow the prediction of diffusion dynamics of other drugs with similar diffusion coefficients.

Development of the Inhibitory Effects of
CNQX/dAPV Through the Cortical Depth

In our experiments, the concentration of CNQX/dAPV (170–500 µM CNQX and 0.7–2 mM dAPV) that fully blocked SEP in L2/3 and L4 was more than 10 times the glutamate receptor antagonist concentration generally used for bath application
A

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

B

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

C

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

D

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

E

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

**FIGURE 3** | Sensory evoked activity resistant to ionotropic glutamate receptor antagonists in the thalamorecipient cortical layers. (A) Top, Stimulus-triggered LFP averages (n = 30 responses; black traces) of sensory responses evoked by principal whisker deflection across cortical depths of the corresponding barrel column overlaid on color-coded current source density plot (CSD) in control, 40 min after 500 µM CNQX and 2 mM dAPV application, and 10 min after consecutive 2% lidocaine application. (B) Corresponding stimulus-triggered averages for MUA (n = 30 responses). Stimulus onset on (A,B) is indicated by the black arrowhead and vertical magenta line. The cortical layer borders are shown left of the CSD and MUA density plots, and recording electrode depths are shown on the right. (C) The time course of early evoked L4 MUA suppression by 2% lidocaine epipially applied (gray bar above the plot) after the full blockade of SEP by CNQX/dAPV. The shaded area around the curve show SE bands (n = 4 animals). The white circle on the curve indicates a half-maximal effect time (T1/2). (D) Averaged L4 LFP and MUA in control, in the presence of CNQX and dAPV, and after lidocaine application. Before averaging, individual LFP and MUA curves were aligned by the onset of SEP. Black arrowhead with an error bar above the dashed magenta line shows a mean time with SE between the stimulus and SEP onset. (E) Normalized mean values of SEP amplitude and early MUA peak in cortical L4 in control, in the presence of CNQX and dAPV, and after lidocaine application. (D,E) Show averaged data from 11 animals and two concentrations, 170 µM CNQX/700 µM dAPV and 0.5 mM CNQX/2 mM dAPV.
<figure><img src="image_12.png" /></figure>

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

**FIGURE 4** | The CNQX and dAPV effect on sensory evoked fast oscillations (FO) in cortical layer 4. (A) Three representative traces (after high-pass filtering, 200–1000 Hz) show fast L4 LFP and MUA oscillations evoked by sensory stimulation. Red bars overlaid on the LFP traces are detected spikes (MUA). Stimulus onset is indicated by the arrowhead. Corresponding PSTH of evoked MUA averaged over 30 responses is shown under the LFP traces. (B) Top, averaged power spectrum of evoked L4 MUA before and during CNQX/dAPV application (application time is shown with a white bar). Bottom, the time-course of sensory-evoked (early + late) MUA suppression by CNQX/dAPV. Shaded area shows SE bands. Data were averaged over 11 animals and two concentrations, 170 µM CNQX/700 µM dAPV and 0.5 mM CNQX/2 mM dAPV.

to suppress thalamocortical excitatory synaptic responses in slices *in vitro* (10 µM CNQX and 50 µM dAPV) (Feldman et al., 1998; Laurent et al., 2002), in intact preparations of superfused hippocampus of P15-25 rats *in vivo* (50 µM CNQX) (K hazipov and Holmes, 2003) and in superfusion of the neonatal rat barrel cortex *in vivo* (20 µM CNQX and 80 µM dAPV) (Minlebaev et al., 2007, 2009). The time course of inhibition produced by epipially applied glutamate receptor antagonists was comparable with bath application *in vitro* only in superficial L2/3, in deeper layers the inhibitory effects developed much slower following a surface-to-depth gradient. This likely reflects passive drug diffusion across the cortical depth counterbalanced by their removal from the extracellular space through blood microcirculation. Both of these factors are age-dependent with larger extracellular space volume (Bondareff and Pysh, 1968; Bondareff and Narotzky, 1972; Lehmenkuhler et al., 1993; Sykova et al., 1998, Sykova et al., 2002; Kilb et al., 2006) and less developed vascularization (Keep and Jones, 1990; Yu et al., 1994; Risser et al., 2009) in the immature animals. Together with a thinner cortex in the immature animals, this likely explains the much higher effective antagonist concentrations of the antagonists during their epipial application in older animals compared to neonates.

Various parameters of cortical activity displayed different sensitivity to the glutamate receptor antagonists. At a given cortical depth, spontaneous activity and the late component of sensory-evoked responses were inhibited earlier and stronger than the early component of sensory-evoked responses including SEP slope and amplitude, and early evoked MUA response. This difference likely involves the monosynaptic nature of SEP, an initial part of which is generated by the activation of

thalamocortical synapses whereas the late part of the sensory-evoked response and spontaneous activity are polysynaptic events supported by the intracortical connections and therefore are more sensitive to the attenuation of excitatory synapses (Johnston and Brown, 1981). Interestingly, the early evoked MUA response in L4 and the L5/6 border was the most resistant to the glutamate receptor antagonists which is likely due to a contribution of the action potentials in the thalamocortical axon terminals to MUA in the thalamorecipient layers (see below). Also, we observed that sensory-evoked responses persisted in the infragranular layers despite complete suppression of sensory responses in the supragranular layers (e.g., see **Figure 7A**), which is consistent with the direct thalamic input to the infragranular layers (Bureau et al., 2006; Brecht, 2007; Meyer et al., 2010; Oberlaender et al., 2012; Constantinople and Bruno, 2013; Crocker-Buque et al., 2015).

## Pharmacokinetics of the Glutamate Receptor Antagonists After Epipial Application

We further attempted to determine the cross-layer pharmacokinetics of glutamate receptor antagonists after epipial application. With this aim, we assessed SEP amplitude as the most robust parameter of glutamatergic synaptic transmission (Figure 7A). The level of SEP inhibition in different layers after epipial application of the glutamate receptor antagonists was further used for an estimation of the instant antagonist concentration in this layer at a given time point. To perform this effect-to-concentration conversion we based our calculations on the concentration dependence of the inhibitory effects of CNQX on evoked fast EPSCs with a Kd of 0.4 µM and Hill coefficient
<figure><img src="image_14.png" /></figure>

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

**FIGURE 5** | The CNQX and dAPV action on spontaneous activity through the depth of the barrel column. (A) Example traces of 16-channel recording of barrel cortex activity in control (top), 15 min after epipal application of 500 µM CNQX and 2 mM dAPV (middle), and 60 min after drug application (bottom). Red bars indicate detected spikes (MUA). (B) Corresponding current source density (CSD, upper panels) and MUA density (lower panels) plots of averaged up-states recorded in control and after drug application (150 events were averaged for each condition). LFP traces (black) are overlaid on color-coded CSD plots. The cortical layer borders on (A,B) are shown left of the CSD and MUA density plots, and recording electrode depths are shown on the right.
<figure><img src="image_16.png" /></figure>

of 1 (Honore et al., 1988; Llano et al., 1991) according to the formula:

$$A = Kd^h/(c^h + Kd^h),$$

where A is a normalized value of the determined parameter (SEP amplitude), h – Hill coefficient and c – drug concentration.

The results of conversion of the inhibitory effects of CNQX on SEP amplitude to the actual concentration of the drug at different cortical layers are shown on **Figure 7B**. According to this analysis, CNQX concentration progressively increased most rapidly at the superficial layers attaining saturation in L2/3 and L4. The CNQX concentration likely further continued increasing in these layers during the drug application but these values above saturating effects are hard to estimate in our assay. In deep layers, CNQX concentration increased much slower and attained only sub-saturating levels of a few µM that is two orders lower than the drug concentration in the epipial solution. We also estimated the time points after the onset of antagonist application at which CNQX concentration attained the Kd value of 0.4 and 10 µM

corresponding to the inhibition of the SEP amplitude to 50 and 3.8% of the control values, respectively (Figure 7C). These estimations also indicated that the effective concentration levels of the antagonist were slowly achieved through the cortical depth during epipial drug application, and the 10 µM concentration level was attained during 1 h of application at only about one half of the cortical depth. Given the remarkably similar diffusion parameters observed in different cortical areas (Vorisek et al., 2002; Slais et al., 2008) of the adult brain, as well as between different brain structures (Nicholson and Phillips, 1981; Svoboda and Sykova, 1991; Lehmenkuhler et al., 1993; Sykova et al., 1994, 1998; Simonova et al., 1996; Vorisek and Sykova, 1997b; Mazel et al., 1998, 2002; Reum et al., 2002) and species (Vargova et al., 2003; Sykova et al., 2005; Homola et al., 2006) the diffusion dynamics of glutamate antagonists and other molecules with similar diffusion properties through the rat somatosensory cortex described in the present study are likely to be relevant for the majority of brain regions. Our results could be of interest for the estimations of pharmacodynamics in other methods of drug delivery including injection into the
<figure><img src="image_17.png" /><figcaption>**FIGURE 7** | Pharmacokinetics of CNQX and dAPV after epipial application. (A) Development of the inhibitory effect of CNQX and dAPV on SEP amplitude across cortical layers over 1 h of application (averaged data from 11 animals). (B) Estimation of CNQX concentration dynamics in different cortical layers after epipial drug application. Instant values of CNQX concentration were calculated based on the concentration dependence of CNQX effect on the fast EP SCs with a Kd of 0.4 µM and Hill coefficient of 1 (Honore et al., 1988; Llano et al., 1991). (C) The time points where epipially applied CNQX (170–500 µM) attains Kd values of 0.4 and 10 µM at different depths of the cortical column.</figcaption></figure>

brain tissue or intraventricular application. It is worthy to note, that drug diffusion dynamics may differ between horizontal and vertical directions in the brain regions exhibiting anisotropic tissue properties (Vorisek and Sykova, 1997a; Mazel et al., 1998; Sykova et al., 1998, Sykova et al., 2002; Vorisek et al., 2002).

## Contribution of Thalamic Axonal Spikes to Cortical MUA

Multiple unit activity is assumed to represent action potentials from multiple neurons whose soma are located near the extracellular recording electrode (Henze et al., 2000; Bartho et al., 2004; Buzsaki et al., 2012). Here, we provide evidence that the action potentials, not only in cortical neurons, but also in thalamocortical axons contribute to MUA in cortical L4. Indeed, we found that spikes with the shortest latency during sensory-evoked responses persist after full blockade of glutamatergic transmission in L4 as evidenced by complete suppression of L4 SEPs. This is in keeping with the findings obtained in the primary somatosensory cortex of piglets using high-resolution magnetoencephalography, where sensory evoked responses consisted of the initial component which was localized in cortical L4 and which was insensitive to glutamate receptor blockade by kynurenic acid, and the subsequent kynurenate - sensitive component (Ikeda et al., 2002). The glutamate receptor antagonist insensitive MUA component which persisted in the sensory-evoked response in L4 has also been reported in the barrel cortex of neonatal rat pups (Minlebaev et al., 2007). This glutamate receptor antagonist resistant MUA was suppressed by voltage-gated sodium channel blockers, lidocaine in the present study and tetrodotoxin in a previous study in rat pups (Minlebaev et al., 2007). Several observations further suggest that this glutamate receptor antagonists resistant MUA represents spikes in thalamocortical axon terminals. Firstly, the mean latency between the stimulus and the peak of early evoked MUA in the presence of CNQX/dAPV was 7.9 ± 0.4 ms, while the latency of spikes recorded in whole-cell configuration from L4 neurons in response to whisker deflection is more than 10 ms (Constantinople and Bruno, 2013). Also, blockade of glutamatergic ionotropic receptors completely suppressed the sensory-evoked postsynaptic response and action potentials in L4 neurons during whole-cell recordings whereas the extracellular sensory-evoked MUA response persisted in the neonatal rat barrel cortex (Minlebaev et al., 2007). Secondly, in the present study, the short latency spikes were mainly restricted to the L4 and L5/6 border where the thalamocortical axon terminals extensively ramify and one could expect efficient summation of the transmembrane currents during action potentials in these compactly organized multiple axonal terminals and generation of large LFP signals (Chmielowska et al., 1989; Lu and Lin, 1993; Bruno and Sakmann, 2006; Bureau et al., 2006; Meyer et al., 2010). Finally, LFP signals generated in L4 by the activity of single thalamocortical neurons are characterized by the short-latency spikes representing afferent volleys followed by glutamatergic postsynaptic current, the latter being selectively suppressed by the intracortical infusion of the AMPA/kainate receptor blocker (Swadlow and Gusev, 2002; Swadlow et al., 2002; Hagen et al., 2017). Taken together, these findings indicate that spikes in thalamocortical axons significantly (~20%) contribute to the initial part of the sensory-evoked MUA response during extracellular intracortical recordings in L4 (but also likely at the L5/6 border) and this contribution should be considered during MUA analysis. Our findings also suggest that thalamic
activation evoked by rapid whisker deflection lasts for only ~5 ms and the rest of the MUA response is mainly generated by cortical neurons. This is different from the neonatal rat pups, where sensory-evoked MUA response after suppression of the glutamatergic ionotropic receptors lasts for several hundreds of milliseconds in keeping with the long-lasting activity bursts that are evoked in the thalamus by brief sensory stimuli (Minlebaev et al., 2011; Yang et al., 2013). Short duration of the CNQX/APV resistant, presumably thalamocortical axonal activity at the onset of sensory-evoked response in adult animals also suggests that thalamic input only generates the initial part of SEP whereas the late part of the response is supported by the local circuits, and that it only triggers but does not pace FOs that is in keeping with the pivotal roles of the intracortical connections, notably GABAergic synapses by the fast spiking interneurons, in the genesis of FOs (Jones et al., 2000).

## DATA AVAILABILITY

Original and processed data, and signal processing and analysis routines are available on request from the authors.

## REFERENCES

Andreasen, M., Lambert, J. D. C., and Jensen, M. S. (1989). Effects of new non-N-methyl-D-aspartate antagonists on synaptic transmission in the in vitro rat hippocampus. *J. Physiol.* 414, 317–336. doi: 10.1113/jphysiol.1989.sp017690

Barth, D. S. (2003). Submillisecond synchronization of fast electrical oscillations in neocortex. *J. Neurosci.* 23, 2502–2510. doi: 10.1523/JNEUROSCI.23-06-02502. 2003

Bartho, P., Hirase, H., Monconduit, L., Zugaro, M., Harris, K. D., and Buzsaki, G. (2004). Characterization of neocortical principal cells and interneurons by network interactions and extracellular features. *J. Neurophysiol.* 92, 600–608. doi: 10.1152/jn.01170.2003

Bondareff, W., and Narotzky, R. (1972). Age changes in the neuronal microenvironment. *Science* 176, 1135-1136. doi: 10.1126/science.176.4039. 1135

Bondareff, W., and Pysh, J. J. (1968). Distribution of the extracellular space during postnatal maturation of rat cerebral cortex. *Anat. Rec.* 160, 773–780. doi: 10. 1002/ar.1091600412

Brecht, M. (2007). Barrel cortex and whisker-mediated behaviors. *Curr. Opin. Neurobiol.* 17, 408–416. doi: 10.1016/j.conb.2007.07.008

Bruno, R. M., and Sakmann, B. (2006). Cortex is driven by weak but synchronously active thalamocortical synapses. *Science* 312, 1622–1627. doi: 10.1126/science. 1124593

Bureau, I., von Saint, P. F., and Svoboda, K. (2006). Interdigitated paralemniscal and lemniscal pathways in the mouse barrel cortex. *PLoS Biol.* 4:e382. doi: 10.1371/journal.pbio.0040382

Buzsaki, G. (2006). *Rhythms of the Brain*. Oxford: Oxford University Press. doi: 10.1093/acprof:oso/9780195301069.001.0001

Buzsaki, 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

Castro-Alamancos, M. A., and Oldford, E. (2002). Cortical sensory suppression during arousal is due to the activity-dependent depression of thalamocortical synapses. *J. Physiol.* 541, 319–331. doi: 10.1113/jphysiol.2002.016857

Chmielowska, J., Carvell, G. E., and Simons, D. J. (1989). Spatial organization of thalamocortical and corticothalamic projection systems in the rat Sml barrel cortex. *J. Comp. Neurol.* 285, 325–338. doi: 10.1002/cne.902850304

Constantinople, M. C., and Bruno, R. M. (2013). Deep cortical layers are activated directly by thalamus. *Science* 340, 1591–1594. doi: 10.1126/science.1236425

## AUTHOR CONTRIBUTIONS

RK conceived the project. DV, JL, NL, GB, and KC performed the experiments. AZ and GV analyzed the data. GV and RK wrote the paper.

## FUNDING

This work was supported by RSF # 17-15-01271 (data acquisition and analysis of the pharmacodynamics *in vivo*) and subsidy allocated to Kazan Federal University for the state assignment in the sphere of scientific activities # 6.5364.2017/9.10 (pharmacokinetics modeling *in silico*) and performed in the frame of the Program of Competitive Growth of Kazan University.

## ACKNOWLEDGMENTS

We thank Drs. A. Sirota and M. Minlebaev for their helpful and critical comments.

Conti, F., and Minelli, A. (1994). Glutamate immunoreactivity in rat cerebral cortex is reversibly abolished by 6-diazo-5-oxo-L-norleucine (DON), an inhibitor of phosphate-activated glutaminase. *J. Histochem. Cytochem.* 42, 717-726. doi: 10.1177/42.6.7910617

Crocker-Buque, A., Brown, S. M., Kind, P. C., Isaac, J. T., and Daw, M. I. (2015). Experience-dependent, layer-specific development of divergent thalamocortical connectivity. *Cereb. Cortex* 25, 2255-2266. doi: 10.1093/cercor/bhu031

Davies, J., Francis, A. A., Jones, A. W., and Watkins, J. C. (1981). 2-Amino-5-phosphonovalerate (2APV), a potent and selective antagonist of amino acid-induced and synaptic excitation. *Neurosci. Lett.* 21, 77-81. doi: 10.1016/0304-3940(81)90061-6

Di, S., Baumgartner, C., and Barth, D. S. (1990). Laminar analysis of extracellular field potentials in rat vibrissa/barrel cortex. *J. Neurophysiol.* 63, 832-840. doi: 10.1152/jn.1990.63.4.832

Douglas, R. J., and Martin, K. A. (2004). Neuronal circuits of the neocortex. *Annu. Rev. Neurosci.* 27, 419-451. doi: 10.1146/annurev.neuro.27.070203.144152

Feldman, D. E., Nicoll, R. A., Malenka, R. C., and Isaac, J. T. (1998). Long-term depression at thalamocortical synapses in developing rat somatosensory cortex. *Neuron* 21, 347-357. doi: 10.1016/S0896-6273(00)80544-9

Feldmeyer, D. (2012). Excitatory neuronal connectivity in the barrel cortex. *Front. Neuroanat.* 6:24. doi: 10.3389/fnana.2012.00024

Feldmeyer, D., Brecht, M., Helmchen, F., Petersen, C. C., Poulet, J. F., Staiger, J. F., et al. (2013). Barrel cortex function. *Prog. Neurobiol.* 103, 3-27. doi: 10.1016/j. pneurobio.2012.11.002

Groothuis, D. R., Vavra, M. W., Schlageter, K. E., Kang, E. W., Itskovich, A. C., and Hertzler, S. (2007). Efflux of drugs and solutes from brain: the interactive roles of diffusional transcapillary transport, bulk flow and capillary transporters. *J. Cereb. Blood Flow Metab.* 27, 43-56. doi: 10.1038/sj.jcbfm.9600315

Hablitz, J. J., and Sutor, B. (1990). Excitatory postsynaptic potentials in rat neocortical neurons in vitro. III. Effects of a quinoxalinedione non-NMDA receptor antagonist. *J. Neurophysiol.* 64, 1282-1290. doi: 10.1152/jn.1990.64.4. 1282

Hagen, E., Fossum, J. C., Pettersen, K. H., Alonso, J. M., Swadlow, H. A., and Einevoll, G. T. (2017). Focal local field potential signature of the single-axon monosynaptic thalamocortical connection. *J. Neurosci.* 37, 5123-5143. doi: 10.1523/JNEUROSCI.2715-16.2017

Henze, D. A., Borhegyi, Z., Csicsvari, J., Mamiya, A., Harris, K. D., and Buzsaki, G. (2000). Intracellular features predicted by extracellular recordings in the
hippocampus in vivo. *J. Neurophysiol.* 84, 390-400. doi: 10.1152/jn.2000.84.1.390

Homola, A., Zoremba, N., Slais, K., Kuhlen, R., and Sykova, E. (2006). Changes in diffusion parameters, energy-related metabolites and glutamate in the rat cortex after transient hypoxia/ischemia. *Neurosci. Lett.* 404, 137-142. doi: 10.1016/j.neulet.2006.05.028

Honore, T., Davies, S. N., Drejer, J., Fletcher, E. J., Jacobson, P., Lodge, D., et al. (1988). Quinoxalinediones: potent competitive non-NMDA glutamate receptor antagonists. *Science* 241, 701-703. doi: 10.1126/science.2899909

Ikeda, H., Leyba, L., Bartolo, A., Wang, Y., and Okada, Y. C. (2002). Synchronized spikes of thalamocortical axonal terminals and cortical neurons are detectable outside the pig brain with MEG. *J. Neurophysiol.* 87, 626-630. doi: 10.1152/jn.00332.2001

Johnston, D., and Brown, T. H. (1981). Giant synaptic potential hypothesis for epileptiform activity. *Science* 211, 294-297. doi: 10.1126/science.7444469

Jones, M. S., and Barth, D. S. (1999). Spatiotemporal organization of fast (>200 Hz) electrical oscillations in rat Vibrissa/Barrel cortex. *J. Neurophysiol.* 82, 1599-1609. doi: 10.1152/jn.1999.82.3.1599

Jones, M. S., MacDonald, K. D., Choi, B., Dudek, F. E., and Barth, D. S. (2000). Intracellular correlates of fast (>200 Hz) electrical oscillations in rat somatosensory cortex. *J. Neurophysiol.* 84, 1505-1518. doi: 10.1152/jn.2000.84.3.1505

Kandel, A., and Buzsaki, G. (1997). Cellular-synaptic generation of sleep spindles, spike-and-wave discharges, and evoked thalamocortical responses in the neocortex of the rat. *J. Neurosci.* 17, 6783-6797. doi: 10.1523/JNEUROSCI.17-06783.1997

Keep, R. F., and Jones, H. C. (1990). Cortical microvessels during brain development: a morphometric study in the rat. *Microvasc. Res.* 40, 412-426. doi: 10.1016/0026-2862(90)90036-Q

Khatipov, R., and Holmes, G. L. (2003). Synchronization of kainate-induced epileptic activity via GABAergic inhibition in the superfused rat hippocampus in vivo. *J. Neurosci.* 23, 5337-5341. doi: 10.1523/JNEUROSCI.23-12-05337. 2003

Khatipov, R., Zaynutdinova, D., Ogievetsky, E., Valeeva, G., Mitrukhina, O., Manent, J. B., et al. (2015). Atlas of the postnatal rat brain in stereotaxic coordinates. *Front. Neuroanat.* 9:161. doi: 10.3389/fnana.2015.00161

Kilb, W., Dierkes, P. W., Sykova, E., Vargova, L., and Luhmann, H. J. (2006). Hypoosmolar conditions reduce extracellular volume fraction and enhance epileptiform activity in the CA3 region of the immature rat hippocampus. *J. Neurosci.* Res. 84, 119-129. doi: 10.1002/jnr.20871

Kohling, R., Lehmenkukhler, A., Nicholson, C., and Speckmann, E. J. (1993). Superfusion of verapamil on the cerebral cortex does not suppress epileptic discharges due to restricted diffusion (rats, in vivo). *Brain Res.* 626, 149-155. doi: 10.1016/0006-8993(93)90574-7

Laurent, A., Goaillard, J. M., Cases, O., Lebrand, C., Gaspar, P., and Ropert, N. (2002). Activity-dependent presynaptic effect of serotonin 1B receptors on the somatosensory thalamocortical transmission in neonatal mice. *J. Neurosci.* 22, 886-900. doi: 10.1523/JNEUROSCI.22-03-00886.2002

Lehmenkuhler, A., Sykova, E., Svoboda, J., Zilles, K., and Nicholson, C. (1993). Extracellular space parameters in the rat neocortex and subcortical white matter during postnatal development determined by diffusion analysis. *Neuroscience* 55, 339-351. doi: 10.1016/0306-4522(93)90503-8

Llano, I., Marty, A., Armstrong, C. M., and Konnerth, A. (1991). Synaptic- and agonist-induced excitatory currents of Purkinje cells in rat cerebellar slices. *J. Physiol.* 434, 183-213. doi: 10.1113/jphysiol.1991.sp018465

Long, S. K., Smith, D. A., Siarey, R. J., and Evans, R. H. (1990). Effect of 6-cyano-2,3-dihydroxy-7-nitroquinoxaline (CNQX) on dorsal root-, NMDA-kainate- and quisqualate-mediated depolarization of rat motoneurones in vitro. *Br. J. Pharmacol.* 100, 850-854. doi: 10.1111/j.1476-5381.1990.tb14103.x

Lu, S. M., and Lin, R. C. (1993). Thalamic afferents of the rat barrel cortex: a light- and electron-microscopic study using Phaseolus vulgaris leucoagglutinin as an anterograde tracer. *Somatosens. Mot. Res.* 10, 1-16. doi: 10.3109/08990229309028819

Mazel, T., Richter, F., Vargova, L., and Sykova, E. (2002). Changes in extracellular space volume and geometry induced by cortical spreading depression in immature and adult rats. *Physiol. Res.* 51(Suppl. 1), S85-S93.

Mazel, T., Simonova, Z., and Sykova, E. (1998). Diffusion heterogeneity and anisotropy in rat hippocampus. *Neuroreport* 9, 1299-1304. doi: 10.1097/00001756-199805110-00008

Meyer, H. S., Wimmer, V. C., Oberlaender, M., de Kock, C. P., Sakmann, B., and Helmstaedter, M. (2010). Number and laminar distribution of neurons in a thalamocortical projection column of rat vibrissal cortex. *Cereb. Cortex* 20, 2277-2286. doi: 10.1093/cercor/bhq067

Minlebaev, M., Ben Ari, Y., and Khazipov, R. (2009). NMDA receptors pattern early activity in the developing barrel cortex in vivo. *Cereb. Cortex* 19, 688-696. doi: 10.1093/cercor/bhn115

Minlebaev, M., Ben-Ari, Y., and Khazipov, R. (2007). Network mechanisms of spindle-burst oscillations in the neonatal rat barrel cortex in vivo. *J. Neurophysiol.* 97, 692-700. doi: 10.1152/jn.00759.2006

Minlebaev, M., Colonnese, M., Tsintsadze, T., Sirota, A., and Khazipov, R. (2011). Early gamma oscillations synchronize developing thalamus and cortex. *Science* 334, 226-229. doi: 10.1126/science.1210574

Neuman, R. S., Ben-Ari, Y., Gho, M., and Cherubini, E. (1988). Blockade of excitatory synaptic transmission by 6-cyano-7- nitroquinoxaline-2,3-dione (CNQX) in the hippocampus in vitro. *Neurosci. Lett.* 92, 64-68. doi: 10.1016/0304-3940(88)90743-4

Nicholson, C., and Phillips, J. M. (1981). Ion diffusion modified by tortuosity and volume fraction in the extracellular microenvironment of the rat cerebellum. *J. Physiol.* 321, 225-257. doi: 10.1113/jphysiol.1981.sp013981

Oberlaender, M., de Kock, C. P., Bruno, R. M., Ramirez, A., Meyer, H. S., Dercksen, V. J., et al. (2012). Cell type-specific three-dimensional structure of thalamocortical circuits in a column of rat vibrissal cortex. *Cereb. Cortex* 22, 2375-2391. doi: 10.1093/cercor/bhr317

Ohata, K., and Marmarou, A. (1992). Clearance of brain edema and macromolecules through the cortical extracellular space. *J. Neurosurg.* 77, 387-396. doi: 10.3171/jns.1992.77.3.0387

Paxinos, G., and Watson, C. (2007). *The Rat Brain in Stereotaxic Coordinates*, 6th Edn. San Diego, CA: Academic Press.

Petersen, C. C. (2007). The functional organization of the barrel cortex. *Neuron* 56, 339-355. doi: 10.1016/j.neuron.2007.09.017

Pockberger, H., Rappelsberger, P., and Petsche, H. (1984). Penicillin-induced epileptic phenomena in the rabbit's neocortex I. The development of interictal spikes after epicortical application of penicillin. *Brain Res.* 309, 247-260. doi: 10.1016/0006-8993(84)90591-2

Reum, T., Olshausen, F., Mazel, T., Vorisek, I., Morgenstern, R., and Sykova, E. (2002). Diffusion parameters in the striatum of rats with 6-hydroxydopamine-induced lesions and with fetal mesencephalic grafts. *Neurosci. Res.* 70, 680-693. doi: 10.1002/jnr.10332

Reyes-Puerta, V., Sun, J. J., Kim, S., and Luhmann, H. J. (2015). Laminar and columnar structure of sensory-evoked multineuronal spike sequences in adult rat barrel cortex in vivo. *Cereb. Cortex* 25, 2001-2021. doi: 10.1093/cercor/bhu007

Risser, L., Plouraboue, F., Cloetens, P., and Fonta, C. (2009). 3D-investigation shows that angiogenesis in primate cerebral cortex mainly occurs at capillary level. *Int. J. Dev. Neurosci.* 27, 185-196. doi: 10.1016/j.ijdevneu.2008.10.006

Roy, N. C., Bessaih, T., and Contreras, D. (2011). Comprehensive mapping of whisker-evoked responses reveals broad, sharply tuned thalamocortical input to layer 4 of barrel cortex. *J. Neurophysiol.* 105, 2421-2437. doi: 10.1152/jn.00939. 2010

Sakata, S., and Harris, K. D. (2009). Laminar structure of spontaneous and sensory-evoked population activity in auditory cortex. *Neuron* 64, 404-418. doi: 10.1016/j.neuron.2009.09.020

Sanchez-Vives, M. V., and McCormick, D. A. (2000). Cellular and network mechanisms of rhythmic recurrent activity in neocortex. *Nat. Neurosci.* 3, 1027-1034. doi: 10.1038/79848

Silva, L. R., Amitai, Y., and Connors, B. W. (1991). Intrinsic oscillations of neocortex generated by layer five pyramidal neurons. *Science* 251, 432-435. doi: 10.1126/science.1824881

Simonova, Z., Svoboda, J., Orkand, P., Bernard, C. C. A., Lassmann, H., and Sykova, E. (1996). Changes of extracellular space volume and tortuosity in the spinal cord of Lewis rats with experimental autoimmune encephalomyelitis. *Physiol. Res.* 45, 11-22.
Slais, K., Vorisek, I., Zoremba, N., Homola, A., Dmytrenko, L., and Sykova, E. (2008). Brain metabolism and diffusion in the rat cerebral cortex during pilocarpine-induced status epilepticus. *Exp. Neurol.* 209, 145–154. doi: 10.1016/j.expneurol.2007.09.008

Steriade, M. (2001). Impact of network activities on neuronal properties in corticothalamic systems. *J. Neurophysiol.* 86, 1–39. doi: 10.1152/jn.2001.86.1.1

Svoboda, J., and Sykova, E. (1991). Extracellular space volume changes in the rat spinal cord produced by nerve stimulation and peripheral injury. *Brain Res.* 560, 216–224. doi: 10.1016/0006-8993(91)91235-S

Swadlow, H. A., and Gusev, A. G. (2002). Receptive-field construction in cortical inhibitory interneurons. *Nat. Neurosci.* 5, 403–404. doi: 10.1038/nn847

Swadlow, H. A., Gusev, A. G., and Bezdudnaya, T. (2002). Activation of a cortical column by a thalamocortical impulse. *J. Neurosci.* 22, 7766–7773. doi: 10.1523/JNEUROSCI.22-17-07766.2002

Sykova, E., Mazel, T., Hasenohrl, R. U., Harvey, A. R., Simonova, Z., Mulders, W., et al. (2002). Learning deficits in aged rats related to decrease in extracellular volume and loss of diffusion anisotropy in hippocampus. *Hippocampus* 12, 269–279. doi: 10.1002/hipo.1101

Sykova, E., Mazel, T., and Simonova, Z. (1998). Diffusion constraints and neuron-glia interaction during aging. *Exp. Gerontol.* 33, 837–851. doi: 10.1016/S0531-5565(98)00038-2

Sykova, E., and Nicholson, C. (2008). Diffusion in brain extracellular space. *Physiol. Rev.* 88, 1277–1340. doi: 10.1152/physrev.00027.2007

Sykova, E., Svoboda, J., Polak, J., and Vatal, A. (1994). Extracellular volume fraction and diffusion characteristics during progressive ischemia and terminal anoxia in the spinal cord of the rat. *Cereb. Blood Flow Metab.* 14, 301–311. doi: 10.1038/jcbfm.1994.37

Sykova, E., Vorisek, I., Mazel, T., Antonova, T., and Schachner, M. (2005). Reduced extracellular space in the brain of tenascin-R- and HNK-1-sulphotransferase deficient mice. *Eur. J. Neurosci.* 22, 1873–1880. doi: 10.1111/j.1460-9568.2005.04375.x

Timofeev, I., Grenier, F., Bazhenov, M., Sejnowski, T. J., and Steriade, M. (2000). Origin of slow cortical oscillations in deafferented cortical slabs. *Cereb. Cortex* 10, 1185–1199. doi: 10.1093/cercor/10.12.1185

Traynelis, S. F., Wollmuth, L. P., McBain, C. J., Menniti, F. S., Vance, K. M., Ogden, K. K., et al. (2010). Glutamate receptor ion channels: structure, regulation, and function. *Pharmacol. Rev.* 62, 405–496. doi: 10.1124/pr.109.002451

Vargova, L., Homola, A., Amecnik, J., Tichy, M., Benes, V., and Sykova, E. (2003). Diffusion parameters of the extracellular space in human gliomas. *Glia* 42, 47–88. doi: 10.1002/glia.10204

Vorisek, I., Hajek, M., Tintera, J., Nicolay, K., and Sykova, E. (2002). Water ADC, extracellular space volume, and tortuosity in the rat cortex after traumatic injury. *Magn. Res. Med.* 48, 994–1003. doi: 10.1002/mrm.10305

Vorisek, I., and Sykova, E. (1997a). Evolution of anisotropic diffusion in the developing rat corpus callosum. *J. Neurophysiol.* 78, 912–919.

Vorisek, I., and Sykova, E. (1997b). Ischemia-induced changes in the extracellular space diffusion parameters, K+, and pH in the developing rat cortex and corpus callosum. *J. Cereb. Blood Flow Metab.* 17, 191–203.

Yang, J. W., An, S., Sun, J. J., Reyes-Puerta, V., Kindler, J., Berger, T., et al. (2013). Thalamic network oscillations synchronize ontogenetic columns in the newborn rat barrel cortex. *Cereb. Cortex* 23, 1299–1316. doi: 10.1093/cercor/bhs103

Yu, B. P., Yu, C. C., and Robertson, R. T. (1994). Patterns of capillaries in developing cerebral and cerebellar cortices of rats. *Acta Anat.* 149, 128–133. doi: 10.1159/000147567

**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.

Copyright © 2018 Vinokurova, Zakharov, Lebedeva, Burkhanova, Chernova, Lotfullina, Khazipov and Valeeva. 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) and the copyright owner(s) 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.