Ionotropic glutamate receptors (iGluRs) are ligand-gated ion channels that are responsible for the majority of excitatory transmission at the synaptic cleft. Mechanically speaking, agonist binding to the ligand binding domain (LBD) activates the receptor by triggering a conformational change that is transmitted to the transmembrane region, opening the ion channel pore. We use fully atomistic molecular dynamics simulations to investigate the binding process in the α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor, an iGluR subtype. The string method with swarms of trajectories was applied to calculate the possible pathways glutamate traverses during ligand binding. Residues peripheral to the binding cleft are found to metastably bind the ligand prior to ligand entry into the binding pocket. Umbrella sampling simulations were performed to compute the free energy barriers along the binding pathways. The calculated free energy profiles demonstrate that metastable interactions contribute substantially to the energetics of ligand binding and form local minima in the overall free energy landscape. Protein-ligand interactions at sites outside of the orthosteric agonist-binding site may serve to lower the transition barriers of the binding process.

At central nervous system synapses, agonist binding to postsynaptic ionotropic glutamate receptors (iGluRs) results in signaling between neurons. N-Methyl-D-aspartic acid (NMDA) receptors are a unique family of iGluRs that activate in response to the concurrent binding of glutamate and glycine. Here, we investigate the process of agonist binding to the GluN2A (glutamate binding) and GluN1 (glycine binding) NMDA receptor subtypes using long-timescale unbiased molecular dynamics simulations. We find that positively charged residues on the surface of the GluN2A ligand-binding domain (LBD) assist glutamate binding via a "guided-diffusion" mechanism, similar in fashion to glutamate binding to the GluA2 LBD of AMPA receptors. Glutamate can also bind in an inverted orientation. Glycine, on the other hand, binds to the GluN1 LBD via an "unguided-diffusion" mechanism, whereby glycine finds its binding site primarily by random thermal fluctuations. Free energy calculations quantify the glutamate- and glycine-binding processes.

The earliest metazoan ancestors of humans include the ctenophore Mnemiopsis leidyi The genome of this comb jelly encodes homologs of vertebrate ionotropic glutamate receptors (iGluRs) that are distantly related to glycine-activated NMDA receptors and that bind glycine with unusually high affinity. Using ligand-binding domain (LBD) mutants for electrophysiological analysis, we demonstrate that perturbing a ctenophore-specific interdomain Arg-Glu salt bridge that is notably absent from vertebrate AMPA, kainate, and NMDA iGluRs greatly increases the rate of recovery from desensitization, while biochemical analysis reveals a large decrease in affinity for glycine. X-ray crystallographic analysis details rearrangements in the binding pocket stemming from the mutations, and molecular dynamics simulations suggest that the interdomain salt bridge acts as a steric barrier regulating ligand binding and that the free energy required to access open conformations in the glycine-bound LBD is largely responsible for differences in ligand affinity among the LBD variants.

K(2P)Ø, the two-pore domain potassium background channel that determines cardiac rhythm in Drosophila melanogaster, and its homologues that establish excitable membrane activity in mammals are of unknown structure. K(2P) subunits have two pore domains flanked by transmembrane (TM) spans: TM1-P1-TM2-TM3-P2-TM4. To establish spatial relationships in K(2P)Ø, we identified pairs of sites that display electrostatic compensation. Channels silenced by the addition of a charge in pore loop 1 (P1) or P2 were restored to function by countercharges at specific second sites. A three-dimensional homology model was determined using the crystal structure of K(V)1.2, effects of K(2P)Ø mutations to establish alignment, and compensatory charge-charge pairs. The model was refined and validated by continuum electrostatic free energy calculations and covalent linkage of introduced cysteines. K(2P) channels use two subunits arranged so that the P1 and P2 loops contribute to one pore, identical P loops face each other diagonally across the pore, and the channel complex has bilateral symmetry with a fourfold symmetric selectivity filter.