Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A3 Receptor.

A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.

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.

Discovery and Structure-Activity Relationship Studies of Novel Adenosine A1 Receptor-Selective Agonists.

A series of benzyloxy and phenoxy derivatives of the adenosine receptor agonists N6-cyclopentyl adenosine (CPA) and N6-cyclopentyl 5'-N-ethylcarboxamidoadenosine (CP-NECA) were synthesized, and their potency and selectivity were assessed. We observed that the most potent were the compounds with a halogen in the meta position on the aromatic ring of the benzyloxy- or phenoxycyclopentyl substituent. In general, the NECA-based compounds displayed greater A1R selectivity than the adenosine-based compounds, with N6-2-(3-bromobenzyloxy)cyclopentyl-NECA and N6-2-(3-methoxyphenoxy)cyclopentyl-NECA showing ∼1500-fold improved A1R selectivity compared to NECA. In addition, we quantified the compounds' affinity and kinetics of binding at both human and rat A1R using a NanoBRET binding assay and found that the halogen substituent in the benzyloxy- or phenoxycyclopentyl moiety seems to confer high affinity for the A1R. Molecular modeling studies suggested a hydrophobic subpocket as contributing to the A1R selectivity displayed. We believe that the identified selective potent A1R agonists are valuable tool compounds for adenosine receptor research.

Dual A1/A3 Adenosine Receptor Antagonists: Binding Kinetics and Structure-Activity Relationship Studies Using Mutagenesis and Alchemical Binding Free Energy Calculations.

Drugs targeting adenosine receptors (AR) can provide treatment for diseases. We report the identification of 7-(phenylamino)-pyrazolo[3,4-c]pyridines L2-L10, A15, and A17 as low-micromolar to low-nanomolar A1R/A3R dual antagonists, with 3-phenyl-5-cyano-7-(trimethoxyphenylamino)-pyrazolo[3,4-c]pyridine (A17) displaying the highest affinity at both receptors with a long residence time of binding, as determined using a NanoBRET-based assay. Two binding orientations of A17 produce stable complexes inside the orthosteric binding area of A1R in molecular dynamics (MD) simulations, and we selected the most plausible orientation based on the agreement with alanine mutagenesis supported by affinity experiments. Interestingly, for drug design purposes, the mutation of L2506.51 to alanine increased the binding affinity of A17 at A1R. We explored the structure-activity relationships against A1R using alchemical binding free energy calculations with the thermodynamic integration coupled with the MD simulation (TI/MD) method, applied on the whole G-protein-coupled receptor-membrane system, which showed a good agreement (r = 0.73) between calculated and experimental relative binding free energies.

Selective activation of Gαob by an adenosine A1 receptor agonist elicits analgesia without cardiorespiratory depression.

The development of therapeutic agonists for G protein-coupled receptors (GPCRs) is hampered by the propensity of GPCRs to couple to multiple intracellular signalling pathways. This promiscuous coupling leads to numerous downstream cellular effects, some of which are therapeutically undesirable. This is especially the case for adenosine A1 receptors (A1Rs) whose clinical potential is undermined by the sedation and cardiorespiratory depression caused by conventional agonists. We have discovered that the A1R-selective agonist, benzyloxy-cyclopentyladenosine (BnOCPA), is a potent and powerful analgesic but does not cause sedation, bradycardia, hypotension or respiratory depression. This unprecedented discrimination between native A1Rs arises from BnOCPA's unique and exquisitely selective activation of Gob among the six Gαi/o subtypes, and in the absence of ß-arrestin recruitment. BnOCPA thus demonstrates a highly-specific Gα-selective activation of the native A1R, sheds new light on GPCR signalling, and reveals new possibilities for the development of novel therapeutics based on the far-reaching concept of selective Gα agonism.

Discovery of a High Affinity Adenosine A1/A3 Receptor Antagonist with a Novel 7-Amino-pyrazolo[3,4-d]pyridazine Scaffold.

