Full and partial agonism of ionotropic glutamate receptors indicated by molecular dynamics simulations

Unlocking the Reverse Targeting Mechanisms of Cannabidiol: Unveiling New Therapeutic Avenues

Cannabidiol (CBD) and Δ9-tetrahydrocannabinol (THC), the main components of Cannabis sativa plants, have attracted a significant amount of attention due to their biological activities. This study identified GPR18 as the target of partial agonist CBD activating the p42/p44 MAPK pathway leading to migration of endometrial epithelial cells. Induced fit docking (IFD) showed that the affinity of THC for GPR18 is higher than that of CBD, and molecular dynamics (MD) simulations showed that CBD-GPR18 complexes at 130/200 ns might have stable conformations, potentially activating GPR18 by changing the distances of key residues in its active pocket. In contrast, THC maintains "metastable" conformations, generating a "shrinking space" leading to full agonism of THC by adding mechanical constraints in GPR18's active pocket. Steered molecular dynamics (SMD) revealed GPR18's active pocket was influenced more by CBD's partial agonism compared with THC. This combined IFD-MD-SMD method may be used to explain the mechanism of activation of partial or full agonists of GPR18.

Willardiine and Its Synthetic Analogues: Biological Aspects and Implications in Peptide Chemistry of This Nucleobase Amino Acid

Willardiine is a nonprotein amino acid containing uracil, and thus classified as nucleobase amino acid or nucleoamino acid, that together with isowillardiine forms the family of uracilylalanines isolated more than six decades ago in higher plants. Willardiine acts as a partial agonist of ionotropic glutamate receptors and more in particular it agonizes the non-N-methyl-D-aspartate (non-NMDA) receptors of L-glutamate: ie. the α-amino-3-hydroxy-5-methyl-4-isoxazole-propionic acid (AMPA) and kainate receptors. Several analogues and derivatives of willardiine have been synthesised in the laboratory in the last decades and these compounds show different binding affinities for the non-NMDA receptors. More in detail, the willardiine analogues have been employed not only in the investigation of the structure of AMPA and kainate receptors, but also to evaluate the effects of receptor activation in the various brain regions. Remarkably, there are a number of neurological diseases determined by alterations in glutamate signaling, and thus, ligands for AMPA and kainate receptors deserve attention as potential neurodrugs. In fact, similar to willardiine its analogues often act as agonists of AMPA and kainate receptors. A particular importance should be recognized to willardiine and its thymine-based analogue AlaT also in the peptide chemistry field. In fact, besides the naturally-occurring short nucleopeptides isolated from plant sources, there are different examples in which this class of nucleoamino acids was investigated for nucleopeptide development. The applications are various ranging from the realization of nucleopeptide/DNA chimeras for diagnostic applications, and nucleoamino acid derivatization of proteins for facilitating protein-nucleic acid interaction, to nucleopeptide-nucleopeptide molecular recognition for nanotechnological applications. All the above aspects on both chemistry and biotechnological applications of willardine/willardine-analogues and nucleopeptide will be reviewed in this work.

Computational Molecular Characterization of the Interaction of Acetylcholine and the NMDA Receptor to Explain the Direct Glycine-Competitive Potentiation of NMDA-Mediated Neuronal Currents

The activation of N-methyl-d-aspartate receptor (NMDAR) is triggered by the closure of bilobed (D1 and D2) clamshell-like clefts upon binding glycine (Gly) and glutamate. There is evidence that cholinergic compounds modulate NMDAR-mediated currents via direct receptor-ligand interactions; however, molecular bases are unknown. Here, we first propose a mechanistic structure-based explanation for the observed ACh-induced submaximal potentiation of NMDA-elicited currents in striatal neurons by predicting competitive inhibition with Gly. Then, the model was validated, in principle, by confirming that the coapplication of Gly and ACh significantly reduces these neuronal currents. Finally, we delineate the interplay of ACh with the NMDAR by a combination of computational strategies. Crystallographic ACh-bound complexes were studied, revealing a similar ACh binding environment on the GluN1 subunit of the NMDAR. We illustrate how ACh can occupy X-ray monomeric open, dimeric "semiopen" cleft conformations obtained by molecular dynamics and a full-active cryo-EM NMDAR structure, explaining the suboptimal NMDAR electrophysiological activity under the "Venus Flytrap model". At an evolutionary biology level, the binding mode of ACh coincides with that of the homologous ornithine-bound periplasmic LAO binding protein complex. Our computed results indicate an analogous mechanism of action, inasmuch as ACh may stabilize the GluN1 subunit "semiclosed" conformations by inducing direct and indirect D1-to-D2 interdomain bonds. Additionally, an alternative binding site was detected, shared by the known NMDAR allosteric modulators. Experimental and computed results strongly suggest that ACh acts as a Gly-competitive, submaximal potentiating agent of the NMDAR, possibly constituting a novel chemotype for multitarget-directed drug development, e.g., to treat Alzheimer's, and it may lead to a new understanding of glutamatergic neurotransmission.

