Front. Comput. Neurosci. Frontiers in Computational Neuroscience Front. Comput. Neurosci. 1662-5188 Frontiers Media S.A. 10.3389/fncom.2013.00145 Neuroscience Original Research Article A possible role of astrocytes in contextual memory retrieval: An analysis obtained using a quantitative framework Tewari Shivendra 1 2 * Parpura Vladimir 3 4 * 1Biotechnology and Bioengineering Center, Medical College of Wisconsin Milwaukee, WI, USA 2Department of Physiology, Medical College of Wisconsin Milwaukee, WI, USA 3Department of Neurobiology, University of Alabama at Birmingham Birmingham, AL, USA 4Department of Biotechnology, University of Rijeka Rijeka, Croatia

Edited by: Gabriel A. Silva, University of California, San Diego, USA

Reviewed by: Gabriel A. Silva, University of California, San Diego, USA; Massoud L. Khraiche, University of California San Diego, USA

*Correspondence: Shivendra Tewari, Molecular & Integrative Physiology, University of Michigan, NCRC Building 10 #A122, 2800 Plymouth Rd, Ann Arbor, MI 48109, USA e-mail: tewarisg@gmail.com;

Vladimir Parpura, Department of Neurobiology, University of Alabama at Birmingham, 1719 6th Avenue South, CIRC 429, Birmingham, AL 35294, USA e-mail: vlad@uab.edu

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

 05 11 2013 2013 7 145 09 07 2013 03 10 2013 Copyright © 2013 Tewari and Parpura. 2013

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

The hippocampus is central to our understanding of memory formation and retrieval. Its CA1 region is known for encoding contextual memory. Here, using a computational approach, which embeds existing physiological data, we propose a particular role of astrocytes in contextual memory retrieval. We provide a quantitative framework under which the astrocyte modulates the firing of a context-associated CA1 pyramidal neurons, resulting in a prominent tuning of neurons to a delta rhythm. Using the very framework, we further studied astrocytic function in the modulation of neuronal firing under pathological conditions, i.e., during astrocytic induction of epileptiform discharge in CA1 pyramidal neurons. Thus, we provide a quantitative framework that would aid understanding of the Schaffer collateral-CA1 tripartite synapse in health and disease.

 hippocampus CA3 region CA1 region contextual memory astrocytes computational biology Introduction

One of the most intriguing properties of the brain is to learn and form new memories. Memories are essential for guiding effective actions, for understanding the words we read, for recognizing the objects we see, etc. It is well established that retrieval of a particular memory is facilitated when the context at the time of memory retrieval is similar to that of the initial learning experience (Maren and Holt, 2000). This phenomenon, termed as context-dependent memory, also provides an explanation as to why memory recall is impaired when people, place or things are experienced “out of context.” The process by which the context is embedded into learning of information, and can later facilitate memory retrieval is known as contextual memory retrieval (Maren and Holt, 2000). Pavlovian conditioning, i.e., classical conditioning, is an exemplar form of contextual fear memory learning and retrieval (Maren, 2001).

Ever since memory deficits, initially characterized in patients PB and FC (Milner and Penfield, 1955) and later HM (Scoville and Milner, 1957), were associated with hippocampus damage, it has become clear that the hippocampal region plays a crucial role in the acquisition of new semantic and episodic memory (Andersen et al., 2007); but see (Corkin, 2002). Apart from storing memories, the appropriate retrieval of stored memories is also of paramount importance for maintaining relationships, achieving personal and professional goals, and carrying out even simple day to day activities. The critical role of the hippocampus in memory retrieval has gone largely uncontested and has been supported by various studies including lesion studies (Holt and Maren, 1999), unit studies (Paller and McCarthy, 2002), as well as functional neuroimaging (Eldridge et al., 2000).

Lee and Kesner (2004) reported a sub–region specific contribution of the hippocampus in memory acquisition and retrieval. Interpretation of the result suggests that the CA1 region is uniquely associated with a reliable retrieval of the contextual memory. Such a role for the CA1 region in contextual memory retrieval has also been highlighted in several studies (Hall et al., 2001; Khuchua et al., 2003; Murchison et al., 2011). Previous work suggests that the firing synchrony in neighboring CA1 pyramidal neurons might be a putative mechanism of contextual memory retrieval in the CA1 region (Hall et al., 2001; Siekmeier et al., 2007; Takahashi and Sakurai, 2009). However, data are accumulating suggesting that the coordination between neurons and astrocytes may play a role in information processing and memory retrieval.

