Abstract
NMDA receptors (NMDARs) may be crucial to working memory (WM). Computational models predict that they sustain neural firing and produce associative memory, which may underpin maintaining and binding information, respectively. We test this in patients with antibodies to NMDAR (n = 10, female) and compare them with healthy control participants (n = 55, 20 male, 35 female). Patients were tested after recovery with a task that separates two aspects of WM: sustaining attention and feature binding. Participants had to remember two colored arrows. Then attention was directed to one of them. After a variable delay, they reported the direction of either the same arrow (congruent cue) or of the other arrow (incongruent cue). We asked how congruency affected recall precision and measured types of error. Patients had difficulty in both sustaining attention to an item over time and feature binding. Controls were less precise after longer delays and incongruent cues. In contrast, patients did not benefit from congruent cues at longer delays [group × congruency (long condition); p = 0.041], indicating they could not sustain attention. Additionally, patients reported the wrong item (misbinding errors) more than controls after congruent cues [group × delay (congruent condition), main effect of group; p ≤ 0.001]. Our results suggest NMDARs are critical for both maintaining attention and feature binding.
Significance Statement
Computational theories suggest NMDA receptors (NMDARs) are critical for actively maintaining information, while other theories propose they allow us to associate or “bind” objects features together. This is the first causal test in humans of the role of NMDARs in actively maintaining attention in working memory and feature binding. We found patients have difficulty with both these processes in support of computational models. Notably, we demonstrate that patients with NMDA receptor-antibody encephalitis are an ideal model condition to study roles of receptors in human cognition. Secondly, few studies follow these patients long after treatment. Our findings demonstrate a specific long-term neuropsychological deficit, previously unreported to our knowledge, that highlights the need for greater focus on neurocognitive rehabilitation with these patients.
Introduction
Working memory (WM) is a limited capacity system for temporary storage of information used for processing (Baddeley, 2003). The neural mechanism of WM storage and the role of attention remains unclear (Oberauer, 2019) but animal and human studies suggest N-methyl-D-aspartate receptors (NMDARs) may be essential for encoding and storage (M. Wang et al., 2013; Koychev et al., 2017).
In WM, an item in the focus of attention is maintained in a prioritized state such that it is more readily available and more accurately maintained (Zokaei et al., 2014a). Prioritization can be studied using a retrospective cue (“retrocue”) task. Here, participants see an array of items to remember. A cue is presented during the retention interval, indicating for example the location or color of one of the items, which guides the focus of attention to that item within WM. Subsequently, they are asked to recall a feature of one item. Recall of cued items is typically more accurate and faster than uncued items (Souza and Oberauer, 2016).
There may be different reasons why retrocues improve performance. Most commonly, cues indicate the probability an item will be probed (Griffin and Nobre, 2003). However, benefits are seen even if the cue causes an item to be merely accessed during the delay (Zokaei et al., 2014b). This is termed an “incidental” retrocue, as it provides no information about which item will be tested. This allows us to measure the benefit of shifting the focus of attention. What drives improvements? To answer this, we can ask what kinds of error people make on this task. Participants can be asked to reconstruct a feature from memory using a continuous scale, allowing response errors to be modeled (Bays et al., 2009; Pertzov et al., 2017) and inform how representations in WM are corrupted.
One proposed mechanism of WM retention involves NMDARs (Murray et al., 2014). Computational models of WM using artificial spiking neural networks demonstrated that without NMDARs, activity representing a stimulus dies away within 200 ms, but rapid transient increases in synaptic strength could support sustained activation during a memory delay (Compte et al., 2000; M. Wang et al., 2013). Simulations with and without NMDARs, with a decay time of 80 ms, demonstrate that NMDARs permit sustained activity. Similarly, models that test the ability of networks to bind features together into objects demonstrate that in rate-coded neurons, removing rapid Hebbian plasticity prevents the network from remembering associations but preserves feature memory (Manohar et al., 2019). Spiking versions of these networks demonstrate NMDARs could underlie this computational property (Fiebig and Lansner, 2017; Fiebig et al., 2020). In summary, computational models have demonstrated both feature binding and active maintenance over time may be crucially underpinned by NMDARs.
In patients with NMDAR-antibody encephalitis (NMDAR-Ab-E), autoantibodies found in patient CSF target the NR1 subunit of the NMDAR (Dalmau et al., 2008). It largely affects young women but is a disease of all ages and both sexes (Irani et al., 2010). Patients present with prominent mixed psychiatric symptoms, amnesia, and a movement disorder, among other features (Dalmau et al., 2008; Irani et al., 2010; Al-Diwani et al., 2019). Longer-term sequelae such as reduced executive function and memory impairments are likely influenced by time to diagnosis and treatment aggression but can persist years posttreatment (McKeon et al., 2018). NMDAR IgG antibodies may induce receptor endocytosis (Hughes et al., 2010) and selectively reduce postsynaptic available NMDARs (Jézéquel et al., 2017; Galovic et al., 2023). Accordingly, reduced NMDARs have been measured both in postmortem brain tissue and living patients with lower open NMDAR channel density on PET imaging, even long into recovery (Galovic et al., 2023). Hence, patients with NMDAR-Ab-E provide an ideal human model to ask whether selective NMDAR modulation impacts WM and attention (Finke et al., 2012; McKeon et al., 2018; Stein et al., 2020). In this study, we explored shifts of attention within WM using a task with an incidental cue and a variable delay. We measure effects of cueing and delay on memory precision and binding.
Materials and Methods
Participants
We recruited 58 healthy volunteers (HC, healthy controls) and 10 patients with NMDAR-Ab-E who provided written consent under ethics approved by the National Research Ethics Service (18/LO/2152). The mean number of years from acute presentation and diagnosis for patients was 4.02. They had all been successfully treated at initial presentation with first line immunotherapy (steroids, intravenous immunoglobulin, and plasma exchange). Three patients had required second-line immunotherapy much later (rituximab and/or cyclophosphamide), two of whom necessitated this for a relapse. Data on specific treatment was unavailable for two patients. Healthy people were recruited via an online sign-up system or a participant database, and patients were recruited from the Oxford Autoimmune Neurology Clinic. We tested 32 HCs and five patients in-person, and following COVID-19 restrictions, 26 HCs and five patients were tested remotely via an adapted online platform participants could access on their own devices. The in-person testing took place individually in a quiet room. Stimuli were displayed on a 20 inch computer monitor at an approximate viewing distance of 70 cm. To screen for mood, behavioral, and major cognitive confounders, each participant answered the HADS (Hospital Anxiety and Depression Scale; Zigmond and Snaith, 1983) and BIS-11 (Barratt Impulsiveness Scale; Patton et al., 1995) questionnaires and performed a digit span test. Both online and in-person sessions lasted ∼1.5 h, and participants were compensated for their time (£10/h).
Experimental protocol
In-person experiments were programmed using MATLAB R2018b and the Psychophysics Toolbox extension (Brainard, 1997; Pelli, 1997). Online testing was programmed with PsychoPy hosted by Pavlovia (Open Science Tools; Peirce et al., 2019). Online surveys and consent forms were generated using Qualtrics software version August 2021 of Qualtrics (Copyright © 2023 Qualtrics, Qualtrics and all other Qualtrics product or service names are registered trademarks or trademarks of Qualtrics; https://www.qualtrics.com).
Trials commenced with a display of a central fixation cross for 500 ms (Fig. 1A). After this, participants saw a memory array (1,000 ms) of two arrows (each 56 pixels long, 250 pixels either side of a white fixation cross on a black background), randomly oriented and differently colored. The colors were selected randomly from eight distinct colors (red, green, blue, yellow, cyan, magenta, orange, and white). After a 500 ms delay, participants saw the name of one of these colors (incidental retrocue, Ariel, 32 pt., text color matched name) and clicked the left or right mouse key to report which side of the screen this arrow appeared on. After a delay of 1,000 or 3,000 ms (50% probability), they were probed with either the cued arrow (congruent condition) or the uncued arrow (incongruent condition), again with 50% probability. Therefore, the retrocue was uninformative about which arrow would be probed (target arrow). The probe stimulus, an arrow, could not be the same orientation as the target arrow and always differed by >20° from the true orientation, so that adjustment was always needed. Participants responded by using the mouse to orient the probe to match the target arrow.
Participants performed four slowed practice trials, where they received feedback about their precision, and then three blocks of 84 trials of normal delay schedule as described above (63 trials per condition). The first three controls tested in person performed four blocks of 100 trials. In the online format, participants completed four blocks of 80 trials (80 trials per condition).
Analysis
We employed a probabilistic three-factor mixture model (Bays et al., 2009; Fougnie et al., 2012) to characterize the source of errors made by participants. In this model, errors can arise from the following: imprecision (variability in memory of arrow orientation, standard deviation of the von Mises distribution), misbinding (misreporting errors where responses are centered on the nontarget arrow's orientation), or random guessing (a fixed probability; Fig. 1B). This can be mathematically described as Pertzov et al. (2017):
To analyze results, we used 2 × 2 × 2 mixed repeated-measures ANOVAs (patients vs HC). Where there were significant main effects or interactions of group, 2 × 2 repeated-measures ANOVAs were conducted. Statistical analyses were calculated with SPSS (v.25). Normality of data was analyzed with the Shapiro–Wilk (SW) test and homogeneity of variance with Levene's test to assess the assumptions of the ANOVA. Misbinding, guess, and target data were arcsine transformed for statistical analyses to meet the assumptions of normality for ANOVA (please see Table 1 for Shapiro–Wilk results). Data and scripts are available at OSF (https://osf.io/rcb3m/).
Results
Trials with an incorrect retrocue response were excluded as well as trials where the time taken to initially move the mouse was >5 s, the time taken to position the arrow after first moving the mouse was <0.2 s, or total reaction time was >10 s (data removed: HC in-person = 3.74%, HC online 4.71%; NMDA in-person = 3.02%, NMDA online = 1.44%); the proportion of trials removed per group did not differ (Table 2). Participants were excluded entirely if they had retrocue accuracy lower than chance (N = 1 in-person in HC) or if the estimated precision was >2SD greater than the group mean (N = 2 online HC). One patient misunderstood instructions and repeated the study; only their second session was analyzed. Online and in-person mean data was not significantly different (Table 3) therefore, we combined the data in our analyses.
For overall error (raw error), defined as the absolute angular difference in degrees between the reported orientation and the actual orientation of the target arrow, patients were worse overall with no interactions with delay or congruency (Table 4). Reaction times, measured from the start of probe presentation on screen to the start of mouse movement, were correspondingly longer in patients (main effect of group, p = 0.004; Fig. 2, Table 5). Congruence sped up reaction times (main effect of congruency p < 0.001), but patients benefited less from congruency specifically when the delay was long (congruency × delay × group; p = 0.016).
We applied a probabilistic mixture model to the orientation recall responses (Fig. 1B) to analyze the contribution of different error sources to responses (Bays et al., 2009). The model yields parameters for three different kinds of error: (1) imprecision indicates the overall degradation in the quality of a memory (standard deviation calculated as 1/κ of the von Mises distribution, in degrees, around the target); (2) misbinding is the probability of reporting the orientation of the wrong item; (3) guessing is the probability of recall being completely at chance, i.e., a uniformly distributed response.
Data and scripts are available at OSF (https://osf.io/rcb3m/). Demographics and questionnaire data are presented in Table 6.
Recovered NMDAR-antibody encephalitis patients show paradoxical congruency effect with long delays
A 2 × 2 × 2 ANOVA was performed on imprecision (1/κ), comparing HCs and patients (Fig. 3A). There was a three-way group × delay × congruency interaction (p = 0.016; Table 7) which showed that while HCs showed the expected benefit for congruent cues and short delays, patients demonstrated a different congruency effect at long delays. To understand this three-way interaction, we examined the congruency effect for long and short delay trials separately (2 × 2 ANOVAs of group × congruency; Table 7). At short delays, there was the expected congruency benefit [main effect of congruency (short condition); p = 0.003] with no group differences. At long delays, however, the congruency benefit was greater in HCs than in the patient groups, who may even show a congruency cost [group × congruency (long condition); p = 0.041]. In other words, patients were unable to sustain the attentional benefit at longer delays. To quantify the degree to which these results were driven by individual patients, we plotted absolute error by trial type for each patient (Fig. 3B). More than half the patients displayed the paradoxical improved precision with incongruency in the long condition compared with 33% HCs. Figure 3C illustrates the negative congruency effect for patients in the long condition compared with the short condition (cluster in the right lower quadrant) while HCs have a positive congruency effect for both delays (for F statistics and effect sizes, see Table 7).
Patients misbind more with congruency
Next, we asked what proportion of errors were due to reporting the wrong arrow (swap errors or “misbinding”). While HCs again show a standard congruency benefit, with less misbinding with congruent cues and with longer delay, unusually patients misbind more in the congruent condition rather than the incongruent condition [Fig. 4A; (group × congruency × delay ANOVA), congruency × group; p = 0.014]. We asked whether this was because patients misbound more than controls on congruent trials, or less on incongruent trials, by performing separate group × delay ANOVAs for the congruent and incongruent trials. Patients misbound more than controls on congruent trials, regardless of delay [group × delay (congruent condition); main effect of group, p ≤ 0.001], demonstrating again a congruency cost. On incongruent trials, there was no difference in misbinding between patients and HCs (Table 8).
Patients guess more than controls and have fewer target responses
Guessing errors occur when participants select responses randomly from a uniform distribution. Patients guessed more overall than HCs (group main effect, p = 0.031), and both groups guessed more in the incongruent condition (congruency main effect, p ≤ 0.001). There was also a congruency × delay interaction (p = 0.040) but no group interactions (Tables 9, 10).
Discussion
NMDA receptors are important in WM storage (Lisman et al., 1998), but few causal studies in humans have asked what component mechanisms they support. Patients with NMDAR-antibody encephalitis have a selective reduction in postsynaptic NMDAR availability, including long into recovery (Galovic et al., 2023). We asked whether sustaining attention to items in WM and feature binding are affected by NMDAR hypofunction using a task that allowed modeling of error types. A striking finding was that in the long-delay condition, patients lacked the typical congruency benefit. In other words, patients had difficulty sustaining attention to prioritized items. Furthermore, patients make more misbinding errors than controls in congruent trials, suggesting they may struggle to remember associations, especially for prioritized information.
In a recent study of one-item WM, patients lacked interference from previous trials, which was explained by modeling in terms of diminished short-term potentiation caused by NMDAR hypofunction (Stein et al., 2020). Previous studies illustrated that incorporating prefrontal synaptic disinhibition (NMDAR hypofunction decreasing activation of GABA interneurons) into spatial WM modeling leads to reduced precision and greater distraction from near distractors (Murray et al., 2014). This may be one explanation for patients’ misbinding errors, as similar representations may be more difficult to differentiate. NMDAR antibodies lower receptor-mediated potentiation and currents (Hughes et al., 2010; X. Wang et al., 2019) which is evident in different cortical areas including the hippocampus (X. Wang et al., 2019) a region of high NMDAR density and a region often linked to memory (McHugh et al., 2007) particularly associative feature binding (Zokaei et al., 2019). It is possible that patients misbind more in the congruent condition as aberrant NMDAR-mediated synaptic function involved in paying attention to an item could make the representation more difficult to differentiate from the other item.
An important point to note is the possibility of potential hippocampal damage in patients which is difficult to completely distinguish from the direct effects of antibodies on NMDARs. There is a predominantly high expression of NMDARs in the hippocampus (Monaghan and Cotman, 1985) and, while we are interested in a receptor level deficit, it is likely that the deficit is partially localized to the hippocampus. MRI evidence suggests only a minority of patients may have long-term microstructural damage to the hippocampus (Finke et al., 2016: imaging mean 26.6 months after disease onset). Consistent with this, in our patient group the most recent MRI imaging of 8/10 patients posttreatment (average 26.1 months postdiagnosis) found the following: 6/10 had largely normal hippocampi; 2/10 had unreported scans—on assessment by the authors; one showed no hippocampal atrophy; and the other demonstrated mild hippocampal atrophy. Two out of 10 had no MRI imaging available. This indicates most of our patients had no substantial residual structural hippocampal damage following the acute phase of the disease.
Hebbian models of WM implicate NMDARs in binding features together to form neural “conjunctions,” potentially in the prefrontal cortex (Manohar et al., 2019). Sustained activity of neurons encoding these conjunctions may constitute attention to a specific item and is also dependent on NMDARs. Within such models, switching attention away from an item disrupts this persistent firing, but the synaptic weight changes remain. This allows unattended information to be reactivated, shifting attention back to that item. In our study, patients lost the congruency effect at long delays, with even decreased precision for the attended item. This indicates patients struggle to sustain attention on prioritized items over time. In the model this would translate to faster decay of persistent firing in conjunctive neurons. It is interesting that patients do not seem to show a decay in precision over time in the incongruent condition. One possibility is that their attention is more likely to switch between items during the delay. Alternatively, decaying persistent firing may weaken feature associations, such that the established synaptic weights are not strong enough to maintain the features of a representation (the features are bound too homogenously). As such, over time the attended representation is less precisely maintained, but the unattended representation is less affected. This may also account for patient misbinding errors in the congruent condition. If the attended item has both faster decay of persistent firing and weaker synaptic weights over time compared with the unattended item, this could lead to more overlap and interference with other features, rendering pattern separation more difficult to achieve and misbinding more likely to occur. A rate-coded model simulating this task has been previously described and demonstrates the need for Hebbian plasticity for both sustained activity and binding (Manohar et al., 2019). To directly simulate the role of NMDARs, the model would need to simulate spiking neurons. Spiking simulations have indeed shown the importance of NMDARs but have not simulated this task specifically (Fiebig et al., 2020). Future simulation work could make more quantitative predictions about how receptor density affects neural activity, memory recall, the effect of delay, and reaction times.
NMDARs are ubiquitous within the brain and this is reflected in the complex neuropsychiatric presentation of the disease, with affective, psychotic, motor and cognitive components. Similarly, WM models require NMDA plasticity in multiple locations and interaction across multiple brain areas. Neural models suggest that NMDAR-related plasticity may be required in both the cortical areas that maintain features (e.g., parietal cortex) and also cortical areas that hold the binding “pointers” (e.g., hippocampus or prefrontal cortex) (Manohar et al., 2019). This makes NMDAR-Ab-E a particularly useful model as it offers a view of specificity in mechanism at the level of cellular computations, even for cognitive functions that, like WM, appear to have little specificity to brain area. This study therefore complements previous work in patients with spatially localized lesions (Aggleton et al., 1986; Pertzov et al., 2013; Yonelinas, 2013; Allen et al., 2014).
The cohort of 10 patients we tested were mostly young women (mean age, 29.2), who had recovered from the acute illness and with good functional outcomes including return to work or full-time education. Some did report residual cognitive difficulty including anxiety, poor sleep, word finding difficulties, and “slowed reading and writing.” Potential confounding factors include the use of psychotropic and anticonvulsant medications. Two patients had been diagnosed with comorbid psychological illnesses which were being treated; two patients reported taking levetiracetam, one sertraline and one benzodiazepine. Medications or untreated symptoms may have interfered with our results and affected performance. Furthermore, some HCs were taking psychotropic medications, two of whom confirmed they were taking antidepressants. Four patients had suffered two clinical episodes of depression, one of whom was diagnosed during an encephalitis relapse. The main limitations of our study are the small sample size and limited diversity of our cohort. This is because NMDAR-E is a rare disease which more commonly occurs in young women though it is a disease of all ages and both sexes. It would therefore be beneficial to replicate the results in a larger and more diverse group of patients. Finally, a few of our patients were not clinically diagnosed with NMDAR-Ab-E until months after symptoms initially started, making the exact time of disease onset difficult to determine. Here, years from diagnosis was based on date of antibody diagnosis where the data was available, and where unavailable, based on patient reported date of diagnosis.
For 7/10 of our patients, serum NMDAR-IgG titers were available [3/7 at diagnosis (1:320, 1:320, 1:160); 4/7 from later timepoints (average of 625 d from diagnosis—1:80 d424; 1:20 d275; 1:40 d1736; 1:500 d74)]. While the presence of current antibodies could theoretically impact function, PET imaging evidence demonstrates that receptor hypofunction continues posttreatment (Galovic et al., 2023) and serum titers do not correlate perfectly with disease course (Gresa-Arribas et al., 2014). Furthermore, the effect we were particularly interested in examining here was not specifically based on the presence of antibodies in these patients but on the NMDA receptor hypofunction induced by the antibodies and the effects of this receptor deficit on WM. Therefore, antibody status/titer or their treatment-related improvements are not directly relevant to our specific objectives.
It should be highlighted this is a disease where treatment largely focuses on acute control with immunomodulation, while longer-term neurological rehabilitation to compensate for NMDAR abnormalities is hugely underexplored. These findings suggest a specific long-term neuropsychological deficit in this patient group which has not previously been reported in any other studies to our knowledge. Our results indicate patients experience WM deficits years postdiagnosis (4.02 mean years since diagnosis) highlighting an important focus for future research. Targeting NMDAR after initial recovery may provide benefit for longer-term neuropsychological sequelae. Agents could include positive allosteric modulators such as D-cycloserine which has previously been trialed with schizophrenia though effects have been poor thus far (Kuppili et al., 2021). Targeting associated proteins such as EphB2 receptors could also prove beneficial. These receptors are known to modulate NMDARs, and evidence has demonstrated that autoantibody-mediated NMDAR loss is prevented by EphB2 receptor activation (Mikasova et al., 2012).
Furthermore, poor attention and WM are important subjective findings but can be difficult to quantify especially in the context of a patient attending a medical clinic. Current clinical tools for assessing NMDAR function are minimal. There may be a clinical utility in this retrocue task, easily completed on a computer, as a neuropsychological assay to help quantify recovery posttreatment in both a clinical setting and outcome studies particularly when assessing symptoms >12 months posttreatment where specific neuropsychological tasks may prove more sensitive (Guasp et al., 2022).
Conclusion
Our data indicate that patients with anti-NMDA-R encephalitis not only have difficulty sustaining attention on prioritized items in WM but also struggle with feature binding of attended items. This implicates a role for NMDAR in maintaining attention on information in WM and in feature binding, in line with some computational models of WM. Moreover, these results highlight an important long-term neuropsychological deficit in this group of patients and suggest greater treatment focus should be directed toward neuropsychological rehabilitation.
Footnotes
S.M.— Medical Research Council (MRC) clinician scientist fellowship MR/P00878/X, Leverhulme RPG2018-310. A.A.-D.—British Medical Association (BMA) Foundation for Medical Research Margaret Temple Grant 2017. National Institute for Health Research (NIHR) Oxford Biomedical Research Centre (BRC) and NIHR Oxford Health BRC. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR, or the Department of Health. S.R.I. and this research were funded in whole or in part by a senior clinical fellowship from the MRC (MR/V007173/1), Wellcome Trust Fellowship (104079/Z/14/Z); BMA Research Grants, Vera Down Grant (2013) and Margaret Temple (2017); Epilepsy Research UK (P1201); the Fulbright UK-US commission (MS-Society research award); and by the NIHR BRC. For the purpose of Open Access, the author has applied a CC BY public copyright license to any Author Accepted Manuscript (AAM) version arising from this submission. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR, or the Department of Health.
A.D., S.M., A.A.-D., J.G., C.H.—none declared. S.R.I. has received honoraria/research support from UCB, Immunovant, MedImmune, Roche, Janssen, Cerebral Therapeutics, ADC Therapeutics, Brain, CSL Behring, and ONO Pharma. S.R.I. is an inventor on “Diagnostic Strategy to improve specificity of CASPR2 antibody detection” (PCT/G82019/051257) and receives royalties on a licensed patent application for LGI1/CASPR2 testing as coapplicant (PCT/GB2009/051441) entitled “Neurological Autoimmune Disorders”.
- Correspondence should be addressed to Afrose Dor at afrose.dor{at}gmail.com or Sanjay Manohar at sanjay.manohar{at}ndcn.ox.ac.uk.