Structure, Function, and Pharmacology of Glutamate Receptor Ion Channels

Many physiologic effects of l-glutamate, the major excitatory neurotransmitter in the mammalian central nervous system, are mediated via signaling by ionotropic glutamate receptors (iGluRs). These ligand-gated ion channels are critical to brain function and are centrally implicated in numerous psychiatric and neurologic disorders. There are different classes of iGluRs with a variety of receptor subtypes in each class that play distinct roles in neuronal functions. The diversity in iGluR subtypes, with their unique functional properties and physiologic roles, has motivated a large number of studies. Our understanding of receptor subtypes has advanced considerably since the first iGluR subunit gene was cloned in 1989, and the research focus has expanded to encompass facets of biology that have been recently discovered and to exploit experimental paradigms made possible by technological advances. Here, we review insights from more than 3 decades of iGluR studies with an emphasis on the progress that has occurred in the past decade. We cover structure, function, pharmacology, roles in neurophysiology, and therapeutic implications for all classes of receptors assembled from the subunits encoded by the 18 ionotropic glutamate receptor genes. SIGNIFICANCE STATEMENT: Glutamate receptors play important roles in virtually all aspects of brain function and are either involved in mediating some clinical features of neurological disease or represent a therapeutic target for treatment. Therefore, understanding the structure, function, and pharmacology of this class of receptors will advance our understanding of many aspects of brain function at molecular, cellular, and system levels and provide new opportunities to treat patients.

Getting Docking into Shape Using Negative Image-Based Rescoring

The failure of default scoring functions to ensure virtual screening enrichment is a persistent problem for the molecular docking algorithms used in structure-based drug discovery. To remedy this problem, elaborate rescoring and postprocessing schemes have been developed with a varying degree of success, specificity, and cost. The negative image-based rescoring (R-NiB) has been shown to improve the flexible docking performance markedly with a variety of drug targets. The yield improvement is achieved by comparing the alternative docking poses against the negative image of the target protein's ligand-binding cavity. In other words, the shape and electrostatics of the binding pocket is directly used in the similarity comparison to rank the explicit docking poses. Here, the PANTHER/ShaEP-based R-NiB methodology is tested with six popular docking softwares, including GLIDE, PLANTS, GOLD, DOCK, AUTODOCK, and AUTODOCK VINA, using five validated benchmark sets. Overall, the results indicate that R-NiB outperforms the default docking scoring consistently and inexpensively, demonstrating that the methodology is ready for wide-scale virtual screening usage.

Discovery of Retinoic Acid-Related Orphan Receptor γt Inverse Agonists via Docking and Negative Image-Based Screening

Retinoic acid-related orphan receptor γt (RORγt) has a vital role in the differentiation of T-helper 17 (TH17) cells. Potent and specific RORγt inverse agonists are sought for treating TH17-related diseases such as psoriasis, rheumatoid arthritis, and type 1 diabetes. Here, the aim was to discover novel RORγt ligands using both standard molecular docking and negative image-based screening. Interestingly, both of these in silico techniques put forward mostly the same compounds for experimental testing. In total, 11 of the 34 molecules purchased for testing were verified as RORγt inverse agonists, thus making the effective hit rate 32%. The pIC₅₀ values for the compounds varied from 4.9 (11 μM) to 6.2 (590 nM). Importantly, the fact that the verified hits represent four different cores highlights the structural diversity of the RORγt inverse agonism and the ability of the applied screening methodologies to facilitate much-desired scaffold hopping for drug design.

