Molecular dynamics study-based mechanism of nefiracetam-induced NMDA receptor potentiation

Benzimidazole compound abrogates SARS-COV-2 receptor-binding domain (RBD)/ACE2 interaction In vitro

The development of clinically actionable pharmaceuticals against coronavirus disease (COVID-19); an infectious disease caused by the SARS-CoV-2 virus is very important for ending the pandemic. Coronavirus spike glycoprotein (GP)-Receptor Binding Domain (RBD) and its interaction with host receptor angiotensin converting enzyme 2 (ACE2) is one of the most structurally understood but therapeutically untapped aspect of COVID-19 pathogenesis. Binding interface based on previous x-ray structure of RBD/ACE2 were virtually screened to identify fragments with high-binding score from 12,000 chemical building blocks. The hit compound was subjected to fingerprint-based similarity search to identify compounds within the FDA-approved drug library containing the same core scaffold. Identified compounds were then re-docked into of RBD/ACE2. The best ranked compound was validated for RBD/ACE2 inhibition using commercial kit. Molecular dynamics simulation was conducted to provide further insight into the mechanism of inhibition. From the original 12000 chemical building blocks, benzimidazole (BAZ) scaffold was identified. Fingerprint-based similarity search of the FDA-approved drug library for BAZ-containing compounds identified 12 drugs with the benzimidazole-like substructure. When these compounds were re-docked into GP/ACE2 interface, the consensus docking identified bazedoxifene as the hit. In vitro RBD/ACE2 inhibition kinetics showed micromolar IC50 value (1.237 μM) in the presence of bazedoxifene. Molecular dynamics simulation of RBD/ACE2 in the presence BAZ resulted in loss of contact and specific hydrogen-bond interaction required for RBD/ACE2 stability. Taken together, these findings identified benzimidazole scaffold as a building block for developing novel RBD/ACE2 complex inhibitor and provided mechanistic basis for the use of bazedoxifene as a repurposable drug for the treatment of COVID-19 acting at RBD/ACE2 interface.

Repurposing of cefpodoxime proxetil as potent neuroprotective agent through computational prediction and in vitro validation

In recent reports, NR2B-NMDA receptor antagonists showed more research value because of its strong targeting ability and less side effects potential. In 2016, EVT-101 was reported to bind in an almost entirely new binding region of this target. Whether strikingly different binding modes can improve targeting and reduce side effects is worth studying. In our preliminary work, we explored the binding patterns of ifenprodil and EVT-101, found the key amino acids and summarized the pharmacophores, hoping to find such antagonists that target the two binding modes simultaneously. In this study, we developed a scalable virtual screening workflow in the FDA-approved drugs library to identify novel NR2B-NMDAR antagonists based on the combination of two pharmacophores. Cefpodoxime proxetil (5) was identified as the hit compound, and it was found for the first time that 5 might have neuroprotective activity as a NR2B-NMDAR antagonist. This result interested us to make further study, the ligand-receptor interactions modeled by molecular docking studies showed that the compound could perfectly merge both the pharmacophore characteristics of ifenprodil and EVT-101 at the binding cavity between the ATDs of GluN1 and GluN2B. The accuracy of molecular docking results and binding stability of ligand-receptor complexes were validated through 100 ns molecular dynamics simulation and binding free energy calculation. Afterwards, MTT assay (49.8%±0.1%, 5 μM) on NMDA injured SH-SY5Y cells and evidence of the effect on attenuating Ca2+ influx induced by NMDA were applied to validate the computational results, further investigation showed that 5 could suppress the NR2B upregulation induced by NMDA. [Formula: see text] Communicated by Ramaswamy H. Sarma.

In silico and in vivo study of radio-iodinated nefiracetam as a radiotracer for brain imaging in mice

This study presents development and characterization of a radiotracer, [125I]iodonefiracetam ([125I]iodoNEF). Labeling with high yield and radiochemical purity was achieved through the formation of a [125I]iodoNEF radiotracer after investigating many factors like oxidizing agent content (chloramines-T (Ch-T)), substrate amount (Nefiracetam (NEF)), pH of reaction mixture, reaction time and temperature. Nefiracetam (NEF) is known as nootropic agent, acting as N-methyl-d-aspartic acid receptor ligand (NMDA). The radiolabeled compound was stable, and exhibited the logarithm of the partition coefficient (log p) value of [125I]iodonefiracetam as 1.85 (lipophilic). Biodistribution studies in normal mice confirmed the suitability of the [125I]iodoNEF radiotracer as a novel tracer for brain imaging. High uptake of 8.61 ± 0.14 percent injected dose/g organ was observed in mice

NMDA Receptor Opening and Closing-Transitions of a Molecular Machine Revealed by Molecular Dynamics

