Structural studies, homology modeling and molecular docking of novel non-competitive antagonists of GluK1/GluK2 receptors

Kainate Receptor Antagonists: Recent Advances and Therapeutic Perspective

Since the 1990s, ionotropic glutamate receptors have served as an outstanding target for drug discovery research aimed at the discovery of new neurotherapeutic agents. With the recent approval of perampanel, the first marketed non-competitive antagonist of AMPA receptors, particular interest has been directed toward 'non-NMDA' (AMPA and kainate) receptor inhibitors. Although the role of AMPA receptors in the development of neurological or psychiatric disorders has been well recognized and characterized, progress in understanding the function of kainate receptors (KARs) has been hampered, mainly due to the lack of specific and selective pharmacological tools. The latest findings in the biology of KA receptors indicate that they are involved in neurophysiological activity and play an important role in both health and disease, including conditions such as anxiety, schizophrenia, epilepsy, neuropathic pain, and migraine. Therefore, we reviewed recent advances in the field of competitive and non-competitive kainate receptor antagonists and their potential therapeutic applications. Due to the high level of structural divergence among the compounds described here, we decided to divide them into seven groups according to their overall structure, presenting a total of 72 active compounds.

Experimental and Computational Structural Studies of 2,3,5-Trisubstituted and 1,2,3,5-Tetrasubstituted Indoles as Non-Competitive Antagonists of GluK1/GluK2 Receptors

The blockade of kainate receptors, in particular with non-competitive antagonists, has—due to their anticonvulsant and neuroprotective properties—therapeutic potential in many central nervous system (CNS) diseases. Deciphering the structural properties of kainate receptor ligands is crucial to designing medicinal compounds that better fit the receptor binding pockets. In light of that fact, here, we report experimental and computational structural studies of four indole derivatives that are non-competitive antagonists of GluK1/GluK2 receptors. We used X-ray studies and Hirshfeld surface analysis to determine the structure of the compounds in the solid state and quantum chemical calculations to compute HOMO and LUMO orbitals and the electrostatic potential. Moreover, non-covalent interaction maps were also calculated. It is worth emphasizing that compounds 3 and 4 are achiral molecules crystallising in non-centrosymmetric space groups, which is a relatively rare phenomenon.

Pharmacoproteomics of Brain Barrier Transporters and Substrate Design for the Brain Targeted Drug Delivery

One of the major reasons why central nervous system (CNS)-drug development has been challenging in the past, is the barriers that prevent substances entering from the blood circulation into the brain. These barriers include the blood-brain barrier (BBB), blood-spinal cord barrier (BSCB), blood-cerebrospinal fluid barrier (BCSFB), and blood-arachnoid barrier (BAB), and they differ from each other in their transporter protein expression and function as well as among the species. The quantitative expression profiles of the transporters in the CNS-barriers have been recently revealed, and in this review, it is described how they affect the pharmacokinetics of compounds and how these expression differences can be taken into account in the prediction of brain drug disposition in humans, an approach called pharmacoproteomics. In recent years, also structural biology and computational resources have progressed remarkably, enabling a detailed understanding of the dynamic processes of transporters. Molecular dynamics simulations (MDS) are currently used commonly to reveal the conformational changes of the transporters and to find the interactions between the substrates and the protein during the binding, translocation in the transporter cavity, and release of the substrate on the other side of the membrane. The computational advancements have also aided in the rational design of transporter-utilizing compounds, including prodrugs that can be actively transported without losing potency towards the pharmacological target. In this review, the state-of-art of these approaches will be also discussed to give insights into the transporter-mediated drug delivery to the CNS.

Assembly and Trafficking of Homomeric and Heteromeric Kainate Receptors with Impaired Ligand Binding Sites

Kainate receptors (KARs) are a subfamily of ionotropic glutamate receptors (iGluRs) mediating excitatory synaptic transmission. Cell surface expressed KARs modulate the excitability of neuronal networks. The transfer of iGluRs from the endoplasmic reticulum (ER) to the cell surface requires occupation of the agonist binding sites. Here we used molecular modelling to produce a range of ligand binding domain (LBD) point mutants of GluK1-3 KAR subunits with and without altered agonist efficacy to further investigate the role of glutamate binding in surface trafficking and activation of homomeric and heteromeric KARs using endoglycosidase digestion, cell surface biotinylation and imaging of changes in intracellular Ca²⁺ concentration [Ca²⁺]_i. Mutations of conserved amino acid residues in the LBD that disrupt agonist binding to GluK1-3 (GluK1-T675V, GluK2-A487L, GluK2-T659V and GluK3-T661V) reduced both the total expression levels and cell surface delivery of all of these mutant subunits compared to the corresponding wild type in transiently transfected human embryonic kidney 293 (HEK293) cells. In contrast, the exchange of non-conserved residues in the LBD that convert antagonist selectivity of GluK1-3 (GluK1-T503A, GluK2-A487T, GluK3-T489A, GluK1-N705S/S706N, GluK2-S689N/N690S, GluK3-N691S) did not alter the biosynthesis and trafficking of subunit proteins. Co-assembly of mutant GluK2 with an impaired LBD and wild type GluK5 subunits enables the cell surface expression of both subunits. However, [Ca²⁺]_i imaging indicates that the occupancy of both GluK2 and GluK5 LBDs is required for the full activation of GluK2/GluK5 heteromeric KAR channels.

