Identification of potential Gly/NMDA receptor antagonists by cheminformatics approach: a combination of pharmacophore modelling, virtual screening and molecular docking studies

Pharmacophore-Based Screening, Molecular Docking, and Dynamic Simulation of Fungal Metabolites as Inhibitors of Multi-Targets in Neurodegenerative Disorders

Neurodegenerative disorders, such as Alzheimer's disease (AD), negatively affect the economic and psychological system. For AD, there is still a lack of disease-altering treatments and promising cures due to its complex pathophysiology. In this study, we computationally screened the natural database of fungal metabolites against three known therapeutic target proteins of AD. Initially, a pharmacophore-based, drug-likeness category was employed for screening, and it filtered the 14 (A-N) best hits out of 17,544 fungal metabolites. The 14 best hits were docked individually against GSK-3β, the NMDA receptor, and BACE-1 to investigate the potential of finding a multitarget inhibitor. We found that compounds B, F, and L were immuno-toxic, whereas E, H, I, and J had a higher LD50 dose (5000 mg/kg). Among the examined metabolites, the Bisacremine-C (compound I) was found to be the most active molecule against GSK-3β (ΔG: -8.7 ± 0.2 Kcal/mol, Ki: 2.4 × 106 M-1), NMDA (ΔG: -9.5 ± 0.1 Kcal/mol, Ki: 9.2 × 106 M-1), and BACE-1 (ΔG: -9.1 ± 0.2 Kcal/mol, Ki: 4.7 × 106 M-1). It showed a 25-fold higher affinity with GSK-3β, 6.3-fold higher affinity with NMDA, and 9.04-fold higher affinity with BACE-1 than their native ligands, respectively. Molecular dynamic simulation parameters, such as RMSD, RMSF, Rg, and SASA, all confirmed that the overall structures of the targeted enzymes did not change significantly after binding with Bisacremine-C, and the ligand remained inside the binding cavity in a stable conformation for most of the simulation time. The most significant hydrophobic contacts for the GSK-3β-Bisacremine-C complex are with ILE62, VAL70, ALA83, and LEU188, whereas GLN185 is significant for H-bonds. In terms of hydrophobic contacts, TYR184 and PHE246 are the most important, while SER180 is vital for H-bonds in NMDA-Bisacremine-C. THR232 is the most crucial for H-bonds in BACE-1-Bisacremine-C and ILE110-produced hydrophobic contacts. This study laid a foundation for further experimental validation and clinical trials regarding the biopotency of Bisacremine-C.

Novel alkyl-substituted 4-methoxy benzaldehyde thiosemicarbazones: Multi-target directed ligands for the treatment of Alzheimer's disease

Alzheimer's disease (AD) is a devastating neurodegenerative disorder affecting mental ability and interrupts neurocognitive functions. Treating multifactorial conditions of AD with a single-target-directed drug is highly difficult. Thus, a multi-target-directed ligand (MTDL) development strategy has been developed as a promising approach for the treatment of AD. Herein, we have synthesized two novel thiosemicarbazones as MTDLs and reported their bioactivities against diverse neuropathological events involved in AD. In vitro studies revealed that both compounds exhibited promising anticholinesterase activity (AChE, IC50 = 15.98 μM, MZET and IC50 = 30.23 μM, MZMT), well supported by a detailed computational study. Both analogs have shown good thermodynamic behaviour and stability through interactions with characteristic amino acid residues throughout simulation of 100 ns against acetylcholinesterase enzyme. In an electrophysiology assay, these analogs have shown a characteristic inhibitory response against the GluN1-1a + GluN2B subunit of N-methyl-D-aspartate receptors. Pre-treatment of BV-2 microglial cells with MZET effectively decreased nitrite production compared to nitrite produced by lipopolysaccharide-treated cells alone. Further, the effect of MZMT and MZET on autophagy regulation was determined using stably transfected SH-SY5Y neuroblastoma cells. MZET significantly enhanced the autophagy flux in neuroblastoma cells. A significant decrease in copper-catalysed oxidation of amyloid-β in presence of synthesized thiosemicarbazones was also observed. Collectively, our findings indicated that these analogs have potential as effective anti-AD candidates and can be used as a prototype to develop more safer multi-targeted anti-AD drugs.