Fellin et al. (2004) demonstrated that the activation of Ca2+ elevation in astrocytes, by various stimuli, produces slow inward currents (SICs) in 2–12 CA1 pyramidal neurons, synchronously. They found that repetitive stimulation of astrocytes leads to synchronous firing of the same CA1 pyramidal neurons, mainly due to the activation of neuronal NR1/NR2B subunits of the extra–synaptic N–Methyl-D–aspartate receptors (eNMDARs). It should be noted that NMDAR NR2B subtype plays a crucial role in encoding and retrieval of contextual fear memory (Zhao et al., 2005; Zhang et al., 2008; Zhuo, 2009); but see Nakazawa et al. (2004). Thus, it appears that the interaction/signaling between neurons and astrocytes plays a role in information processing and memory retrieval.

Astrocytes in the hippocampus have numerous highly branched fine processes that extend to distances up to 100 μm and enwrap asymmetrical (mainly excitatory/glutamatergic) synapses of a relatively small group of neurons (Fellin et al., 2006); of note, individual astrocytic domains do not overlap (Bushong et al., 2002, 2004; Ogata and Kosaka, 2002; Nedergaard et al., 2003; Oberheim et al., 2009). Thus, strategically, astrocytes are ideally positioned to sample electrical activity patterns of the Schaffer collateral input from CA3 pyramidal neurons and retrieve their associated context in neighboring CA1 pyramidal neurons. During the retrieval, the firing pattern of the pyramidal neurons, in the CA3–CA1 region, is quite similar to that recorded during memory formation (Small et al., 2001). As astrocytes play an important role in short-term and long-term potentiation of individual synapses (Martin et al., 2007; Henneberger et al., 2010), it is possible that they also contribute to various aspects of memory. We propose that astrocytes facilitate the memory storage and retrieval by sensing encoding/retrieval patterns of neurons that reside within individual astrocytic domains and, in turn, astrocytes signal the associated CA1 pyramidal neurons.

Based on these recent experimental observations, we developed a computational framework/model, where an astrocyte modulates synchrony among CA1 pyramidal neurons within a domain of 100 μm in diameter, which qualitatively reproduces the biological experimental findings (Fellin et al., 2004). Power spectral density analysis of the computed local field potentials indicates a shift in neuronal firing rhythm from broad-frequency range in the absence of astrocytic gliotransmission to a delta-dominant rhythm when gliotransmission is turned on. As delta rhythm has been implicated in cognitive domain (Knyazev, 2012) and consolidation of declarative memories (Ruch et al., 2012), this might be a putative mechanism of memory retrieval tuned by astrocytes, which will require further experimental verification in biological preparations. Apart from a possible role in memory retrieval, using the computational model, we propose that astrocytes may lead to epileptiform discharges in CA1 pyramidal neurons, as has been previously demonstrated in biological preparations (Kang et al., 2005; Tian et al., 2005). Such pathological discharges occuring under computational conditions warrant an increase in eNMDAR plasmalemmal activity that could arise from an overexpression of such receptors and/or positive modulation of their conductance (Omkumar et al., 1996; Liao et al., 2001).

 Materials and methods Overview of the model architecture

The complete network (shown in Figure 1) comprises of three major components: a CA3 pyramidal neuron, four CA1 pyramidal neurons and an astrocyte in the stratum oriens (s.o.). A CA3 pyramidal neuron is computationally described by Traub's branching dendrite model of the CA3 pyramidal neuron comprising of 21 compartments (Traub et al., 1994). Each of the four CA1 pyramidal neurons is described by the 19-compartment CA1 pyramidal neuron model (Traub et al., 1991), which was modified to incorporate four additional compartments accounting for the axon initial segment (IS), axon proper, and two spines (placed at the 6th and 12th compartment of the original CA1 pyramidal neuron model) to yield a 23-compartment CA1 pyramidal neuron model. Axon and spine compartments were added to form CA1–CA1 synapses in the s.o. at the basal dendrites. The dendritic spine model utilized here is from Tewari and Majumdar (Tewari and Majumdar, 2012a), but was modified to include synaptic NMDARs and eNMDARs. The astrocyte model is also adopted from Tewari and Majumdar (Tewari and Majumdar, 2012a,b), and it accounts for gliotransmitter (glutamate) release/astrocyte-dependent SICs through the activation of NR2B-containing receptors on the CA1 neurons. It is assumed that the CA3 pyramidal neuron makes 25 synapses with the CA1 neuron 1 (Figure 1). The number of synapses at the CA1–CA1 synapses is assumed to be 20. These numbers were obtained by a hit and trial approach, and reflect the smallest numbers possible to have a successful signal transduction of the proposed network. All associated model equations and parameters are provided in the supplementary materials.

