Targeting hPKM2 in cancer: A bio isosteric approach for ligand design.

The term cancer refers to a plethora of diseases characterized by the development of abnormal cells that divide uncontrollably and can infiltrate further proximal or distal body tissues. Each type of cancer can be defined by aggressiveness, localization, metabolism, and response to available treatments. Among the most common hallmarks of cancer is a more acidic intracellular microenvironment. Offset pH values are due to an excess of lactate and an increased hypoxia-inducible factor (HIF) expression, which leads to a hypoxic state and a metabolic shift towards glycolysis to produce adenosine-5'-triphosphate (ATP) necessary for cellular metabolism. Warburg's hypothesis underpins this concept, making glycolysis and its central enzyme pyruvate kinase (hPKM2), an ideal target for drug development. Using molecular docking and extensive molecular dynamics (MD) simulations we investigated the binding mode of phosphoenolpyruvate (PEP) inside the hPKM2 active site, and then evaluated a set of known bio-isosteric inhibitors to understand the differences caused by their substitutions on their binding mode. Ultimately, we propose a new molecular entity to hamper hPKM2, unbalance cellular energy, and possibly trigger autophagic mechanisms.

Understanding VPAC receptor family peptide binding and selectivity.

The vasoactive intestinal peptide (VIP) and pituitary adenylate cyclase-activating polypeptide (PACAP) receptors are key regulators of neurological processes. Despite recent structural data, a comprehensive understanding of peptide binding and selectivity among different subfamily receptors is lacking. Here, we determine structures of active, Gs-coupled, VIP-VPAC1R, PACAP27-VPAC1R, and PACAP27-PAC1R complexes. Cryo-EM structural analyses and molecular dynamics simulations (MDSs) reveal fewer stable interactions between VPAC1R and VIP than for PACAP27, more extensive dynamics of VIP interaction with extracellular loop 3, and receptor-dependent differences in interactions of conserved N-terminal peptide residues with the receptor core. MD of VIP modelled into PAC1R predicts more transient VIP-PAC1R interactions in the receptor core, compared to VIP-VPAC1R, which may underlie the selectivity of VIP for VPAC1R over PAC1R. Collectively, our work improves molecular understanding of peptide engagement with the PAC1R and VPAC1R that may benefit the development of novel selective agonists.

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.

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.

Structural and functional diversity among agonist-bound states of the GLP-1 receptor.

Recent advances in G-protein-coupled receptor (GPCR) structural elucidation have strengthened previous hypotheses that multidimensional signal propagation mediated by these receptors depends, in part, on their conformational mobility; however, the relationship between receptor function and static structures is inherently uncertain. Here, we examine the contribution of peptide agonist conformational plasticity to activation of the glucagon-like peptide 1 receptor (GLP-1R), an important clinical target. We use variants of the peptides GLP-1 and exendin-4 (Ex4) to explore the interplay between helical propensity near the agonist N terminus and the ability to bind to and activate the receptor. Cryo-EM analysis of a complex involving an Ex4 analog, the GLP-1R and Gs heterotrimer revealed two receptor conformers with distinct modes of peptide-receptor engagement. Our functional and structural data, along with molecular dynamics (MD) simulations, suggest that receptor conformational dynamics associated with flexibility of the peptide N-terminal activation domain may be a key determinant of agonist efficacy.

Exploring Ligand Binding to Calcitonin Gene-Related Peptide Receptors.

Class B1 G protein-coupled receptors (GPCRs) are important targets for many diseases, including cancer, diabetes, and heart disease. All the approved drugs for this receptor family are peptides that mimic the endogenous activating hormones. An understanding of how agonists bind and activate class B1 GPCRs is fundamental for the development of therapeutic small molecules. We combined supervised molecular dynamics (SuMD) and classic molecular dynamics (cMD) simulations to study the binding of the calcitonin gene-related peptide (CGRP) to the CGRP receptor (CGRPR). We also evaluated the association and dissociation of the antagonist telcagepant from the extracellular domain (ECD) of CGRPR and the water network perturbation upon binding. This study, which represents the first example of dynamic docking of a class B1 GPCR peptide, delivers insights on several aspects of ligand binding to CGRPR, expanding understanding of the role of the ECD and the receptor-activity modifying protein 1 (RAMP1) on agonist selectivity.

Insights into agonist-elicited activation of the human glucose-dependent insulinotropic polypeptide receptor.

