Development and Characterization of Novel Selective, Non-Basic Dopamine D2 Receptor Antagonists for the Treatment of Schizophrenia.

The dopamine D2 receptor, which belongs to the family of G protein-coupled receptors (GPCR), is an important and well-validated drug target in the field of medicinal chemistry due to its wide distribution, particularly in the central nervous system, and involvement in the pathomechanism of many disorders thereof. Schizophrenia is one of the most frequent diseases associated with disorders in dopaminergic neurotransmission, and in which the D2 receptor is the main target for the drugs used. In this work, we aimed at discovering new selective D2 receptor antagonists with potential antipsychotic activity. Twenty-three compounds were synthesized, based on the scaffold represented by the D2AAK2 compound, which was discovered by our group. This compound is an interesting example of a D2 receptor ligand because of its non-classical binding to this target. Radioligand binding assays and SAR analysis indicated structural modifications of D2AAK2 that are possible to maintain its activity. These findings were further rationalized using molecular modeling. Three active derivatives were identified as D2 receptor antagonists in cAMP signaling assays, and the selected most active compound 17 was subjected to X-ray studies to investigate its stable conformation in the solid state. Finally, effects of 17 assessed in animal models confirmed its antipsychotic activity in vivo.

Allosteric Modulators of Dopamine D2 Receptors for Fine-Tuning of Dopaminergic Neurotransmission in CNS Diseases: Overview, Pharmacology, Structural Aspects and Synthesis.

Allosteric modulation of G protein-coupled receptors (GPCRs) is nowadays a hot topic in medicinal chemistry. Allosteric modulators, i.e., compounds which bind in a receptor site topologically distinct from orthosteric sites, exhibit a number of advantages. They are more selective, safer and display a ceiling effect which prevents overdosing. Allosteric modulators of dopamine D2 receptor are potential drugs against a number of psychiatric and neurological diseases, such as schizophrenia and Parkinson's disease. In this review, an insightful summary of current research on D2 receptor modulators is presented, ranging from their pharmacology and structural aspects of ligand-receptor interactions to their synthesis.

The Role of Lipids in Allosteric Modulation of Dopamine D2 Receptor-In Silico Study.

The dopamine D2 receptor, belonging to the class A G protein-coupled receptors (GPCRs), is an important drug target for several diseases, including schizophrenia and Parkinson's disease. The D2 receptor can be activated by the natural neurotransmitter dopamine or by synthetic ligands, which in both cases leads to the receptor coupling with a G protein. In addition to receptor modulation by orthosteric or allosteric ligands, it has been shown that lipids may affect the behaviour of membrane proteins. We constructed a model of a D2 receptor with a long intracellular loop (ICL3) coupled with Giα1 or Giα2 proteins, embedded in a complex asymmetric membrane, and simulated it in complex with positive, negative or neutral allosteric ligands. In this study, we focused on the influence of ligand binding and G protein coupling on the membrane-receptor interactions. We show that there is a noticeable interplay between the cell membrane, G proteins, D2 receptor and its modulators.

Some Dietary Phenolic Compounds Can Activate Thyroid Peroxidase and Inhibit Lipoxygenase-Preliminary Study in the Model Systems.

The presented research concerns the triple activity of trans-cinnamic (tCA), ferulic (FA) and syringic acids (SA). They act as thyroid peroxidase (TPO) activators, lipoxygenase (LOX) inhibitors and show antiradical activity. All compounds showed a dose-dependent TPO activatory effect, thus the AC50 value (the concentration resulting in 50% activation) was determined. The tested compounds can be ranked as follows: tCA > FA > SA with AC50 = 0.10, 0.39, 0.69 mM, respectively. Strong synergism was found between FA and SA. The activatory effects of all tested compounds may result from interaction with the TPO allosteric site. It was proposed that conformational change resulting from activator binding to TPO allosteric pocket results from the flexibility of a nearby loop formed by residues Val352-Tyr363. All compounds act as uncompetitive LOX inhibitors. The most effective were tCA and SA, whereas the weakest was FA (IC50 = 0.009 mM and IC50 0.027 mM, respectively). In all cases, an interaction between the inhibitors carboxylic groups and side-chain atoms of Arg102 and Arg139 in an allosteric pocket of LOX was suggested. FA/tCA and FA/SA acted synergistically, whereas tCA/SA demonstrated antagonism. The highest antiradical activity was found in the case of SA (IC50 = 0.22 mM). FA/tCA and tCA/SA acted synergistically, whereas antagonism was found for the SA/FA mixture.

