The Molecular Mechanism of P2Y1 Receptor Activation

Role of astrocytes in Alzheimer's disease pathogenesis and the impact of exercise-induced remodeling

Alzheimer's disease (AD) is a prevalent and debilitating brain disorder that worsens progressively with age, characterized by cognitive decline and memory impairment. The accumulation of amyloid-beta (Aβ) leading to amyloid plaques and hyperphosphorylation of Tau, resulting in intracellular neurofibrillary tangles (NFTs), are primary pathological features of AD. Despite significant research investment and effort, therapies targeting Aβ and NFTs have proven limited in efficacy for treating or slowing AD progression. Consequently, there is a growing interest in non-invasive therapeutic strategies for AD prevention. Exercise, a low-cost and non-invasive intervention, has demonstrated promising neuroprotective potential in AD prevention. Astrocytes, among the most abundant glial cells in the brain, play essential roles in various physiological processes and are implicated in AD initiation and progression. Exercise delays pathological progression and mitigates cognitive dysfunction in AD by modulating astrocyte morphological and phenotypic changes and fostering crosstalk with other glial cells. This review aims to consolidate the current understanding of how exercise influences astrocyte dynamics in AD, with a focus on elucidating the molecular and cellular mechanisms underlying astrocyte remodeling. The review begins with an overview of the neuropathological changes observed in AD, followed by an examination of astrocyte dysfunction as a feature of the disease. Lastly, the review explores the potential therapeutic implications of exercise-induced astrocyte remodeling in the context of AD.

Genetics of ischemic stroke functional outcome

Ischemic stroke, which accounts for 87% of cerebrovascular accidents, is responsible for massive global burden both in terms of economic cost and personal hardship. Many stroke survivors face long-term disability-a phenotype associated with an increasing number of genetic variants. While clinical variables such as stroke severity greatly impact recovery, genetic polymorphisms linked to functional outcome may offer physicians a unique opportunity to deliver personalized care based on their patient's genetic makeup, leading to improved outcomes. A comprehensive catalogue of the variants at play is required for such an approach. In this review, we compile and describe the polymorphisms associated with outcome scores such as modified Rankin Scale and Barthel Index. Our search identified 74 known genetic polymorphisms spread across 48 features associated with various poststroke disability metrics. The known variants span diverse biological systems and are related to inflammation, vascular homeostasis, growth factors, metabolism, the p53 regulatory pathway, and mitochondrial variation. Understanding how these variants influence functional outcome may be helpful in maximizing poststroke recovery.

G protein-coupled receptors (GPCRs): advances in structures, mechanisms, and drug discovery

G protein-coupled receptors (GPCRs), the largest family of human membrane proteins and an important class of drug targets, play a role in maintaining numerous physiological processes. Agonist or antagonist, orthosteric effects or allosteric effects, and biased signaling or balanced signaling, characterize the complexity of GPCR dynamic features. In this study, we first review the structural advancements, activation mechanisms, and functional diversity of GPCRs. We then focus on GPCR drug discovery by revealing the detailed drug-target interactions and the underlying mechanisms of orthosteric drugs approved by the US Food and Drug Administration in the past five years. Particularly, an up-to-date analysis is performed on available GPCR structures complexed with synthetic small-molecule allosteric modulators to elucidate key receptor-ligand interactions and allosteric mechanisms. Finally, we highlight how the widespread GPCR-druggable allosteric sites can guide structure- or mechanism-based drug design and propose prospects of designing bitopic ligands for the future therapeutic potential of targeting this receptor family.

Exploring bias in platelet P2Y1 signalling: Host defence versus haemostasis

Platelets are necessary for maintaining haemostasis. Separately, platelets are important for the propagation of inflammation during the host immune response against infection. The activation of platelets also causes inappropriate inflammation in various disease pathologies, often in the absence of changes to haemostasis. The separate functions of platelets during inflammation compared with haemostasis are therefore varied and this will be reflected in distinct pathways of activation. The activation of platelets by the nucleotide adenosine diphosphate (ADP) acting on P2Y1 and P2Y12 receptors is important for the development of platelet thrombi during haemostasis. However, P2Y1 stimulation of platelets is also important during the inflammatory response and paradoxically in scenarios where no changes to haemostasis and platelet aggregation occur. In these events, Rho-GTPase signalling, rather than the canonical phospholipase Cβ (PLCβ) signalling pathway, is necessary. We describe our current understanding of these differences, reflecting on recent advances in knowledge of P2Y1 structure, and the possibility of biased agonism occurring from activation via other endogenous nucleotides compared with ADP. Knowledge arising from these different pathways of P2Y1 stimulation of platelets during inflammation compared with haemostasis may help therapeutic control of platelet function during inflammation or infection, while preserving essential haemostasis. LINKED ARTICLES: This article is part of a themed issue on Platelet purinergic receptor and non-thrombotic disease. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.4/issuetoc.

The P2Y1 receptor in the colonic myenteric plexus of rats and its correlation with opioid-induced constipation

This study was designed to explore the expression changes of P2Y1 receptors in the distal colonic myenteric layer of rats. An opioid induced constipation(OIC) rat model was generated by intraperitoneal (i.p) injection of loperamide. At 7 days post-treatment, the model rats were assessed by calculating the fecal water content and the gastrointestinal transit ratio. The immunofluorescence (IF)-based histochemical study was used to observe the distribution of P2Y1 receptors in the distal colonic myenteric plexus. Western blotting (WB) was performed to evaluate the expression changes of P2Y1 proteins in the myenteric layer, and the electrophysiological approaches were carried out to determine the regulatory roles of P2Y1 receptors on distal colonic motor function. IF showed that P2Y1 receptors are co-expressed MOR in the enteric nerve cells of the distal colonic myenteric plexus. Moreover, the WB revealed that the protein levels of P2Y1 were significantly decreased in the distal colonic myenteric layer of OIC rats. In vitro tension experiments exhibited that the P2Y1 receptor antagonist MRS2500 enhanced the spontaneous contraction amplitude, adding EM2 and β-FNA did not have any effect on MRS2500. Therefore, P2Y1 receptor expression could be associated with the occurrence of OIC in this rat model and the regulation of colonic motility by MOR may be related to the release of purine neurotransmitters such as ATP in the colonic nervous system.

A comprehensive review of the ten main platelet receptors involved in platelet activity and cardiovascular disease