Glucose-dependent insulinotropic polypeptide (GIP) and its receptor (GIPR) are part of the incretin system that regulates glucose homeostasis. A series of GIPR residues putatively important for ligand binding and receptor activation were mutated and pharmacologically evaluated using GIPR selective agonists in cAMP accumulation, ERK1/2 phosphorylation (pERK1/2) and ß-arrestin 2 recruitment assays. The impact of mutation on ligand efficacy was determined by operational modelling of experimental data for each mutant, with results mapped onto the full-length, active-state GIPR structure. Two interaction networks, comprising transmembrane helix (TM) 7, TM1 and TM2, and extracellular loop (ECL) 2, TM5 and ECL3 were revealed, respectively. Both networks were critical for Gαs-mediated cAMP accumulation and the recruitment of ß-arrestin 2, however, cAMP response was more sensitive to alanine substitution, with most mutated residues displaying reduced signaling. Unlike the other two assays, activation of ERK1/2 was largely independent of the network involving ECL2, TM5 and ECL3, indicating that pERK1/2 is at least partially distinct from Gαs or ß-arrestin pathways and this network is also crucial for potential biased agonism at GIPR. Collectively, our work advances understanding of the structure-function relationship of GIPR and provides a framework for the design and/or interpretation of GIP analogues with unique signaling profiles.

Partial agonism improves the anti-hyperglycaemic efficacy of an oxyntomodulin-derived GLP-1R/GCGR co-agonist.

OBJECTIVE: Glucagon-like peptide-1 and glucagon receptor (GLP-1R/GCGR) co-agonism can maximise weight loss and improve glycaemic control in type 2 diabetes and obesity. In this study, we investigated the cellular and metabolic effects of modulating the balance between G protein and ß-arrestin-2 recruitment at GLP-1R and GCGR using oxyntomodulin (OXM)-derived co-agonists. This strategy has been previously shown to improve the duration of action of GLP-1R mono-agonists by reducing target desensitisation and downregulation. METHODS: Dipeptidyl dipeptidase-4 (DPP-4)-resistant OXM analogues were generated and assessed for a variety of cellular readouts. Molecular dynamic simulations were used to gain insights into the molecular interactions involved. In vivo studies were performed in mice to identify the effects on glucose homeostasis and weight loss. RESULTS: Ligand-specific reductions in ß-arrestin-2 recruitment were associated with slower GLP-1R internalisation and prolonged glucose-lowering action in vivo. The putative benefits of GCGR agonism were retained, with equivalent weight loss compared to the GLP-1R mono-agonist liraglutide despite a lesser degree of food intake suppression. The compounds tested showed only a minor degree of biased agonism between G protein and ß-arrestin-2 recruitment at both receptors and were best classified as partial agonists for the two pathways measured. CONCLUSIONS: Diminishing ß-arrestin-2 recruitment may be an effective way to increase the therapeutic efficacy of GLP-1R/GCGR co-agonists. These benefits can be achieved by partial rather than biased agonism.

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.

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.

Structure and dynamics of the active Gs-coupled human secretin receptor.

The class B secretin GPCR (SecR) has broad physiological effects, with target potential for treatment of metabolic and cardiovascular disease. Molecular understanding of SecR binding and activation is important for its therapeutic exploitation. We combined cryo-electron microscopy, molecular dynamics, and biochemical cross-linking to determine a 2.3 Å structure, and interrogate dynamics, of secretin bound to the SecR:Gs complex. SecR exhibited a unique organization of its extracellular domain (ECD) relative to its 7-transmembrane (TM) core, forming more extended interactions than other family members. Numerous polar interactions formed between secretin and the receptor extracellular loops (ECLs) and TM helices. Cysteine-cross-linking, cryo-electron microscopy multivariate analysis and molecular dynamics simulations revealed that interactions between peptide and receptor were dynamic, and suggested a model for initial peptide engagement where early interactions between the far N-terminus of the peptide and SecR ECL2 likely occur following initial binding of the peptide C-terminus to the ECD.

New Insights into Key Determinants for Adenosine 1 Receptor Antagonists Selectivity Using Supervised Molecular Dynamics Simulations.

Adenosine receptors (ARs), like many otherGprotein-coupledreceptors (GPCRs), are targets of primary interest indrug design. However, one of the main limits for the development of drugs for this class of GPCRs is the complex selectivity profile usually displayed by ligands. Numerous efforts have been madefor clarifying the selectivity of ARs, leading to the development of many ligand-based models. The structure of the AR subtype A1 (A1AR) has been recently solved,providing important structural insights. In the present work, we rationalized the selectivity profile of two selective A1AR and A2AAR antagonists, investigating their recognition trajectories obtained by Supervised Molecular Dynamics from an unbound state and monitoring the role of the water molecules in the binding site.