Here we describe the design and synthesis of pyrazolo[3,4-d]pyridazines as adenosine receptor (AR) ligands. We demonstrate that the introduction of a 3-phenyl group, together with a 7-benzylamino and 1-methyl group at the pyrazolopyridazine scaffold, generated the antagonist compound 10b, which displayed 21 nM affinity and a residence time of ∼60 min, for the human A1R, 55 nM affinity and a residence time of ∼73 min, for the human A3R and 1.7 µΜ affinity for the human A2BR while not being toxic. Strikingly, the 2-methyl analog of 10b, 15b, had no significant affinity. Docking calculations and molecular dynamics simulations of the ligands inside the orthosteric binding area suggested that the 2-methyl group in 15b hinders the formation of hydrogen bonding interactions with N6.55 which are considered critical for the stabilization inside the orthosteric binding cavity. We, therefore, demonstrate that 10a is a novel scaffold for the development of high affinity AR ligands. From the mutagenesis experiments the biggest effect was observed for the Y2717.46A mutation which caused an ∼10-fold reduction in the binding affinity of 10b.

Suppression of Proliferation of Human Glioblastoma Cells by Combined Phosphodiesterase and Multidrug Resistance-Associated Protein 1 Inhibition.

The paucity of currently available therapies for glioblastoma multiforme requires novel approaches to the treatment of this brain tumour. Disrupting cyclic nucleotide-signalling through phosphodiesterase (PDE) inhibition may be a promising way of suppressing glioblastoma growth. Here, we examined the effects of 28 PDE inhibitors, covering all the major PDE classes, on the proliferation of the human U87MG, A172 and T98G glioblastoma cells. The PDE10A inhibitors PF-2545920, PQ10 and papaverine, the PDE3/4 inhibitor trequinsin and the putative PDE5 inhibitor MY-5445 potently decreased glioblastoma cell proliferation. The synergistic suppression of glioblastoma cell proliferation was achieved by combining PF-2545920 and MY-5445. Furthermore, a co-incubation with drugs that block the activity of the multidrug resistance-associated protein 1 (MRP1) augmented these effects. In particular, a combination comprising the MRP1 inhibitor reversan, PF-2545920 and MY-5445, all at low micromolar concentrations, afforded nearly complete inhibition of glioblastoma cell growth. Thus, the potent suppression of glioblastoma cell viability may be achieved by combining MRP1 inhibitors with PDE inhibitors at a lower toxicity than that of the standard chemotherapeutic agents, thereby providing a new combination therapy for this challenging malignancy.

Transcriptomics-Based Phenotypic Screening Supports Drug Discovery in Human Glioblastoma Cells.

We have used three established human glioblastoma (GBM) cell lines-U87MG, A172, and T98G-as cellular systems to examine the plasticity of the drug-induced GBM cell phenotype, focusing on two clinical drugs, the phosphodiesterase PDE10A inhibitor Mardepodect and the multi-kinase inhibitor Regorafenib, using genome-wide drug-induced gene expression (DIGEX) to examine the drug response. Both drugs upregulate genes encoding specific growth factors, transcription factors, cellular signaling molecules, and cell surface proteins, while downregulating a broad range of targetable cell cycle and apoptosis-associated genes. A few upregulated genes encode therapeutic targets already addressed by FDA approved drugs, but the majority encode targets for which there are no approved drugs. Amongst the latter, we identify many novel druggable targets that could qualify for chemistry-led drug discovery campaigns. We also observe several highly upregulated transmembrane proteins suitable for combined drug, immunotherapy, and RNA vaccine approaches. DIGEX is a powerful way of visualizing the complex drug response networks emerging during GBM drug treatment, defining a phenotypic landscape which offers many new diagnostic and therapeutic opportunities. Nevertheless, the extreme heterogeneity we observe within drug-treated cells using this technique suggests that effective pan-GBM drug treatment will remain a significant challenge for many years to come.

Emerging roles of adhesion G protein-coupled receptors.