Design, synthesis, and pharmacological evaluation of [1, 3] dioxolo-chromeno[2,3-b]pyridines as anti-seizure agents

An efficient one-pot three-component reaction for the synthesis of [1,3]dioxolo[4',5':6,7]chromeno[2,3-b]pyridines 4(a-i) has been developed. Synthesis was achieved by reacting sesamol (1), aromatic aldehydes 2(a-i), and 2-aminopropene-1,1,3-tricarbonitrile (3) in the presence of triethylamine at 100 °C under neat reaction condition. Simple operational procedure, broad substrate scope, column chromatography free separations, and high yield of products make it an efficient and largely acceptable synthetic strategy. Synthesized compounds 4(a-i) were further screened for preliminary anticonvulsant activity using MES and scPTZ tests. These analogs were also checked for neurotoxicity and hepatotoxicity. Selected active compounds have been then screened quantitatively to determine ED50 and TD50 values. Analog 4h was found effective in both preclinical seizure models with significant therapeutic/toxicity profile (4h: ED50 = 34.7 mg/kg, MES test; ED50 = 37.9 mg/kg, scPTZ test; TD50 = 308.7 mg/kg). Molecular dynamic simulation for 100 ns of compound 4h-complexed with GABAA receptor revealed good thermodynamic behavior and fairly stable interactions (4h, Docking score =  - 10.94). In conclusion, effective synthetic strategy, significant anticonvulsant activity with good toxicity profile and detailed molecular modeling studies led us to anticipate the emergence of these analogs as valid leads for the development of future effective neurotherapeutic agents.

Multitargeted Virtual Screening and Molecular Simulation of Natural Product-like Compounds against GSK3β, NMDA-Receptor, and BACE-1 for the Management of Alzheimer's Disease

The complexity of Alzheimer's disease (AD) and several side effects of currently available medication inclined us to search for a novel natural cure by targeting multiple key regulatory proteins. We initially virtually screened the natural product-like compounds against GSK3β, NMDA receptor, and BACE-1 and thereafter validated the best hit through molecular dynamics simulation (MDS). The results demonstrated that out of 2029 compounds, only 51 compounds exhibited better binding interactions than native ligands, with all three protein targets (NMDA, GSK3β, and BACE) considered multitarget inhibitors. Among them, F1094-0201 is the most potent inhibitor against multiple targets with binding energy -11.7, -10.6, and -12 kcal/mol, respectively. ADME-T analysis results showed that F1094-0201 was found to be suitable for CNS drug-likeness in addition to their other drug-likeness properties. The MDS results of RMSD, RMSF, Rg, SASA, SSE and residue interactions indicated the formation of a strong and stable association in the complex of ligands (F1094-0201) and proteins. These findings confirm the F1094-0201's ability to remain inside target proteins' binding pockets while forming a stable complex of protein-ligand. The free energies (MM/GBSA) of BACE-F1094-0201, GSK3β-F1094-0201, and NMDA-F1094-0201 complex formation were -73.78 ± 4.31 kcal mol^-1, -72.77 ± 3.43 kcal mol^-1, and -52.51 ± 2.85 kcal mol^-1, respectively. Amongst the target proteins, F1094-0201 have a more stable association with BACE, followed by NMDA and GSK3β. These attributes of F1094-0201 indicate it as a possible option for the management of pathophysiological pathways associated with AD.

QSAR, simulation techniques, and ADMET/pharmacokinetics assessment of a set of compounds that target MAO-B as anti-Alzheimer agent

Background

Alzheimer’s disease (AD), the most common cause of dementia in the elderly, is a progressive neurodegenerative disorder that gradually affects cognitive function and eventually causes death. Most approved drugs can only treat the disease alleviating the disease symptoms; therefore, there is a need to develop drugs that can treat this illness holistically. The medical community is searching for new drugs and new drug targets to cure this disease. In this study, QSAR, molecular docking evaluation, and ADMET/pharmacokinetics assessment were used as modeling methods to identify the compounds with outstanding physicochemical properties.

Results