Novel Positive Allosteric Modulators of µ Opioid Receptor-Insight from In Silico and In Vivo Studies.

Opioids are the drugs of choice in severe pain management. Unfortunately, their use involves serious, potentially lethal side effects. Therefore, efforts in opioid drug design turn toward safer and more effective mechanisms, including allosteric modulation. In this study, molecular dynamics simulations in silico and 'writhing' tests in vivo were used to characterize potential allosteric mechanism of two previously reported compounds. The results suggest that investigated compounds bind to µ opioid receptor in an allosteric site, augmenting action of morphine at subeffective doses, and exerting antinociceptive effect alone at higher doses. Detailed analysis of in silico calculations suggests that first of the compounds behaves more like allosteric agonist, while the second compound acts mainly as a positive allosteric modulator.

Preferential Coupling of Dopamine D2S and D2L Receptor Isoforms with Gi1 and Gi2 Proteins-In Silico Study.

The dopamine D2 receptor belongs to rhodopsin-like G protein-coupled receptors (GPCRs) and it is an important molecular target for the treatment of many disorders, including schizophrenia and Parkinson's disease. Here, computational methods were used to construct the full models of the dopamine D2 receptor short (D2S) and long (D2L) isoforms (differing with 29 amino acids insertion in the third intracellular loop, ICL3) and to study their coupling with Gi1 and Gi2 proteins. It was found that the D2L isoform preferentially couples with the Gi2 protein and D2S isoform with the Gi1 protein, which is in accordance with experimental data. Our findings give mechanistic insight into the interplay between isoforms of dopamine D2 receptors and Gi proteins subtypes, which is important to understand signaling by these receptors and their mediation by pharmaceuticals, in particular psychotic and antipsychotic agents.

The Universal 3D QSAR Model for Dopamine D2 Receptor Antagonists.

In order to search for novel antipsychotics acting through the D2 receptor, it is necessary to know the structure-activity relationships for dopamine D2 receptor antagonists. In this context, we constructed the universal three-dimensional quantitative structure-activity relationship (3D- QSAR) model for competitive dopamine D2 receptor antagonists. We took 176 compounds from chemically different groups characterized by the half maximal inhibitory concentration (IC50)from the CHEMBL database and docked them to the X-ray structure of the human D2 receptor in the inactive state. Selected docking poses were applied for Comparative Molecular Field Analysis (CoMFA) alignment. The obtained CoMFA model is characterized by a cross-validated coefficient Q2 of 0.76 with an optimal component of 5, R2 of 0.92, and an F value of 338.9. The steric and electrostatic field contributions are 67.4% and 32.6%, respectively. The statistics obtained prove that the CoMFA model is significant. Next, the IC50 of the 16 compounds from the test set was predicted with R2 of 0.95. Finally, a progressive scrambling test was carried out for additional validation. The CoMFA fields were mapped onto the dopamine D2 receptor binding site, which enabled a discussion of the structure-activity relationship based on ligand-receptor interactions. In particular, it was found that one of the desired steric interactions covers the area of a putative common allosteric pocket suggested for some other G protein-coupled receptors (GPCRs), which would suggest that some of the known dopamine receptor antagonists are bitopic in their essence. The CoMFA model can be applied to predict the potential activity of novel dopamine D2 receptor antagonists.

Signaling within Allosteric Machines: Signal Transmission Pathways Inside G Protein-Coupled Receptors.

In recent years, our understanding of function of G protein-coupled receptors (GPCRs) has changed from a picture of simple signal relays, transmitting only a particular signal to a particular G protein heterotrimer, to versatile machines, capable of various responses to different stimuli and being modulated by various factors. Some recent reports provide not only the data on ligands/modulators and resultant signals induced by them, but also deeper insights into exact pathways of signal migration and mechanisms of signal transmission through receptor structure. Combination of these computational and experimental data sheds more light on underlying mechanisms of signal transmission and signaling bias in GPCRs. In this review we focus on available clues on allosteric pathways responsible for complex signal processing within GPCRs structures, with particular emphasis on linking compatible in silico- and in vitro-derived data on the most probable allosteric connections.