Adhesion G protein-coupled receptors (aGPCRs) form a sub-group within the GPCR superfamily. Their distinctive structure contains an abnormally large N-terminal, extracellular region with a GPCR autoproteolysis-inducing (GAIN) domain. In most aGPCRs, the GAIN domain constitutively cleaves the receptor into two fragments. This process is often required for aGPCR signalling. Over the last two decades, much research has focussed on aGPCR-ligand interactions, in an attempt to deorphanize the family. Most ligands have been found to bind to regions N-terminal to the GAIN domain. These receptors may bind a variety of ligands, ranging across membrane-bound proteins and extracellular matrix components. Recent advancements have revealed a conserved method of aGPCR activation involving a tethered ligand within the GAIN domain. Evidence for this comes from increased activity in receptor mutants exposing the tethered ligand. As a result, G protein-coupling partners of aGPCRs have been more extensively characterised, making use of their tethered ligand to create constitutively active mutants. This has led to demonstrations of aGPCR function in, for example, neurodevelopment and tumour growth. However, questions remain around the ligands that may bind many aGPCRs, how this binding is translated into changes in the GAIN domain, and the exact mechanism of aGPCR activation following GAIN domain conformational changes. This review aims to examine the current knowledge around aGPCR activation, including ligand binding sites, the mechanism of GAIN domain-mediated receptor activation and how aGPCR transmembrane domains may relate to activation. Other aspects of aGPCR signalling will be touched upon, such as downstream effectors and physiological roles.

CGRP, adrenomedullin and adrenomedullin 2 display endogenous GPCR agonist bias in primary human cardiovascular cells.

Agonist bias occurs when different ligands produce distinct signalling outputs when acting at the same receptor. However, its physiological relevance is not always clear. Using primary human cells and gene editing techniques, we demonstrate endogenous agonist bias with physiological consequences for the calcitonin receptor-like receptor, CLR. By switching the receptor-activity modifying protein (RAMP) associated with CLR we can "re-route" the physiological pathways activated by endogenous agonists calcitonin gene-related peptide (CGRP), adrenomedullin (AM) and adrenomedullin 2 (AM2). AM2 promotes calcium-mediated nitric oxide signalling whereas CGRP and AM show pro-proliferative effects in cardiovascular cells, thus providing a rationale for the expression of the three peptides. CLR-based agonist bias occurs naturally in human cells and has a fundamental purpose for its existence. We anticipate this will be a starting point for more studies into RAMP function in native environments and their importance in endogenous GPCR signalling.

Multisite Model of Allosterism for the Adenosine A1 Receptor.

Despite being a target for about one-third of approved drugs, G protein-coupled receptors (GPCRs) still represent a tremendous reservoir for therapeutic strategies against diseases. For example, several cardiovascular and central nervous system conditions could benefit from clinical agents that activate the adenosine 1 receptor (A1R); however, the pursuit of A1R agonists for clinical use is usually impeded by both on- and off-target side effects. One of the possible strategies to overcome this issue is the development of positive allosteric modulators (PAMs) capable of selectively enhancing the effect of a specific receptor subtype and triggering functional selectivity (a phenomenon also referred to as bias). Intriguingly, besides enforcing the effect of agonists upon binding to an allosteric site, most of the A1R PAMs display intrinsic partial agonism and orthosteric competition with antagonists. To rationalize this behavior, we simulated the binding of the prototypical PAMs PD81723 and VCP171, the full-agonist NECA, the antagonist 13B, and the bitopic agonist VCP746. We propose that a single PAM can bind several A1R sites rather than a unique allosteric pocket, reconciling the structure-activity relationship and the mutagenesis results.

Structure-based identification of dual ligands at the A2AR and PDE10A with anti-proliferative effects in lung cancer cell-lines.

Enhanced/prolonged cAMP signalling has been suggested as a suppressor of cancer proliferation. Interestingly, two key modulators that elevate cAMP, the A2A receptor (A2AR) and phosphodiesterase 10A (PDE10A), are differentially co-expressed in various types of non-small lung cancer (NSCLC) cell-lines. Thus, finding dual-target compounds, which are simultaneously agonists at the A2AR whilst also inhibiting PDE10A, could be a novel anti-proliferative approach. Using ligand- and structure-based modelling combined with MD simulations (which identified Val84 displacement as a novel conformational descriptor of A2AR activation), a series of known PDE10A inhibitors were shown to dock to the orthosteric site of the A2AR. Subsequent in-vitro analysis confirmed that these compounds bind to the A2AR and exhibit dual-activity at both the A2AR and PDE10A. Furthermore, many of the compounds exhibited promising anti-proliferative effects upon NSCLC cell-lines, which directly correlated with the expression of both PDE10A and the A2AR. Thus, we propose a structure-based methodology, which has been validated in in-vitro binding and functional assays, and demonstrated a promising therapeutic value.