The 37 MAO-B compounds were screened using the aforementioned methods and yielded a model with the following molecular properties: AATS1v, AATS3v, GATS4m, and GATS6e. Good statistical values were R2train = 0.69, R2adj = 0.63, R2pred = 0.57, LOF = 0.23, and RMSE = 0.38. The model was validated using an evaluation set that confirmed its robustness. The molecular docking was also utilized using crystal structure of human monoamine oxidase B in complex with chlorophenylchromone-carboxamide with ID code of 6FW0, and three compounds were identified with outstanding high binding affinity (13 = − 30.51 kcal mol−1, 31 = − 31.85 kcal mol−1, and 33 = − 33.70 kcal mol−1), and better than the Eldepryl (referenced) drug (− 11.40 kcal mol−1).

Conclusions

These three compounds (13, 31, and 33) were analyzed for ADMET/pharmacokinetics evaluation and found worthy of further analysis as promising drug candidates to cure AD and could also serve as a template to design several monoamine oxidase B inhibitors in the future to cure AD.

Computational identification of 2,4-disubstituted amino-pyrimidines as L858R/T790M-EGFR double mutant inhibitors using pharmacophore mapping, molecular docking, binding free energy calculation, DFT study and molecular dynamic simulation

Pharmacophore modelling studies have been performed for a series of 2,4-disubstituted-pyrimidines derivatives as EGFR L858R/T790M tyrosine kinase inhibitors. The high scoring AARR.15 hypothesis was selected as the best pharmacophore model with the highest survival score of 3.436 having two hydrogen bond acceptors and two aromatic ring features. Pharmacophore-based virtual screening followed by structure-based yielded the six molecules (ZINC17013227, ZINC17013215, ZINC9573324, ZINC9573445, ZINC24023331 and ZINC17013503) from the ZINC database with significant in silico predicted activity and strong binding affinity towords the EGFR L858R/T790M tyrosine kinase. In silico toxicity and cytochrome profiling indicates that all the 06 virtually screened compounds were substrate/inhibitors of the CYP-3A4 metabolizing enzyme and were non-carcinogenic and devoid of Ames mutagenesis. Density functional theory (DFT) and molecular dynamic (MD) simulation further validated the obtained hits.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-021-00113-x.

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.

GluN2B/N-methyl-D-aspartate Receptor Antagonists: Advances in Design, Synthesis, and Pharmacological Evaluation Studies

Selective GluN2B/N-methyl-D-aspartate receptor (NMDAR) antagonists have exposed their clinical effectiveness in a cluster of neurodegenerative diseases, such as epilepsy, Alzheimer's disease, Parkinson's disease, pain, and depression. Hence, GluN2B/NMDARs are considered to be a prospective target for the management of neurodegenerative diseases. Here, we have discussed the current results and significance of subunit selective GluN2B/NMDAR antagonists to pave the way for the establishment of new, safe, and economical drug candidates in the near future. By using summarized data of selective GluN2B/NMDAR antagonists, medicinal chemists are certainly a step closer to the goal of improving the therapeutic and side effect profile of selective antagonists. Outlined summary of designing strategies, synthetic schemes, and pharmacological evaluation studies reinvigorate efforts to identify, modify, and synthesize novel GluN2B/NMDAR antagonists for treating neurodegenerative diseases.

Molecular docking utilising the OliveNet™ library reveals novel phenolic compounds which may potentially target key proteins associated with major depressive disorder

The antidepressant medications that are currently prescribed to patients suffering from major depressive disorder (MDD) have limitations and as a result, there is an urgent need to increase the options that are available. A number of studies have found that natural polyphenols have neuroprotective properties and there is evidence to suggest that they modulate neurotransmitter systems. There are more than 200 phenolic compounds that have been identified in Olea europaea, many of which have not yet been investigated for their potential biological effects. In this study, in silico methods were used to screen the phenolic library from the OliveNet™ database and identify novel lead compounds for proteins implicated in the pathophysiology of MDD. The molecular docking results revealed that the monoamine oxidase enzyme isoforms (MAO-A/MAO-B) had binding specificities for certain phenolic subclasses. The lead ligands that were identified from these subclasses were positioned near the flavin adenine dinucleotide (FAD) cofactor, interacting in a similar manner as known inhibitors. In addition to the MAO enzymes, several phenolic compounds were docked to neurotransmitter transporters and postsynaptic receptors, as well as proteins involved in neuroinflammation, oxidative stress and the endocannabinoid system. Based on the binding affinity, position, orientation and interactions of the lead phenolic compounds identified in this study, it is predicted that they may have antidepressant properties. The results should be validated further using molecular dynamics (MD) simulations, as well as in vivo and in vitro techniques.