Recent Advances and Applications of Molecular Docking to G Protein-Coupled Receptors.

The growing number of studies on G protein-coupled receptors (GPCRs) family are a source of noticeable improvement in our understanding of the functioning of these proteins. GPCRs are responsible for a vast part of signaling in vertebrates and, as such, invariably remain in the spotlight of medicinal chemistry. A deeper insight into the underlying mechanisms of interesting phenomena observed in GPCRs, such as biased signaling or allosteric modulation, can be gained with experimental and computational studies. The latter play an important role in this process, since they allow for observations on scales inaccessible for most other methods. One of the key steps in such studies is proper computational reconstruction of actual ligand-receptor or protein-protein interactions, a process called molecular docking. A number of improvements and innovative applications of this method were documented recently. In this review, we focus particularly on innovations in docking to GPCRs.

Interplay between Two Allosteric Sites and Their Influence on Agonist Binding in Human µ Opioid Receptor.

Allostery is a widespread mechanism that allows for precise protein tuning. Its underlying mechanisms are elusive, particularly when there are multiple allosteric sites at the protein. This concerns also G-protein-coupled receptors (GPCRs), which are targets for a vast part of currently used drugs. To address this issue, we performed molecular dynamics simulations of a GPCR-human µ opioid receptor (MOR) in a native-like environment, with full agonist (R)-methadone, Na(+) ions, and a positive modulator BMS986122 in various configurations. We found that MOR's seventh transmembrane helix (TM VII) is central for allosteric signal transmission, and modulators affect its bending and rotation. The PAM stabilizes favorable agonist interactions, while Na(+) tends to disrupt agonist binding. We identified two residues involved in allosteric signal transmission: Trp 7.35 at the top and Tyr 7.53 at the bottom of TM VII.

Activation and Allosteric Modulation of Human µ Opioid Receptor in Molecular Dynamics.

Allosteric protein modulation has gained increasing attention in drug design. Its application as a mechanism of action could bring forth safer and more effective medicines. Targeting opioid receptors with allosteric modulators can result in better treatment of pain, depression, and respiratory and immune disorders. In this work we use recent reports on negative modulators of µ opioid receptor as a starting point for identification of allosteric sites and mechanisms of opioid receptor modulation using homology modeling and docking and molecular dynamics studies. An allosteric binding site description is presented. Results suggest a shared binding region for lipophilic allosteric ligands, reveal possible differences in the modulation mechanism between cannabinoids and salvinorin A, and show ambiguous properties of the latter. Also, they emphasize the importance of native-like environment in molecular dynamics simulations and uncover relationships between modulator and orthosteric ligand binding and receptor behavior. Relationships between ligands, transmission switch, and hydrophobic lock are analyzed.

Discovery and in†vitro Evaluation of Novel Serotonin 5-HT2A Receptor Ligands Identified Through Virtual Screening.

The 5-HT2A receptor is a molecular target of high pharmacological importance. Ligands of this protein, particularly atypical antipsychotics, are useful in the treatment of numerous mental disorders, including schizophrenia and major depressive disorder. Structure-based virtual screening using a 5-HT2A receptor complex was performed to identify novel ligands for the 5-HT2A receptor, serving as potential antidepressants. From the Enamine screening library, containing over 4 million compounds, 48 molecules were selected for subsequent experimental validation. These compounds were tested against the 5-HT2A receptor in radioligand binding assays. From the tested batch, six molecules were identified as ligands of the main molecular target and were forwarded to a more detailed in vitro profiling. This included radioligand binding assays at 5-HT1A, 5-HT7, and D2 receptors and functional studies at 5-HT2A receptors. These compounds were confirmed to show a binding affinity for at least one of the targets tested in vitro. The success rate for the inactive template-based screening reached 17 %, while it was 9 % for the active template-based screening. Similarity and fragment analysis indicated the structural novelty of the identified compounds. Pharmacokinetics for these molecules was determined using in silico approaches.

Illuminating the "Twilight Zone": Advances in Difficult Protein Modeling.