Synthesis, Structural and Thermal Studies of 3-(1-Benzyl-1,2,3,6-tetrahydropyridin-4-yl)-5-ethoxy-1H-indole (D2AAK1_3) as Dopamine D₂ Receptor Ligand

Compound D2AAK1_3 was designed as a modification of the lead structure D2AAK1 (an in vivo active multi-target compound with nanomolar affinity to a number of aminergic GPCRs) and synthesized in the reaction of 5-ethoxyindole and 1-benzyl-4-piperidone. This compound has an affinity to the human dopamine D₂ receptor with K_i of 151 nM. The aim of these studies was the structural and thermal characterization of the compound D2AAK1_3. In particular; X-ray studies; molecular docking and molecular dynamics as well as thermal analysis were performed. The studied compound crystallizes in orthorhombic system; in chiral space group P2₁2₁2₁. The compound has a non-planar conformation. The studied compound was docked to the novel X-ray structure of the human dopamine D₂ receptor in the inactive state (PDB ID: 6CM4) and established the main contact between its protonatable nitrogen atom and Asp (3.32) of the receptor. The obtained binding pose was stable in molecular dynamics simulations. Thermal stability of the compound was investigated using the TG-DSC technique in the air atmosphere, while TG-FTIR analyses in air and nitrogen atmospheres were also performed. The studied compound is characterized by good thermal stability. The main volatile products of combustion are the following gases: CO₂; H₂O toluene and CO while in the case of pyrolysis process in the FTIR spectra; the characteristic bands of NH₃; piperidine and indole are additionally observed.

Comparative molecular field analysis and molecular dynamics studies of the dopamine D2 receptor antagonists without a protonatable nitrogen atom

The dopaminergic hypothesis of schizophrenia is the main concept explaining the direct reasons of schizophrenia and the effectiveness of current antipsychotics. All antipsychotics present on the market are potent dopamine D₂ receptor antagonists or partial agonists. In this work we investigate a series of dopamine D₂ receptor antagonists which do not fulfill the criteria of the classical pharmacophore model as they do not possess a protonatable nitrogen atom necessary to interact with the conserved Asp(3.32). Such compounds are interesting, inter alia, due to possible better pharmacokinetic profile when compared to basic, ionizable molecules. By means of homology modeling, molecular docking and molecular dynamics we determined that the compounds investigated interact with Asp(3.32) via their amide nitrogen atom. It was found that the studied compounds stabilize the receptor inactive conformation through the effect on the ionic lock, which is typical for GPCR antagonists. We constructed a CoMFA model for the studied compounds with the following statistics: R² = 0.95, Q² = 0.63. The quality of the CoMFA model was confirmed by high value of R² of the test set, equal 0.96. The CoMFA model indicated two regions where bulky substituents are favored and two regions where bulky substituents are not beneficial. Two red contour regions near carbonyl groups were identified meaning that negative charge would be favored here. Furthermore, the S-oxide group is connected with blue contour region meaning that positive charge is favored in this position. These findings may be applied for further optimization of the studied compound series.

The dopamine D2 receptor dimer and its interaction with homobivalent antagonists: homology modeling, docking and molecular dynamics

In order to apply structure-based drug design techniques to G protein-coupled receptor complexes, it is essential to model their 3D structure and to identify regions that are suitable for selective drug binding. For this purpose, we have developed and tested a multi-component protocol to model the inactive conformation of the dopamine D2 receptor dimer, suitable for interaction with homobivalent antagonists. Our approach was based on protein-protein docking, applying the Rosetta software to obtain populations of dimers as present in membranes with all the main possible interfaces. Consensus scoring based on the values and frequencies of best interfaces regarding four scoring parameters, Rosetta interface score, interface area, free energy of binding and energy of hydrogen bond interactions indicated that the best scored dimer model possesses a TM4-TM5-TM7-TM1 interface, which is in agreement with experimental data. This model was used to study interactions of the previously published dopamine D2 receptor homobivalent antagonists based on clozapine,1,4-disubstituted aromatic piperidines/piperazines and arylamidoalkyl substituted phenylpiperazine pharmacophores. It was found that the homobivalent antagonists stabilize the receptor-inactive conformation by maintaining the ionic lock interaction, and change the dimer interface by disrupting a set of hydrogen bonds and maintaining water- and ligand-mediated hydrogen bonds in the extracellular and intracellular part of the interface. Graphical Abstract Structure of the final model of the dopamine D2 receptor homodimer, indicating the distancebetween Tyr37 and Tyr 5.42 in the apo form (left) and in the complex with the ligand (right).

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

