Нелинейная
динамика и нейронаука
УДК 530.182,519.876.5
https://doi.org/10.18500/0869-6632-2020-28-l-90-97
Взаимодействие между ультрамедленными флуктуациями нейронных сетей префронтальной коры и колебаниями мозга
М. Рой1'2, H.A. Новиков1, Д. Г. Захаров1 и Б. С. Гуткин1,2
1 Национальный исследовательский университет «Высшая школа экономики» Россия, Москва, ул. Мясницкая 20 2 Group of Neural Theory, LNC INSERM U960, École Normale Supérieure PSL* University
29 rue d'Ulm, Paris, France E-mail: n.novikov@hse.ru, dgzakharov@hse.ru, boris.gutkin@ens.fr Автор для переписки Денис Геннадьевич Захаров, dgzakharov@hse.ru Поступила в редакцию 15.10.2019, принята к публикации 3.12.2019, опубликована 26.02.2020
Целью нашей работы было изучение влияния различных ритмов головного мозга (тета, бета, гамма ритмы в диапазоне частот от 5 до 80 Гц) на ультрамедленные колебания (с частотой 0.5 Гц и ниже), проявляющиеся в чередовании состояний с высокой и низкой активностью. Эти ультрамедленные колебания обычно наблюдаются при нервной деятельности в человеческом мозге и, в частности, в префронтальной коре во время отдыха. Считается, что они генерируются локальными кортикальными сетями при наличии импульсных входов и нейронного шума. Структура этих колебаний имеет специфическую статистику, а их характеристики связаны с когнитивными способностями, такими как, например, эффективность и емкость рабочей памяти. Методы. В нашем исследовании мы использовали ранее построенную математическую модель, описывающую активность кортикальной сети, состоящее из популяций пирамидных клеток и интернейронов. Эта модель была разработана для описания глобального входного воздействия на локальные сети префронтальной коры из других кортикальных областей или подкорковых структур. Динамика модели исследовалась численно. Результаты. Мы обнаружили, что увеличение частоты существенно увеличивает время пребывания в состоянии с высокой активностью и, следовательно, повышает устойчивость самоподдерживающейся колебательной активности в гамма-диапазоне. Обсуждение. Мы считаем, что такие эффекты были бы полезны для обработки и передачи информации в кортикальных сетях с иерархическим торможением.
Ключевые слова: нейронные сети, частотные модели, префронтальная кора, ультрамедленные колебания, ритмы мозга.
Образец цитирования: Рой М., Новиков H.A., Захаров Д.Г., Гуткин Б.С. Взаимодействие между ультрамедленными флуктуациями нейронных сетей префронтальной коры и колебаниями мозга//Известия вузов. ПНД. 2020. Т. 28, № 1. С. 90-97. https://doi.org/10.18500/0869-6632-2020-28-l-90-97
Статья опубликована на условиях лицензии Creative Commons Attribution License (CC-BY 4.0).
Финансовая поддержка. Работа выполнена при поддержке Российского научного фонда, грант № 17-11-01273.
90
© Рой М., Новиков H.A., Захаров Д.Г., Гуткин Б.С., 2020
https://doi.org/10.18500/0869-6632-2020-28-1-90-97
Interaction between PFC neural networks ultraslow fluctuations
and brain oscillations
M. Rooy1,2, N.A. Novikov1, D. G. Zakharov1 and B. S. Gutkin1,2
1National Research University Higher School of Economics Myasnitskaya St. 20, Moscow, Russia 2Group of Neural Theory, LNC INSERM U960, Ecole Normale Superieure PSL* University
29 rue d'Ulm, Paris, France E-mail: n.novikov@hse.ru, dgzakharov@hse.ru, boris.gutkin@ens.f Correspondence should be addressed to Denis Zakharov, dgzakharov@hse.ru Received 15.10.2019, accepted 3.12.2019, published 26.02.2020
Aim of the work was to study the influence of different brain rhythms (i.e. theta, beta, gamma ranges with frequencies from 5 to 80 Hz) on the ultraslow oscillations with frequency of 0.5 Hz and below, where high and low activity states alternate. Ultraslow oscillations are usually observed within neural activity in the human brain and in the prefrontal cortex in particular during rest. Ultraslow oscillations are considered to be generated by local cortical circuitry together with pulse-like inputs and neuronal noise. Structure of ultraslow oscillations shows specific statistics and their characteristics has been connected with cognitive abilities, such as working memory performance and capacity. Methods. In the study we used previously constructed computational model describing activity of a cortical circuit consisting of the populations of pyramidal cells and interneurons. This model was developed to mimic global input impinging on the local prefrontal cortex circuit from other cortical areas or subcortical structures. The model dynamics was studied numerically. Results. We found that frequency increase deferentially lengthens the up states and therefore increases stability of self-sustained activity with oscillations in the gamma band. Discussion. We argue that such effects would be beneficial to information processing and transfer in cortical networks with hierarchical inhibition.
Key words: neuronal networks, rate models, prefrontal cortex, ultraslow oscillations, brain oscillations.
Reference: Rooy M., Novikov N.A., Zakharov D.G. and Gutkin B.S. Interaction between PFC neural networks ultraslow fluctuations and brain oscillations. Izvestiya VUZ. Applied Nonlinear Dynamics, 2020, vol. 28, no. 1, pp. 90-97. https:// doi.org/10.18500/0869-6632-2020-28-1-90-97
This is an open access article distributed under the terms of Creative Commons Attribution License (CC-BY 4.0). Acknowledgements. This work was supported by Russian Science Foundation, grant no. 17-11-01273.
Introduction
In the prefrontal cortex a canonical local circuit consists of pyramidal (PYR) neurons coupled together with three types of inhibitory interneurons: parvalbumin (PV), somatostatin (SOM) and vasoactive intestinal polypeptide (VIP) neurons (see schema in Fig. B1). Experimental results showed in simultaneously recorded neural populations of the prefrontal cortex (PFC) superficial layers that activity transits synchronously between high and low activity states lasting several to tens of seconds [1]. This constitutes resting state of activity and has been shown to be important for working memory performance. Interestingly, in the context of low neural firing rates activities, crucial for cellular regeneration mechanisms because of reduced metabolic costs [2], the slower is the switching between high and low activity states, the higher is the rate of information transmission [3]. Therefore, ultraslow activity fluctuations are computationally interesting in the context of quiet wakefulness, notably for efficient detection of important signals in the environment, or for correct transmission of internal signals mediating the consolidation of memories. Furthermore, the mechanisms involved in the generation of ultraslow network activity fluctuations may have important implications for the mechanisms of the generation of persistent activity during working memory tasks and for gating of sensory information
during these tasks. Since specific oscillatory frequency bands modulate brain activity during working memory tasks [4], we tested effects of an oscillatory background on PFC network bistable dynamics in order to study the relative roles of different interneuronal types on persistent activity and gating.
1. Methods
In order to better understand mechanisms underlying the changes in endogenous activity under oscillations, we modeled the layer II/III PFC local circuit. PV interneurons, target the axosomatic region of pyramidal cells, exerting divisive effects on their output activities, while SOM interneurons, target the dendrites of pyramidal cells, exerting substractive type of inhibition by increasing the spiking threshold [5]. VIP interneurons in layers II/III exert inhibitory control of SOM and PV interneurons in the PFC, with a stronger effect on SOM interneurons [6]. Model variables described the firing rates dynamics of various neuronal populations (rpyr, rsom, rpv, and rvip) in a local PFC circuit [7], based on [8] theoretical work for substractive vs divisive inhibition of pyramidal activity by SOM and PV interneuron populations respectively. We further incorporated substractive inhibition between SOM, VIP and PV populations [7]. Full model is given by following set of differential equation:
drpyr _
Ts ~dT — - rpyr+ _ _
I k ( , > r _r )F ( m pyr-pyrrpyr_ ft) som-pyrrsom+1ext-pyr_1adapt+1oscilA + _
+ ke("Jpyr-pyr'pyr ' pyr) F e I 1+ f r l+osS(t),
\ 1 + pv-pyr' pv /
drpv
— rpv + (ki rpv )Fi( ®pyr-pv'pyr ®pv-pv rpv ®vip-pv'vip + lext-pv) + os'%(t) (1)
S dt dr.
som
dt dr у]
~dt
Ts — rsom + rsom)Fi( ®pyr-somrpyr ®vip-somrvip + ^ext-som) + ®s?(t))
Ts уу — rvip + (ki rvip )Fi ( ® pyr-viprpyr ® som-viprsom + ^ext-vip) + ^s^(t),
with /ext-pyr, /ext-pv, /ext-som, and /ext-vip external constant inputs to each neural population type. 1oscill is the sinusoidal component to the external input to PYR population, centered around 0, such that: /oscill — /max sin(2nf). /max is the oscillations amplitude, set at 10% of /ext-pyr, while f is the oscillation frequency. We set ts — 0.02 ms, close to each populations' type membrane time constant [9]. Fe(Fi) is a sigmoid response function characteristic of an excitatory (inhibitory) population, which gives a nonlinear relationship between input currents to a population, and its output firing rate.
Parameters ke(ki) modulate the amplitude of firing rate response to input current for pyramidals (inhibitory) neurons. Here ki are constants and Fe, Fi5 ke we taken as they were developed in [7]:
Fe'i(X, Qe'i)— !+ e(U,i) , (2)
e(x6e)
ke(x) —-r-r—. (3)
ew 1 + e(x0e) v y
The wx-x (with x — pyr, pv or vip) are the self excitatory (or self inhibitory) synaptic coupling of the excitatory (inhibitory) neural populations. The wx-y (with x — y, x — pyr, pv, som or vip and y — pyr, pv, som or vip) are the excitatory (or inhibitory) synaptic coupling from one population to another. We did not consider self-inhibition in SOM and VIP interneuron populations, since inhibitory chemical synapses between those neurons are rarely observed [9, 10], and did not study PV and SOM direct connections. The parameter os controls the strength of the random fluctuations of neural
populations' firing rate, and ^(t) is a white noise process. In simulations we used the Euler method and this terms were produced by appropriate sampling of a normal distribution. In order to model spike frequency adaptation from PYR neurons in the network model, we used the following equation, adapted from [11]:
dIadapt _ I . t (4)
Tadapt dt — 1 adapt + rpyrJadap^ (4)
where Tadapt and Jadapt are the adaptation time constant and the adaptation strength of PYR population. We chose xadapt _ 600 ms. We used an optimization algorithm, to select parameters (mx — x, mi — y, Jadapt and I0-x) reproducing network bistability of pyramidal neurons and interneurons in awake wild-type (WT) mice, with mean activities and durations close to the ones found experimentally. This algorithm created random associations of values between parameters and computed an error for activity properties between simulations and experiments. Details of the methodology are given in [11]. We selected the association of parameter values minimizing the error. The high activity/low activity state durations are computed as for the experimental results (see Methods in [1]).
2. Results
We developed a computational model of the local PFC circuit that reproduces the dynamics of resting activity seen in imaging experiments (see Figure 1) in both qualitative appearance of spiking and in statistics of two activity states. This provided a basis for further investigation of the influence of oscillatory inputs on transitions structure of the resting state. The model has following dynamical structure (see [11] for full analysis): there is a bistability between a stable «low activity» where both excitatory and inhibitory neuronal populations are nearly quiescent state and a «high activity state», where both populations produce significant level of activity. Both of these states are nodes. There is a third fixed point of the system that is a saddle. The stable manifold of the saddle form the separatrix between two stable states. For purposes of present study, the system is subjected to an additive noise that produces random transitions between two states. We have previously shown [11] that an appropriate choice of the noise strength can produce transition statistics (frequency and the life-times of two states) that reproduce the experimental data [1]. Figure 1 shows an exemplary experimental data obtained from calcium imaging (panels A1,2,3) using the model structured as shown in panel B1. Panel B3 shows an example of the pyramidal population firing rate dynamics as a function of time. Please note the transitions between two stable states of activity; these transitions are due to additive noise and the transition from the high-to-low states is also partially caused by the adaptation of the PYR firing rate (see equation 2 and the Iadapt terms in system (1)). In panel B2 we visualize a simulated modulation of spiking in a single pyramidal neuron that is produced by generating spike time using an Poisson process with a time-dependent (instantaneous) rate \(t) equal to firing rate of the pyramidal population rpyr (t) seen in the dynamics of the model (1) as a function of time. We emphasize that the dynamics are generated by the population model and the spike raster is a visualization of spiking of a potential neuron constituent of this population.
We stimulated our optimized model with sinusoid inputs of moderate amplitude and variable frequency. This stimulation was designed to mimic global input impinging on the local PFC circuit from other cortical areas or subcortical structures. We wanted to know how the statistics of the activity states are affected by the sinusoid input. Using computational modeling, we found that externally driven oscillations influence the duration of high and low activity states and that this influence depends on the frequency of the input. In figure 2, we can see that high frequencies (40 Hz, gamma oscillations) increase the durations of both high and low activity states (H-states and L-states, respectively), while low frequencies (5 Hz, theta oscillations) induce shorter H-states and L-states. Futhermore, we notice
bioRxiv preprint first posted online Apr. 13, 2018; doi: http://dx.doi.org/10.1101/301051. The copyright holder for this preprint (which was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. All rights reserved. No reuse allowed without permission.
A1
cholinergic input
A2
A3 E 80
CO JD
9-60
50 100 150 time (seconds)
200
> 40
9-20
H-states L-states
l/Li
1
50 100 150 time (seconds)
200
A4
100 80
о 60
40 20
I i
Ш
Q л i i i
PYR PYR PV PV SOM SOM VIP VIP HLHLHLHL
B2 30
^ 25
q3
| 20 J 15
8 io
5
1 и.. T. 1
г':;}''л
in. ju .. |i nj
■ »w
i
0
50 100 150 200 time (seconds)
B4
100
180 I 60
CT3
I 40
(/)
В 20
со
" 0
50 100 150 200 time (seconds)
PYR PYR PV PV SOM SOM VIP VIP HLHLHLHL
Fig. 1. Bistable firing rate dynamics of interconnected neural populations replicates ultraslow fluctuations recorded in the PFC of WT mice. (Al) Two-photon image of GCaMP6f expressing neurons, modified with permission from [1]. (Scale bar: 20 ¡im). (A2) Spike trains of a population of simultaneously recorded cells in a WT mouse, obtained through deconvolution of spontaneous Ca2+ transients. 80% of the recorded cells are PYR neurons. Modified with permission from [1]. (A3) Time varying population mean activity of the neurons shown in A2. The dashed red line delineates the threshold between high and low activity states (H-states and L-states, respectively). Red periods correspond to H-states and blue periods to L-states. Modified with permission from [1]. (See [1] for more info on the methods). (Bl) Schematic of the studied circuitry. (B2) Visualization of single pyramidal neurons activity from the population rate model. To produce this we use the mean population rate produced by the simulated model and creating spike times according to the Poisson process with the rate \(t) = mean population rate as a function of time. Note that since this rate fluctuates between the low and high activity states (see panel B3), we see the same in the single cell firing rastergram. (B3) Time varying mean population activity of pyramidal neurons, computed from the network model. We use the same method as [1] to delineate H-states (in red) and L-states (in blue)
o
<D
3-
c
0
"S
1_
"O
£ rtj a) E
£ to
ti
0_i_i_i_i_
5 10 15 20 25 30 35 40
oscillation frequency (Hz)
Fig. 2. Active state duration increases with frequency of input oscillations. State mean duration function of oscillation frequency (Hz) of external inputs to PYR population. The blue line shows the evolution of L-state mean durations (sec), while the red line shows H-state mean durations (sec). The shaded areas are ± sem
that as frequency increases, the relative effect on the duration of the up-states becomes stronger relative to the duration of the down-states. This means that as input frequency increases into the gamma range, the network becomes more active (mean firing rates increase). During working memory tasks, gamma oscillations correlates positively with memory maintenance, consistent with increased H-state durations found in our results. We further predict higher rates of information transmission associated to gamma bands, since our study associate them with slow switching between high and low activity states.
3. Conclusions
In this work we designed a minimal prefrontal cortical circuit model to analyze the dynamics of firing rates of the multiple neuronal population circuit. We developed a dynamical parameter optimization procedure for the circuit model. The optimization allowed us to identify model types that reproduced the experimentally observed data. We previously conducted analysis of the resulting model dynamics and found that the model was bistable with a slow alternation between high and low activity states, just as in the data. We made a stochastic version of the model that with the optimized parameters reproduced the statistics of the experimentally observed upper and lower sate durations. We then focused on these statistics and designed oscillatory inputs to the model. Notably we projected the oscillations to the pyramidal neuron population in the model. We then performed analysis of the effects of the oscillation frequency on the state statistics. We found that frequency increase differentially lengthens the up states and therefore increases stability of self-sustained activity with oscillations in the gamma band. We argue that such effects would be beneficial to information processing and transfer in cortical networks with hierarchal inhibition.
References
1. Koukouli F., Rooy M., Changeux J.-P., and Maskos U. Nicotinic receptors in mouse prefrontal cortex modulate ultraslow fluctuations related to conscious processing. PNAS, 2016, vol. 113, no. 51, pp. 14823-14828.
2. Vyazovskiy V.V. and Harris K.D. Sleep and the single neuron: The role of global slow oscillations in individual cell rest. Nat Rev Neurosci, 2013, vol. 14, pp. 443-451.
3. Droste F. and Lindner B. Up-down-like background spiking can enhance neural information transmission. eNeuro, 2017, vol. 4. ENEUR0.0282-17.2017.
4. Fell J., Axmacher N. The role of phase synchronization in memory processes. Nat Rev Neurosci, 2011, vol. 12, no. 2, pp. 105-118.
5. Jadi M., Polsky A., Schiller J., Mel B.W. Location-dependent effects of inhibition on local spiking in pyramidal neuron dendrites. PLoS Comput Biol, 2012, vol. 8, no. 6, el002550.
6. Pi H.-J. et al. Cortical interneurons that specialize in disinhibitory control. Nature, 2013, vol. 503, pp. 521-52.
7. Papasavvas C.A., Wang Y., Trevelyan A.J., Kaiser M. Gain control through divisive inhibition prevents abrupt transition to chaos in a neural mass model. Phys Rev E Stat Nonlin Soft Matter Phys., 2015, vol. 9, no. 3, p. 032723.
8. Chance F.S., Abbott L.F. Divisive inhibition in recurrent networks. Network, 2000, vol. 11, no. 2, pp. 119-129.
9. Beierlein M., Gibson J.R., Connors B.W. Two Dynamically Distinct Inhibitory Networks in Layer 4 of the Neocortex. J NeumphysioL, 2003, vol. 90, pp. 2987-3000.
10. Pfeffer C.K., Xue M., He V., Huang Z.J., Scanziani M. Inhibition of inhibition in visual cortex: The logic of connections between molecularly distinct interneurons. Nature Neuroscience, 2013, vol. 16, pp. 1068-1076.
11. Rooy M., Koukouli F., Maskos U. and Gutkin B. Nicotinic modulation of hierarchal inhibitory control over prefrontal cortex resting state dynamics: Modeling of genetic modification and schizophrenia-related pathology, 2018. bioRxiv. 301051.
Мари Рой - родилась (1991) в Ля Гарен Коломб, Франция. Училась на факультете физики университета Пьера и Мари Кюри (Université Pierre et Marie Curie, Париж, Франция). Закончила магистратуру по когнитивным наукам в Высшей нормальной школе (École Normale Supérieure, Париж, Франция). Там же защитила диссертацию (PhD) по теоретической нейронауке (2019).
29 rue d'Ulm, Paris, France
Group of Neural Theory, LNC INSERM U960, École Normale Supérieure PSL* University
Новиков Никита Александрович - родился в Москве (1986). Окончил Российский технологический университет МИРЭА (2009). Окончил аспирантуру НИУ ВШЭ по специальности «Психофизиология» (2018), на данный момент готовится к защите кандидатской диссертации. Работал инженером-программистом в НИИ «Аргон» (2007-2013), младшим научным сотрудником Вычислительного центра им. Дородницына РАН (2012-2013), младшим научным сотрудником в группе когнитивной психофизиологии НИУ ВШЭ (2013-2015). В настоящее время является младшим научным сотрудником Центра нейроэкономики и когнитивных исследований, ведет научную деятельность под руководством Бориса Гутки-на. Научные интересы - изучение влияния нейрональных осцилляций на вычислительные модели рабочей памяти. Автор 10 публикаций в рецензируемых научных изданиях.
Россия, 101000 Москва, ул. Мясницкая, 20
Центр нейроэкономики и когнитивных исследований, Национальный исследовательский университет «Высшая школа экономики» E-mail: n.novikov@hse.ru
Захаров Денис Геннадьевич - родился в Нижнем Новгороде (1976), окончил радиофизический факультет Нижегородского государственного университета им. Н.И. Лобачевского (1999). Защитил диссертацию на соискание ученой степени кандидата физико-математических наук по направлению радиофизика (2005). Более 20 лет работал в Институте прикладной физики РАН. В настоящее время является старшим научным сотрудником Центра нейроэкономики и когнитивных исследований Национального исследовательского университета «Высшая школа экономики». Область научных интересов - теория нелинейных колебаний и волн, нелинейная динамика, моделирование активности нейронных сетей. Имеет более 50 научных публикаций.
Россия, 101000 Москва, ул. Мясницкая, 20
Центр нейроэкономики и когнитивных исследований, Национальный исследовательский университет «Высшая школа экономики» E-mail: dgzakharov@hse.ru
ffi
Гуткин Борис Самуилович - родился в Санкт Петербурге (1966). В 1985 году закончил Университет штата Северная Каролина в Роли (North Carolina State University at Raleigh, США). Учился в аспирантуре факультета математики и центра когнитивных нейрона-ук Питтсбургского Университета (University of Pittsburgh, Пенсильвания, США). Защитил PhD по нейронауке (1999). Автор более 100 статей по вычислительному и математическому моделированию нейробиологических процессов. Главный редактор журнала In Silico Pharmacology; член редколлегии журнала PLOS Computational Biology.
29 rue d'Ulm, Paris, France
Group of Neural Theory, LNC INSERM U960, École Normale Supérieure PSL* University Россия, 101000 Москва, ул. Мясницкая, 20
Центр нейроэкономики и когнитивных исследований, Национальный исследовательский университет «Высшая школа экономики» E-mail: boris.gutkin@ens.fr