Current and emerging therapeutic targets of alzheimer's disease for the design of multi-target directed ligands

Alzheimer's disease (AD) is the most prevalent neurodegenerative disease, and a major cause of death worldwide. The number of people suffering from this debilitating disorder is rising at an unprecedented rate, with a subsequent surge in healthcare costs. Only four drugs are clinically available for the treatment of AD symptoms, but they are not disease-modifying. Consequently, there is an urgent need for a cure. Although the cause of this debilitating condition remains poorly understood, it is believed that several factors may be involved in combination - including, health and lifestyle, environmental, and genetic factors. In recent years, a number of hallmarks of the disease have also been discovered, and it is believed that these factors may play an important role in the development of AD. Amyloid aggregation is one such factor which has been highly investigated, in addition to cholinesterase enzymes and tau aggregation. In the last decade, multi-target drugs have been increasingly investigated for their application to AD treatment. By combining two or more pharmacophores in a single compound, it is possible to synthesise a drug which can target several factors that are involved in AD development. This is a particularly attractive approach as it would avoid the use of combination therapies. As a result, it could reduce the burden on carers and families, and decrease healthcare and social care costs. Many active pharmacophores have been employed for the development of hybrid drugs, due to their abilities to inhibit the factors currently widely recognised to be involved in AD. These compounds have demonstrated promising results; however, research is still required to optimise the pharmacological profiles of the drugs, in addition to their potencies. Meanwhile, extensive research is continuously being performed into other potential targets for the treatment of AD. Based on the results obtained thus far, it is likely that multi-target compounds will continue to be increasingly studied in the future as potential treatments for AD.

Structure based virtual screening of novel noncompetitive antagonist of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor

The α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) subtype ionotropic glutamate receptors are attractive antiepileptic targets responsible for mediating the majority of excitatory neurotransmission and plasticity. The noncompetitive antagonists obtain more and more attention as drug candidates for treatment of the neurological diseases involving excessive activity of AMPARs, due to they regulate AMPA receptors (AMPARs) activity independently of endogenous glutamate levels unlike the competitive antagonists. Development of novel AMPAR noncompetitive antagonists, which are safer and more efficacious than competitive antagonists, is highly under demand. Here, we present the discovery of novel antagonists against AMPAR through Structure-Based Virtual Screening (SBVS). Three compounds were successfully distinguished by several different filtering strategies, namely STOCK6S-10902, STOCK1N-49134 and STOCK5S-68665. The interaction mode of these compounds was further explored through molecular dynamics simulation, binding free energy calculation and the binding free energy decomposition. It is demonstrated that some residues within the binding pocket, which have been proved their importance in antagonist binding and gating, form strong hydrogen bond interactions with these three molecules. In particular, H-bond interactions with high occupancies between Ser516, Ser788 and STOCK6S-10902 and Ser516, Asn791 and STOCK1N-49134 were observed. The three hit compounds with new scaffolds and the detailed binding modes could potentially serve as a starting point for further optimization and development.

Xylazine Regulates the Release of Glycine and Aspartic Acid in Rat Brain

Introduction

Xylazine, a type of α₂-adrenoceptors, is a commonly used drug in veterinary medicine. Xylazine-induced changes in the content of amino acid neurotransmitters - glycine (Gly) and aspartic acid (Asp), in different brain regions and neurons were studied.

Material and Methods

Wistar rats were administered 50 mg/kg or 70 mg/kg of xylazine by intraperitoneal injection. In addition, in vitro experiments were conducted, in which neurons were treated with 15 μg/mL, 25 μg/mL, 35μg/mL, and 45 μg/mL of xylazine. Test methods were based on the enzyme-linked immunosorbent assays (ELISA).

Results

During anaesthesia, Asp levels in each brain area were significantly lower compared to the control group. Except for the cerebrum, levels of Gly in other brain areas were significantly increased during the anaesthesia period. In vitro, xylazine-related neuron secretion of Gly increased significantly compared to the control group at 60 min and 90 min. Moreover, xylazine caused a significant decrease in the levels of Asp secreted by neurons at 20 min, but gradually returned to the level of the control group.

Conclusion