Schematic diagram of the proposed network model of contextual memory retrieval within the hippocampus where an astrocyte modulates synchronous firing of a pack of four CA1 pyramidal neurons during context-dependent memory retrieval. The pack of neurons is aligned within 100 μm of the stratum pyramidale along the input, i.e., Schaffer collateral (SC) axons of the CA3 pyramidal neuron, which makes synapses with the spines of the proximal apical dendrites of CA1 pyramidal cells in the stratum radiatum. The field potential recording electrode is placed next to the CA1 neuron 4, at the most distal end to the input of this pack. The fine processes of the astrocyte, allied with the neuronal pack and located in the stratum oriens, listen to the synaptic transmission which takes place in the stratum radiatum, at the synapse made by the CA3 pyramidal neuron onto the lead CA1 pyramidal neuron 1 (the utmost left, proximal to the input). The astrocyte integrates received information and relays the signal to the following CA1 pyramidal neurons within the pack (in this case the three remaining CA1 pyramidal neurons 2–4, distal to the input) associated with the context of the input signal from the CA3 pyramidal neuron concurrently, the astrocyte also potentiates the SC-CA1 synapse strengthening, the small network activity.

As shown in Figure 1, the CA3 neuron forms synapses with the CA1 neuron 1, at the spines in the stratum radiatum (s.r.). To model glutamate release from the CA3 and CA1 neurons, we use a continuous function that transforms axon voltage into glutamate concentration in the cleft (Destexhe et al., 1994): [g]=gmax1+exp​[ −(Vaxon−Vp)/Kp ], where gmax is the maximum concentration of glutamate in the synaptic cleft, Vaxon is the membrane potential in the axon compartment of the glutamate releasing neuron, Kp gives the steepness, and Vp is the axon voltage at which glutamate concentration is half of the maximum (see Table 1). The glutamate released from CA3 neuron is sensed by the astrocyte in the s.o. (Perea and Araque, 2005) leading to astrocytic Ca2+ elevation. Simultaneously, glutamate released from the CA3 neuron leads to an inward current in the CA1 neuron 1. The SICs evoked by the astrocyte through eNMDARs are: ISIC=gNR2B·h·vs, where gNR2B is the conductance through NR1/NR2B-containing receptors, vs is the spine membrane potential and h is the gating variable of eNMDAR that depends on the astrocyte as described by the following equation: dhdt=στ1−hτ2.

List of parameters involved at the synapses.

Symbol Value Reference
gmax 3.5–4 mM Danbolt, 2001
Vp 72 mV† Destexhe et al., 1994
Kp 5 mV Destexhe et al., 1994
gNR2B 0.6 nS Gasparini et al., 2004
τ1 92.3 ms Fellin et al., 2004
τ2 568.5 ms Fellin et al., 2004

This value is reported with respect to the resting membrane potential.

The above relation is regulated by astrocytic Ca2+ by σ, which has a value of 0.4 (representing the maximum open probability of eNMDAR) whenever astrocytic Ca2+ goes above 200 nM (Parpura and Haydon, 2000; Nadkarni and Jung, 2003; Tewari and Majumdar, 2012b), otherwise set to zero. The time constants τ1 and τ2 govern the rise-time and decay-time of the eNMDAR currents, and are based on the experimental observation (Fellin et al., 2004); parameter values are listed in Table 1.

A given CA1 neuron forms synapses with the other CA1 neurons in the basal dendrite region (Andersen et al., 2007); in the proposed network model, the synapses are formed at the 6th compartment of the CA1 neuron model. The sequence in which the CA1 neurons make synapses is as follows: 1→2, 2→3, 3→4 (Figure 1, left to right); the CA1 neuron 4 is not connected to any other neuron. The astrocyte senses the signal (glutamate) released by the CA3 neuron and potentiates the above CA1 neuron synapses by releasing gliotransmitter (glutamate) in the vicinity of the eNMDARs on the spines, which leads to synchronized firing of the CA1 neurons. The astrocyte releases glutamate from glutamate secretory vesicles utilizing Ca2+-dependent exocytosis as soon as astrocytic Ca2+ concentration exceeds the threshold Ca2+ concentration for triggering secretion [see (Tewari and Majumdar, 2012a) for details].

 Computer simulations

The proposed network model invoked a total of 1247 ordinary differential equations which were solved using a quasi-constant step sized backward-difference numerical solver implemented in MATLAB (The MathWorks, Natick, MA). The system of ordinary differential equations was solved for a duration of 50 s (model time) which took approximately 8 h (real time) to compute.

 Local field potential calculations