Computer Simulations Predict High Structural Heterogeneity of Functional State of NMDA Receptors

N-methyl-D-aspartate receptors (NMDARs)-i.e., transmembrane proteins expressed in neurons-play a central role in the molecular mechanisms of learning and memory formation. It is unclear how the known atomic structures of NMDARs determined by x-ray crystallography and electron cryomicroscopy (18 published Protein Data Bank entries) relate to the functional states of NMDARs inferred from electrophysiological recordings (multiple closed, open, preopen, etc. states). We address this problem by using molecular dynamics simulations at atomic resolution, a method successfully applied in the past to much smaller biomolecules. Our simulations predict that several conformations of NMDARs with experimentally determined geometries, including four "nonactive" electron cryomicroscopy structures, rapidly interconvert on submicrosecond timescales and therefore may correspond to the same functional state of the receptor (specifically, one of the closed states). This conclusion is not trivial because these conformational transitions involve changes in certain interatomic distances as large as tens of Å. The simulations also predict differences in the conformational dynamics of the apo and holo (i.e., agonist and coagonist bound) forms of the receptor on the microsecond timescale. To our knowledge, five new conformations of NMDARs, with geometries joining various features from different known experimental structures, are also predicted by the model. The main limitation of this work stems from its limited sampling (30 μs of aggregate length of molecular dynamics trajectories). Though this level significantly exceeds the sampling in previous simulations of parts of NMDARs, it is still much lower than the sampling recently achieved for smaller biomolecules (up to a few milliseconds), thus precluding, in particular, the observation of transitions between different functional states of NMDARs. Despite this limitation, such computational predictions may guide further experimental studies on the structure, dynamics, and function of NMDARs, for example by suggesting optimal locations of spectroscopic probes. Overall, atomic resolution simulations provide, to our knowledge, a novel perspective on the structure and dynamics of NMDARs, complementing information obtained by experimental methods.

Simulated Dynamics of Glycans on Ligand-Binding Domain of NMDA Receptors Reveals Strong Dynamic Coupling between Glycans and Protein Core

N-Methyl-d-aspartate (NMDA) receptors, key neuronal receptors playing the central role in learning and memory, are heavily glycosylated in vivo. Astonishingly little is known about the structure, dynamics, and physiological relevance of glycans attached to them. We recently demonstrated that certain glycans on the ligand binding domain (LBD) of NMDA receptors (NMDARs) can serve as intramolecular potentiators, changing EC50 of NMDAR coagonists. In this work, we use molecular dynamics trajectories, in aggregate 86.5 μs long, of the glycosylated LBD of the GluN1 subunit of the NMDAR to investigate the behavior of glycans on NMDARs. Though all glycans in our simulations were structurally the same (Man5), the dynamics of glycans at different locations on NMDARs was surprisingly different. The slowest-time scale motions that we detected in various glycans in some cases corresponded to a flipping of parts of glycans relative to each other, while in other cases they reduced to a head-to-tail bending of a glycan. We predict that time scales of conformational changes in glycans on the GluN1 LBD of NMDARs range from nanoseconds to at least hundreds of microseconds. Some of the conformational changes in the glycans correlate with the physiologically important clamshell-like opening and closing of the GluN1 LBD domain. Thus, glycans are an integral part of NMDARs, and computational models of NMDARs should include glycans to faithfully represent the structure and the dynamics of these receptors.

Role of the Ion Channel Extracellular Collar in AMPA Receptor Gating