The data showed that during anaesthesia the overall levels of Asp decreased and overall levels of Gly increased. In addition, the inhibitory effect of xylazine on Asp and the promotion of Gly were dose-dependent. Our data showed that different effects of xylazine on excitatory and inhibitory neurotransmitters provided a theoretical basis for the mechanism of xylazine activity in clinical anaesthesia.

Xylazine regulates the release of glycine and aspartic acid in rat brain

Introduction

Xylazine, a type of α2-adrenoceptors, is a commonly used drug in veterinary medicine. Xylazine-induced changes in the content of amino acid neurotransmitters – glycine (Gly) and aspartic acid (Asp), in different brain regions and neurons were studied.

Material and Methods

Wistar rats were administered 50 mg/kg or 70 mg/kg of xylazine by intraperitoneal injection. In addition, in vitro experiments were conducted, in which neurons were treated with 15 μg/mL, 25 μg/mL, 35μg/mL, and 45 μg/mL of xylazine. Test methods were based on the enzyme-linked immunosorbent assays (ELISA).

Results

During anaesthesia, Asp levels in each brain area were significantly lower compared to the control group. Except for the cerebrum, levels of Gly in other brain areas were significantly increased during the anaesthesia period. In vitro, xylazine-related neuron secretion of Gly increased significantly compared to the control group at 60 min and 90 min. Moreover, xylazine caused a significant decrease in the levels of Asp secreted by neurons at 20 min, but gradually returned to the level of the control group.

Conclusion

The data showed that during anaesthesia the overall levels of Asp decreased and overall levels of Gly increased. In addition, the inhibitory effect of xylazine on Asp and the promotion of Gly were dose-dependent. Our data showed that different effects of xylazine on excitatory and inhibitory neurotransmitters provided a theoretical basis for the mechanism of xylazine activity in clinical anaesthesia.

Development of Pharmacophore Model for Indeno[1,2-b]indoles as Human Protein Kinase CK2 Inhibitors and Database Mining

Protein kinase CK2, initially designated as casein kinase 2, is an ubiquitously expressed serine/threonine kinase. This enzyme, implicated in many cellular processes, is highly expressed and active in many tumor cells. A large number of compounds has been developed as inhibitors comprising different backbones. Beside others, structures with an indeno[1,2-b]indole scaffold turned out to be potent new leads. With the aim of developing new inhibitors of human protein kinase CK2, we report here on the generation of common feature pharmacophore model to further explain the binding requirements for human CK2 inhibitors. Nine common chemical features of indeno[1,2-b]indole-type CK2 inhibitors were determined using MOE software (Chemical Computing Group, Montreal, Canada). This pharmacophore model was used for database mining with the aim to identify novel scaffolds for developing new potent and selective CK2 inhibitors. Using this strategy several structures were selected by searching inside the ZINC compound database. One of the selected compounds was bikaverin (6,11-dihydroxy-3,8-dimethoxy-1-methylbenzo[b]xanthene-7,10,12-trione), a natural compound which is produced by several kinds of fungi. This compound was tested on human recombinant CK2 and turned out to be an active inhibitor with an IC₅₀ value of 1.24 µM.

Identification of novel PfDHODH inhibitors as antimalarial agents via pharmacophore-based virtual screening followed by molecular docking and in vivo antimalarial activity

Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH) catalyses the fourth reaction of de novo pyrimidine biosynthesis in parasites, and represents an important target for the treatment of malaria. In this study, we describe pharmacophore-based virtual screening combined with docking study and biological evaluation as a rational strategy for identification of novel hits as antimalarial agents. Pharmacophore models were established from known PfDHODH inhibitors using the GALAHAD module with IC50 values ranging from 0.033 μM to 142 μM. The best pharmacophore model consisted of three hydrogen bond acceptor, one hydrogen bond donor and one hydrophobic features. The pharmacophore models were validated through receiver operating characteristic and Günere-Henry scoring methods. The best pharmacophore model as a 3D search query was searched against the IBS database. Several compounds with different structures (scaffolds) were retrieved as hit molecules. Among these compounds, those with a QFIT value of more than 81 were docked in the PfDHODH enzyme to further explore the binding modes of these compounds. In silico pharmacokinetic and toxicities were predicted for the best docked molecules. Finally, the identified hits were evaluated in vivo for their antimalarial activity in a parasite inhibition assay. The hits reported here showed good potential to become novel antimalarial agents.