Excitatory and inhibitory D-serine binding to the NMDA receptor

  1. Remy A Yovanno
  2. Tsung Han Chou
  3. Sarah J Brantley
  4. Hiro Furukawa
  5. Albert Y Lau  Is a corresponding author
  1. Department of Biophysics and Biophysical Chemistry, Johns Hopkins University School of Medicine, United States
  2. W.M. Keck Structural Biology Laboratory, Cold Spring Harbor Laboratory, United States
  3. Department of Biology, Johns Hopkins University, United States
7 figures, 2 videos and 1 additional file

Figures

Figure 1 with 4 supplements
Identifying D-serine binding pathways for GluN2A using pathway similarity analysis (PSA).

(A) Overview of the PSA workflow for quantifying similarity between D-serine binding pathways. (B) 2-dimensional order parameter (ξ1,ξ2) that describes the degree of GluN2A LBD closure. For each of the …

Figure 1—source data 1

Simulation summary: overview of simulation systems.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig1-data1-v1.xlsx
Figure 1—source data 2

Record of all successful binding pathways in each simulation system for D-serine binding to GluN2A.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig1-data2-v1.xlsx
Figure 1—source data 3

Per-residue contact frequency analysis for D-serine binding to GluN2A by cluster identified with PSA.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig1-data3-v1.xlsx
Figure 1—source data 4

GluN2A residues most frequently contacted by D-serine given that the pathway results in successful binding – listed for each simulation system.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig1-data4-v1.xlsx
Figure 1—figure supplement 1
Overlap coefficient analysis for GluN2A and GluN1 binding pathways.

This figure supplement relates to Figures 1 and 2. The overlap coefficient was computed for each pair of paths within each cluster for (A) GluN2A and (B) GluN1. The dotted line in each plot …

Figure 1—figure supplement 2
D-serine pathway residues mapped onto the intact GluN2A NMDAR (PDB ID: 6MMM Hansen et al., 2021).

This figure supplement relates to Figures 1 and 2. Sidechain atoms are shown as spheres for (A) GluN2A in teal and (B) GluN1 in magenta and the adjacent domains labeled.

Figure 1—figure supplement 3
Degree of LBD closure during D-serine binding pathways to (A) GluN2A and (B) GluN1.

This figure supplement relates to Figures 1 and 2. The two-dimensional (ξ1,ξ2) order parameter was computed for all frames in binding/unbinding pathways and was plotted as a function of relative …

Figure 1—figure supplement 4
Dendrograms for hierarchical clustering of weighted average Hausdorff distances for D-serine binding pathways for (A) GluN2A and (B) GluN1 according to the Ward linkage criterion.

This figure supplement relates to Figures 1 and 2.

Identifying D-serine binding pathways for GluN1 using pathway similarity analysis (PSA).

(A) 2-dimensional order parameter (ξ1,ξ2) that describes the degree of GluN1 LBD closure. For each of the above (B–E), the left image shows D-serine density, while the right image shows the residues …

Figure 2—source data 1

Record of all successful binding pathways in each simulation system for D-serine binding to GluN1.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig2-data1-v1.xlsx
Figure 2—source data 2

Per residue contact frequency analysis for D-serine binding to GluN1 by cluster identified with PSA.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig2-data2-v1.xlsx
Figure 2—source data 3

GluN1 residues most frequently contacted by D-serine given that the pathway results in successful binding – listed for each simulation system.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig2-data3-v1.xlsx
Figure 3 with 2 supplements
Conformational free energy landscapes for GluN2A and GluN1 LBDs.

Umbrella sampling molecular dynamics simulations were used to compute the potential of mean force (PMF) along the (ξ1,ξ2) order parameter for (A) D-serine bound to GluN2A, (B) glycine bound to GluN2A, …

Figure 3—source data 1

Per-residue contact frequency analysis of the bound state for each agonist computed from lowest-energy conformers extracted from umbrella sampling simulations.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig3-data1-v1.xlsx
Figure 3—figure supplement 1
GluN2A residues contacting D-serine in lowest-energy conformers.

(A) Binding-site residues for D-serine bound to the GluN2A LBD computed from lowest-energy conformers (1 kcal mol–1) from umbrella sampling simulations. Residues with >50% contact frequency in the …

Figure 3—figure supplement 2
Error of umbrella sampling PMFs computed by block averaging for (A) D-serine bound to GluN2A, (B) glycine bound to GluN2A, (C) glutamate bound to GluN2A in the inverted pose, and (D) D-serine bound to GluN1.

The colors of the colorbar correspond to the standard deviation in kcal mol–1, which was computed for each window over 10 blocks.

D-serine competes glutamate binding as an antagonist at high concentration.

(A) Representative Two-electrode voltage clamp (TEVC) recording on GluN1/GluN2A NMDARs expressing oocytes. The traces show inhibition of the NMDAR current by the GluN1 agonists D-serine (left) and …

Comparison of D-serine and glutamate binding to GluN2A.

(A) Overlay of D-serine (teal) and glutamate (gray) density. (B) Residues that distinguish D-serine (teal) from glutamate (gray) binding pathways (see Figure 5—source data 2).

Figure 5—source data 1

Record of all successful binding pathways in each simulation system for glutamate binding to GluN2A.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig5-data1-v1.xlsx
Figure 5—source data 2

Comparison of relative residue contact frequency for D-serine and glutamate.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig5-data2-v1.xlsx
Figure 6 with 2 supplements
Conformational dynamics of near-pocket glycans.

N-linked Man5GlcNAc2 (Man5) glycans (A) N443-Man5 and N444-Man5 for GluN2A and (B) N491-Man5 for GluN1. Glycan conformational energy landscapes for (C) GluN2A N443-Man5, (D) GluN2A N444-Man5, and (E)…

Figure 6—source data 1

Comparison of relative residue contact frequency during GluN2A and GluN1 binding pathways for glycosylated and non-glycosylated simulations.

https://cdn.elifesciences.org/articles/77645/elife-77645-fig6-data1-v1.xlsx
Figure 6—figure supplement 1
N-linked Man5GlcNAc2 (Man5) glycans interacts with D-serine as it binds.

(A) Schematic of the Man5 glycan. (B) Contact network involving the GluN2A N443-Man5 glycan and residues Glu-412 and, Lys-438, (C) Lys-738 and Glu-413, (D) Glu-413,Tyr-730, and Ser-511. (E) Contact …

Figure 6—figure supplement 2
Glycan-D2 distance dependence on agonist binding for the (A) GluN2A N443-Man5 glycan, (B) GluN2A N444-Man5 glycan, and (C) GluN1 N491-Man5 glycan.
Glycan-D2 distance dependence on LBD closure for the (A) GluN2A N443-Man5 glycan, (B) GluN2A N444-Man5 glycan, and (C) GluN1 N491-Man5 glycan.

Videos

Video 1
Process of D-serine binding to the GluN2A LBD.
Video 2
Process of D-serine binding to the GluN1 LBD.

Additional files

Download links