The somata of the CA1 neurons were spaced 25 μm apart. The recording electrode was placed 25 μm away from the soma of the CA1 neuron 4 (see Figure 1). The local field potential was calculated using line source approximation (Gold et al., 2006). Based on the ionic currents generated from the network simulation (in the presence and absence of the astrocyte), trans-membrane currents across each compartment were calculated for each of the CA1 neurons. Please note that the contribution from IS and axon compartments was neglected because myelinated axon and nodes of Ranvier make no contribution to the local field potential (Gold et al., 2006). Afterwards, linear source approximation was used to compute local potentials at the recording electrode using a custom-written MATLAB program. Mathematically, field potential is calculated using the following expression: Φ(t)=∑∀cρcIc4πrc, where Ic is the trans-membrane current across compartment c, rc is the distance of compartment c from the recording electrode, and ρc is the extracellular resistivity at compartment c. The value of ρc for soma and proximal apical dendritic compartments is 643 Ω cm (Lopez-Aguado et al., 2001), while for the basal and distal apical dendritic compartments it is 260 and 287 Ω cm, respectively (Lopez-Aguado et al., 2001).

 Results Network simulations

The network is simulated in the absence and in the presence of astrocytic activity, i.e., a feedback to neurons in respect to synaptic transmission. Figure 2 shows representative action potential (AP) bursts in the axon compartment of the CA3 neuron that lead to subsequent firing in the CA1 neurons 1–4, with a slight delay, as expected from the flow of the input signal in the CA1 region. The network was then simulated again, while incorporating the astrocytic response to neuronal synaptic transition, i.e., gliotransmission, that leads to an induction of SICs in the CA1 neurons. We show a representative astrocytic response and its associated effects on the firing of the CA1 neurons (Figure 3). It is apparent by comparing the two different simulations (compare Figures 2 and 3; and Figures 4A and 4B) that the astrocyte makes the CA1 neurons to fire in synchrony, and in a different pattern. Note that the astrocyte induces SIC twice (data not shown; this effect is captured by the gating variable h in Equation 3) while the astrocytic Ca2+ concentration is above the threshold for triggering exocytotic release of glutamate (see the shaded area in Figure 3).

Model simulations of the CA3–CA1 pyramidal neuron network (shown in Figure 1), in the absence of s.o. astrocytic activity, showing action potential generation in the axon compartment of the CA3 and CA1 pyramidal neuron models. The CA3 neuron was given a constant current of 0.6 nA, in the soma compartment, which evokes a high-frequency burst of ~100 Hz and of various durations, which repeats itself at ~0.66 Hz. There were corresponding responses in the post-synaptic CA1 neurons.

Model simulations of the CA3–CA1 pyramidal neuron network (shown in Figure 1), in the presence of s.o. astrocytic activity i.e., gliotransmission. The CA3 neuron was stimulated again in the soma with an input current of 0.6 nA (therefore the CA3 neuron activity is not shown). Representative astrocytic Ca2+ cycle and its effect on the CA1 neurons are shown. The shaded area of the astrocytic Ca2+ indicates the time during which SICs are induced in the basal dendrites of the CA1 neurons 2–4 and the apical dendrites of the CA1 neuron 1.

Calculated field potential based on model simulations of the network in the absence (A) and in the presence (B) of the astrocyte. The simulation omits the first 10 s of simulation as, in the model, astrocytic Ca2+ oscillations begin after 10 s of the CA3 neuronal activity.

 Local field potentials

Local field potential (LFP) calculation at the point of interest is based on linear source approximation in the absence and presence of the astrocyte. It is apparent from the peaks of LFP that the CA1 neurons fire synchronously under the astrocyte-dependent pathway (Figures 4 and 5). Hence, after the addition of the astrocyte-induced SICs, the CA1 neuron pack starts firing at a different rate, which is dependent on astrocyte Ca2+ oscillations (Figure 4). The LFP duration is much shorter in the absence of astrocytic activity (Figure 5A) compared to the LFP duration in the presence of an active astrocyte (Figure 5B); the input signal, from the CA3 neuron, was same for both pathways. To assess a possible difference in the dominant frequency between the two pathways, we used Fast Fourier Transform (FFT) spectrum analysis (Figures 5D,E). Based on this analysis, the astrocyte-independent pathway has wide-spread activity with multiple peaks, a pattern likely resulting due to the differences in the firing rate of each CA1 neuron within the pack; the main proportion of the activity forms an envelope, around 4–11 Hz, lacking any typical rhythm (see Figure 5D). The FFT analysis of the astrocyte-dependent pathway shows a shift towards lower frequency with sharpened tuning, resulting in an occurrence of a delta rhythm (Figure 5E, asterisk), possibly due to synchronous firing of the CA1 neurons in tune with the astrocyte gliotransmission.

Astrocytic contribution to field potential recordings. Duration of evoked field potentials in the absence of astrocyte (A), in the presence of astrocyte activity (B), and in the presence of astrocytic activity along with an enhanced eNMDAR conductance in CA1 neurons (C). Corresponding FFT spectra are shown in (D,E, and F), respectively. Note that FFT spectra are calculated based on complete simulated data of 50 s. Asterisks denote the effect of astrocyte-induced SICs.

 Epileptiform simulations

