Computational study of synthetic agonist ligands of ionotropic glutamate receptors

Neurotransmitter Funneling Optimizes Glutamate Receptor Kinetics

Ionotropic glutamate receptors (iGluRs) mediate neurotransmission at the majority of excitatory synapses in the brain. Little is known, however, about how glutamate reaches the recessed binding pocket in iGluR ligand-binding domains (LBDs). Here we report the process of glutamate binding to a prototypical iGluR, GluA2, in atomistic detail using unbiased molecular simulations. Charged residues on the LBD surface form pathways that facilitate glutamate binding by effectively reducing a three-dimensional diffusion process to a spatially constrained, two-dimensional one. Free energy calculations identify residues that metastably bind glutamate and help guide it into the binding pocket. These simulations also reveal that glutamate can bind in an inverted conformation and also reorient while in its pocket. Electrophysiological recordings demonstrate that eliminating these transient binding sites slows activation and deactivation, consistent with slower glutamate binding and unbinding. These results suggest that binding pathways have evolved to optimize rapid responses of AMPA-type iGluRs at synapses.

An improved quantum biochemistry description of the glutamate–GluA2 receptor binding within an inhomogeneous dielectric function framework

A new methodology to define the inhomogeneous dielectric constant of protein residues, to apply to the calculation of protein–ligand properties such as the electrostatic interaction.

Molecular Dynamics Investigation of gluazo, a Photo-Switchable Ligand for the Glutamate Receptor GluK2

Photochromic ligands (PCLs), defined as photoswitchable molecules that are able to endow native receptors with a sensitivity towards light, have become a promising photopharmacological tool for various applications in biology. In general, PCLs consist of a ligand of the target receptor covalently linked to an azobenzene, which can be reversibly switched between two configurations upon light illumination. Gluazo, as a PCL that targets excitatory amino acid receptors, in its dark-adapted trans iso-form was characterized to be a partial agonist of the kainate glutamate receptor GluK2. Application of UV light leads to the formation of the cis form, with remarkedly reduced affinity towards GluK2. The mechanism of the change of ligand affinity induced by the photoisomerization was unresolved. The presented computational study explains how the isomerization of such a PCL affects the structural changes in the target receptor that lead to its activation.

Ligand photo-isomerization triggers conformational changes in iGluR2 ligand binding domain

Neurological glutamate receptors bind a variety of artificial ligands, both agonistic and antagonistic, in addition to glutamate. Studying their small molecule binding properties increases our understanding of the central nervous system and a variety of associated pathologies. The large, oligomeric multidomain membrane protein contains a large and flexible ligand binding domains which undergoes large conformational changes upon binding different ligands. A recent application of glutamate receptors is their activation or inhibition via photo-switchable ligands, making them key systems in the emerging field of optochemical genetics. In this work, we present a theoretical study on the binding mode and complex stability of a novel photo-switchable ligand, ATA-3, which reversibly binds to glutamate receptors ligand binding domains (LBDs). We propose two possible binding modes for this ligand based on flexible ligand docking calculations and show one of them to be analogues to the binding mode of a similar ligand, 2-BnTetAMPA. In long MD simulations, it was observed that transitions between both binding poses involve breaking and reforming the T686-E402 protein hydrogen bond. Simulating the ligand photo-isomerization process shows that the two possible configurations of the ligand azo-group have markedly different complex stabilities and equilibrium binding modes. A strong but slow protein response is observed after ligand configuration changes. This provides a microscopic foundation for the observed difference in ligand activity upon light-switching.

Computation of standard binding free energies of polar and charged ligands to the glutamate receptor GluA2

Accurate calculation of the binding affinity of small molecules to proteins has the potential to become an important tool in rational drug design. In this study, we use the free energy perturbation (FEP) method with restraints to calculate the standard binding free energy of five ligands (ACPA, AMPA, CNQX, DNQX, and glutamate) to the glutamate receptor GluA2, which plays an essential role in synaptic transmission. To deal with the convergence problem in FEP calculations with charged ligands, we use a protocol where the ligand is coupled in the binding site while it is decoupled in bulk solution simultaneously. The contributions from the conformational, rotational, and translational entropies to the standard binding free energy are determined by applying/releasing respective restraints to the ligand in bulk/binding site. We also employ the confine-and-release approach, which helps to resolve convergence problems in FEP calculations. Our results are in good agreement with the experimental values for all five ligands, including the charged ones which are often problematic in FEP calculations. We also analyze the different contributions to the binding free energy of each ligand to GluA2 and discuss the nature of these interactions.

Conformational analysis of NMDA receptor GluN1, GluN2, and GluN3 ligand-binding domains reveals subtype-specific characteristics

The NMDA receptor family of glutamate receptor ion channels is formed by obligate heteromeric assemblies of GluN1, GluN2, and GluN3 subunits. GluN1 and GluN3 bind glycine, whereas GluN2 binds glutamate. Crystal structures of the GluN1 and GluN3A ligand-binding domains (LBDs) in their apo states unexpectedly reveal open- and closed-cleft conformations, respectively, with water molecules filling the binding pockets. Computed conformational free energy landscapes for GluN1, GluN2A, and GluN3A LBDs reveal that the apo-state LBDs sample closed-cleft conformations, suggesting that their agonists bind via a conformational selection mechanism. By contrast, free energy landscapes for the AMPA receptor GluA2 LBD suggest binding of glutamate via an induced-fit mechanism. Principal component analysis reveals a rich spectrum of hinge bending, rocking, twisting, and sweeping motions that are different for the GluN1, GluN2A, GluN3A, and GluA2 LBDs. This variation highlights the structural complexity of signaling by glutamate receptor ion channels.