The mechanism of partial agonism at N-methyl-D-aspartate receptors is an unresolved issue, especially with respect to the role of protein dynamics. We have performed multiple molecular dynamics simulations (7 ؋ 20 ns) to examine the behavior of the ligand-binding core of the NR1 subunit with a series of ligands. Our results show that water plays an important role in stabilizing different conformations of the core and how a closed cleft conformation of the protein might be stabilized in the absence of ligands. In the case of ligand-bound simulations with both full and partial agonists, we observed that ligands within the binding cleft may undergo distinct conformational changes, without grossly influencing the degree of cleft closure within the ligand-binding domain. In agreement with recently published crystallographic data, we also observe similar changes in backbone torsions corresponding to the hinge region between the two lobes for the partial agonist, D-cycloserine. This observation rationalizes the classification of D-cycloserine as a partial agonist and should provide a basis with which to predict partial agonism in this class of receptor by analyzing the behavior of these torsions with other potential ligands.
N-Methyl-D-aspartate (NMDA)4 receptors are ligand-gated ion channels that play a major role in learning and memory (1). They have also been implicated in a number of injury or disease states including for example, schizophrenia (2-5). NMDA receptors are a distinct subtype of ionotropic glutamate receptor (6) in that they require both glutamate and glycine for activation and membrane depolarization (7). The receptor is a heterotetrameric cation channel typically comprised of two NR1 subunits that contain the glycine-binding site along with two NR2A-D glutamate-binding sites. The heterogeneity of the complex is further extended if the third subtype (NR3), which also binds glycine, is also taken into account.The architecture of each subunit of NMDA receptors is similar to that of non-NMDA receptors. Each subunit is comprised of an aminoterminal binding domain that shows homology to the LIVBP (leucineisoleucine-valine-binding protein) (8), three transmembrane helices (M1, M2, and M3), a re-entrant P-loop, and a ligand-binding domain that is comprised of two discontinuous polypeptide chains, S1 and S2, that form two distinct lobes or subdomains, D1 and D2 (Fig. 1A). An isolated D1D2 ligand-binding core has been shown to bind ligands with similar affinity as wild-type receptors (9, 10).Recently, the crystal structure of the ligand-binding domain of the NR1 subunit has been solved by x-ray diffraction with both full (9) and partial (10) agonists. The overall fold and structure (Fig. 1B) is very similar to those observed for crystal structures of the ligand-binding domains for the (non-NMDA) GluR2 (11-13), GluR5 (14,15), and GluR6 receptors (Refs. 15 and 16; for recent reviews see Refs. 17-19). The major difference is an insertion of 30 residues into loop 1, which also contains two disulfide bridges. Inter...