Cardiovascular disease (CVD) is a major cause of death worldwide. Although there are many variables that contribute to the development of this disease, it is predominantly the activity of platelets that provides the mechanisms by which this disease prevails. While there are numerous platelet receptors expressed on the surface of platelets, it is largely the consensus that there are 10 main platelet receptors that contribute to a majority of platelet function. Understanding these key platelet receptors is vitally important for patients suffering from myocardial infarction, CVD, and many other diseases that arise due to overactivation or mutations of these receptors. The goal of this manuscript is to review the main platelet receptors that contribute most to platelet activity.

The full activation mechanism of the adenosine A1 receptor revealed by GaMD and Su-GaMD simulations

The full activation process of G protein-coupled receptor (GPCR) plays an important role in cellular signal transduction. However, it remains challenging to simulate the whole process in which the GPCR is recognized and activated by a ligand and then couples to the G protein on a reasonable simulation timescale. Here, we developed a molecular dynamics (MD) approach named supervised (Su) Gaussian accelerated MD (GaMD) by incorporating a tabu-like supervision algorithm into a standard GaMD simulation. By using this Su-GaMD method, from the active and inactive structure of adenosine A₁ receptor (A₁R), we successfully revealed the full activation mechanism of A₁R, including adenosine (Ado)-A₁R recognition, preactivation of A₁R, and A₁R-G protein recognition, in hundreds of nanoseconds of simulations. The binding of Ado to the extracellular side of A₁R initiates conformational changes and the preactivation of A₁R. In turn, the binding of G_i2 to the intracellular side of A₁R causes a decrease in the volume of the extracellular orthosteric site and stabilizes the binding of Ado to A₁R. Su-GaMD could be a useful tool to reconstruct or even predict ligand-protein and protein-protein recognition pathways on a short timescale. The intermediate states revealed in this study could provide more detailed complementary structural characterizations to facilitate the drug design of A₁R in the future.

Extracellular loop 2 of G protein-coupled olfactory receptors is critical for odorant recognition

G protein–coupled olfactory receptors (ORs) enable us to detect innumerous odorants. They are also ectopically expressed in nonolfactory tissues and emerging as attractive drug targets. ORs can be promiscuous or highly specific, which is part of a larger mechanism for odor discrimination. Here, we demonstrate that the OR extracellular loop 2 (ECL2) plays critical roles in OR promiscuity and specificity. Using site-directed mutagenesis and molecular modeling, we constructed 3D OR models in which ECL2 forms a lid over the orthosteric pocket. We demonstrate using molecular dynamics simulations that ECL2 controls the shape and volume of the odorant-binding pocket, maintains the pocket hydrophobicity, and acts as a gatekeeper of odorant binding. Therefore, we propose the interplay between the specific orthosteric pocket and the variable, less specific ECL2 controls OR specificity and promiscuity. Furthermore, the 3D models created here enabled virtual screening of new OR agonists and antagonists, which exhibited a 70% hit rate in cell assays. Our approach can potentially be generalized to structure-based ligand screening for other G protein–coupled receptors that lack high-resolution 3D structures.

Interaction With the Lipid Membrane Influences Fentanyl Pharmacology

Overdose deaths from fentanyl have reached epidemic proportions in the USA and are increasing worldwide. Fentanyl is a potent opioid agonist that is less well reversed by naloxone than morphine. Due to fentanyl’s high lipophilicity and elongated structure we hypothesised that its unusual pharmacology may be explained by its interactions with the lipid membrane on route to binding to the μ-opioid receptor (MOPr). Through coarse-grained molecular dynamics simulations, electrophysiological recordings and cell signalling assays, we determined how fentanyl and morphine access the orthosteric pocket of MOPr. Morphine accesses MOPr via the aqueous pathway; first binding to an extracellular vestibule, then diffusing into the orthosteric pocket. In contrast, fentanyl may take a novel route; first partitioning into the membrane, before accessing the orthosteric site by diffusing through a ligand-induced gap between the transmembrane helices. In electrophysiological recordings fentanyl-induced currents returned after washout, suggesting fentanyl deposits in the lipid membrane. However, mutation of residues forming the potential MOPr transmembrane access site did not alter fentanyl’s pharmacological profile in vitro. A high local concentration of fentanyl in the lipid membrane, possibly in combination with a novel lipophilic binding route, may explain the high potency and lower susceptibility of fentanyl to reversal by naloxone.

Mapping the binding sites of UDP and prostaglandin E2 glyceryl ester in the nucleotide receptor P2Y6

Cyclooxygenase‐2 catalyzes the biosynthesis of prostaglandins from arachidonic acid and the biosynthesis of prostaglandin glycerol esters (PG‐Gs) from 2‐arachidonoylglycerol. PG‐Gs are mediators of several biological actions such as macrophage activation, hyperalgesia, synaptic plasticity, and intraocular pressure. Recently, the human UDP receptor P2Y₆ was identified as a target for the prostaglandin E2 glycerol ester (PGE₂‐G). Here, we show that UDP and PGE₂‐G are evolutionary conserved endogenous agonists at vertebrate P2Y₆ orthologs. Using sequence comparison of P2Y₆ orthologs, homology modeling, and ligand docking studies, we proposed several receptor positions participating in agonist binding. Site‐directed mutagenesis and functional analysis of these P2Y₆ mutants revealed that both UDP and PGE₂‐G share in parts one ligand‐binding site. Thus, the convergent signaling of these two chemically very different agonists has already been manifested in the evolutionary design of the ligand‐binding pocket.

Structural basis for molecular recognition of G protein-coupled estrogen receptor by selected bisphenols