Deciphering the Agonist Binding Mechanism to the Adenosine A1 Receptor.

Despite being among the most characterized G protein-coupled receptors (GPCRs), adenosine receptors (ARs) have always been a difficult target in drug design. To date, no agonist other than the natural effector and the diagnostic regadenoson has been approved for human use. Recently, the structure of the adenosine A1 receptor (A1R) was determined in the active, Gi protein complexed state; this has important repercussions for structure-based drug design. Here, we employed supervised molecular dynamics simulations and mutagenesis experiments to extend the structural knowledge of the binding of selective agonists to A1R. Our results identify new residues involved in the association and dissociation pathway, they suggest the binding mode of N6-cyclopentyladenosine (CPA) related ligands, and they highlight the dramatic effect that chemical modifications can have on the overall binding mechanism, paving the way for the rational development of a structure-kinetics relationship of A1R agonists.

The Role of ICL1 and H8 in Class B1 GPCRs; Implications for Receptor Activation.

The first intracellular loop (ICL1) of G protein-coupled receptors (GPCRs) has received little attention, although there is evidence that, with the 8th helix (H8), it is involved in early conformational changes following receptor activation as well as contacting the G protein ß subunit. In class B1 GPCRs, the distal part of ICL1 contains a conserved R12.48KLRCxR2.46b motif that extends into the base of the second transmembrane helix; this is weakly conserved as a [R/H]12.48KL[R/H] motif in class A GPCRs. In the current study, the role of ICL1 and H8 in signaling through cAMP, iCa2+ and ERK1/2 has been examined in two class B1 GPCRs, using mutagenesis and molecular dynamics. Mutations throughout ICL1 can either enhance or disrupt cAMP production by CGRP at the CGRP receptor. Alanine mutagenesis identified subtle differences with regard elevation of iCa2+, with the distal end of the loop being particularly sensitive. ERK1/2 activation displayed little sensitivity to ICL1 mutation. A broadly similar pattern was observed with the glucagon receptor, although there were differences in significance of individual residues. Extending the study revealed that at the CRF1 receptor, an insertion in ICL1 switched signaling bias between iCa2+ and cAMP. Molecular dynamics suggested that changes in ICL1 altered the conformation of ICL2 and the H8/TM7 junction (ICL4). For H8, alanine mutagenesis showed the importance of E3908.49b for all three signal transduction pathways, for the CGRP receptor, but mutations of other residues largely just altered ERK1/2 activation. Thus, ICL1 may modulate GPCR bias via interactions with ICL2, ICL4 and the Gß subunit.

Pharmacological characterisation of novel adenosine A3 receptor antagonists.

The adenosine A3 receptor (A3R) belongs to a family of four adenosine receptor (AR) subtypes which all play distinct roles throughout the body. A3R antagonists have been described as potential treatments for numerous diseases including asthma. Given the similarity between (adenosine receptors) orthosteric binding sites, obtaining highly selective antagonists is a challenging but critical task. Here we screen 39 potential A3R, antagonists using agonist-induced inhibition of cAMP. Positive hits were assessed for AR subtype selectivity through cAMP accumulation assays. The antagonist affinity was determined using Schild analysis (pA2 values) and fluorescent ligand binding. Structure-activity relationship investigations revealed that loss of the 3-(dichlorophenyl)-isoxazolyl moiety or the aromatic nitrogen heterocycle with nitrogen at α-position to the carbon of carboximidamide group significantly attenuated K18 antagonistic potency. Mutagenic studies supported by molecular dynamic simulations combined with Molecular Mechanics-Poisson Boltzmann Surface Area calculations identified the residues important for binding in the A3R orthosteric site. We demonstrate that K18, which contains a 3-(dichlorophenyl)-isoxazole group connected through carbonyloxycarboximidamide fragment with a 1,3-thiazole ring, is a specific A3R (< 1 µM) competitive antagonist. Finally, we introduce a model that enables estimates of the equilibrium binding affinity for rapidly disassociating compounds from real-time fluorescent ligand-binding studies. These results demonstrate the pharmacological characterisation of a selective competitive A3R antagonist and the description of its orthosteric binding mode. Our findings may provide new insights for drug discovery.