In vitro, molecular modeling and behavioral studies of 3-{[4-(5-methoxy-1H-indol-3-yl)-1,2,3,6-tetrahydropyridin-1-yl]methyl}-1,2-dihydroquinolin-2-one (D2AAK1) as a potential antipsychotic

Antipsychotics currently available to treat schizophrenia suffer several limitations: (1) they are efficient against positive but not negative and cognitive symptoms of the disease; (2) they help only a half of patients; (3) they have severe side effects including neurological and metabolic side effects. Thus, novel drugs to treat schizophrenia are highly demanded. We identified a novel dopamine D2 receptor antagonist, D2AAK1, with Ki of 58 nM using structure-based virtual screening. D2AAK1 possesses additional nanomolar or low micromolar affinity to D1, D3, 5-HT1A and 5-HT2A receptors, making it an ideal candidate for a multi-target drug. Here we present homology modeling, molecular docking and molecular dynamics of D2AAK1 and its molecular targets and animal studies of D2AAK1 as a potential antipsychotic. The main contact of D2AAK1 and all the receptors studied is the electrostatic interaction between the protonable nitrogen atom of the ligand and the conserved Asp(3.32) as typical for orthosteric ligands of aminergic GPCRs. We confirmed antagonistic/partial agonistic properties of D2AAK1 towards all the receptors in in vitro essays and in in silico studies as the ligand stabilizes the ionic lock interaction. We also demonstrated neuroleptic, anxiolytic and, importantly, procognitive properties of D2AAK1 in mouse models.

Synthesis and molecular docking of novel non-competitive antagonists of GluK2 receptor

Here we present the synthesis, pharmacological activity, and molecular docking of novel non-competitive antagonists of GluK2 receptor. The compounds concerned are derivatives of indole and carbazole and are the second reported series of non-competitive antagonists of the GluK2 receptor (the first one was also published by our group). The activity of the indole derivatives is in the micromolar range, as in the case of the first series of non-competitive GluK2 receptor antagonists. We have found that designed carbazole derivatives are devoid of activity. Active indole derivatives interact with the transduction domain of the GluK2 receptor, i.e., the domain which links the transmembrane region of the receptor with the agonist-binding domain. The binding pocket is situated within one receptor subunit.

Pharmacological and molecular studies on the interaction of varenicline with different nicotinic acetylcholine receptor subtypes. Potential mechanism underlying partial agonism at human α4β2 and α3β4 subtypes

To determine the structural components underlying differences in affinity, potency, and selectivity of varenicline for several human (h) nicotinic acetylcholine receptors (nAChRs), functional and structural experiments were performed. The Ca2+ influx results established that: (a) varenicline activates (μM range) nAChR subtypes with the following rank sequence: hα7>hα4β4>hα4β2>hα3β4>>>hα1β1γδ; (b) varenicline binds to nAChR subtypes with the following affinity order (nM range): hα4β2~hα4β4>hα3β4>hα7>>>Torpedo α1β1γδ. The molecular docking results indicating that more hydrogen bond interactions are apparent for α4-containing nAChRs in comparison to other nAChRs may explain the observed higher affinity; and that (c) varenicline is a full agonist at hα7 (101%) and hα4β4 (93%), and a partial agonist at hα4β2 (20%) and hα3β4 (45%), relative to (±)-epibatidine. The allosteric sites found at the extracellular domain (EXD) of hα3β4 and hα4β2 nAChRs could explain the partial agonistic activity of varenicline on these nAChR subtypes. Molecular dynamics simulations show that the interaction of varenicline to each allosteric site decreases the capping of Loop C at the hα4β2 nAChR, suggesting that these allosteric interactions limit the initial step in the gating process. In conclusion, we propose that in addition to hα4β2 nAChRs, hα4β4 nAChRs can be considered as potential targets for the clinical activity of varenicline, and that the allosteric interactions at the hα3β4- and hα4β2-EXDs are alternative mechanisms underlying partial agonism at these nAChRs.