Complicated ligand-dependent signaling pathways of bisphenol A (BPA) and its analogues involve not only intranuclear estrogen receptor but also membrane receptor G protein-coupled estrogen receptor (GPER). However, the structural basis for molecular recognition of GPER by the environmental chemicals remains unknown. To reveal the structural dependence of GPER recognition by bisphenols, a systematic molecular dynamics simulation study was performed for selected bisphenols with different electron hybrid orbitals and substituents on their C atoms connecting two phenol rings. BPA was used as a control, bisphenol C(BPC) as an example for a connecting C with sp² hybrid orbitals to provide more ligand rigidity, bisphenol E(BPE) and bisphenol F(BPF) for decreased steric hindrance and hydrophobicity around the connecting C, and bisphenol B(BPB) and bisphenol AF(BPAF) for increased hydrophobicity and steric hindrance. All the tested bisphenols can bind with GPER at its classic orthosteric site to obtain GPER-ligand complexes, while van der Waals interactions and direct inter-molecular electrostatic energies provide the driving forces for ligand binding. Bulky substituents and structural rigidity of the connecting C dramatically impair hydrogen bonding between GPER and the bisphenols, which results in decreased contribution of both favorable intermolecular hydrogen bonds and unfavorable polar solvation effect to complex stability of BPB and BPC since decreased number of key residues is expected. Increase in substituent lipophilicity enhances the van der Waals interactions and favorable non-polar solvation effect. The six bisphenols of high structural similarity shared two key recognition residues, Leu137^TM3 and Trp272^TM6, the latter of which was in the highly conserved CWxP motif of TM6 and has been reported as key residue for G protein-coupled receptor activation. Based on the obtained knowledge, GPER affinity and relevant toxicity of BPA alternatives can be easily predicted, and the calculated binding free energies are consistent with the available experimental observations.

Neuroprotective Roles of the Adenosine A3 Receptor Agonist AST-004 in Mouse Model of Traumatic Brain Injury

Traumatic brain injury (TBI) remains one of the greatest public health concerns with increasing morbidity and mortality rates worldwide. Our group reported that stimulation of astrocyte mitochondrial metabolism by P2Y1 receptor agonists significantly reduced cerebral edema and reactive gliosis in a TBI model. Subsequent data on the pharmacokinetics (PK) and rapid metabolism of these compounds suggested that neuroprotection was likely mediated by a metabolite, AST-004, which binding data indicated was an adenosine A3 receptor (A3R) agonist. The neuroprotective efficacy of AST-004 was tested in a control closed cortical injury (CCCI) model of TBI in mice. Twenty-four (24) hours post-injury, mice subjected to CCCI and treated with AST-004 (0.22 mg/kg, injected 30 min post-trauma) exhibited significantly less secondary brain injury. These effects were quantified with less cell death (PSVue794 fluorescence) and loss of blood brain barrier breakdown (Evans blue extravasation assay), compared to vehicle-treated TBI mice. TBI-treated mice also exhibited significantly reduced neuroinflammatory markers, glial-fibrillary acidic protein (GFAP, astrogliosis) and ionized Ca2+-binding adaptor molecule 1 (Iba1, microgliosis), both at the mRNA (qRT-PCR) and protein (Western blot and immunofluorescence) levels, respectively. Four (4) weeks post-injury, both male and female TBI mice presented a significant reduction in freezing behavior during contextual fear conditioning (after foot shock). AST-004 treatment prevented this TBI-induced impairment in male mice, but did not significantly affect impairment in female mice. Impairment of spatial memory, assessed 24 and 48 h after the initial fear conditioning, was also reduced in AST-004-treated TBI-male mice. Female TBI mice did not exhibit memory impairment 24 and 48 h after contextual fear conditioning and similarly, AST-004-treated female TBI mice were comparable to sham mice. Finally, AST-004 treatments were found to increase in vivo ATP production in astrocytes (GFAP-targeted luciferase activity), consistent with the proposed mechanism of action. These data reveal AST-004 as a novel A3R agonist that increases astrocyte energy production and enhances their neuroprotective efficacy after brain injury.

Different responses of the blockade of the P2Y1 receptor with BPTU in human and porcine intestinal tissues and in cell cultures

BACKGROUND

Gastrointestinal smooth muscle relaxation is accomplished by activation of P2Y₁ receptors, therefore this receptor plays an important role in regulation of gut motility. Recently, BPTU was developed as a negative allosteric modulator of the P2Y₁ receptor. Accordingly, the aim of this study was to assess the effect of BPTU on purinergic neurotransmission in pig and human gastrointestinal tissues.

METHODS

Ca²⁺ imaging in tSA201 cells that express the human P2Y₁ receptor, organ bath and microelectrodes in tissues were used to evaluate the effects of BPTU on purinergic responses.

KEY RESULTS

BPTU concentration dependently (0.1 and 1 µmol L^-1 ) inhibited the rise in intracellular Ca²⁺ evoked by ADP in tSA201 cells. In the pig small intestine, 30 µmol L^-1 BPTU reduced the fast inhibitory junction potential by 80%. Smooth muscle relaxations induced by electrical field stimulation were reduced both in pig ileum (EC₅₀  = 6 µmol L^-1 ) and colon (EC₅₀  = 35 µmol L^-1 ), but high concentrations of BPTU (up to 100 µmol L^-1 ) had no effect on human colonic muscle. MRS2500 (1 µmol L^-1 ) abolished all responses. Finally, 10 µmol L^-1 ADPβS inhibited spontaneous motility and this was partially reversed by 30 µmol L^-1 BPTU in pig, but not human colonic tissue and abolished by MRS2500 (1 µmol L^-1 ).

CONCLUSIONS & INFERENCES

BPTU blocks purinergic responses elicited via P2Y₁ receptors in cell cultures and in pig gastrointestinal tissue. However, the concentrations needed are higher in pig tissue compared to cell cultures and BPTU was ineffective in human colonic tissue.

In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors

Molecular modeling has contributed to drug discovery for purinergic GPCRs, including adenosine receptors (ARs) and P2Y receptors (P2YRs). Experimental structures and homology modeling have proven to be useful in understanding and predicting structure activity relationships (SAR) of agonists and antagonists. This review provides an excursus on molecular dynamics (MD) simulations applied to ARs and P2YRs. The binding modes of newly synthesized A1AR- and A3AR-selective nucleoside derivatives, potentially of use against depression and inflammation, respectively, have been predicted to recapitulate their SAR and the species dependence of A3AR affinity. P2Y12R and P2Y1R crystallographic structures, respectively, have provided a detailed understanding of the recognition of anti-inflammatory P2Y14R antagonists and a large group of allosteric and orthosteric antagonists of P2Y1R, an antithrombotic and neuroprotective target. MD of A2AAR (an anticancer and neuroprotective target), A3AR, and P2Y1R has identified microswitches that are putatively involved in receptor activation. The approach pathways of different ligands toward A2AAR and P2Y1R binding sites have also been explored. A1AR, A2AAR, and A3AR were utilizes to study allosteric phenomena, but locating the binding site of structurally diverse allosteric modulators, such as an A3AR enhancer LUF6000, is challenging. Ligand residence time, a predictor of in vivo efficacy, and the structural role of water were investigated through A2AAR MD simulations. Thus, new MD and other modeling algorithms have contributed to purinergic GPCR drug discovery.