NMDARs are known to play a role in the initiation and propagation of epileptiform discharges (Auzmendi et al., 2009). To study the possible effect of glutamatergic gliotransmission on the eNMDARs, as demonstrated in biological preparations [e.g., (Fellin et al., 2004)], we simulated the astrocyte-dependent pathway after doubling the eNMDAR conductance in the CA1 neurons (Table 1). This computational alteration emulates a biological condition of an increased presence of plasmalemmal eNMDARs (due to an overexpression/increased trafficking) and/or flux through an unaltered number of channels. As a control, we simulated the same network under conditions incapacitating gliotransmission of the s.o. astrocyte, caused by clamping astrocytic Ca2+ at 100 nM, which is below the threshold for glutamate release. The simulation data suggest that the presence of the astrocyte lacking gliotransmission does not affect the firing of the CA1 neurons overexpressing eNMDARs upon their stimulation by the CA3 neuron (compare Figure S1 and Figure 2). However, the astrocytic activity/gliotransmission may induce epileptiform discharges in the CA1 neurons via an increased eNMDAR activity, as evidenced by the grossly exaggerated duration of LFPs (Figure 5C). FFT spectrum analysis of the epileptiform activity further indicates that the astrocyte activity tunes the CA1 neurons to fire in a delta rhythm, while other surrounding frequencies appeared suppressed (Figure 5F, asterisk).

 Discussion

Herein we presented a biophysically detailed mathematical framework of the CA3/Schaffer collateral-CA1 pyramidal neuron signal transduction that is based on experimental reports suggesting a possible role of astrocytes in contextual memory retrieval by synchronizing CA1 pyramidal neuronal network and achieving a delta rhythm. We were interested in studying a small neuronal network for two reasons: (i) a small network model can be kept biophysically detailed, and (ii) to overcome the experimental limitations in biological preparations (i.e., a lack of the spatial resolution to record field potentials from such a small network, as opposed to recording from a much larger area in reality). We did analysis of the simulated local field potentials generated in the absence and in the presence of astrocytic activity. The results suggest that astrocytic gliotransmission potentiates the synapses through eNMDARs, synchronizes CA1 neurons and modulates their overall firing rate. Frequency domain analysis of the simulated LFP suggests that astrocytic Ca2+ signaling shifts the neuronal activity in the CA1 region from a broad spectrum frequency lacking any typical rhythm towards a delta frequency. Interestingly, a recent study (Grover et al., 2009) described the induction of long-term potentiation in the CA1 region due to a burst stimulation over a broad frequency range centered around delta frequency, albeit the underlying mechanism was not investigated. Thus, our in silico model offers a possible computational explanation to this delta frequency and it should be additionally explored in biological preparations.

Delta rhythm is a low frequency brain activity which is implicated in cognitive domain such as attention, salience detection, and subliminal perception (Knyazev, 2012). Moreover, it is known that sensory cues innervate the CA1 pyramidal neurons associated with the context. Under such circumstances repetition of an item stimulus will recover the contextual state, in the entorhinal cortex, in response to firing of the associated pyramidal cells in the CA1 region. In our circuitry model, we observe that astrocytic input, which makes CA1 pyramidal neurons to fire in a synchronous manner, eliminates group firing activity at other frequencies except delta, suggesting a possible role of astrocytes in memory retrieval. From a functional point of view, this finding proposes a cognitive role of astrocytes where an environmental cue demands attention, e.g., delta rhythm generation in visual cortex depends on the activity to which attention is paid (Lakatos et al., 2008). The model suggests that similar activity of delta rhythm generation might exist in the CA1 region of hippocampus during events of contextual memory retrieval.

NR2B subunit overexpression enhances chronic pain in mice (Zhuo, 2009). Additionally, NR2B-containing NMDARs are linked to cell death and seizures (Ding et al., 2007). To test the functional consequence of NR2B overexpression on CA1 neuronal activity, we increased the conductance of NR2B-containing NMDARs which is equivalent to an insertion of eNMDARs in the plasmalemma or enhanced unitary conductance of the eNMDARs. Such a condition simulates an increased access of astrocytic transmission to the NR2B-containing NMDARs. The results suggest that enhanced conductance of NR2B-containing NMDARs leads to a high-frequency oscillation of field potentials. The duration of these field potentials is ~600 ms which is equal to the duration of SIC via NR2B-containing NMDARs. This activity could be a molecular correlate of seizures seen in the hippocampus, which are sensitive to NR2B subunit antagonists (Kolarova et al., 2003; Butler et al., 2010).

In our present work, we have not adjusted our model for micro-anatomical restructuring at the tripartite synapse, which would almost inevitably occur due to the dynamics at astrocytic filopodia/perisynaptic processes (Hirrlinger et al., 2004; Lavialle et al., 2011). Such restructuring could lead to changes in the proximity of neuronal and astrocytic release and/or input sites. Future dynamic modeling to simulate micro-anatomical arrangements would be needed to computationally address the role that astrocytic structural plasticity may play in tuning neuronal firing rhythm.

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.