A Supervised Molecular Dynamics Approach to Unbiased Ligand-Protein Unbinding.

The recent paradigm shift toward the use of the kinetics parameters in place of thermodynamic constants is leading the computational chemistry community to develop methods for studying the mechanisms of drug binding and unbinding. From this standpoint, molecular dynamics (MD) plays an important role in delivering insight at the molecular scale. However, a known limitation of MD is that the time scales are usually far from those involved in ligand-receptor unbinding events. Here, we show that the algorithm behind supervised MD (SuMD) can simulate the dissociation mechanism of druglike small molecules while avoiding the input of any energy bias to facilitate the transition. SuMD was tested on seven different intermolecular complexes, covering four G protein-coupled receptors: the A2A and A1 adenosine receptors, the orexin 2 and the muscarinic 2 receptors, and the soluble globular enzyme epoxide hydrolase. SuMD well-described the multistep nature of ligand-receptor dissociation, rationalized previous experimental data and produced valuable working hypotheses for structure-kinetics relationships.

A2A and A2B adenosine receptors: The extracellular loop 2 determines high (A2A) or low affinity (A2B) for adenosine.

A2A and A2B adenosine receptors (ARs) are closely related G protein-coupled receptor subtypes, which represent important (potential) drug targets. Despite their almost identical binding sites for adenosine, A2AARs are activated by low (nanomolar) adenosine concentrations, while A2BARs require micromolar concentrations. In the present study, we exchanged the extracellular loop 2 (ECL2) of the human A2AAR for that of the A2BAR. The resulting chimeric A2A(ECL2-A2B)AR was investigated in radioligand binding and cAMP accumulation assays in comparison to the wildtype A2AAR. While the ribose-modified adenosine analog N-ethylcarboxamidoadenosine (NECA) and its 2-substituted derivative CGS-21680 did not exhibit significant changes, adenosine showed dramatically reduced potency and affinity for the A2A(ECL2-A2B)AR mutant displaying similarly low potency as for the wt A2BAR. Supervised molecular dynamics simulation studies predicted a meta-binding site with high affinity for adenosine, but not for NECA, which may contribute to the observed effects.

Revisiting the Allosteric Regulation of Sodium Cation on the Binding of Adenosine at the Human A2A Adenosine Receptor: Insights from Supervised Molecular Dynamics (SuMD) Simulations.

One of the most intriguing findings highlighted from G protein-coupled receptor (GPCR) crystallography is the presence, in many members of class A, of a partially hydrated sodium ion in the middle of the seven transmembrane helices (7TM) bundle. In particular, the human adenosine A2A receptor (A2A AR) is the first GPCR in which a monovalent sodium ion was crystallized in a distal site from the canonical orthosteric one, corroborating, from a structural point of view, its role as a negative allosteric modulator. However, the molecular mechanism by which the sodium ion influences the recognition of the A2A AR agonists is not yet fully understood. In this study, the supervised molecular dynamics (SuMD) technique was exploited to analyse the sodium ion recognition mechanism and how its presence influences the binding of the endogenous agonist adenosine. Due to a higher degree of flexibility of the receptor extracellular (EC) vestibule, we propose the sodium-bound A2A AR as less efficient in stabilizing the adenosine during the different steps of binding.

Could the presence of sodium ion influence the accuracy and precision of the ligand-posing in the human A2A adenosine receptor orthosteric binding site using a molecular docking approach? Insights from Dockbench.

The allosteric modulation of G protein-coupled receptors (GPCRs) by sodium ions has received considerable attention as crystal structures of several receptors, in their inactive conformation, show a Na+ ion bound to specific residues which, in the human A2A adenosine receptor (hA2A AR), are Ser913.39, Trp2466.48, Asn2807.45, and Asn2847.49. A cluster of water molecules completes the coordination of the sodium ion in the putative allosteric site. It is absolutely consolidated that the progress made in the field of GPCRs structural determination has increased the adoption of docking-driven approaches for the identification or the optimization of novel potent and selective ligands. Despite the extensive use of docking protocols in virtual screening approaches, to date, almost any of these studies have been carried out without taking into account the presence of the sodium cation and its first solvation shell in the putative allosteric binding site. In this study, we have focused our attention on determining how the presence of sodium ion binding and additionally its first hydration sphere, in hA2AAR could influence the ligand positioning accuracy during molecular docking simulations for most of the available resting and activated hA2A AR crystal structures, using DockBench as a comparative benchmarking tool and implementing a new correlation coefficient (EM). This work provides indications on the evidence that the posing performance (accuracy and/or precision) of the docking protocols in reproducing the crystallographic poses of different hA2A AR antagonists is generally increased in the presence of the sodium cation and its first solvation shell, in agreement with experimental observations. Consequently, the inclusion of sodium ion and its first solvation shell should be considered in order to facilitate the selection of new potential ligands in all molecular docking-based virtual screening protocols that aim to find novel GPCRs antagonists and inverse agonists.