Computational Investigations on the Binding Mode of Ligands for the Cannabinoid-Activated G Protein-Coupled Receptor GPR18

GPR18 is an orphan G protein-coupled receptor (GPCR) expressed in cells of the immune system. It is activated by the cannabinoid receptor (CB) agonist ∆⁹-tetrahydrocannabinol (THC). Several further lipids have been proposed to act as GPR18 agonists, but these results still require unambiguous confirmation. In the present study, we constructed a homology model of the human GPR18 based on an ensemble of three GPCR crystal structures to investigate the binding modes of the agonist THC and the recently reported antagonists which feature an imidazothiazinone core to which a (substituted) phenyl ring is connected via a lipophilic linker. Docking and molecular dynamics simulation studies were performed. As a result, a hydrophobic binding pocket is predicted to accommodate the imidazothiazinone core, while the terminal phenyl ring projects towards an aromatic pocket. Hydrophobic interaction of Cys251 with substituents on the phenyl ring could explain the high potency of the most potent derivatives. Molecular dynamics simulation studies suggest that the binding of imidazothiazinone antagonists stabilizes transmembrane regions TM1, TM6 and TM7 of the receptor through a salt bridge between Asp118 and Lys133. The agonist THC is presumed to bind differently to GPR18 than to the distantly related CB receptors. This study provides insights into the binding mode of GPR18 agonists and antagonists which will facilitate future drug design for this promising potential drug target.

Update of P2Y receptor pharmacology: IUPHAR Review 27

Eight G protein-coupled P2Y receptor subtypes respond to extracellular adenine and uracil mononucleotides and dinucleotides. P2Y receptors belong to the δ group of rhodopsin-like GPCRs and contain two structurally distinct subfamilies: P2Y₁ , P2Y₂ , P2Y₄ , P2Y₆ , and P2Y₁₁ (principally G_q protein-coupled P2Y₁ -like) and P2Y_12-14 (principally G_i protein-coupled P2Y₁₂ -like) receptors. Brain P2Y receptors occur in neurons, glial cells, and vasculature. Endothelial P2Y₁ , P2Y₂ , P2Y₄ , and P2Y₆ receptors induce vasodilation, while smooth muscle P2Y₂ , P2Y₄ , and P2Y₆ receptor activation leads to vasoconstriction. Pancreatic P2Y₁ and P2Y₆ receptors stimulate while P2Y₁₃ receptors inhibits insulin secretion. Antagonists of P2Y₁₂ receptors, and potentially P2Y₁ receptors, are anti-thrombotic agents, and a P2Y₂ /P2Y₄ receptor agonist treats dry eye syndrome in Asia. P2Y receptor agonists are generally pro-inflammatory, and antagonists may eventually treat inflammatory conditions. This article reviews recent developments in P2Y receptor pharmacology (using synthetic agonists and antagonists), structure and biophysical properties (using X-ray crystallography, mutagenesis and modelling), physiological and pathophysiological roles, and present and potentially future therapeutic targeting.

Ligand binding and activation of UTP-activated G protein-coupled P2Y2 and P2Y4 receptors elucidated by mutagenesis, pharmacological and computational studies

The nucleotide receptors P2Y₂ and P2Y₄ are the most closely related G protein-coupled receptors (GPCRs) of the P2Y receptor (P2YR) family. Both subtypes couple to G_q proteins and are activated by the pyrimidine nucleotide UTP, but only P2Y₂R is also activated by the purine nucleotide ATP. Agonists and antagonists of both receptor subtypes have potential as drugs e.g. for neurodegenerative and inflammatory diseases. So far, potent and selective, "drug-like" ligands for both receptors are scarce, but would be required for target validation and as lead structures for drug development. Structural information on the receptors is lacking since no X-ray structures or cryo-electron microscopy images are available. Thus, we performed receptor homology modeling and docking studies combined with mutagenesis experiments on both receptors to address the question how ligand binding selectivity for these closely related P2YR subtypes can be achieved. The orthosteric binding site of P2Y₂R appeared to be more spacious than that of P2Y₄R. Mutation of Y197 to alanine in P2Y₄R resulted in a gain of ATP sensitivity. Anthraquinone-derived antagonists are likely to bind to the orthosteric or an allosteric site depending on their substitution pattern and the nature of the orthosteric binding site of the respective P2YR subtype. These insights into the architecture of P2Y₂- and P2Y₄Rs and their interactions with structurally diverse agonists and antagonist provide a solid basis for the future design of potent and selective ligands.

The effect of dimerization on the activation and conformational dynamics of adenosine A1 receptor

The adenosine A1 receptor (A1R) is one of four adenosine receptors in humans, which are involved in the function of the cardiovascular, respiratory and central nervous systems. Experimental results indicate that A1R can form a homodimer and that the protomer-protomer interaction in the A1R dimer is related to certain pharmacological characteristics of A1R activation. In this work, we performed docking, metadynamics simulation, conventional molecular dynamics simulations, Gaussian-accelerated molecular dynamics simulations, potential of mean force calculations, dynamic cross-correlation motions analysis and community network analysis to study the binding mode of 5'-N-ethylcarboxamidoadenosine (NECA) to A1R and the effect of dimerization on the activation of A1R. Our results show that NECA binds to A1R in a similar mode to adenosine in the A1R crystal structure and NECA in the A2AR crystal structure. The A1R homodimer can be activated by one or two agonists with NECA occupying its orthosteric pockets in one (which we call the NECA-A1R system) or both protomers (which we call the dNECA-A1R system). In the NECA-A1R system, activation is predicated in the protomer without NECA bound. In the dNECA-A1R system, only one protomer achieves the active state. These findings suggest an asymmetrical activation mechanism of the homodimer and a negative cooperativity between the two protomers. We envision that our results may further facilitate the drug development of A1R.

Molecular Modeling of an Orphan GPR18 Receptor

Background: :

