- 1Melbourne Brain Centre, The University of Melbourne, Parkville, VIC, Australia
- 2Department of Biomedical Engineering, The University of Melbourne, Parkville, VIC, Australia
- 3Department of Medicine, The Royal Melbourne Hospital, The University of Melbourne, Parkville, VIC, Australia
Gamma and theta oscillations have been functionally associated with cognitive processes, such as learning and memory. Synaptic conductances play an important role in the generation of intrinsic network rhythmicity, but few studies have examined the effects of voltage-gated ion channels (VGICs) on these rhythms. In this report, we have used a pyramidal-interneuron-gamma (PING) network consisting of excitatory pyramidal cells and two types of inhibitory interneurons. We have constructed a conductance-based neural network incorporating a persistent sodium current (INaP), a delayed rectifier potassium current (IKDR), a inactivating potassium current (IA) and a hyperpolarization-activated current (IH). We have investigated the effects of several conductances on network theta and gamma frequency oscillations. Variation of all conductances of interest changed network rhythmicity. Theta power was altered by all conductances tested. Gamma rhythmogenesis was dependent on IA and IH. The IKDR currents in excitatory pyramidal cells as well as both types of inhibitory interneurons were essential for theta rhythmogenesis and altered gamma rhythm properties. Increasing INaP suppressed both gamma and theta rhythms. Addition of noise did not alter these patterns. Our findings suggest that VGICs strongly affect brain network rhythms. Further investigations in vivo will be of great interest, including potential effects on neural function and cognition.
1. Introduction
Synchronous oscillations in the brain, including hippocampal networks, have attracted much attention in recent years. These oscillations can be investigated in vivo and in vitro, and more recently can be modeled in silico due to the wide availability of advanced desktop processors and simulation software. These rhythmic activities are believed to represent temporal links between compartments of larger neural ensembles (Colgin and Moser, 2010). Increasing evidence suggests that brain oscillations have a key function in mechanisms of sensory-cognitive processes and that some forms of cognitive neuropathology may represent “rhythmopathies,” i.e., disorders of normal rhythmicity (Mably and Colgin, 2018). Abnormal neural oscillations have been reported in schizophrenia (Uhlhaas and Singer, 2015), cognitive and motor phenotypes of Parkinson's disease (Oswal et al., 2013), as well as mild cognitive impairment and Alzheimer's disease (Başar, 2013; Başar-Eroğlu et al., 2013; Vecchio et al., 2013; Yener and Başar, 2013). Thus, understanding the basis of these oscillations may improve the functional deficits of these diseases as well as potential therapeutic targets.
The hippocampus plays a significant role in learning and memory encoding and retrieval (Nadel and Moscovitch, 1997; Douchamps et al., 2013). These functions are related to hippocampal oscillatory activities (Xu et al., 2004). Hippocampal gamma (40–80 Hz) and theta (3–8 Hz) rhythms have been specifically implicated in cognitive functions (Rutishauser et al., 2010; Wang, 2010; Lisman and Jensen, 2013). Theta rhythm is also believed to be significant in exploratory movement in rats (Buzsáki, 2002). Synaptic currents and network connectivity have significant roles in network rhythms (Wang, 2010). In addition, voltage-gated conductances play an important role in network rhythmogenesis which has been less extensively explored. Voltage-gated conductances, such as the hyperpolarization activated non-specific cation current (IH) and the persistent sodium current (INaP), have been shown to affect theta (Xu et al., 2004) and gamma rhythms (Wang, 1993), respectively. This implies a potent role of IH and INaP in rhythmogenesis and potential modulatory targets. To further explore the effects of voltage-gated ion channels (VGICs) on rhythmicity, we constructed a pyramidal-interneuron-gamma (PING) network model containing several voltage-gated conductances and three distinct populations of neurons: excitatory pyramidal cells (E-cells), oriens lacunosum moleculare interneurons (O-cells), and inhibitory interneurons (I-cells). In this form of gamma, individual E-cells fire at or near gamma frequency, and their active participation is crucial; that is, the E-cells drive and synchronize the I-cells, and the I-cells gate and synchronize the E-cells (Traub et al., 1997; Börgers et al., 2005). In this network, there is a very strong connectivity between excitatory and inhibitory sub-populations, with such short rise times that an increase in E-cell spiking quickly evokes a surge in spiking of the I-cells (Rich et al., 2017). The resulting pulse of inhibitory input to the E-cells leads them toward synchrony. When inhibition wears off, the E-cells resume spiking, causing the cycle to repeat. Such strong connectivities can then create a 1:1 bursting ratio between inhibitory and excitatory cells that is a hallmark of classic PING theory (Kopell et al., 2010). A similar ratio in excitatory and inhibitory firing rates is supported by in vitro experimental evidence (Whittington et al., 2000). In such a model, the external drives play a critical role in generating the PING activity as well (Börgers and Kopell, 2005). A very weak drive current to the E-cells, similar to a very strong drive to I-cells, can abolish the PING rhythms. In these scenarios, I-cells can synchronize at frequencies too high to be entrained by the E-cells, or act asynchronously, which inhibits E-cell activity and hence, the synchronous rhythms vanish. Nevertheless, there is a setting in this parameter space where both suppression of the E-cells by asynchronous activity of the I-cells and PING may occur. In a setting where E-cells and I-cells produce oscillations of roughly the same frequency, in the absence of O-cells, the system is pushed toward instability in the presence of noise (see Börgers and Kopell, 2005). The O-cell spikes affect both E-cell and I-cell populations. O-cells are a type of inhibitory neurons in the hippocampus that gate information flow, while firing phase-locked to and taking part in generating theta rhythms (Sekulić and Skinner, 2017). Their firing activity during the trough of theta rhythms is recorded in the local field potentials in the pyramidal layer of the CA1 region of the hippocampus (Klausberger et al., 2003; Varga et al., 2012). These cells target distal dendrites of pyramidal cells in CA1 and are known to express hyperpolarization-activated inward channels (Maccaferri and McBain, 1996) that enables them to express post-inhibitory rebound spiking and contribute to in vivo theta rhythms (Sekulić and Skinner, 2017).
By utilizing this model in the appropriate parameter regime, we report effects on theta and gamma oscillations in response to changes in a range of VGIC conductance amplitudes. Finally, we examine whether the nature and degree of alterations in VGIC conductances affect rhythmogenesis of theta and gamma oscillations.
2. Methods
2.1. Noiseless Simulations
A model developed by Kopell et al. (2010) (see also Gloveli et al., 2005; Tort et al., 2007) was implemented in NEURON. This model has been used in recent studies (Scheffer-Teixeira and Tort, 2016; Keeley et al., 2017; Jansen et al., 2021), and in this work was augmented with a persistent sodium conductance (gNaP) added to the excitatory pyramidal cells (E-cells). We then constructed a network consisting of 40 E-cells, five oriens lacunosum moleculare interneurons (O-cells), and five inhibitory interneurons (I-cells) within CA1 hippocampus. A persistent sodium conductance (gNaP) was incorporated in the E-cells, and was utilized along with a delayed rectifier potassium conductance in each cell type (gKDRe for E-cells, gKDRi for I-cells, and gKDRo for O-cells), a inactivating potassium conductance (gA) and a hyperpolarization-activated current conductance (gH) in the O-cells. Initial simulations were performed without added noise, with baseline values of gA= 16, gNa = 100, gH = 12, gKDRe = 80, gKDRi = 9, and gKDRo = 23 (). The gNaP was set to 2e−5 (Hsu et al., 2018). The network was driven by currents applied somatically to the E, I, and O-cell populations (Kopell et al., 2010) with values of IE = 0.8, II = 0.8, and IO = 3 in . The details of the implemented model and the parameter settings can be found in the Supplementary Material of this article. Local field potentials (LFPs) were measured by using an E cell with no drive currents as a probe (see Kopell et al., 2010). Power spectra calculations on LFPs were performed by high and low pass filtering at 0.5 and 100 Hz, respectively. The spectrum.periodogram function from the Signal Processing Toolbox in MATLAB (Chambers et al., 2012) was utilized to measure the power spectrum. The activity of the network with the baseline settings is shown in Figure 1 and serves as control data.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/630271/fncom-15-630271-HTML/image_m/fncom-15-630271-g001.jpg)
Figure 1. The activity of the network in baseline conditions. (A) The raster plot indicates the spike activity of individual neurons from different populations in the last 2 s of the simulation. (B) (Top) For better readability, the LFP signal is shown only for the duration of 50–500 ms of the simulation after the system reaches steady-state. This period is long enough to reach the steady state. (Bottom) The theta rhythm and the gamma nested within it are noticeable in the power spectrum for the signal. Note the different scales used for each frequency band.
Simulations of 10,000 ms were performed on an Acer Nitro 5 computer with a quadcore Intel Core i5-8300H CPU and 8GB RAM, running Windows 10. In each run, one conductance value was modified and all remaining conductance values were kept constant to observe shifts in spectra after the model reached steady-state once again. For gA, gH, gKDRe, gKDRi, and gKDRo trials, conductances were set to ratios of 0, 0.05, 0.5, 0.8, 1.5, 1.8, and 2 g, whereby g represents baseline conductance. In gNaP trials, larger incremental steps of conductance values were required to observe network changes. Thus, the set of 0, 0.05, 20, 35, 50, 65, and 80 g was employed for gNaP trials.
2.2. Noisy Simulations
Additional simulations were performed with noise incorporated. The applied current was set to IE = 0.8 +W, whereby W depicts a white noise process with an SEM (standard error of the mean) of 1.35e−3 , so the external drive has deterministic and stochastic components. In order to determine the effects of conductance change on power spectra, 11 simulations were performed at each conductance level. For each set, runs were averaged to report peak power and frequency. Confidence intervals were generated for both mean peak power and frequency. For statistical analysis we conducted a two-tailed Student's t-test with 10 degrees of freedom and a significance level of p < 0.05 (i.e., 95% confidence levels). The results are reported as mean ± standard error of the mean (SEM).
3. Results
3.1. Baseline Network Rhythmicity
Control simulations with no added noise are shown in Figure 1. The resulting network rhythmicity is consistent with previous studies of this PING network (Kopell et al., 2010), with clearly visible theta and gamma activity peaking at 3.70 and 59.89 Hz, respectively. Theta band power is more prominent than the gamma band (0.575 vs. 0.288 ), which is also within expectations for gamma nested in theta in vivo (Colgin and Moser, 2010; Jacobson et al., 2013; Butler et al., 2016). There is a lesser peak present in higher frequency activity (e.g., fast gamma) around 120 Hz. However, this peak is not representing a true oscillation but a harmonic of the fundamental gamma peak at 60 Hz.
Baseline simulations with noise are displayed in Figure 2. Figures 2A,B depict the baseline network activity. Figure 2C shows the mean power spectrum for all 11 trials. Only significant findings from the noisy experiments are reported. highlighting differences between noisy and noiseless trials. Theta power peaked with a value of 0.356 ± 0.111 (Figure 2B), lower than the noiseless peak (0.5752). The gamma peak showed no significant shift in power and remained lower than theta (0.285 ± 0.010 , p < 0.05). Mean gamma and theta peak frequencies were also consistent with noise-free values, at 59.910 ± 0.033 and 3.766 ± 0.198 Hz, respectively. Figure 2C indicates that the SEM was greater in theta compared to gamma. Thus, the theta band appeared more sensitive to changes in noisy conditions.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/630271/fncom-15-630271-HTML/image_m/fncom-15-630271-g002.jpg)
Figure 2. The activity of the network in baseline conditions after applying the white noise to the E-cell population. (A) The raster plot indicating the spike activity of individual neurons. (B) (Top) For better readability, the LFP signal is shown only for the duration of 50–500 ms of the simulation for a sample trial run. (Bottom) The theta rhythm and the gamma nested within it are noticeable in the power spectrum for the signal. Note the different scales used for each frequency band. (C) The mean power spectrum using 11 trial runs and the corresponding 95% confidence intervals around it using a Student's t-test.
3.2. Network Rhythmicity in Response to Variation of VGIC Conductance Amplitudes
The following sections report the results on LFP spectral properties resulting from adjusting the amplitude of each conductance separately. Figure 3 provides an overview of the variations in power spectrum in response to change of each of the 6 VGIC conductance values. In Figures 4–9, the effects of altering each conductance value is displayed separately. Figure 10 shows the shifts in theta and gamma peak frequencies and power amplitudes, in response to modification of VGIC conductances with noise in the system.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/630271/fncom-15-630271-HTML/image_m/fncom-15-630271-g003.jpg)
Figure 3. The behavior of the system with changing VGIC conductances, with closeups of the theta (3–12 Hz) and gamma (40–80 Hz) band frequency ranges. (A) Any change in the gA value compared to baseline greatly reduces the power of oscillations in the theta range. It also causes the gamma rhythms to lower and vanish eventually. (B) Any variation from the baseline setting in the gH reduces the power of both theta and gamma rhythms. (C,D) Deviations from baseline in the gKDR value in both groups of inhibitory interneurons reduces the theta band power. Gamma rhythms also decrease but less significantly. (E) Conductance shifts in the potassium delayed rectifier in the E-cells weaken the theta band significantly. This current also weakens the gamma rhythm. (F) Lowering the value of gNaP does not have any appreciable effect on the network rhythmicity. Increasing conductances, however, greatly lower the theta band power.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/630271/fncom-15-630271-HTML/image_m/fncom-15-630271-g004.jpg)
Figure 4. Impact of gA on network rhythmicity. (A–D) Plots showing the effects of selected gA values on the overall LFP signal (top) and power spectrum (bottom) of the network. (E,F) Graphical representation of how conductance values affect the (E) power and (F) frequency of peak rhythms in the theta (top) and gamma (bottom) bands. The red bar with the asterisk sign above indicates the results for the baseline value of the conductance under investigation. The results demonstrate that the A current has a very significant effect on both theta and gamma rhythm generation. It may be worth noting that decreasing this conductance seems to push the gamma rhythm toward higher frequencies before dropping off. Theta rhythms are weakened in all cases compared to the baseline in Figure 1.
3.2.1. Effects of VGIC Modulation on Theta Rhythmicity
Theta rhythmicity was found to be affected significantly by changes in all VGIC conductances from baseline values. Generally, a decrease in the power of theta was observed and in some cases was completely abolished. Peak frequency shifts were also often observed.
Effects of changing gA. Single trial runs modifying gA showed that this conductance strongly influenced theta band power (see Figure 4E). Theta rapidly declined in power with decreasing gA, reducing by 98.2% when gA was set to zero. The power was reduced to 4.3% of baseline value at 1.5 g, rose to 8.8% at 1.8 g, before being lost once again at 2 g with a reduction of 99.2%.
Theta peak frequency increased at 0.05, 1.5, and 2 g, rising by 21.7–24.7%, and remained stable at the remaining conductances (Figure 4E).
Effects of changing gH. Theta power decreased steeply when gH was reduced and was almost completely eliminated with a 98.3% reduction at 0.5 g. Further decreases in conductance caused a slight rebound of the peak's power to 16.4% of baseline. Raising gH also caused theta power reduction, with a 30.5% decrease at 1.5 g, then being eliminated (99.5 and 99.4% reduction, respectively) at 1.5 and 2 g, its lowest points for this conductance (see Figure 5E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/630271/fncom-15-630271-HTML/image_m/fncom-15-630271-g005.jpg)
Figure 5. Impact of gH on network rhythmicity. (A–D) Network LFP signal (top) and power spectra (bottom) at selected gH values. (E,F) Changes in (E) power and (F) frequency in response to gH for both theta (top) and gamma (bottom) peaks. The red bar with the asterisk sign indicates the results for the baseline value of gH. The results illustrate that the presence of H current is necessary for the generation of the gamma oscillations, while theta rhythms conversely seem to be lost at high conductance. Increased gH values will initially shift the gamma oscillations more toward fast gamma. Although further increases appear to shift the gamma peak toward higher values, it is vital to note that the gamma power is vanished at these values; Hence, the detected peaks are not of any significance. Lowering gH values also decreases the gamma power. Theta peaks are all greatly reduced compared to the power seen in Figure 1.
The frequency of the theta peak remained largely stable throughout variations. At 0.8 and 1.8 g, where the peak frequencies shifted higher, the corresponding power was nearly abolished and hence, no significant peak was detected in the theta range (see Figures 5E,F). In these cases, the highest power content was identified in the delta range (1–3 Hz).
Effects of changing gKDR in excitatory neurons (gKDRe). Figure 6E depicts that shifting gKDRe in either direction from baseline causes a drastic decrease in theta band power. gKDRe = 0.064 (0.8 g) and gKDRe = 0.04 (0.5 g) displayed the greatest theta power reductions (99.5 and 99.1%, respectively).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/630271/fncom-15-630271-HTML/image_m/fncom-15-630271-g006.jpg)
Figure 6. Impact of gKDRe on network rhythmicity. (A–D) Graphs showing the network's LFP signal (top) and power spectrum (bottom) at selected conductances. (E,F) Bar plot showing changes in theta (top) and gamma (bottom) peaks in terms of (E) power and (F) frequency at varying values of gKDRe. Baseline value of gKDRe and the corresponding results are represented with the red bars and the asterisk sign above them. Effects of varying gKDRe show that the delayed rectifier potassium current in the E-cells have a very significant effect on theta rhythm generation, with its power greatly weakened in most cases. Gamma also appears to be reduced at very small conductance values for this current, but when the current is nullified altogether, the power rises above that seen in Figure 1.
At 0.8, 1.8, and 2 g, the peak frequency of theta increased by 18.6–21.7% (see Figure 6F). A large shift in peak frequency (by 197.9%) was recorded at 0.5 g, but theta power appeared totally suppressed at this conductance, so this was likely a composite peak between theta and gamma range rather than a peak shift.
Effects of changing gKDR in interneurons (gKDRi). As seen in Figure 7E, theta power consistently declined with delayed rectifier conductance (gKDRi) reduction in interneurons. A large initial decrease of 65.8% was observed at 0.8 g, and theta was almost completely abolished as gKDRi approached 0, with a 98.5% reduction in power. Initial increases in this conductance caused an even sharper reduction (93.0%), but at 2 g the power recovered to 14.7% of the baseline value.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/630271/fncom-15-630271-HTML/image_m/fncom-15-630271-g007.jpg)
Figure 7. Impact of gKDRi on network rhythmicity. (A–D) Graphic representation of network LFP signal (top) and power spectrum (bottom) for several selected gKDRi values. (E,F) Fluctuations in (E) power and (F) frequency as seen in the theta (top) and gamma (bottom) bands in response to changing conductance. As mentioned before, the red bars with the asterisk sign denote the results at the baseline value of gKDRi. Similar to previous cases, the theta oscillations show a large decrease in power with any variation of gKDRi when compared to the baseline conditions in Figure 1, and theta appears to be almost completely lost when this current is removed. gKDRi has a more modest effect on gamma, but still notable; very small conductance values seem to reduce gamma power as well as push peak frequency higher.
The peak in theta frequency range shifted higher by 156.7% at 0.5 g but was associated with a large amplitude reduction of 97.7% (see Figure 7F). Other variations in gKDRi had smaller effects causing either no shift (0.8, 2 g) or between 10.3 and 21.6% shift toward higher frequencies (0, 0.05, and 1.5 g).
Effects of changing gKDR in oriens lacunosum interneurons (gKDRo). gKDRo modulation in either direction from the baseline value significantly reduced theta power (Figure 8E). Decreases in conductance lowered theta power, with a 98.2% reduction at gKDRo = 0.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/630271/fncom-15-630271-HTML/image_m/fncom-15-630271-g008.jpg)
Figure 8. Impact of gKDRo on network rhythmicity. (A–D) Changes in rhythmicity of model network shown through the LFP signal (top) and power spectra (bottom) at selected values of gKDRo. (E,F) Effects of modifying conductance on theta (top) and gamma (bottom) rhythms, in terms of (E) power and (F) frequency. Red bars with the asterisk sign are correspondent to the baseline values of gKDRo. The results suggest that changes in gKDRo causes large variations in the theta range, with the power being reduced greatly in all cases compared to Figure 1. Gamma seems to be reduced at low values of this conductance, to a more moderate degree. In addition, it appears that the gamma peak starts being pushed toward higher frequencies at large increases in this current's conductance.
Some variation in peak frequency was observed at 0, 0.05, and 1.5 g, but generally with reduced power (see Figure 8F).
Effects of changing gNaP. For gNaP, there was no significant change in network activity in the theta band when conductance was decreased. This conductance was very robust to variations of a similar scale to those tested with the other conductances, so larger incremental steps were used. These larger changes reduced theta power, by 77.7% at gNaP = 20 g, with even larger reductions at higher amplitudes (see Figures 9C,E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/630271/fncom-15-630271-HTML/image_m/fncom-15-630271-g009.jpg)
Figure 9. Impact of gNaP on network rhythmicity. (A–D) Network LFP signal (top) and power spectra (bottom) observed at selected gNaP values. (E,F) The values for (E) power and (F) frequency of peak theta (top) and gamma (bottom) network rhythms at varying gNaP. We have again used the red bars to highlight the results of the study with baseline gNaP values. It is illustrated here that the increase in the persistent sodium current causes a significant reduction in theta band power. Gamma oscillations are weakened as well, though less substantially. Decreases in gNaP do not appear to have any meaningful effect.
Theta peak frequency was not strongly affected by gNaP change (see Figure 9F) with a maximal increase of 25.7% observed at gNaP = 80 g.
3.2.2. Effects of VGIC Modulation on Gamma Rhythmicity
Gamma rhythmicity was also found to be significantly modulated by alteration of VGIC conductances, though it was generally less sensitive than the theta band. In some cases, larger changes in conductance were required before the rhythm's stability was affected or more gradual changes were observed. Most effects observed were in the form of power decreases, sometimes leading to elimination of the rhythm, but at least one case was observed where power increased. Both increases and decreases in peak frequency were observed.
Effects of changing gA. Gamma power tended to decrease with lowering gA. Notably, both gamma and higher frequency power were abolished when gA was set to 0 (see Figure 4E). Like the theta band, gamma power is decreased but fluctuated as gA increased (see Figure 4E).
Peak frequency shifted 19.8% higher at gA = 0.8 g, and higher still (by 80.5%) at 0.5 g. However, it diminished greatly below 0.5 g and resulted in a decrease as high as 64.9% at 0.05 g, as seen in Figure 4F.
Effects of changing gH. Gamma power consistently decreased with reduced gH, as seen in Figure 5E. The gamma power is almost abolished as gH approaches 0 with 99% reduction from baseline power (see Figure 5A). Conversely, at gH = 1.5 g, power increased by 38.1%, but the activity also shifted outside the typical gamma range to over 90 Hz (a 57.5% increase). Gamma power and frequency had opposing trends as gH increased further. At 2 g, the power dropped significantly to 97.8% less than baseline (Figure 5D), but the peak continued rising to higher frequencies, reaching over 130 Hz, as shown in Figures 5E,F. While the detected peak frequencies at 1.8 and 2 g are outside the desirable range for gamma, the peak power is nearly abolished at these points which determines the insignificance of these disappearing peaks.
Effects of changing gKDR in excitatory neurons (gKDRe). Gamma power displayed much less variation in gKDRe trials, decreasing maximally by 8.4% in the trials ranging from 0.5 to 2 g. However, at 0.05 g, gamma power dropped abruptly by 78.8% of baseline value. At gKDRe = 0, the power then increased to 14.2% above baseline value (see Figures 6A,E, respectively).
The gamma peak frequency was generally stable with a maximal 6.2% increase at gKDRe = 0 and shifts of <4% otherwise (see Figure 6F).
Effects of changing gKDR in interneurons (gKDRi). Like in the E-cells, gamma power was mostly stable in gKDRi trials as seen in Figure 7E. The peak power was reduced by 10.3% of the baseline value at 1.5 g, and less variation was observed in most other instances. This deviated at gKDRi = 0 and gKDRi = 0.05 g, where the power dropped by 50.5%, indicating a significant role of this conductance in gamma rhythmogenesis.
In addition, gamma peak frequency was increased 30.5% to nearly 80 Hz at gKDRi = 0 and gKDRi = 0.05 g (Figures 7A,F). The high frequency peak normally present at 120 Hz was abolished in these trials.
Effects of changing gKDR in oriens lacunosum interneurons (gKDRo). Gamma behaved similarly in gKDRo and gKDRi, though a steady downward trend in power was observed both when increasing and decreasing gKDRo. The power reduction did not exceed 20% of the baseline until 0.05 g, where the power dropped by 30.9%.
Gamma peak frequency only shifted at 1.8 and 2 g, where it was increased by 18.6 and 19.4%, respectively. 1.8 and 2 g were also the two conductance values where loss of the high-frequency peak at 120 Hz occurred (see Figure 3D).
Effects of changing gNaP. As with theta, reduced amplitude of gNaP compared to baseline had no appreciable effect on gamma rhythmicity. Increased gNaP gradually decreased gamma power, with the greatest effect at gNaP = 50 g (see Figure 9D), where the gamma peak's power abruptly dipped to 38.0% of baseline conditions.
Gamma peak frequency remained stable with variations in this conductance (Figure 9F).
3.3. Effects of Noise on Network Response to VGIC Modulation
In this section, we employ the simulated model with added noise to depict how VGIC conductance variations can affect this system.
3.3.1. Effects of Noise With VGIC Modulation on Theta Rhythmicity
Though many of the rhythmicity trends observed with noise were comparable to our noiseless trials, theta displayed significant sensitivity to this additional factor being introduced into the system. Theta power was more affected by the addition of noise than gamma and continued to be more affected by conductance modulation as well, consistent with the noiseless trials. The results in this section are obtained using 11 trials and reported as (mean ± SEM) in a Student's t-test.
Effects of changing gA. Theta power reduction was consistent with noiseless runs with respect to gA, with a key difference at 2 g. This gA value had the highest mean power for theta among gA trials, not counting the baseline case (0.043 ± 0.036 ); theta was functionally lost at this conductance without noise (see Figure 10A). Theta power still became negligibly small (0.015 ± 0.009 ) with gA = 0, reinforcing the significant role of this current in theta rhythmogenesis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/630271/fncom-15-630271-HTML/image_m/fncom-15-630271-g010.jpg)
Figure 10. Effects of changing conductance values in a network containing noise. Graphs display the baseline conductance (red bars) and four adjusted conductance values for each current, (A) gA, (B) gH, (C) gKDRe, (D) gKDRi, (E) gKDRo, and (F) gNaP. Fluctuations in peak frequency (top) and power (bottom) are given at each conductance value, using 11 trials in a system containing noise to calculate average values. Bars display the mean, whereas vertical error lines are used to indicate the 95% confidence interval around the means in a Student's t-test. Notably, theta power and frequency appear to be much more sensitive to the existence of noise in the network than the gamma band, showing much greater variance in general.
Theta peak frequency was shifted upwards for all variations of gA moving from 3.70 Hz at baseline to its highest point of 4.907 ± 0.932 Hz at 0.5 g.
Effects of changing gH. When examining gH modulation with noise (see Figure 10B), the most significant differences were in theta power. While the theta peak was almost completely nullified at 0.5 and 2 g without noise, with added noise the theta power followed an increasing trend in moving from 0 until 2 g. Nevertheless, compared to the baseline power observed with noise, 0.356 ± 0.111 , it dropped significantly for all these variations of gH (compare Figures 5E, 10B). The lowest mean power was at gH = 0 (0.015 ± 0.009 ) but great variability was observed in this simulation, with the SEM amplitude comparable to the mean value (±0.035 ).
Mean theta peak frequency was elevated (13.5–22.4%) compared to the baseline in all four gH modulations tested.
Effects of changing gKDR in excitatory neurons (gKDRe). For gKDRe (see Figure 10C), theta power trends remained similar to those in the noiseless system; theta was greatly reduced from the baseline condition in all cases, down to 0.029 ± 0.012 at 2 g. However, the SEM values at 0, 0.5, 1.5, and 2 g are all large compared to the mean (±0.052, ±0.028, ±0.040, and ±0.012 , respectively).
Theta frequency remained relatively stable among gKDRe values, lacking the large shift seen in the noiseless case at 0.5 g (compare Figures 6F, 10C).
Effects of changing gKDR in interneurons (gKDRi). Mean theta power remained consistent for most gKDRi cases with noise compared to noiseless runs (compare Figures 7E, 10D), showing clear decreased values for all variations of gKDRi. The notable exception was 0.5 g where without noise, the theta peak was lost. With added noise, it remained present, while still greatly reduced from 0.356 ± 0.111 at baseline to 0.074 ± 0.040 . Although the mean theta power at 2 g was similar to the noiseless equivalent, it displayed significant run to run variability with a large SEM (0.082 ± 0.072 ).
Theta frequency showed deviations at 0.5 and 2 g. The mean frequency did not increase as markedly as noiseless experiments at 0.5 g and had a very large SEM (5.687 ± 1.957 ). An increase in peak theta frequency from 3.766 ± 0.198 Hz at baseline to 4.286 ± 0.711 Hz at 2 g was also observed, unlike noiseless trials.
Effects of changing gKDR in oriens lacunosum interneurons (gKDRo). The addition of noise led to significant deviation in the theta power band with gKDRo variations (see Figure 10E). Though the mean power at 2 g was consistent with noiseless runs, trials at 0, 0.5, and 1.5 g had mean power significantly above their equivalents in Figure 8E. This was particularly apparent at gKDRo = 0 and 1.5 g, where the theta peak was abolished in noiseless runs. With noise, gKDRo = 0 had the highest power in this band (0.192 ± 0.058 ) outside of the baseline conductance. Theta power trends seemed to lower with increasing gKDRo; conversely, after an initial reduction, power increased with decreasing gKDRo. However, the reported standard error values were very large, particularly at 0.5 and 1.5 g, where theta power was 0.130 ± 0.143 and 0.067 ± 0.085, respectively, vs. 0.356 ± 0.111 at baseline. This indicates the broad disparity among the trials leading to occasional elimination of theta peak at some trials of 0.5 and 1.5 g.
Theta frequency mostly followed a similar trend in the noiseless and noisy cases. However, 0.5 g had an increased mean frequency (4.768 Hz) in this band compared to the baseline value, and gKDRo = 0 was unchanged from baseline frequency. These both differ from the noiseless results in Figure 8F. Most theta frequency values had a wide confidence interval between ±0.776 and ±1.178 , except gKDRo = 0 (3.707 ± 0.015 ).
Effects of changing gNaP. gNaP effects with noise are displayed in Figure 10F. As with other conductances, the largest deviations were in theta power (compare with Figure 9E). Interestingly, though no change was observed in noiseless trials when gNaP was decreased, the mean power decreased when doing the same with noise, going to 0.235 ± 0.112 at 20 g and as low as 0.048 ± 0.054 at 50 g. The latter implies that peak theta power would be completely diminished on the lower end of the confidence interval and the disparity levels are significantly high in the gNaP = 50 g trials. Further increasing gNaP still decreased the overall theta power. Nevertheless, the initial decrease observed at gNaP = 20 g was no longer as steep, the mean theta power only being reduced by 34.0% of baseline power. The SEM is considerably large compared to the mean (0.235 ± 0.112 ).
The peak frequency at theta range stayed highly stable, similar to noiseless runs, with its largest change being from 3.766 ± 0.198 Hz at baseline to 4.026 ± 0.240 Hz at gNaP = 80 g.
3.3.2. Effects of Noise With VGIC Modulation on Gamma Rhythmicity
Introducing noise to the system had significant effects on several cases of VGIC modulation for the gamma band. However, many trends remained very similar or identical to those observed in the absence of noise. In most cases, the gamma power and frequency displayed smaller and more uniform SEM values than theta. Overall, this rhythm appeared very robust to the addition of noise into the network.
Effects of changing gA. The introduction of noise interacted significantly with gA amplitude effects on gamma power as shown in Figure 10A. For example, at 0.5 g, gamma power was reduced by 76.2% without noise, but with noise there was an increase of 20.7% of the mean power. At gA = 0, gamma power was absent, while after that increasing gA values up to 2 g reduced gamma power.
Gamma peak frequency showed very little variance across trials and was consistent with noiseless trials; the only exception was for gA= 0, but this was more likely an effect of the peak being negligibly small rather than a significant shift. Interestingly, 0.5 g resulted in a greatly increased peak frequency (95.797 ± 2.436 vs. 59.910 ± 0.033 Hz at baseline), which matched the noiseless run.
Effects of changing gH. Gamma peak power followed the same trend as seen in standard conditions for gH (compare Figure 10B with Figure 5E). This includes the power increase at 1.5 g, and the loss of gamma at gH = 0. However, the power at 0.5 g was greater than corresponding noiseless conditions (0.092 ± 0.015 ). Conversely, power at 2 g was considerably smaller than the noiseless run, measuring 0.015 ± 0.020 —the peak would be completely lost on the lower end of the confidence interval.
The peak frequency for 2 g showed a very wide error margin, and the mean was lower than baseline value (56.080 ± 29.190 vs. 59.910 ± 0.033 Hz). This was in contrary to the standard run, where this conductance value measured the highest frequency gamma peak. The other conductances tested had similar peaks to the noiseless run, with 1.5 g being particularly resistant to variability. However, the mean frequency at gH = 0 was higher than that of 0.5 g under noisy conditions (48.740 vs. 44.031 Hz), opposite to the noiseless network. gH = 0 also showed a large SEM in frequency (±17.432 Hz), but as the peak was functionally absent at this conductance, this variability is less likely to be an effect of network noise.
Effects of changing gKDR in excitatory neurons (gKDRe). Both gamma power and frequency trends remained nearly identical to standard noiseless simulations for gKDRe. The gamma band showed much less fluctuation across noise trials than theta and seemed especially resistant to variation in peak frequency. A notable difference is at gKDRe = 0 where despite gamma power still increasing, the mean (0.383 ± 0.076 ) was higher than the noiseless value (compare with Figure 6E).
Effects of changing gKDR in interneurons (gKDRi). For gKDRi, gamma had mean power and peak frequency values and trends identical to those seen without noise, with only one significant deviation at gKDRi = 0. Rather than being reduced, the mean gamma power observed was significantly higher than the baseline value (0.438 vs. 0.285 ), while the peak frequency remained consistent.
Effects of changing gKDR in oriens lacunosum interneurons (gKDRo). Like the I-cell counterpart, gamma band frequency and power under gKDRo modulation were nearly identical to noiseless trials and followed the same trends.
Effects of changing gNaP. The same general trends were observed in gamma power as the noiseless runs for gNaP, but a significant deviation occurred at gNaP = 50 g. The sharp decrease in peak gamma power at this conductance was no longer observed, as the mean power instead decreased more gradually, by 10.5% of the baseline power.
Gamma peak frequencies were mostly consistent with the noise-free simulations, except for the gamma peak at gNaP = 50 g, which no longer decreased from baseline frequency.
4. Discussion
We have simulated a PING network of pyramidal excitatory neurons, inhibitory interneurons, and oriens lacunosum moleculare interneurons. Subsequently, we have examined the effects of varying VGIC conductance through multiple simulation runs. These experiments were then repeated with noise applied to the E-drive used in the model in order to examine the effect of random noise on these rhythms and to simulate more realistic network conditions. Apart from improving understanding of the fundamental mechanisms of rhythmogenesis in neural networks, these observations are also of interest as many commonly used psychoactive drugs are known to modulate voltage-gated ion channels either directly (e.g., phenytoin and lamotrigine on Nav channels) or indirectly (e.g., SSRIs on Nav and Kv channels and antipsychotic drugs on Kv class channels).
The results in section 3 show that theta rhythmicity is very responsive to changes in VGIC conductance values. The most widely observed effects of conductance modulation were substantial decreases in the power of the theta rhythm and, in many cases, resulting in loss of the peak. There were no observed cases where the power of the theta peak increased compared to baseline. Under noise-free conditions, absence of IA or any of the IKDR currents lead to loss of theta peak, suggesting a likely dependence of the theta peak on these currents. Conversely, relatively high conductance in IA, IH, IKDRe, or INaP resulted in loss of this peak. In light of these findings, introduction of noise appeared to prevent theta loss at high IH, low IKDRo, and both high and low IKDRe. Additionally, removing IH in a noisy network resulted in theta loss instead. The generally large amounts of variance observed across all the trials containing noise suggest that this rhythm may also be highly sensitive to other network conditions as well as VGIC conductance shifts.
Gamma band oscillations were also significantly affected by variation in VGIC conductances. Gamma rhythmicity showed less susceptibility to changing conductance compared to the theta band. Most notably, the gamma band disappeared when the IH and IA currents were set to 0. This implies a likely dependence on these currents for the gamma rhythm to be generated. Specifically, IH and IA seemed to have inverse effects on gamma. Increasing the H-current increased gamma power at first, but this trend was reversed at higher conductances. Conversely, decreasing the A-current had little effect on gamma rhythm at low amplitude changes, but weakened it steeply as the conductance declined further. IKDRi and IKDRe currents demonstrated similar trends whereby gamma power declined with decreasing g values. This trend was not seen in the E-cells, where the absence of IKDRe enhanced gamma. This could imply that the current poses a regulatory effect in maintaining gamma within a certain range. INaP generally decreased the gamma peak as it was elevated, though the extent of the decrease fluctuated. The addition of noise demonstrated less impact on gamma rhythmicity compared to the theta band. Nevertheless, some noteworthy shifts still occurred. Primarily, the IKDRi current showed the same trend as its E-cell counterpart with noise, whereas removal of the current resulted in gamma being amplified. Additionally, decreasing gA under noisy conditions increased gamma power at first. Further decreases in conductance followed the same trend as before, with the peak being lost following removal of the current. Hence, the inverse correlation of IA and IH seems even more apparent under noisy conditions.
Our findings demonstrate a strong association between voltage-gated ion channel conductance and network rhythmicity. Few studies have examined this phenomenon and covered the wide range of VGICs as accomplished in this work. Wang (1993) investigated persistent sodium and slowly inactivating potassium currents and examined how the presence of their interplay produces oscillations in the 10–50 Hz range. Our findings differ from this, as we did not observe any significant decrease in gamma activity with lowered gNaP. Instead, we saw a decrease in theta power following introduction of noise. The Wang model contained fewer currents and observed a single cell environment, so it is unclear if this deviation can simply be attributed to differences in simulation conditions. Another study by Xu et al. (2004) showed that the use of a targeted H-current blocker on septo-hippocampal neurons in vivo resulted in a decrease in hippocampal theta rhythm. Though the range at which they observed their theta peak (5–8 Hz) differs from that seen in our model, we also observed a decrease in theta when lowering gH. In the absence of other clear analogs, it otherwise appears that the observations seen in our simulations are novel findings requiring further investigation. Abnormalities of neural network rhythmicity have been associated with several disorders of the brain. Alzheimer's disease has been associated with a decrease in gamma activity (Herrmann and Demiralp, 2005; Mably and Colgin, 2018) and an increase in theta (Bhattacharya et al., 2011; Vecchio et al., 2013); both phenomena are associated with cognitive deficits. Gamma band irregularities have also been implicated in epilepsy (Herrmann and Demiralp, 2005), schizophrenia (Herrmann and Demiralp, 2005; Roopun et al., 2008; Uhlhaas and Singer, 2015) and depression (Fitzgerald and Watson, 2018). The relevance of network oscillations is further corroborated by research connecting them to memory and learning (Yamamoto et al., 2014). Theta and gamma oscillations, as well as theta–gamma coupling, are believed to be differentially modulated during the learning phase (Priori et al., 2004).
Furthermore, many pharmacological treatments for neurodegenerative disorders have known relevant effects on network rhythmicity. Dopaminergic drugs are found to enhance theta rhythms in Parkinson's disease (Priori et al., 2004). Antipsychotic drugs reduce abnormally high gamma in schizophrenia (Ahnaou et al., 2017). Acetylcholinesterase inhibitors and the NMDA antagonist memantine are cognition-enhancing drugs used in Alzheimer's treatment; they have been found to increase both theta and gamma rhythmicity in rats (Guadagna et al., 2012). Interestingly, some treatment methods have shown effects on VGICs but have not been linked with any observed abnormal rhythms. Inhibition of INa and INaP are well-established as mechanisms of epilepsy treatment (Köhling, 2002; Ahnaou et al., 2014). Potassium channels are also speculated to have a significant role in epilepsy (Errington et al., 2005), as some anti-convulsives seem to work primarily through potassium channel modulation (Köhling, 2002; Wickenden, 2002). Acetylcholinesterase inhibitors have also been found to inhibit both Na+ and K+ VGICs (Gunthorpe et al., 2012).
This study reveals strong modulatory effects of several classes of VGICs on functionally significant neural network oscillations in the utilized PING network model. It will be of interest to explore pharmacological blockade of these channels in vivo, in association with cognitive performance. These results also raise the question whether some commonly used drugs with VGIC targets, such as antipsychotics and antiepileptics, may also target network frequency modulation. Further computational or in vivo studies are required to better characterize and corroborate the effects IA, IH, IKDRe, IKDRi, IKDRo, and INaP currents have on neural oscillations. Investigating all these parameters in the hippocampal neural network model leads to a very broad exploring space. While we demonstrate in Supplementary Material that manipulating the drive currents (specifically IE and II) visibly affects the operating regime and the oscillations of the network, we leave further investigation of these effects to future works. Another consideration for future studies would be to experiment with co-varying VGIC conductances and observing their effects on VGICs, as co-variance has been explored successfully in other work (Ratté et al., 2014). This would further characterize the role of particular VGICs in network rhythmicity as well as their relationship to each other.
The current model successfully reproduces theta and gamma oscillations and is based on plausible physiological data. To the best of our knowledge, the Kopell model is among the very few studies which enable us to study the effects of several voltage gated ion channels. It is also shown that there is concordance between biological phenomena and the results of this model (Scheffer-Teixeira and Tort, 2016; Jansen et al., 2021). Nevertheless, there are also limitations to this model. For instance, alternate regimes can be considered for future iterations of this model or creation of new models which benefit from more realistic pyramidal neuron interconnectivity (Thomson and Bannister, 2003; Douglas and Martin, 2004) and alternative firing patterns (Moca et al., 2014). These are among the elements which were not included in this study for simplicity. Moreover, to avoid excess complexity of the model, this work does not incorporate other important components, such as the three types of glial cells. It has been well-established that all three types of glial cells are important for memory formation and learning (Tewari and Parpura, 2013, 2016; Sibille et al., 2014; Verkhratsky et al., 2015; Gibbs, 2016; Tewari et al., 2016). Therefore, an augmented version of the present model which encompasses other components of the hippocampal region remains to be investigated in the future works. This could lead to a better understanding of all the important and role-playing elements in modulating the local field potential rhythms. Expanding our understanding of all these interactions will likely lead to novel treatments of disorders characterized by neural rhythmic abnormalities.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
CF defined the project. CO and FH conceived and designed the experiments, performed the experiments, analyzed the data, contributed the reagents, materials, analysis tools, and wrote the paper. All authors reviewed the manuscript.
Funding
This work was supported by the RMH Neuroscience Foundation.
Conflict of Interest
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.
Acknowledgments
The authors would like to thank the reviewers for their helpful comments.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncom.2021.630271/full#supplementary-material
References
Ahnaou, A., Huysmans, H., Jacobs, T., and Drinkenburg, W. (2014). Cortical EEG oscillations and network connectivity as efficacy indices for assessing drugs with cognition enhancing potential. Neuropharmacology 86, 362–377. doi: 10.1016/j.neuropharm.2014.08.015
Ahnaou, A., Huysmans, H., Van de Casteele, T., and Drinkenburg, W. (2017). Cortical high gamma network oscillations and connectivity: a translational index for antipsychotics to normalize aberrant neurophysiological activity. Transl. Psychiatry 7, 1–14. doi: 10.1038/s41398-017-0002-9
Başar, E. (2013). Brain oscillations in neuropsychiatric disease. Dialog. Clin. Neurosci. 15:291. doi: 10.31887/DCNS.2013.15.3/ebasar
Başar-Eroğlu, C., Schmiedt-Fehr, C., and Mathes, B. (2013). Auditory-evoked alpha oscillations imply reduced anterior and increased posterior amplitudes in schizophrenia. Suppl. Clin. Neurophysiol. 62, 121–129. doi: 10.1016/B978-0-7020-5307-8.00008-9
Bhattacharya, B. S., Coyle, D., and Maguire, L. P. (2011). Alpha and theta rhythm abnormality in Alzheimer's disease: a study using a computational model. Adv. Exp. Med. Biol. 718, 57–73. doi: 10.1007/978-1-4614-0164-3_6
Börgers, C., Epstein, S., and Kopell, N. J. (2005). Background gamma rhythmicity and attention in cortical local circuits: a computational study. Proc. Natl. Acad. Sci. U.S.A. 102, 7002–7007. doi: 10.1073/pnas.0502366102
Börgers, C., and Kopell, N. (2005). Effects of noisy drive on rhythms in networks of excitatory and inhibitory neurons. Neural Comput. 17, 557–608. doi: 10.1162/0899766053019908
Butler, J. L., Mendonça, P. R., Robinson, H. P., and Paulsen, O. (2016). Intrinsic cornu ammonis area 1 theta-nested gamma oscillations induced by optogenetic theta frequency stimulation. J. Neurosci. 36, 4155–4169. doi: 10.1523/JNEUROSCI.3150-15.2016
Buzsáki, G. (2002). Theta oscillations in the hippocampus. Neuron 33, 325–340. doi: 10.1016/S0896-6273(02)00586-X
Chambers, J., Bethwaite, B., Diamond, N., Peachey, T., Abramson, D., Petrou, S., et al. (2012). Parametric computation predicts a multiplicative interaction between synaptic strength parameters that control gamma oscillations. Front. Comput. Neurosci. 6:53. doi: 10.3389/fncom.2012.00053
Colgin, L. L., and Moser, E. I. (2010). Gamma oscillations in the hippocampus. Physiology 25, 319–329. doi: 10.1152/physiol.00021.2010
Douchamps, V., Jeewajee, A., Blundell, P., Burgess, N., and Lever, C. (2013). Evidence for encoding versus retrieval scheduling in the hippocampus by theta phase and acetylcholine. J. Neurosci. 33, 8689–8704. doi: 10.1523/JNEUROSCI.4483-12.2013
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
Errington, A. C., Stohr, T., and Lees, G. (2005). Voltage gated ion channels: targets for anticonvulsant drugs. Curr. Top. Med. Chem. 5, 15–30. doi: 10.2174/1568026053386872
Fitzgerald, P. J., and Watson, B. O. (2018). Gamma oscillations as a biomarker for major depression: an emerging topic. Transl. Psychiatry 8, 1–7. doi: 10.1038/s41398-018-0239-y
Gibbs, M. E. (2016). Role of glycogenolysis in memory and learning: regulation by noradrenaline, serotonin and atp. Front. Integr. Neurosci. 9:70. doi: 10.3389/fnint.2015.00070
Gloveli, T., Dugladze, T., Rotstein, H. G., Traub, R. D., Monyer, H., Heinemann, U., et al. (2005). Orthogonal arrangement of rhythm-generating microcircuits in the hippocampus. Proc. Natl. Acad. Sci. U.S.A. 102, 13295–13300. doi: 10.1073/pnas.0506259102
Guadagna, S., Bundgaard, C., Hovelsø, N., Volbracht, C., Francis, P. T., Egebjerg, J., et al. (2012). Memantine potentiates hippocampal theta oscillations at a therapeutic dose in anesthetized mice: a mechanistic link to its cognitive-enhancing properties. Neuropharmacology 62, 2208–2218. doi: 10.1016/j.neuropharm.2012.01.014
Gunthorpe, M. J., Large, C. H., and Sankar, R. (2012). The mechanism of action of retigabine (ezogabine), a first-in-class K+ channel opener for the treatment of epilepsy. Epilepsia 53, 412–424. doi: 10.1111/j.1528-1167.2011.03365.x
Herrmann, C., and Demiralp, T. (2005). Human eeg gamma oscillations in neuropsychiatric disorders. Clin. Neurophysiol. 116, 2719–2733. doi: 10.1016/j.clinph.2005.07.007
Hsu, C. L., Zhao, X., Milstein, A. D., and Spruston, N. (2018). Persistent sodium current mediates the steep voltage dependence of spatial coding in hippocampal pyramidal neurons. Neuron 99, 147–162. doi: 10.1016/j.neuron.2018.05.025
Jacobson, T. K., Howe, M. D., Schmidt, B., Hinman, J. R., Escabí, M. A., and Markus, E. J. (2013). Hippocampal theta, gamma, and theta-gamma coupling: effects of aging, environmental change, and cholinergic activation. J. Neurophysiol. 109, 1852–1865. doi: 10.1152/jn.00409.2012
Jansen, N. A., Perez, C., Schenke, M., van Beurden, A. W., Dehghani, A., Voskuyl, R. A., et al. (2021). Impaired θ-γ coupling indicates inhibitory dysfunction and seizure risk in a dravet syndrome mouse model. J. Neurosci. 41, 524–537. doi: 10.1523/JNEUROSCI.2132-20.2020
Keeley, S., Fenton, A. A., and Rinzel, J. (2017). Modeling fast and slow gamma oscillations with interneurons of different subtype. J. Neurophysiol. 117, 950–965. doi: 10.1152/jn.00490.2016
Klausberger, T., Magill, P. J., Márton, L. F., Roberts, J. D. B., Cobden, P. M., Buzsáki, G., et al. (2003). Brain-state-and cell-type-specific firing of hippocampal interneurons in vivo. Nature 421, 844–848. doi: 10.1038/nature01374
Köhling, R. (2002). Voltage-gated sodium channels in epilepsy. Epilepsia 43, 1278–1295. doi: 10.1046/j.1528-1157.2002.40501.x
Kopell, N., Börgers, C., Pervouchine, D., Malerba, P., and Tort, A. (2010). “Gamma and theta rhythms in biophysical models of hippocampal circuits,” in Hippocampal Microcircuits. Springer Series in Computational Neuroscience, Vol. 5, eds V. Cutsuridis, B. Graham, S. Cobb, and I. Vida (New York, NY: Springer). doi: 10.1007/978-1-4419-0996-1_15
Lisman, J. E., and Jensen, O. (2013). The theta-gamma neural code. Neuron 77, 1002–1016. doi: 10.1016/j.neuron.2013.03.007
Mably, A. J., and Colgin, L. L. (2018). Gamma oscillations in cognitive disorders. Curr. Opin. Neurobiol. 52, 182–187. doi: 10.1016/j.conb.2018.07.009
Maccaferri, G., and McBain, C. J. (1996). The hyperpolarization-activated current (Ih) and its contribution to pacemaker activity in rat CA1 hippocampal stratum oriens-alveus interneurones. J. Physiol. 497, 119–130. doi: 10.1113/jphysiol.1996.sp021754
Moca, V. V., Nikolić, D., Singer, W., and Mureşan, R. C. (2014). Membrane resonance enables stable and robust gamma oscillations. Cereb. Cortex 24, 119–142. doi: 10.1093/cercor/bhs293
Nadel, L., and Moscovitch, M. (1997). Memory consolidation, retrograde amnesia and the hippocampal complex. Curr. Opin. Neurobiol. 7, 217–227. doi: 10.1016/S0959-4388(97)80010-4
Oswal, A., Brown, P., and Litvak, V. (2013). Synchronized neural oscillations and the pathophysiology of Parkinson's disease. Curr. Opin. Neurol. 26, 662–670. doi: 10.1097/WCO.0000000000000034
Priori, A., Foffani, G., Pesenti, A., Tamma, F., Bianchi, A. M., Pellegrini, M., et al. (2004). Rhythm-specific pharmacological modulation of subthalamic activity in Parkinson's disease. Exp. Neurol. 189, 369–379. doi: 10.1016/j.expneurol.2004.06.001
Ratté, S., Zhu, Y., Lee, K. Y., and Prescott, S. A. (2014). Criticality and degeneracy in injury-induced changes in primary afferent excitability and the implications for neuropathic pain. Elife 3:e02370. doi: 10.7554/eLife.02370
Rich, S., Zochowski, M., and Booth, V. (2017). Dichotomous dynamics in ei networks with strongly and weakly intra-connected inhibitory neurons. Front. Neural Circuits 11:104. doi: 10.3389/fncir.2017.00104
Roopun, A. K., Cunningham, M. O., Racca, C., Alter, K., Traub, R. D., and Whittington, M. A. (2008). Region-specific changes in gamma and beta2 rhythms in nmda receptor dysfunction models of schizophrenia. Schizophr. Bull. 34, 962–973. doi: 10.1093/schbul/sbn059
Rutishauser, U., Ross, I. B., Mamelak, A. N., and Schuman, E. M. (2010). Human memory strength is predicted by theta-frequency phase-locking of single neurons. Nature 464, 903–907. doi: 10.1038/nature08860
Scheffer-Teixeira, R., and Tort, A. B. (2016). On cross-frequency phase-phase coupling between theta and gamma oscillations in the hippocampus. Elife 5:e20515. doi: 10.7554/eLife.20515
Sekulić, V., and Skinner, F. K. (2017). Computational models of o-lm cells are recruited by low or high theta frequency inputs depending on h-channel distributions. eLife 6:e22962. doi: 10.7554/eLife.22962
Sibille, J., Pannasch, U., and Rouach, N. (2014). Astroglial potassium clearance contributes to short-term plasticity of synaptically evoked currents at the tripartite synapse. J. Physiol. 592, 87–102. doi: 10.1113/jphysiol.2013.261735
Tewari, S., and Parpura, V. (2013). A possible role of astrocytes in contextual memory retrieval: an analysis obtained using a quantitative framework. Front. Comput. Neurosci. 7:145. doi: 10.3389/fncom.2013.00145
Tewari, S. G., Gottipati, M. K., and Parpura, V. (2016). Mathematical modeling in neuroscience: neuronal activity and its modulation by astrocytes. Front. Integr. Neurosci. 10:3. doi: 10.3389/fnint.2016.00003
Tewari, S. G., and Parpura, V. (2016). Astrocytes modulate local field potential rhythm. Front. Integr. Neurosci. 9:69. doi: 10.3389/fnint.2015.00069
Thomson, A. M., and Bannister, A. P. (2003). Interlaminar connections in the neocortex. Cereb. Cortex 13, 5–14. doi: 10.1093/cercor/13.1.5
Tort, A. B., Rotstein, H. G., Dugladze, T., Gloveli, T., and Kopell, N. J. (2007). On the formation of gamma-coherent cell assemblies by oriens lacunosum-moleculare interneurons in the hippocampus. Proc. Natl. Acad. Sci. U.S.A. 104, 13490–13495. doi: 10.1073/pnas.0705708104
Traub, R. D., Jefferys, J. G., and Whittington, M. A. (1997). Simulation of gamma rhythms in networks of interneurons and pyramidal cells. J. Comput. Neurosci. 4, 141–150. doi: 10.1023/A:1008839312043
Uhlhaas, P. J., and Singer, W. (2015). Oscillations and neuronal dynamics in schizophrenia: the search for basic symptoms and translational opportunities. Biol. Psychiatry 77, 1001–1009. doi: 10.1016/j.biopsych.2014.11.019
Varga, C., Golshani, P., and Soltesz, I. (2012). Frequency-invariant temporal ordering of interneuronal discharges during hippocampal oscillations in awake mice. Proc. Natl. Acad. Sci. U.S.A. 109, E2726–E2734. doi: 10.1073/pnas.1210929109
Vecchio, F., Babiloni, C., Lizio, R., Fallani, F. D. V., Blinowska, K., Verrienti, G., et al. (2013). Resting state cortical EEG rhythms in Alzheimer's disease: toward EEG markers for clinical applications: a review. Suppl. Clin. Neurophysiol. 62, 223–236. doi: 10.1016/B978-0-7020-5307-8.00015-6
Verkhratsky, A., Horvat, A., and Vardjan, N. (2015). Memory formation shaped by astroglia. Front. Integr. Neurosci. 9:56. doi: 10.3389/fnint.2015.00056
Wang, X. J. (1993). Ionic basis for intrinsic 40 Hz neuronal oscillations. Neuroreport 5, 221–224. doi: 10.1097/00001756-199312000-00008
Wang, X. J. (2010). Neurophysiological and computational principles of cortical rhythms in cognition. Physiol. Rev. 90, 1195–1268. doi: 10.1152/physrev.00035.2008
Whittington, M. A., Traub, R., Kopell, N., Ermentrout, B., and Buhl, E. (2000). Inhibition-based rhythms: experimental and mathematical observations on network dynamics. Int. J. Psychophysiol. 38, 315–336. doi: 10.1016/S0167-8760(00)00173-2
Wickenden, A. D. (2002). Potassium channels as anti-epileptic drug targets. Neuropharmacology 43, 1055–1060. doi: 10.1016/S0028-3908(02)00237-X
Xu, C., Datta, S., Wu, M., and Alreja, M. (2004). Hippocampal theta rhythm is reduced by suppression of the H-current in septohippocampal gabaergic neurons. Eur. J. Neurosci. 19, 2299–2309. doi: 10.1111/j.0953-816X.2004.03316.x
Yamamoto, J., Suh, J., Takeuchi, D., and Tonegawa, S. (2014). Successful execution of working memory linked to synchronized high-frequency gamma oscillations. Cell 157, 845–857. doi: 10.1016/j.cell.2014.04.009
Keywords: neural oscillations, theta rhythmogenesis, gamma rhythmogenesis, voltage-gated ion channels, pyramidal-interneuron-gamma network, hippocampal neural network
Citation: Olteanu C, Habibollahi F and French C (2021) Effects of Several Classes of Voltage-Gated Ion Channel Conductances on Gamma and Theta Oscillations in a Hippocampal Microcircuit Model. Front. Comput. Neurosci. 15:630271. doi: 10.3389/fncom.2021.630271
Received: 17 November 2020; Accepted: 08 March 2021;
Published: 01 April 2021.
Edited by:
Jue Zhang, Peking University, ChinaReviewed by:
Shivendra Gajraj Tewari, Biotechnology HPC Software Applications Institute (BHSAI), United StatesRaul Cristian Muresan, Transylvanian Institute of Neuroscience (TINS), Romania
Copyright © 2021 Olteanu, Habibollahi and French. 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.
*Correspondence: Chris French, frenchc@unimelb.edu.au
†These authors have contributed equally to this work