AMPA subtype ionotropic glutamate receptors mediate fast excitatory neurotransmission and are implicated in numerous neurological diseases. Ionic currents through AMPA receptor channels can be allosterically regulated via different sites on the receptor protein. We used site-directed mutagenesis and patch-clamp recordings to probe the ion channel extracellular collar, the binding region for noncompetitive allosteric inhibitors. We found position and substitution-dependent effects for introduced mutations at this region on AMPA receptor gating. The results of mutagenesis suggested that the transmembrane domains M1, M3 and M4, which contribute to the ion channel extracellular collar, undergo significant relative displacement during gating. We used molecular dynamics simulations to predict an AMPA receptor open state structure and rationalize the results of mutagenesis. We conclude that the ion channel extracellular collar plays a distinct role in gating and represents a hub for powerful allosteric modulation of AMPA receptor function that can be used for developing novel therapeutics.

Computationally Discovered Potentiating Role of Glycans on NMDA Receptors

N-methyl-D-aspartate receptors (NMDARs) are glycoproteins in the brain central to learning and memory. The effects of glycosylation on the structure and dynamics of NMDARs are largely unknown. In this work, we use extensive molecular dynamics simulations of GluN1 and GluN2B ligand binding domains (LBDs) of NMDARs to investigate these effects. Our simulations predict that intra-domain interactions involving the glycan attached to residue GluN1-N440 stabilize closed-clamshell conformations of the GluN1 LBD. The glycan on GluN2B-N688 shows a similar, though weaker, effect. Based on these results, and assuming the transferability of the results of LBD simulations to the full receptor, we predict that glycans at GluN1-N440 might play a potentiator role in NMDARs. To validate this prediction, we perform electrophysiological analysis of full-length NMDARs with a glycosylation-preventing GluN1-N440Q mutation, and demonstrate an increase in the glycine EC50 value. Overall, our results suggest an intramolecular potentiating role of glycans on NMDA receptors.

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.

Conformational transitions in the glycine-bound GluN1 NMDA receptor LBD via single-molecule FRET

The N-methyl-D-aspartate receptor (NMDAR) is a member of the glutamate receptor family of proteins and is responsible for excitatory transmission. Activation of the receptor is thought to be controlled by conformational changes in the ligand binding domain (LBD); however, glutamate receptor LBDs can occupy multiple conformations even in the activated form. This work probes equilibrium transitions among NMDAR LBD conformations by monitoring the distance across the glycine-bound LBD cleft using single-molecule Förster resonance energy transfer (smFRET). Recent improvements in photoprotection solutions allowed us to monitor transitions among the multiple conformations. Also, we applied a recently developed model-free algorithm called "step transition and state identification" to identify the number of states, their smFRET efficiencies, and their interstate kinetics. Reversible interstate conversions, corresponding to transitions among a wide range of cleft widths, were identified in the glycine-bound LBD, on much longer timescales compared to channel opening. These transitions were confirmed to be equilibrium in nature by shifting the distribution reversibly via denaturant. We found that the NMDAR LBD proceeds primarily from one adjacent smFRET state to the next under equilibrium conditions, consistent with a cleft-opening/closing mechanism. Overall, by analyzing the state-to-state transition dynamics and distributions, we achieve insight into specifics of long-lived LBD equilibrium structural dynamics, as well as obtain a more general description of equilibrium folding/unfolding in a conformationally dynamic protein. The relationship between such long-lived LBD dynamics and channel function in the full receptor remains an open and interesting question.

Selective effect of cell membrane on synaptic neurotransmission

Atomistic molecular dynamics simulations were performed with 13 non-peptidic neurotransmitters (NTs) in three different membrane environments. The results provide compelling evidence that NTs are divided into membrane-binding and membrane-nonbinding molecules. NTs adhere to the postsynaptic membrane surface whenever the ligand-binding sites of their synaptic receptors are buried in the lipid bilayer. In contrast, NTs that have extracellular ligand-binding sites do not have a similar tendency to adhere to the membrane surface. This finding is a seemingly simple yet important addition to the paradigm of neurotransmission, essentially dividing it into membrane-independent and membrane-dependent mechanisms. Moreover, the simulations also indicate that the lipid composition especially in terms of charged lipids can affect the membrane partitioning of NTs. The revised paradigm, highlighting the importance of cell membrane and specific lipids for neurotransmission, should to be of interest to neuroscientists, drug industry and the general public alike.

Identification of estrogen receptor α ligands with virtual screening techniques