One of the best known to date GPCR class A (Rhodopsin) includes more than 100 orphan receptors for which the endogenous ligand is not known or is unclear. One of them is N-arachidonyl glycine receptor, named GPR18, a receptor that has been reported to be activated by Δ9-THC, endogenous cannabinoid receptors agonist anandamide and other cannabinoid receptor ligands suggesting it could be considered as third cannabinoid receptor. GPR18 activity, as well as its distribution might suggest usage of GPR18 ligands in treatment of endometriosis, cancer, and neurodegenerative disorders. Yet, so far only few GPR18 antagonists have been described, thus only ligand-based design approaches appear to be most useful to identify new ligands for this orphan receptor.

Methods: :

Main goal of this study, GPR18 inactive form homology model was built on the basis of the evolutionary closest homologous template: Human P2Y1 Receptor crystal structure.

Results: :

Obtained model was further evaluated and showed active/nonactive ligands differentiating properties with acceptable confidence. Moreover, it allowed for preliminary assessment of proteinligand interactions for a set of previously described ligands.

Conclusion::

Thus collected data might serve as a starting point for a discovery of novel, active GPR18 blocking ligands.

Toward Understanding the Impact of Dimerization Interfaces in Angiotensin II Type 1 Receptor

Angiotensin II type 1 receptor (AT1R) is a prototypical class A G protein-coupled receptor (GPCR) that has an important role in cardiovascular pathologies and blood pressure regulation as well as in the central nervous system. GPCRs may exist and function as monomers; however, they can assemble to form higher order structures, and as a result of oligomerization, their function and signaling profiles can be altered. In the case of AT1R, the classical Gα_q/11 pathway is initiated with endogenous agonist angiotensin II binding. A variety of cardiovascular pathologies such as heart failure, diabetic nephropathy, atherosclerosis, and hypertension are associated with this pathway. Recent findings reveal that AT1R can form homodimers and activate the noncanonical (β-arrestin-mediated) pathway. Nevertheless, the exact dimerization interface and atomic details of AT1R homodimerization have not been still elucidated. Here, six different symmetrical dimer interfaces of AT1R are considered, and homodimers were constructed using other published GPCR crystal dimer interfaces as template structures. These AT1R homodimers were then inserted into the model membrane bilayers and subjected to all-atom molecular dynamics simulations. Our simulation results along with the principal component analysis and water pathway analysis suggest four different interfaces as the most plausible: symmetrical transmembrane (TM)1,2,8; TM5; TM4; and TM4,5 AT1R dimer interfaces that consist of one inactive and one active protomer. Moreover, we identified ILE238^6.33 as a hub residue in the stabilization of the inactive state of AT1R.

Computational modeling of the olfactory receptor Olfr73 suggests a molecular basis for low potency of olfactory receptor-activating compounds

The mammalian olfactory system uses hundreds of specialized G-protein-coupled olfactory receptors (ORs) to discriminate a nearly unlimited number of odorants. Cognate agonists of most ORs have not yet been identified and potential non-olfactory processes mediated by ORs are unknown. Here, we used molecular modeling, fingerprint interaction analysis and molecular dynamics simulations to show that the binding pocket of the prototypical olfactory receptor Olfr73 is smaller, but more flexible, than binding pockets of typical non-olfactory G-protein-coupled receptors. We extended our modeling to virtual screening of a library of 1.6 million compounds against Olfr73. Our screen predicted 25 Olfr73 agonists beyond traditional odorants, of which 17 compounds, some with therapeutic potential, were validated in cell-based assays. Our modeling suggests a molecular basis for reduced interaction contacts between an odorant and its OR and thus the typical low potency of OR-activating compounds. These results provide a proof-of-principle for identifying novel therapeutic OR agonists.

Pharmacology of P2Y receptors

P2Y receptors are G-protein-coupled receptors (GPCRs) for extracellular nucleotides. There are eight mammalian P2Y receptor subtypes divided into two subgroups (P2Y₁, P2Y₂, P2Y₄, P2Y₆, and P2Y₁₁) and (P2Y₁₂, P2Y₁₃, and P2Y₁₄). The P2Y receptors are expressed in various cell types and play important roles in physiology and pathophysiology including inflammatory responses and neuropathic pain. The antagonism of P2Y₁₂ receptors is used in pharmacotherapy for the prevention and therapy of cardiovascular events. The nucleoside analogue ticagrelor and active metabolites of the thienopyridine compounds ticlopidine, clopidogrel and prasugrel inhibit platelet P2Y₁₂ receptors and reduce thereby platelet aggregation. The P2Y₂ receptor agonist diquafosol is used for the treatment of the dry eye syndrome. The P2Y receptor subtypes differ in their amino acid sequences, their pharmacological profiles and their signaling transduction pathways. Recently, selective receptor ligands have been developed for all subtypes. The published crystal structures of the human P2Y₁ and P2Y₁₂ receptors as well as receptor models will facilitate the development of novel drugs for pharmacotherapy.

Allosteric Modulation of Class A GPCRs: Targets, Agents, and Emerging Concepts

G-protein-coupled receptors (GPCRs) have been tractable drug targets for decades with over one-third of currently marketed drugs targeting GPCRs. Of these, the class A GPCR superfamily is highly represented, and continued drug discovery for this family of receptors may provide novel therapeutics for a vast range of diseases. GPCR allosteric modulation is an innovative targeting approach that broadens the available small molecule toolbox and is proving to be a viable drug discovery strategy, as evidenced by recent FDA approvals and clinical trials. Numerous class A GPCR allosteric modulators have been discovered recently, and emerging trends such as the availability of GPCR crystal structures, diverse functional assays, and structure-based computational approaches are improving optimization and development. This Perspective provides an update on allosterically targeted class A GPCRs and their disease indications and the medicinal chemistry approaches toward novel allosteric modulators and highlights emerging trends and opportunities in the field.

Studies towards the development of a PET radiotracer for imaging of the P2Y1 receptors in the brain: synthesis, 18F-labeling and preliminary biological evaluation