Could Adenosine Recognize its Receptors with a Stoichiometry Other than 1 : 1?

One of the most largely accepted concepts in the G protein-coupled receptors (GPCRs) field is that the ligand, either agonist or antagonist, recognizes its receptor with a stoichiometry of 1 : 1. Recent experimental evidence, reporting ternary complexes formed by GPCR:orthosteric: allosteric ligands, has complicated the ligand-receptor 1 : 1 binding scenario. Molecular modeling simulations have been used to retrieve insights on the whole ligand-receptor recognition process, beyond information on the final bound state provided by experimental techniques. The simulation of adenosine binding pathways towards the A2A adenosine receptor highlighted the presence of alternative binding sites (meta-binding sites) beside the canonical orthosteric one, mainly in proximity to the extracellular vestibule. In light of all these considerations, we investigated the possibility that a second molecule of adenosine engages its receptor when this is already in the holo form, generating a ternary complex with a stoichiometry of 2 : 1. Unexpectedly, supervised molecular dynamics (SuMD) simulations showed that the A2A adenosine receptor could bind the second molecule of adenosine in one of the possible meta-binding sites as well as into its orthosteric site. The formation of this ternary complex, which favored the formation of the intracellular "ionic lock" between R102 (3.50) and E228 (6.30), could putatively be framed in the context of a negative allosteric regulation.

AquaMMapS: An Alternative Tool to Monitor the Role of Water Molecules During Protein-Ligand Association.

Unquestionably, water appears to be an active player in the noncovalent protein-ligand binding process, as it can either bridge interactions between protein and ligand or can be replaced by the bound ligand. Accordingly, in the last decade, alternative computational methodologies have been sought with the aim of predicting the position and thermodynamic profile of water molecules (i.e., hydration sites) in the binding site using either the ligand-bound or ligand-free protein conformation. Herein, we present an alternative approach, named AquaMMapS, that provides a three-dimensional sampling of putative hydration sites. Interestingly, AquaMMapS can post-inspect molecular dynamics (MD) trajectories obtained from different MD engines using indifferently crystallographic or docking-driven structures as a starting point. Moreover, AquaMMapS is naturally integrated into supervised molecular dynamics (SuMD) simulations, presenting the possibility to inspect hydration sites during the ligand-protein association process. Finally, a penalty scoring method, named AquaMMapScoring(AMS), was developed to evaluate the number and nature of the water molecules displaced by a ligand approaching its binding site during the binding event, guiding a medicinal chemist to explore the most suitable regions of a ligand that can be decorated either with or without interfering with the interaction networks mediated by water molecules with specific recognition regions of the protein.

Comparison of the Human A2A Adenosine Receptor Recognition by Adenosine and Inosine: New Insight from Supervised Molecular Dynamics Simulations.

Adenosine deaminase converts adenosine into inosine. In contrast to adenosine, relatively little attention has been paid to the physiological roles of inosine. Nevertheless, recent studies have demonstrated that inosine has neuroprotective, cardioprotective, immunomodulatory, and antidepressive effects. Inosine was recently shown to be a less potent agonist than adenosine at the A2A adenosine receptor. To better depict the differences in the mechanisms of receptor recognition between adenosine and inosine, we carried out supervised molecular dynamics (SuMD) simulations, and the results are analyzed herein.

Supporting the Identification of Novel Fragment-Based Positive Allosteric Modulators Using a Supervised Molecular Dynamics Approach: A Retrospective Analysis Considering the Human A2A Adenosine Receptor as a Key Example.

Structure-driven fragment-based (SDFB) approaches have provided efficient methods for the identification of novel drug candidates. This strategy has been largely applied in discovering several pharmacological ligand classes, including enzyme inhibitors, receptor antagonists and, more recently, also allosteric (positive and negative) modulators. Recently, Siegal and collaborators reported an interesting study, performed on a detergent-solubilized StaR adenosine A2A receptor, describing the existence of both fragment-like negative allosteric modulators (NAMs), and fragment-like positive allosteric modulators (PAMs). From this retrospective study, our results suggest that Supervised Molecular Dynamics (SuMD) simulations can support, on a reasonable time scale, the identification of fragment-like PAMs following their receptor recognition pathways and characterizing the possible allosteric binding sites.