Vladimir Parpura acknowledges the support of this work by a National Science Foundation Award CBET 0943343; he thanks Manoj K. Gottipati for constructive comments on the manuscript. Shivendra Tewari acknowledges useful comments on a previous version of the article by Drs. Kaushik Majumdar and Catherine Kaczorowski.

 Supplementary material

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

 References Andersen P. Morris R. Amaral D. Bliss T. O'keefe J. (2007). The Hippocampus Book. New York, NY: Oxford University Press. Auzmendi J. Gonzalez N. Girardi E. (2009). The NMDAR subunit NR2B expression is modified in hippocampus after repetitive seizures. Neurochem. Res. 34, 819–826. 10.1007/s11064-008-9828-018751892 Bushong E. A. Martone M. E. Ellisman M. H. (2004). Maturation of astrocyte morphology and the establishment of astrocyte domains during postnatal hippocampal development. Int. J. Dev. Neurosci. 22, 73–86. 10.1016/j.ijdevneu.2003.12.00815036382 Bushong E. A. Martone M. E. Jones Y. Z. Ellisman M. H. (2002). Protoplasmic astrocytes in CA1 stratum radiatum occupy separate anatomical domains. J. Neurosci. 22, 183–192. 11756501 Butler T. R. Self R. L. Smith K. J. Sharrett-Field L. J. Berry J. N. Littleton J. M. . (2010). Selective vulnerability of hippocampal cornu ammonis 1 pyramidal cells to excitotoxic insult is associated with the expression of polyamine-sensitive N-methyl-D-asparate-type glutamate receptors. Neuroscience 165, 525–534. 10.1016/j.neuroscience.2009.10.01819837138 Corkin S. (2002). What's new with the amnesic patient H.M. Nat. Rev. Neurosci. 3, 153–160. 10.1038/nrn72611836523 Danbolt N. C. (2001). Glutamate uptake. Prog. Neurobiol. 65, 1–105. 10.1016/S0301-0082(00)00067-811369436 Destexhe A. Mainen Z. F. Sejnowski T. J. (1994). Synthesis of models for excitable membranes, synaptic transmission and neuromodulation using a common kinetic formalism. J. Comput. Neurosci. 1, 195–230. 10.1007/BF009617348792231 Ding S. Fellin T. Zhu Y. Lee S. Y. Auberson Y. P. Meaney D. F. . (2007). Enhanced astrocytic Ca2+ signals contribute to neuronal excitotoxicity after status epilepticus. J. Neurosci. 27, 10674–10684. 10.1523/JNEUROSCI.2001-07.200717913901 Eldridge L. L. Knowlton B. J. Furmanski C. S. Bookheimer S. Y. Engel S. A. (2000). Remembering episodes: a selective role for the hippocampus during retrieval. Nat. Neurosci. 3, 1149–1152. 10.1038/8067111036273 Fellin T. Pascual O. Gobbo S. Pozzan T. Haydon P. G. Carmignoto G. (2004). Neuronal synchrony mediated by astrocytic glutamate through activation of extrasynaptic NMDA receptors. Neuron 43, 729–743. 10.1016/j.neuron.2004.08.01115339653 Fellin T. Pascual O. Haydon P. G. (2006). Astrocytes coordinate synaptic networks: balanced excitation and inhibition. Physiology 21, 208–215. 10.1152/physiol.00161.200516714479 Gasparini S. Migliore M. Magee J. C. (2004). On the initiation and propagation of dendritic spikes in CA1 pyramidal neurons. J. Neurosci. 24, 11046–11056. 10.1523/JNEUROSCI.2520-04.200415590921 Gold C. Henze D. A. Koch C. Buzsaki G. (2006). On the origin of the extracellular action potential waveform: a modeling study. J. Neurophysiol. 95, 3113–3128. 10.1152/jn.00979.200516467426 Grover L. M. Kim E. Cooke J. D. Holmes W. R. (2009). LTP in hippocampal area CA1 is induced by burst stimulation over a broad frequency range centered around delta. Learn. Mem. 16, 69–81. 10.1101/lm.117910919144965 Hall J. Thomas K. L. Everitt B. J. (2001). Cellular imaging of zif268 expression in the hippocampus and amygdala during contextual and cued fear memory retrieval: selective activation of hippocampal CA1 neurons during the recall of contextual memories. J. Neurosci. 21, 2186–2193. 11245703 Henneberger C. Papouin T. Oliet S. H. Rusakov D. A. (2010). Long-term potentiation depends on release of D-serine from astrocytes. Nature 463, 232–236. 10.1038/nature0867320075918 Hirrlinger J. Hülsmann S. Kirchhoff F. (2004). Astroglial processes show spontaneous motility at active synaptic terminals in situ. Eur. J. Neurosci. 20, 2235–2239. 10.1111/j.1460-9568.2004.03689.x15450103 Holt W. Maren S. (1999). Muscimol inactivation of the dorsal hippocampus impairs contextual retrieval of fear memory. J. Neurosci. 19, 9054–9062. 10516322 Kang N. Xu J. Xu Q. Nedergaard M. Kang J. (2005). Astrocytic glutamate release-induced transient depolarization and epileptiform discharges in hippocampal CA1 pyramidal neurons. J. Neurophysiol. 94, 4121–4130. 10.1152/jn.00448.200516162834 Khuchua Z. Wozniak D. F. Bardgett M. E. Yue Z. McDonald M. Boero J. (2003). Deletion of the N-terminus of Murine MAP2 by gene targeting disrupts Hippocampal CA1 neuron architecture and alters Contextual Memory. Neuroscience 119, 101–111. 10.1016/S0306-4522(03)00094-012763072 Knyazev G. G. (2012). EEG delta oscillations as a correlate of basic homeostatic and motivational processes. Neurosci. Biobehav. Rev. 36, 677–695. 10.1016/j.neubiorev.2011.10.00222020231 Kolarova A. Ringer R. Tauber M. G. Leib S. L. (2003). Blockade of NMDA receptor subtype NR2B prevents seizures but not apoptosis of dentate gyrus neurons in bacterial meningitis in infant rats. BMC Neurosci. 4:21. 10.1186/1471-2202-4-2113129439 Lakatos P. Karmos G. Mehta A. D. Ulbert I. Schroeder C. E. (2008). Entrainment of neuronal oscillations as a mechanism of attentional selection. Science 320, 110–113. 10.1126/science.115473518388295 Lavialle M. Aumann G. Anlauf E. Prols F. Arpin M. Derouiche A. (2011). Structural plasticity of perisynaptic astrocyte processes involves ezrin and metabotropic glutamate receptors. Proc. Natl. Acad. Sci. U.S.A. 108, 12915–12919. 10.1073/pnas.110095710821753079 Lee I. Kesner R. P. (2004). Differential contributions of dorsal hippocampal subregions to memory acquisition and retrieval in contextual fear-conditioning. Hippocampus 14, 301–310. 10.1002/hipo.1017715132429 Liao G. Y. Wagner D. A. Hsu M. H. Leonard J. P. (2001). Evidence for direct protein kinase-C mediated modulation of N-methyl-D-aspartate receptor current. Mol. Pharmacol. 59, 960–964. 11306676 Lopez-Aguado L. Ibarz J. M. Herreras O. (2001). Activity-dependent changes of tissue resistivity in the CA1 region in vivo are layer-specific: modulation of evoked potentials. Neuroscience 108, 249–262. 10.1016/S0306-4522(01)00417-111734358 Maren S. (2001). Neurobiology of pavlovian fear conditioning. Annu. Rev. Neurosci. 24, 897–931. 10.1146/annurev.neuro.24.1.89711520922 Maren S. Holt W. (2000). The hippocampus and contextual memory retrieval in Pavlovian conditioning. Behav. Brain Res. 110, 97–108. 10.1016/S0166-4328(99)00188-610802307 Martin E. D. Fernandez M. Perea G. Pascual O. Haydon P. G. Araque A. . (2007). Adenosine released by astrocytes contributes to hypoxia-induced modulation of synaptic transmission. Glia 55, 36–45. 10.1002/glia.2043117004232 Milner B. Penfield W. (1955). The effect of hippocampal lesions on recent memory. Trans. Am. Neurol. Assoc. 80th mtg, 42–48. 13311995 Murchison C. F. Schutsky K. Jin S.-H. Thomas S. A. (2011). Norepinephrine and ß1-adrenergic signaling facilitate activation of hippocampal CA1 pyramidal neurons during contextual memory retrieval. Neuroscience 181, 109–116. 10.1016/j.neuroscience.2011.02.04921377513 Nadkarni S. Jung P. (2003). Spontaneous oscillations of dressed neurons: a new mechanism for epilepsy. Phys. Rev. Lett. 91, 268101. 10.1103/PhysRevLett.91.26810114754091 Nakazawa K. McHugh T. J. Wilson M. A. Tonegawa S. (2004). NMDA receptors, place cells and hippocampal spatial memory. Nat. Rev. Neurosci. 5, 361–372. 10.1038/nrn138515100719 Nedergaard M. Ransom B. Goldman S. A. (2003). New roles for astrocytes: redefining the functional architecture of the brain. Trends Neurosci. 26, 523–530. 10.1016/j.tins.2003.08.00814522144 Oberheim N. A. Takano T. Han X. He W. Lin J. H. Wang F. . (2009). Uniquely hominid features of adult human astrocytes. J. Neurosci. 29, 3276–3287. 10.1523/JNEUROSCI.4707-08.200919279265 Ogata K. Kosaka T. (2002). Structural and quantitative analysis of astrocytes in the mouse hippocampus. Neuroscience 113, 221–233. 10.1016/S0306-4522(02)00041-612123700 Omkumar R. V. Kiely M. J. Rosenstein A. J. Min K. T. Kennedy M. B. (1996). Identification of a phosphorylation site for calcium/calmodulin dependent protein kinase II in the NR2B subunit of the N-methyl-D-aspartate receptor. J. Biol. Chem. 271, 31670–31678. 10.1074/jbc.271.49.316708940188 Paller K. A. McCarthy G. (2002). Field potentials in the human hippocampus during the encoding and recognition of visual stimuli. Hippocampus 12, 415–420. 10.1002/hipo.1005312099492 Parpura V. Haydon P. G. (2000). Physiological astrocytic calcium levels stimulate glutamate release to modulate adjacent neurons. Proc. Natl. Acad. Sci. U.S.A. 97, 8629–8634. 10.1073/pnas.97.15.862910900020 Perea G. Araque A. (2005). Properties of synaptically evoked astrocyte calcium signal reveal synaptic information processing by astrocytes. J. Neurosci. 25, 2192–2203. 10.1523/JNEUROSCI.3965-04.200515745945 Ruch S. Markes O. Duss S. B. Oppliger D. Reber T. P. Koenig T. . (2012). Sleep stage II contributes to the consolidation of declarative memories. Neuropsychologia 50, 2389–2396. 10.1016/j.neuropsychologia.2012.06.00822750121 Scoville W. B. Milner B. (1957). Loss of recent memory after bilateral hippocampal lesions. J. Neurol. Neurosurg. Psychiatry 20, 11–21. 10.1136/jnnp.20.1.1113406589 Siekmeier P. J. Hasselmo M. E. Howard M. W. Coyle J. (2007). Modeling of context-dependent retrieval in hippocampal region CA1: implications for cognitive function in schizophrenia. Schizophr. Res. 89, 177–190. 10.1016/j.schres.2006.08.00717055702 Small S. A. Nava A. S. Perera G. M. Delapaz R. Mayeux R. (2001). Circuit mechanisms underlying memory encoding and retrieval in the long axis of the hippocampal formation. Nat. Neurosci. 4, 442–449. 10.1038/8611511276237 Takahashi S. Sakurai Y. (2009). Sub-millisecond fi ring synchrony of closely neighboring pyramidal neurons in hippocampal CA1 of rats during delayed non-matching to sample task. Front. Neural Circuits 3:9. 10.3389/neuro.04.009.200919753324 Tewari S. Majumdar K. (2012a). A mathematical model for astrocytes mediated LTP at single hippocampal synapses. J. Comput. Neurosci. 33, 341–370. 10.1007/s10827-012-0389-522454034 Tewari S. G. Majumdar K. K. (2012b). A mathematical model of the tripartite synapse: astrocyte-induced synaptic plasticity. J. Biol. Phys. 38, 465–496. 10.1007/s10867-012-9267-723729909 Tian G. F. Azmi H. Takano T. Xu Q. Peng W. Lin J. . (2005). An astrocytic basis of epilepsy. Nat. Med. 11, 973–981. 10.1038/nm127716116433 Traub R. D. Jefferys J. G. Miles R. Whittington M. A. Toth K. (1994). A branching dendritic model of a rodent CA3 pyramidal neurone. J. Physiol. 481(Pt 1), 79–95. 7853251 Traub R. D. Wong R. K. Miles R. Michelson H. (1991). A model of a CA3 hippocampal pyramidal neuron incorporating voltage-clamp data on intrinsic conductances. J. Neurophysiol. 66, 635–650. 1663538 Zhang X.-H. Wu L.-J. Gong B. Ren M. Li B.-M. Zhuo M. (2008). Induction- and conditioning-protocol dependent involvement of NR2B-containing NMDA receptors in synaptic potentiation and contextual fear memory in the hippocampal CA1 region of rats. Mol. Brain 1:9. 10.1186/1756-6606-1-918826591 Zhao M.-G. Toyoda H. Lee Y.-S. Wu L.-J. Ko S. W. Zhang X.-H. . (2005). Roles of NMDA NR2B subtype receptor in prefrontal long-term potentiation and contextual fear memory. Neuron 47, 859–872. 10.1016/j.neuron.2005.08.01416157280 Zhuo M. (2009). Plasticity of NMDA receptor NR2B subunit in memory and chronic pain. Mol. Brain 2, 4. 10.1186/1756-6606-2-419192303