Purine nucleotides such as ATP and ADP are important extracellular signaling molecules in almost all tissues activating various subtypes of purinoreceptors. In the brain, the P2Y₁ receptor (P2Y₁R) subtype mediates trophic functions like differentiation and proliferation, and modulates fast synaptic transmission, both suggested to be affected in diseases of the central nervous system. Research on P2Y₁R is limited because suitable brain-penetrating P2Y₁R-selective tracers are not yet available. Here, we describe the first efforts to develop an ¹⁸F-labeled PET tracer based on the structure of the highly affine and selective, non-nucleotidic P2Y₁R allosteric modulator 1-(2-[2-(tert-butyl)phenoxy]pyridin-3-yl)-3-[4-(trifluoromethoxy)phenyl]urea (7). A small series of fluorinated compounds was developed by systematic modification of the p-(trifluoromethoxy)phenyl, the urea and the 2-pyridyl subunits of the lead compound 7. Additionally, the p-(trifluoromethoxy)phenyl subunit was substituted by carborane, a boron-rich cluster with potential applicability in boron neutron capture therapy (BNCT). By functional assays, the new fluorinated derivative 1-{2-[2-(tert-butyl)phenoxy]pyridin-3-yl}-3-[4-(2-fluoroethyl)phenyl]urea (18) was identified with a high P2Y₁R antagonistic potency (IC₅₀ = 10 nM). Compound [¹⁸F]18 was radiosynthesized by using tetra-n-butyl ammonium [¹⁸F]fluoride with high radiochemical purity, radiochemical yield and molar activities. Investigation of brain homogenates using hydrophilic interaction chromatography (HILIC) revealed [¹⁸F]fluoride as major radiometabolite. Although [¹⁸F]18 showed fast in vivo metabolization, the high potency and unique allosteric binding mode makes this class of compounds interesting for further optimizations and investigation of the theranostic potential as PET tracer and BNCT agent.

New Binding Sites, New Opportunities for GPCR Drug Discovery

Many central biological events rely on protein-ligand interactions. The identification and characterization of protein-binding sites for ligands are crucial for the understanding of functions of both endogenous ligands and synthetic drug molecules. G protein-coupled receptors (GPCRs) typically detect extracellular signal molecules on the cell surface and transfer these chemical signals across the membrane, inducing downstream cellular responses via G proteins or β-arrestin. GPCRs mediate many central physiological processes, making them important targets for modern drug discovery. Here, we focus on the most recent breakthroughs in finding new binding sites and binding modes of GPCRs and their potentials for the development of new medicines.

Exploring a new ligand binding site of G protein-coupled receptors

Identifying a target ligand binding site is an important step for structure-based rational drug design as shown here for G protein-coupled receptors (GPCRs), which are among the most popular drug targets. We applied long-time scale molecular dynamics simulations, coupled with mutagenesis studies, to two prototypical GPCRs, the M3 and M4 muscarinic acetylcholine receptors. Our results indicate that unlike synthetic antagonists, which bind to the classic orthosteric site, the endogenous agonist acetylcholine is able to diffuse into a much deeper binding pocket. We also discovered that the most recently resolved crystal structure of the LTB4 receptor comprised a bound inverse agonist, which extended its benzamidine moiety to the same binding pocket discovered in this work. Analysis on all resolved GPCR crystal structures indicated that this new pocket could exist in most receptors. Our findings provide new opportunities for GPCR drug discovery.

The Molecular Mechanism Underlying Ligand Binding to the Membrane-Embedded Site of a G-Protein-Coupled Receptor

The crystal structure of P2Y₁ receptor (P2Y₁R), a class A GPCR, revealed a special extra-helical site for its antagonist, BPTU, which locates in-between the membrane and the protein. However, due to the limitation of crystallization experiments, the membrane was mimicked by use of detergents, and the information related to the binding of BPTU to the receptor in the membrane environment is rather limited. In the present work, we conducted a total of ∼7.5 μs all-atom simulations in explicit solvent using conventional molecular dynamics and multiple enhanced sampling methods, with models of BPTU and a POPC bilayer, both in the absence and presence of P2Y₁R. Our simulations revealed that BPTU prefers partitioning into the interface of polar/lipophilic region of the lipid bilayer before associating with the receptor. Then, it interacts with the second extracellular loop of the receptor and reaches the binding site through the lipid-receptor interface. In addition, by use of funnel-metadynamics simulations which efficiently enhance the sampling of bound and unbound states, we provide a statistically accurate description of the underlying binding free energy landscape. The calculated absolute ligand-receptor binding affinity is in excellent agreement with the experimental data (Δ G_{b0_theo} = -11.5 kcal mol^-1, Δ G_{b0_exp}= -11.7 kcal mol^-1). Our study broadens the view of the current experimental/theoretical models and our understanding of the protein-ligand recognition mechanism in the lipid environment. The strategy used in this work is potentially applicable to investigate ligands association/dissociation with other membrane-embedded sites, allowing identification of compounds targeting membrane receptors of pharmacological interest.

Diverse signalling of the platelet P2Y1 receptor leads to a dichotomy in platelet function

Platelet P2Y₁ receptor signalling via RhoGTPases is necessary for platelet-dependent leukocyte recruitment, where no platelet aggregation is observed. We investigated signalling cascades involved in distinct P2Y₁-dependent platelet activities in vitro, using specific inhibitors for phospholipase C (PLC) (U73122, to inhibit the canonical pathway), and RhoGTPases: Rac1 (NSC23766) and RhoA (ROCK inhibitor GSK429286). Human platelet rich plasma (for platelet aggregation) or isolated washed platelets (for chemotaxis assays) was treated with U73122, GSK429286 or NSC23766 prior to stimulation with adenosine diphosphate (ADP) or the P2Y₁ specific agonist MRS2365. Aggregation, chemotaxis (towards f-MLP), or platelet-induced human neutrophil chemotaxis (PINC) towards macrophage derived chemokine (MDC) was assessed. Molecular docking of ADP and MRS2365 to P2Y₁ was analysed using AutoDock Smina followed by GOLD molecular docking in the Accelrys Discovery Studio software. Inhibition of PLC, but not Rac1 or RhoA, suppressed platelet aggregation induced by ADP and MRS2365. In contrast, platelet chemotaxis and PINC, were significantly attenuated by inhibition of platelet Rac1 or RhoA, but not PLC. MRS2365, compared to ADP had a less pronounced effect on P2Y₁-induced aggregation, but a similar efficacy to stimulate platelet chemotaxis and PINC, which might be explained by differences in molecular interaction of ADP compared to MRS2365 with the P2Y₁ receptor. Platelet P2Y₁ receptor activation during inflammation signals through alternate pathways involving Rho GTPases in contrast to canonical P2Y₁ receptor induced PLC signalling. This might be explained by selective molecular interactions of ligands within the orthosteric site of the P2Y₁ receptor.

Conformation and dynamics of the C-terminal region in human phosphoglycerate mutase 1