We report the first complete description of the molecular mechanisms behind the transition of the N-methyl-d-aspartate (NMDA) receptor from the state where the transmembrane domain (TMD) and the ion channel are in the open configuration to the relaxed unliganded state where the channel is closed. Using an aggregate of nearly 1 µs of unbiased all-atom implicit membrane and solvent molecular dynamics (MD) simulations we identified distinct structural states of the NMDA receptor and revealed functionally important residues (GluN1/Glu522, GluN1/Arg695, and GluN2B/Asp786). The role of the "clamshell" motion of the ligand binding domain (LBD) lobes in the structural transition is supplemented by the observed structural similarity at the level of protein domains during the structural transition, combined with the overall large rearrangement necessary for the opening and closing of the receptor. The activated and open states of the receptor are structurally similar to the liganded crystal structure, while in the unliganded receptor the extracellular domains perform rearrangements leading to a clockwise rotation of up to 45 degrees around the longitudinal axis of the receptor, which closes the ion channel. The ligand-induced rotation of extracellular domains transferred by LBD-TMD linkers to the membrane-anchored ion channel is responsible for the opening and closing of the transmembrane ion channel, revealing the properties of NMDA receptor as a finely tuned molecular machine.

Morelloflavone as a novel inhibitor of mitotic kinesin Eg5

Among 40 plant-derived biflavonoids with inhibitory potential against Eg5, morelloflavone from Garcinia dulcis leaves was selected for further testing based on in silico analysis of binding modes, molecular interactions, binding energies and functional groups that interact with Eg5. Computational models predicted that morelloflavone binds the putative allosteric pocket of Eg5, within the cavity surrounded by amino acid residues of Ile-136, Glu-116, Glu-118, Trp-127, Gly-117, Ala-133, Glu-215, Leu-214 and Tyr-211. Binding energy was -8.4 kcal/mol, with a single hydrogen bond formed between morelloflavone and Tyr-211. The binding configuration was comparable to that of a reference inhibitor, S-trityl-L-cysteine. Subsequent biochemical analysis in vitro confirmed that morelloflavone inhibited both the basal and microtubule-activated ATPase activity of Eg5 in a manner that does not compete with ATP binding. Morelloflavone also suppressed Eg5 gliding along microtubules. These results suggest that morelloflavone binds the allosteric binding site in Eg5 and thereby inhibits ATPase activity and motor function of Eg5.

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.

Glutamate and Glycine Binding to the NMDA Receptor

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.

Flavonoid-rich extract of Chromolaena odorata modulate circulating GLP-1 in Wistar rats: computational evaluation of TGR5 involvement

Chromolaena odorata is a major bio-resource in folkloric treatment of diabetes. In the present study, its anti-diabetic component and underlying mechanism were investigated. A library containing 140 phytocompounds previously characterized from C. odorata was generated and docked (Autodock Vina) into homology models of dipeptidyl peptidase (DPP)-4, Takeda-G-protein-receptor-5 (TGR5), glucagon-like peptide 1 (GLP1) receptor, renal sodium dependent glucose transporter (SGLUT)-1/2 and nucleotide-binding oligomerization domain (NOD) proteins 1&2. GLP-1 gene (RT-PCR) modulation and its release (EIA) by C. odorata were confirmed in vivo. From the docking result above, TGR5 was identified as a major target for two key C. odorata flavonoids (5,7-dihydroxy-6-4-dimethoxyflavanone and homoesperetin-7-rutinoside); sodium taurocholate and C. odorata powder included into the diet of the animals both raised the intestinal GLP-1 expression versus control (p < 0.05); When treated with flavonoid-rich extract of C. odorata (CoF) or malvidin, circulating GLP-1 increased by 130.7% in malvidin-treated subjects (0 vs. 45 min). CoF treatment also resulted in 128.5 and 275% increase for 10 and 30 mg/kg b.w., respectively.

CONCLUSIONS

The results of this study support that C. odorata flavonoids may modulate the expression of GLP-1 and its release via TGR5. This finding may underscore its anti-diabetic potency.

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.

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.

NMDA receptor agonists reverse impaired psychomotor and cognitive functions associated with hippocampal Hbegf-deficiency in mice

BACKGROUND

Structural and functional changes of the hippocampus are correlated with psychiatric disorders and cognitive dysfunctions. Genetic deletion of heparin-binding epidermal growth factor-like growth factor (HB-EGF), which is predominantly expressed in cortex and hippocampus, also causes similar psychiatric and cognitive dysfunctions, accompanying down-regulated NMDA receptor signaling. However, little is known of such dysfunctions in hippocampus-specific Hbegf cKO mice.

RESULTS

We successfully developed hippocampus-specific cKO mice by crossbreeding floxed Hbegf and Gng7-Cre knock-in mice, as Gng7 promoter-driven Cre is highly expressed in hippocampal neurons as well as striatal medium spiny neurons. In mice lacking hippocampus Hbegf gene, there was a decreased neurogenesis in the subgranular zone (SGZ) of the dentate gyrus as well as down-regulation of PSD-95/NMDA-receptor-NR1/NR2B subunits and related NMDA receptor signaling. Psychiatric, social-behavioral and cognitive abnormalities were also observed in hippocampal cKO mice. Interestingly, D-cycloserine and nefiracetam, positive allosteric modulators (PAMs) of NMDA receptor reversed the apparent reduction in NMDA receptor signaling and most behavioral abnormalities. Furthermore, decreased SGZ neurogenesis in hippocampal cKO mice was reversed by nefiracetam.

CONCLUSIONS

The present study demonstrates that PAMs of NMDA receptor have pharmacotherapeutic potentials to reverse down-regulated NMDA receptor signaling, neuro-socio-cognitive abnormalities and decreased neurogenesis in hippocampal cKO mice.