Homology modeling was long considered a method of choice in tertiary protein structure prediction. However, it used to provide models of acceptable quality only when templates with appreciable sequence identity with a target could be found. The threshold value was long assumed to be around 20-30%. Below this level, obtained sequence identity was getting dangerously close to values that can be obtained by chance, after aligning any random, unrelated sequences. In these cases, other approaches, including ab initio folding simulations or fragment assembly, were usually employed. The most recent editions of the CASP and CAMEO community-wide modeling methods assessment have brought some surprising outcomes, proving that much more clues can be inferred from protein sequence analyses than previously thought. In this chapter, we focus on recent advances in the field of difficult protein modeling, pushing the threshold deep into the "twilight zone", with particular attention devoted to improvements in applications of machine learning and model evaluation.

Allosteric modulation of dopamine D2L receptor in complex with Gi1 and Gi2 proteins: the effect of subtle structural and stereochemical ligand modifications.

BACKGROUND: Allosteric modulation of G protein-coupled receptors (GPCRs) is nowadays one of the hot topics in drug discovery. In particular, allosteric modulators of D2 receptor have been proposed as potential modern therapeutics to treat schizophrenia and Parkinson's disease. METHODS: To address some subtle structural and stereochemical aspects of allosteric modulation of D2 receptor, we performed extensive in silico studies of both enantiomers of two compounds (compound 1 and compound 2), and one of them (compound 2) was synthesized as a racemate in-house and studied in vitro. RESULTS: Our molecular dynamics simulations confirmed literature reports that the R enantiomer of compound 1 is a positive allosteric modulator of the D2L receptor, while its S enantiomer is a negative allosteric modulator. Moreover, based on the principal component analysis (PCA), we hypothesized that both enantiomers of compound 2 behave as silent allosteric modulators, in line with our in vitro studies. PCA calculations suggest that the most pronounced modulator-induced receptor rearrangements occur at the transmembrane helix 7 (TM7). In particular, TM7 bending at the conserved P7.50 and G7.42 was observed. The latter resides next to the Y7.43, which is a significant part of the orthosteric binding site. Moreover, the W7.40 conformation seems to be affected by the presence of the positive allosteric modulator. CONCLUSIONS: Our work reveals that allosteric modulation of the D2L receptor can be affected by subtle ligand modifications. A change in configuration of a chiral carbon and/or minor structural modulator modifications are solely responsible for the functional outcome of the allosteric modulator.

Funnel metadynamics and behavioral studies reveal complex effect of D2AAK1 ligand on anxiety-like processes.

Anxiety is a troublesome symptom for many patients, especially those suffering from schizophrenia. Its regulation involves serotonin receptors, targeted e.g. by antipsychotics or psychedelics such as LSD. 5-HT2A receptors are known for an extremely long LSD residence time, enabling minute doses to exert a long-lasting effect. In this work, we explore the changes in anxiety-like processes induced by the previously reported antipsychotic, D2AAK1. In vivo studies revealed that the effect of D2AAK1 on the anxiety is mediated through serotonin 5-HT1A and 5-HT2A receptors, and that it is time-dependent (anxiogenic after 30 min, anxiolytic after 60 min) and dose-dependent. The funnel metadynamics simulations suggest complicated ligand-5HT2AR interactions, involving an allosteric site located under the third extracellular loop, which is a possible explanation of the time-dependency. The binding of D2AAK1 at the allosteric site results in a broader opening of the extracellular receptor entry, possibly altering the binding kinetics of orthosteric ligands.

Tuning Down the Pain - An Overview of Allosteric Modulation of Opioid Receptors: Mechanisms of Modulation, Allosteric Sites, Modulator Syntheses.

Opioid signaling plays a central role in pain perception. As such, it remains the main target in the development of antinociceptive agents, despite serious side effects involved. In recent years, hopes for improved opioid painkillers are rising, together with our understanding of allosterism and biased signaling mechanisms. In this review, we focus on recently discovered allosteric modulators of opioid receptors, insights into phenomena underlying their action, as well as on how they extend our understanding of mechanisms of previously known compounds. A brief overlook of their synthesis is also presented.

Role of N-terminus in function and dynamics of sirtuin 7: an in silico study.