Phosphoglycerate mutase 1 (PGAM1), an important enzyme in glycolysis, is overexpressed in a number of human cancers, thus has been proposed as a promising metabolic target for cancer treatments. The C-terminal portion of the available crystal structures of PGAM1 and its homologous proteins is partially disordered, as evidenced by weak electron density. In this study, we identified the conformational behavior of the C-terminal region of PGAM1 as well as its role during the catalytic cycle. Using the PONDR-FIT server, we demonstrated that the C-terminal region was intrinsically disordered. We applied the Monte Carlo (MC) method to explore the conformational space of the C-terminus and conducted a series of explicit-solvent molecular dynamics (MD) simulations, and revealed that the C-terminal region is inherently dynamic; large-scale conformational changes in the C-terminal segment led to the structural transition of PGAM1 from the closed state to the open state. Furthermore, the C-terminal segment influenced 2,3-bisphosphoglycerate (2,3-BPG) binding. The proposed swing model illustrated a critical role of the C-terminus in the catalytic cycle through the conformational changes. In conclusion, the C-terminal region induces large movements of PGAM1 from the closed state to the open state and influences cofactor binding during the catalytic cycle. This report describes the dynamic features of the C-terminal region in detail and should aid in design of novel and efficient inhibitors of PGAM1. A swing mechanism of the C-terminal region is proposed, to facilitate further studies of the catalytic mechanism and the physiological functions of its homologues.

Demystifying P2Y1 Receptor Ligand Recognition through Docking and Molecular Dynamics Analyses

We performed a molecular modeling analysis of 100 nucleotide-like bisphosphates and 46 non-nucleotide arylurea derivatives previously reported as P2Y₁R binders using the recently solved hP2Y₁R structures. We initially docked the compounds at the X-ray structures and identified the binding modes of representative compounds highlighting key patterns in the structure-activity relationship (SAR). We subsequently subjected receptor complexes with selected key agonists (2MeSADP and MRS2268) and antagonists (MRS2500 and BPTU) to membrane molecular dynamics (MD) simulations (at least 200 ns run in triplicate, simulation time 0.6-1.6 μs per ligand system) while considering alternative protonation states of nucleotides. Comparing the temporal evolution of the ligand-protein interaction patterns with available site-directed mutagenesis (SDM) data and P2Y₁R apo state simulation provided further SAR insights and suggested reasonable explanations for loss/gain of binding affinity as well as the most relevant charged species for nucleotide ligands. The MD analysis also predicted local conformational changes required for the receptor inactive state to accommodate nucleotide agonists.

Identification of a Different Agonist-Binding Site and Activation Mechanism of the Human P2Y1 Receptor

The human P2Y1 receptor (P2Y1R) is a purinergic G-protein-coupled receptor (GPCR) that functions as a receptor for adenosine 5'-diphosphate (ADP). An antagonist of P2Y1R might potentially have antithrombotic effects, whereas agonists might serve as antidiabetic agents. On the basis of the antagonist-bound MRS2500-P2Y1R crystal structure, we constructed computational models of apo-P2Y1R and the agonist-receptor complex 2MeSADP-P2Y1R. We then performed conventional molecular dynamics (cMD) and accelerated molecular dynamics (aMD) simulations to study the conformational dynamics after binding with agonist/antagonist as well as the P2Y1R activation mechanism. We identified a new agonist-binding site of P2Y1R that is consistent with previous mutagenesis data. This new site is deeper than those of the agonist ADP in the recently simulated ADP-P2Y1R structure and the antagonist MRS2500 in the MRS2500-P2Y1R crystal structure. During P2Y1R activation, the cytoplasmic end of helix VI shifts outward 9.1 Å, the Ser1463.47-Tyr2375.58 hydrogen bond breaks, a Tyr2375.58-Val2626.37 hydrogen bond forms, and the conformation of the χ1 rotamer of Phe2696.44 changes from parallel to perpendicular to helix VI. The apo-P2Y1R system and the MRS2500-P2Y1R system remain inactive. The newly identified agonist binding site and activation mechanism revealed in this study may aid in the design of P2Y1R antagonists/agonists as antithrombotic/antidiabetic agents, respectively.

Molecular Recognition of Agonists and Antagonists by the Nucleotide-Activated G Protein-Coupled P2Y2 Receptor

A homology model of the nucleotide-activated P2Y₂R was created based on the X-ray structures of the P2Y₁ receptor. Docking studies were performed, and receptor mutants were created to probe the identified binding interactions. Mutation of residues predicted to interact with the ribose (Arg110) and the phosphates of the nucleotide agonists (Arg265, Arg292) or that contribute indirectly to binding (Tyr288) abolished activity. The Y114F, R194A, and F261A mutations led to inactivity of diadenosine tetraphosphate and to a reduced response of UTP. Significant reduction in agonist potency was observed for all other receptor mutants (Phe111, His184, Ser193, Phe261, Tyr268, Tyr269) predicted to be involved in agonist recognition. An ionic lock between Asp185 and Arg292 that is probably involved in receptor activation interacts with the phosphate groups. The antagonist AR-C118925 and anthraquinones likely bind to the orthosteric site. The updated homology models will be useful for virtual screening and drug design.

Activation-Induced Conformational Changes of Dopamine D3 Receptor Promote the Formation of the Internal Water Channel

The atomic-level dopamine activation mechanism for transmitting extracellular ligand binding events through transmembrane helices to the cytoplasmic G protein remains unclear. In the present study, the complete dopamine D3 receptor (D3R), with a homology-modeled N-terminus, was constructed to dock different ligands to simulate conformational alterations in the receptor’s active and inactive forms during microsecond-timescale molecular dynamic simulations. In agonist-bound systems, the D3R N-terminus formed a “lid-like” structure and lay flat on the binding site opening, whereas in antagonist and inverse agonist-bound systems, the N-terminus exposed the binding cavity. Receptor activation was characterized using the different molecular switch residue distances, and G protein-binding site volumes. A continuous water pathway was observed only in the dopamine-G_αi-bound system. In the inactive D3Rs, water entry was hindered by the hydrophobic layers. Finally, a complete activation mechanism of D3R was proposed. Upon agonist binding, the “lid-like” conformation of the N-terminus induces a series of molecular switches to increase the volume of the D3R cytoplasmic binding part for G protein association. Meanwhile, water enters the transmembrane region inducing molecular switches to assist in opening the hydrophobic layers to form a continuous water channel, which is crucial for maintaining a fully active conformation for signal transduction.