Utilization of computer-aided molecular discovery methods in virtual screening (VS) is a cost-effective approach to identify novel bioactive small molecules. Unfortunately, no universal VS strategy can guarantee high hit rates for all biological targets, but each target requires distinct, fine-tuned solutions. Here, we have studied in retrospective manner the effectiveness and usefulness of common pharmacophore hypothesis, molecular docking and negative image-based screening as potential VS tools for a widely applied drug discovery target, estrogen receptor α (ERα). The comparison of the methods helps to demonstrate the differences in their ability to identify active molecules. For example, structure-based methods identified an already known active ligand from the widely-used bechmarking decoy molecule set. Although prospective VS against one commercially available database with around 100,000 drug-like molecules did not retrieve many testworthy hits, one novel hit molecule with pIC50 value of 6.6, was identified. Furthermore, our small in-house compound collection of easy-to-synthesize molecules was virtually screened against ERα, yielding to five hit candidates, which were found to be active in vitro having pIC50 values from 5.5 to 6.5.

Mechanism-Based Mathematical Model for Gating of Ionotropic Glutamate Receptors

We present a mathematical model for ionotropic glutamate receptors (iGluR's) that is built on mechanistic understanding and yields a number of thermodynamic and kinetic properties of channel gating. iGluR's are ligand-gated ion channels responsible for the vast majority of fast excitatory neurotransmission in the central nervous system. The effects of agonist-induced closure of the ligand-binding domain (LBD) are transmitted to the transmembrane channel (TMC) via interdomain linkers. Our model demonstrates that, relative to full agonists, partial agonists may reduce either the degree of LBD closure or the curvature of the LBD free energy basin, leading to less stabilization of the channel open state and hence lower channel open probability. A rigorous relation is derived between the channel closed-to-open free energy difference and the tension within the linker. Finally, by treating LBD closure and TMC opening as diffusive motions, we obtain gating trajectories that resemble stochastic current traces from single-channel recordings and calculate the rate constants for transitions between the channel open and closed states. Our model can be implemented by molecular dynamics simulations to realistically depict iGluR gating and may guide functional experiments in gaining deeper insight into this essential family of channel proteins.

A quantum biochemistry investigation of willardiine partial agonism in AMPA receptors

We employ quantum biochemistry methods based on the Density Functional Theory (DFT) approach to unveil detailed binding energy features of willardiines co-crystallized with the AMPA receptor.

Reduced curvature of ligand-binding domain free-energy surface underlies partial agonism at NMDA receptors

NMDA receptors are ligand-gated ion channels that mediate excitatory synaptic transmission in the central nervous system. Partial agonists elicit submaximal channel activation, but crystal structures of the ligand-binding domains (LBDs) bound with partial and full agonists show little difference. To uncover the molecular mechanism for partial agonism, here we computed the free-energy surfaces of the GluN1 (an obligatory subunit of NMDA receptors) LBD bound with a variety of ligands. The free-energy minima are similarly positioned for full and partial agonists, but the curvatures are significantly reduced in the latter case, indicating higher probabilities for sampling conformations with a not fully closed domain cleft. The free-energy surfaces for antagonists have both shifted minima and further reduced curvatures. Reduced curvature of free-energy surface appears to explain well the partial agonism at NMDA receptors and may present a unique paradigm in producing graded responses for receptors in general.

Structure of an agonist-bound ionotropic glutamate receptor

Ionotropic glutamate receptors (iGluRs) mediate most excitatory neurotransmission in the central nervous system and function by opening their ion channel in response to binding of agonist glutamate. Here, we report a structure of a homotetrameric rat GluA2 receptor in complex with partial agonist (S)-5-nitrowillardiine. Comparison of this structure with the closed-state structure in complex with competitive antagonist ZK 200775 suggests conformational changes that occur during iGluR gating. Guided by the structures, we engineered disulfide cross-links to probe domain interactions that are important for iGluR gating events. The combination of structural information, kinetic modeling, and biochemical and electrophysiological experiments provides insight into the mechanism of iGluR gating.

Structure and gating of tetrameric glutamate receptors

Ionotropic glutamate receptors (iGluRs) are ligand-gated ion channels that open their ion-conducting pores in response to the binding of agonist glutamate. In recent years, significant progress has been achieved in studies of iGluRs by determining numerous structures of isolated water-soluble ligand-binding and amino-terminal domains, as well as solving the first crystal structure of the full-length AMPA receptor in the closed, antagonist-bound state. These structural data combined with electrophysiological and fluorescence recordings, biochemical experiments, mutagenesis and molecular dynamics simulations have greatly improved our understanding of iGluR assembly, activation and desensitization processes. This article reviews the recent structural and functional advances in the iGluR field and summarizes them in a simplified model of full-length iGluR gating.

MMGBSA as a tool to understand the binding affinities of filamin-peptide interactions

Filamins (FLN) are large dimeric proteins that cross-link actin and work as important scaffolds in human cells. FLNs consist of an N-terminal actin-binding domain followed by 24 immunoglobulin-like domains (FLN1-24). FLN domains are divided into four subgroups based on their amino acid sequences. One of these subgroups, including domains 4, 9, 12, 17, 19, 21, and 23, shares a similar ligand-binding site between the β strands C and D. Several proteins, such as integrins β2 and β7, glycoprotein Ibα (GPIbα), and migfilin, have been shown to bind to this site. Here, we computationally estimated the binding free energies of filamin A (FLNa) subunits with bound peptides using the molecular mechanics-generalized Born surface area (MMGBSA) method. The obtained computational results correlated well with the experimental data, and they ranked efficiently both the binding of one ligand to all used FLNa-domains and the binding of all used ligands to FLNa21. Furthermore, the steered molecular dynamics (SMD) simulations pinpointed the binding hot spots for these complexes. These results demonstrate that molecular dynamics combined with free energy calculations are applicable to estimating the energetics of protein-protein interactions and can be used to direct the development of novel FLN function modulators.

Structural mechanism of N-methyl-D-aspartate receptor type 1 partial agonism

N-methyl-D-aspartate (NMDA) receptors belong to a family of ionotropic glutamate receptors that contribute to the signal transmission in the central nervous system. NMDA receptors are heterotetramers that usually consist of two GluN1 and GluN2 monomers. The extracellular ligand-binding domain (LBD) of a monomer is comprised of discontinuous segments that form the functional domains D1 and D2. While the binding of a full agonist glycine to LBD of GluN1 is linked to cleft closure and subsequent ion-channel opening, partial agonists are known to activate the receptor only sub-maximally. Although the crystal structures of the LBD of related GluA2 receptor explain the mechanism for the partial agonism, structures of GluN1-LBD cannot distinguish the difference between full and partial agonists. It is, however, probable that the partial agonists of GluN1 alter the structure of the LBD in order to result in a different pharmacological response than seen with full agonists. In this study, we used molecular dynamics simulations to reveal an intermediate closure-stage for GluN1, which is unseen in crystal structures. According to our calculations, this intermediate closure is not a transient stage but an energetically stable conformation. Our results demonstrate that the partial agonist cannot exert firm GluN1-LBD closure, especially if there is even a small force that disrupts the LBD closure. Accordingly, this result suggests the importance of forces from the ion channel for the relationship between pharmacological response and the structure of the LBD of members of this receptor family.

Functional insights from glutamate receptor ion channel structures

X-ray crystal structures for the soluble amino-terminal and ligand-binding domains of glutamate receptor ion channels, combined with a 3.6-Å-resolution structure of the full-length AMPA receptor GluA2 homotetramer, provide unique insights into the mechanisms of the assembly and function of glutamate receptor ion channels. Increasingly sophisticated biochemical, computational, and electrophysiological experiments are beginning to reveal the mechanism of action of partial agonists and suggest new models for the mechanism of action of allosteric modulators. Newly identified NMDA receptor ligands acting at novel sites offer hope for the development of subtype-selective modulators. The many unresolved issues include the role of the amino-terminal domain in AMPA receptor signaling and the mechanisms by which auxiliary proteins regulate receptor activity. The structural basis for ion permeation and ion channel block also remain areas of uncertainty, and despite substantial progress, molecular dynamics simulations have yet to reveal how glutamate binding opens the ion channel pore.