Proposed model of the Dictyostelium cAMP receptors bound to cAMP.

3',5'-cyclic adenosine monophosphate (cAMP) is well known as a ubiquitous intracellular messenger regulating a diverse array of cellular processes. However, for a group of social amoebae or Dictyostelia undergoing starvation, intracellular cAMP is secreted in a pulsatile manner to their exterior. This then uniquely acts as a first messenger, triggering aggregation of the starving amoebae followed by their developmental progression towards multicellular fruiting bodies formation. Such developmental signalling for extracellularly-acting cAMP is well studied in the popular dictyostelid, Dictyostelium discoideum, and is mediated by a distinct family ('class E') of G protein-coupled receptors (GPCRs) collectively designated as the cAMP receptors (cARs). Whilst the biochemical aspects of these receptors are well characterised, little is known about their overall 3D architecture and structural basis for cAMP recognition and subtype-dependent changes in binding affinity. Using a ligand docking-guided homology modelling approach, we hereby present for the first time, plausible models of active forms of the cARs from D. discoideum. Our models highlight some structural features that may underlie the differential affinities of cAR isoforms for cAMP binding and also suggest few residues that may play important roles for the activation mechanism of this GPCR family.

Elevated intracellular cAMP concentration mediates growth suppression in glioma cells.

Supressed levels of intracellular cAMP have been associated with malignancy. Thus, elevating cAMP through activation of adenylyl cyclase (AC) or by inhibition of phosphodiesterase (PDE) may be therapeutically beneficial. Here, we demonstrate that elevated cAMP levels suppress growth in C6 cells (a model of glioma) through treatment with forskolin, an AC activator, or a range of small molecule PDE inhibitors with differing selectivity profiles. Forskolin suppressed cell growth in a PKA-dependent manner by inducing a G2/M phase cell cycle arrest. In contrast, trequinsin (a non-selective PDE2/3/7 inhibitor), not only inhibited cell growth via PKA, but also stimulated (independent of PKA) caspase-3/-7 and induced an aneuploidy phenotype. Interestingly, a cocktail of individual PDE 2,3,7 inhibitors suppressed cell growth in a manner analogous to forskolin but not trequinsin. Finally, we demonstrate that concomitant targeting of both AC and PDEs synergistically elevated intracellular cAMP levels thereby potentiating their antiproliferative actions.

Receptor component protein, an endogenous allosteric modulator of family B G protein coupled receptors.

Receptor component protein (RCP) is a 148 amino acid intracellular peripheral membrane protein, previously identified as promoting the coupling of CGRP to cAMP production at the CGRP receptor, a heterodimer of calcitonin receptor like-receptor (CLR), a family B G protein-coupled receptor (GPCR) and receptor activity modifying protein 1 (RAMP1). We extend these observations to show that it selectively enhances CGRP receptor coupling to Gs but not Gq or pERK activation. At other family B GPCRs, it enhances cAMP production at the calcitonin, corticotrophin releasing factor type 1a and glucagon-like peptide type 2 receptors with their cognate ligands but not at the adrenomedullin type 1 (AM1), gastric inhibitory peptide and glucagon-like peptide type 1 receptors, all expressed in transfected HEK293S cells. However, there is also cell-line variability as RCP did not enhance cAMP production at the endogenous calcitonin receptor in HEK293T cells and it has previously been reported that it is active on the AM1 receptor expressed on NIH3T3 cells. RCP appears to behave as a positive allosteric modulator at coupling a number of family B GPCRs to Gs, albeit in a manner that is regulated by cell-specific factors. It may exert its effects at the interface between the 2nd intracellular loop of the GPCR and Gs, although there is likely to be some overlap between this location and that occupied by the C-terminus of RAMPs if they bind to the GPCRs.