Designing Safer Analgesics via μ-Opioid Receptor Pathways

Pain is both a major clinical and economic problem, affecting more people than diabetes, heart disease, and cancer combined. While a variety of prescribed or over-the-counter (OTC) medications are available for pain management, opioid medications, especially those acting on the μ-opioid receptor (μOR) and related pathways, have proven to be the most effective, despite some serious side effects including respiration depression, pruritus, dependence, and constipation. It is therefore imperative that both academia and industry develop novel μOR analgesics which retain their opioid analgesic properties but with fewer or no adverse effects. In this review we outline recent progress towards the discovery of safer opioid analgesics.

Conformational Transitions of the Pituitary Adenylate Cyclase-Activating Polypeptide Receptor, a Human Class B GPCR

The G protein-coupled pituitary adenylate cyclase-activating polypeptide receptor (PAC1R) is a potential therapeutic target for endocrine, metabolic and stress-related disorders. However, many questions regarding the protein structure and dynamics of PAC1R remain largely unanswered. Using microsecond-long simulations, we examined the open and closed PAC1R conformations interconnected within an ensemble of transitional states. The open-to-closed transition can be initiated by “unzipping” the extracellular domain and the transmembrane domain, mediated by a unique segment within the β3-β4 loop. Transitions between different conformational states range between microseconds to milliseconds, which clearly implicate allosteric effects propagating from the extracellular face of the receptor to the intracellular G protein-binding site. Such allosteric dynamics provides structural and mechanistic insights for the activation and modulation of PAC1R and related class B receptors.

Ligand binding to a G protein-coupled receptor captured in a mass spectrometer

G protein (heterotrimeric guanine nucleotide-binding protein)-coupled receptors belong to the largest family of membrane-embedded cell surface proteins and are involved in a diverse array of physiological processes. Despite progress in the mass spectrometry of membrane protein complexes, G protein-coupled receptors have remained intractable because of their low yield and instability after extraction from cell membranes. We established conditions in the mass spectrometer that preserve noncovalent ligand binding to the human purinergic receptor P2Y1. Results established differing affinities for nucleotides and the drug MRS2500 and link antagonist binding with the absence of receptor phosphorylation. Overall, therefore, our results are consistent with drug binding, preventing the conformational changes that facilitate downstream signaling. More generally, we highlight opportunities for mass spectrometry to probe effects of ligand binding on G protein-coupled receptors.

Development of Potent and Selective Antagonists for the UTP-Activated P2Y4 Receptor

P2Y₄ is a G_q protein-coupled receptor activated by uridine-5'-triphosphate (UTP), which is widely expressed in the body, e.g., in intestine, heart, and brain. No selective P2Y₄ receptor antagonist has been described so far. Therefore, we developed and optimized P2Y₄ receptor antagonists based on an anthraquinone scaffold. Potency was assessed by a fluorescence-based assay measuring inhibition of UTP-induced intracellular calcium release in 1321N1 astrocytoma cells stably transfected with the human P2Y₄ receptor. The most potent compound of the present series, sodium 1-amino-4-[4-(2,4-dimethylphenylthio)phenylamino]-9,10-dioxo-9,10-dihydroanthracene-2-sulfonate (PSB-16133, 61) exhibited an IC₅₀ value of 233 nM, selectivity versus other P2Y receptor subtypes, and is thought to act as an allosteric antagonist. A receptor homology model was built and docking studies were performed to analyze ligand-receptor interactions. Compound 64 (PSB-1699, sodium 1-amino-4-[4-(3-pyridin-3-ylmethylthio)phenylamino]-9,10-dioxo-9,10-dihydroanthracene-2-sulfonate) represents the most selective P2Y₄ receptor antagonist known to date. Compounds 61 and 64 are therefore anticipated to become useful tools for studying this scarcely investigated receptor.

The molecular mechanism of hPPARα activation

MD simulations were performed to explore the molecular mechanism of hPPARα activation. 11 key residues favouring binding ligands and the movements of helices and loops playing important roles in inducing the active conformation change of hPPARα were discovered.

The Principles of Ligand Specificity on beta-2-adrenergic receptor

G protein-coupled receptors are recognized as one of the largest families of membrane proteins. Despite sharing a characteristic seven-transmembrane topology, G protein-coupled receptors regulate a wide range of cellular signaling pathways in response to various physical and chemical stimuli, and prevail as an important target for drug discovery. Notably, the recent progress in crystallographic methods led to a breakthrough in elucidating the structures of membrane proteins. The structures of β₂-adrenergic receptor bound with a variety of ligands provide atomic details of the binding modes of agonists, antagonists and inverse agonists. In this study, we selected four representative molecules from each functional class of ligands and investigated their impacts on β₂-adrenergic receptor through a total of 12 × 100 ns molecular dynamics simulations. From the obtained trajectories, we generated molecular fingerprints exemplifying propensities of protein-ligand interactions. For each functional class of compounds, we characterized and compared the fluctuation of the protein backbone, the volumes in the intracellular pockets, the water densities in the receptors, the domain interaction networks as well as the movements of transmembrane helices. We discovered that each class of ligands exhibits a distinct mode of interactions with mainly TM5 and TM6, altering the shape and eventually the state of the receptor. Our findings provide insightful prospective into GPCR targeted structure-based drug discoveries.

The Molecular Mechanism of P2Y1 Receptor Activation

Human purinergic G protein-coupled receptor P2Y1 (P2Y1 R) is activated by adenosine 5'-diphosphate (ADP) to induce platelet activation and thereby serves as an important antithrombotic drug target. Crystal structures of P2Y1 R revealed that one ligand (MRS2500) binds to the extracellular vestibule of this GPCR, whereas another (BPTU) occupies the surface between transmembrane (TM) helices TM2 and TM3. We introduced a total of 20 μs all-atom long-timescale molecular dynamic (MD) simulations to inquire why two molecules in completely different locations both serve as antagonists while ADP activates the receptor. Our results indicate that BPTU acts as an antagonist by stabilizing extracellular helix bundles leading to an increase of the lipid order, whereas MRS2500 blocks signaling by occupying the ligand binding site. Both antagonists stabilize an ionic lock within the receptor. However, binding of ADP breaks this ionic lock, forming a continuous water channel that leads to P2Y1 R activation.