The sirtuin family comprises seven NAD+-dependent histone deacetylases named SIRT1 to SIRT7. The least investigated SIRT7 is currently considered as a promising therapeutic target for cardiovascular diseases, diabetes and different types of cancer. So far, its structure was not experimentally resolved, except of a fragment of its N-terminus. The aim of this study was to create in silico model of SIRT7 containing its core together with N-terminus, which is known to affect the enzyme's catalytic activity and to find pockets that could be targeted by structure-based virtual screening. Homology model of SIRT7 was prepared using X-ray structures of other sirtuins and a resolved fragment of the N-terminus of SIRT7 as templates. All atom-unbiased molecular dynamics simulations were performed. It was found that N-terminus of SIRT7 remains in spatial proximity of the catalytic core for considerable fraction of time, and therefore, it may affect its catalytic activity by helping the enzyme to hold the substrate peptide. It may also participate in holding and release of the cofactor. Preferred orientations of NAD+ and acetyl-lysine inside SIRT7 were found, with all components forming a stable complex. Molecular dynamics provided an ensemble of conformations that will be targeted with virtual screening. Reliable in silico structure of SIRT7 will be a useful tool in searching for its inhibitors, which can be potential drugs in cancer treatment.Communicated by Ramaswamy H. Sarma.

Molecular mechanisms of allosteric probe dependence in µ opioid receptor.

Allostery is one of the most important features of proteins. It greatly contributes to the complexity of life, since it enables possibility of precise tuning of protein function, as well as performing more than one function per protein. Probe dependence is one of the unique features of allostery. It allows a protein to respond differently to the same allosteric modulator when different drugs or transmitters are bound. Unfortunately, allosteric mechanisms are difficult to investigate experimentally. Instead, they can be reproduced artificially in simulations. We simulated in silico a native-like cell membrane fragment with an active-state human µ opioid receptor (MOR) in order to investigate diverse effects of a receptor's positive allosteric modulator on various agonists. Particular emphasis on native-likeness of the environment was put. We managed to reproduce the experimentally observed effects, which allowed us to take deeper insight into their underlying mechanisms. We found an allosteric pathway in the receptor, leading from the ligand binding site to the intracellular, effector site. We observed that the modulator affected the pathway, inducing different resultant responses for full and partial agonists.

Opportunities and Challenges in the Discovery of Allosteric Modulators of GPCRs.

From the pharmacological point of view, allosteric modulators may present numerous advantages over orthosteric ligands. Growing availability of novel tools and experimental data provides a tempting opportunity to apply computational methods to improve known modulators and design novel ones. However, recent progress in understanding of complexity of allostery increases awareness of problems involved in design of modulators with desired properties. Deeper insight into phenomena such as probe dependence, altering signaling bias with minor changes in ligand structure, as well as influence of subtle endogenous allosteric factors turns out to be fundamental. These effects make the design of a modulator with precise pharmacological outcome a very challenging task, and need to be taken into consideration throughout the design process. In this chapter, we focus on nuances of targeting GPCR allosteric sites in computational drug design efforts, in particular with application of docking, virtual screening, and molecular dynamics.

Protein-Protein Docking in Drug Design and Discovery.

Protein-protein interactions (PPIs) are responsible for a number of key physiological processes in the living cells and underlie the pathomechanism of many diseases. Nowadays, along with the concept of so-called "hot spots" in protein-protein interactions, which are well-defined interface regions responsible for most of the binding energy, these interfaces can be targeted with modulators. In order to apply structure-based design techniques to design PPIs modulators, a three-dimensional structure of protein complex has to be available. In this context in silico approaches, in particular protein-protein docking, are a valuable complement to experimental methods for elucidating 3D structure of protein complexes. Protein-protein docking is easy to use and does not require significant computer resources and time (in contrast to molecular dynamics) and it results in 3D structure of a protein complex (in contrast to sequence-based methods of predicting binding interfaces). However, protein-protein docking cannot address all the aspects of protein dynamics, in particular the global conformational changes during protein complex formation. In spite of this fact, protein-protein docking is widely used to model complexes of water-soluble proteins and less commonly to predict structures of transmembrane protein assemblies, including dimers and oligomers of G protein-coupled receptors (GPCRs). In this chapter we review the principles of protein-protein docking, available algorithms and software and discuss the recent examples, benefits, and drawbacks of protein-protein docking application to water-soluble proteins, membrane anchoring and transmembrane proteins, including GPCRs.