Allosteric coupling of sub-millisecond clamshell motions in ionotropic glutamate receptor ligand-binding domains

Arthropod repellent interactions with olfactory receptors and ionotropic receptors analyzed by molecular modeling

The main insect chemoreceptors are olfactory receptors (ORs), gustatory receptors (GRs) and ionotropic receptors (IRs). The odorant binding sites of many insect ORs appear to be occluded and inaccessible from the surface of the receptor protein, based on the three-dimensional structure of OR5 from the jumping bristletail Machilis hrabei (MhraOR5) and a survey of a sample of vinegar fly (Drosophila melanogaster) OR structures obtained from artificial intellegence (A.I.) modeling. Molecular dynamics simulations revealed that the occluded site can become accessible through tunnels that transiently open and close. The present study extends this analysis to examine seventeen ORs and one GR docking with ligands that have known valence: nine that signal attraction and nine that signal aversion. All but one of the receptors displayed occluded ligand binding sites analogous to MhraOR5, and docking software predicted the known attractant and repellent ligands will bind to the occluded sites. Docking of the repellent DEET was examined, and more than half of the OR ligand sites were predicted to bind DEET, including receptors that signal aversion as well as those that signal attraction. However, DEET may not actually have access to all the attractant binding sites. The larger size and lower flexibility of repellent molecules may restrict their passage through the tunnel bottlenecks, which could act as filters to select access to the ligand binding sites. In contrast to ORs and GRs, the IR ligand binding site is in an extracellular domain known to undergo a large conformational change from an open to a closed state. A.I. models of two D. melanogaster IRs of known valence and two blacklegged tick (Ixodes scapularis) IRs having unknown ligands were computationally tested for attractant and repellent binding. The ligand-binding sites in the closed state appear inaccessible to the protein surface, so attractants and repellents must bind initially at an accessible site in the open state before triggering the conformational change. In some IRs, repellent binding sites were identified at exterior sites adjacent to the ligand-binding site. These may be allosteric sites that, when occupied by repellents, can stabilize the open state of an attractant IR, or stabilize the closed state of an IR in the absence of its activating ligand. The model of D. melanogaster IR64a suggests a possible molecular mechanism for the activation of this IR by H+. The amino acids involved in this proposed mechanism are conserved in IR64a from several Dipteran pest species and disease vectors, potentially offering a route to discovery of new repellents that act via the allosteric site.

Partial agonism in heteromeric GLUK2/GLUK5 kainate receptor

Kainate receptors are a subtype of ionotropic glutamate receptors that form transmembrane channels upon binding glutamate. Here, we have investigated the mechanism of partial agonism in heteromeric GluK2/K5 receptors, where the GluK2 and GluK5 subunits have distinct agonist binding profiles. Using single-molecule Förster resonance energy transfer, we found that at the bi-lobed agonist-binding domain, the partial agonist AMPA-bound receptor occupied intermediate cleft closure conformational states at the GluK2 cleft, compared to the more open cleft conformations in apo form and more closed cleft conformations in the full agonist glutamate-bound form. In contrast, there is no significant difference in cleft closure states at the GluK5 agonist-binding domain between the partial agonist AMPA- and full agonist glutamate-bound states. Additionally, unlike the glutamate-bound state, the dimer interface at the agonist-binding domain is not decoupled in the AMPA-bound state. Our findings suggest that partial agonism observed with AMPA binding is mediated primarily due to differences in the GluK2 subunit, highlighting the distinct contributions of the subunits towards activation.

Excitatory and inhibitory D-serine binding to the NMDA receptor

N-methyl-D-aspartate receptors (NMDARs) uniquely require binding of two different neurotransmitter agonists for synaptic transmission. D-serine and glycine bind to one subunit, GluN1, while glutamate binds to the other, GluN2. These agonists bind to the receptor's bi-lobed ligand-binding domains (LBDs), which close around the agonist during receptor activation. To better understand the unexplored mechanisms by which D-serine contributes to receptor activation, we performed multi-microsecond molecular dynamics simulations of the GluN1/GluN2A LBD dimer with free D-serine and glutamate agonists. Surprisingly, we observed D-serine binding to both GluN1 and GluN2A LBDs, suggesting that D-serine competes with glutamate for binding to GluN2A. This mechanism is confirmed by our electrophysiology experiments, which show that D-serine is indeed inhibitory at high concentrations. Although free energy calculations indicate that D-serine stabilizes the closed GluN2A LBD, its inhibitory behavior suggests that it either does not remain bound long enough or does not generate sufficient force for ion channel gating. We developed a workflow using pathway similarity analysis to identify groups of residues working together to promote binding. These conformation-dependent pathways were not significantly impacted by the presence of N-linked glycans, which act primarily by interacting with the LBD bottom lobe to stabilize the closed LBD.