From structure to clinic: Design of a muscarinic M1 receptor agonist with potential to treatment of Alzheimer's disease.

Current therapies for Alzheimer's disease seek to correct for defective cholinergic transmission by preventing the breakdown of acetylcholine through inhibition of acetylcholinesterase, these however have limited clinical efficacy. An alternative approach is to directly activate cholinergic receptors responsible for learning and memory. The M1-muscarinic acetylcholine (M1) receptor is the target of choice but has been hampered by adverse effects. Here we aimed to design the drug properties needed for a well-tolerated M1-agonist with the potential to alleviate cognitive loss by taking a stepwise translational approach from atomic structure, cell/tissue-based assays, evaluation in preclinical species, clinical safety testing, and finally establishing activity in memory centers in humans. Through this approach, we rationally designed the optimal properties, including selectivity and partial agonism, into HTL9936-a potential candidate for the treatment of memory loss in Alzheimer's disease. More broadly, this demonstrates a strategy for targeting difficult GPCR targets from structure to clinic.

Crystal Structure and Subsequent Ligand Design of a Nonriboside Partial Agonist Bound to the Adenosine A2A Receptor.

In this study, we determined the crystal structure of an engineered human adenosine A2A receptor bound to a partial agonist and compared it to structures cocrystallized with either a full agonist or an antagonist/inverse agonist. The interaction between the partial agonist, belonging to a class of dicyanopyridines, and amino acids in the ligand binding pocket inspired us to develop a small library of derivatives and assess their affinity in radioligand binding studies and potency and intrinsic activity in a functional, label-free, intact cell assay. It appeared that some of the derivatives retained the partial agonist profile, whereas other ligands turned into inverse agonists. We rationalized this remarkable behavior with additional computational docking studies.

X-Ray Crystallography and Free Energy Calculations Reveal the Binding Mechanism of A2A Adenosine Receptor Antagonists.

We present a robust protocol based on iterations of free energy perturbation (FEP) calculations, chemical synthesis, biophysical mapping and X-ray crystallography to reveal the binding mode of an antagonist series to the A2A adenosine receptor (AR). Eight A2A AR binding site mutations from biophysical mapping experiments were initially analyzed with sidechain FEP simulations, performed on alternate binding modes. The results distinctively supported one binding mode, which was subsequently used to design new chromone derivatives. Their affinities for the A2A AR were experimentally determined and investigated through a cycle of ligand-FEP calculations, validating the binding orientation of the different chemical substituents proposed. Subsequent X-ray crystallography of the A2A AR with a low and a high affinity chromone derivative confirmed the predicted binding orientation. The new molecules and structures here reported were driven by free energy calculations, and provide new insights on antagonist binding to the A2A AR, an emerging target in immuno-oncology.

Comparison of Orexin 1 and Orexin 2 Ligand Binding Modes Using X-ray Crystallography and Computational Analysis.

The orexin system, which consists of the two G protein-coupled receptors OX1 and OX2, activated by the neuropeptides OX-A and OX-B, is firmly established as a key regulator of behavioral arousal, sleep, and wakefulness and has been an area of intense research effort over the past two decades. X-ray structures of the receptors in complex with 10 new antagonist ligands from diverse chemotypes are presented, which complement the existing structural information for the system and highlight the critical importance of lipophilic hotspots and water molecules for these peptidergic GPCR targets. Learnings from the structural information regarding the utility of pharmacophore models and how selectivity between OX1 and OX2 can be achieved are discussed.

Identification of a novel allosteric GLP-1R antagonist HTL26119 using structure- based drug design.

A series of novel allosteric antagonists of the GLP-1 receptor (GLP-1R), exemplified by HTL26119, are described. SBDD approaches were employed to identify HTL26119, exploiting structural understanding of the allosteric binding site of the closely related Glucagon receptor (GCGR) (Jazayeri et al., 2016) and the homology relationships between GCGR and GLP-1R. The region around residue C3476.36b of the GLP-1R receptor represents a key difference from GCGR and was targeted for selectivity for GLP-1R.

Agonists and Antagonists of Protease-Activated Receptor 2 Discovered within a DNA-Encoded Chemical Library Using Mutational Stabilization of the Target.

The discovery of ligands via affinity-mediated selection of DNA-encoded chemical libraries is driven by the quality and concentration of the protein target. G-protein-coupled receptors (GPCRs) and other membrane-bound targets can be difficult to isolate in their functional state and at high concentrations, and therefore have been challenging for affinity-mediated selection. Here, we report a successful selection campaign against protease-activated receptor 2 (PAR2). Using a thermo-stabilized mutant of PAR2, we conducted affinity selection using our >100-billion-compound DNA-encoded library. We observed a number of putative ligands enriched upon selection, and subsequent cellular profiling revealed these ligands to comprise both agonists and antagonists. The agonist series shared structural similarity with known agonists. The antagonists were shown to bind in a novel allosteric binding site on the PAR2 protein. This report serves to demonstrate that cell-free affinity selection against GPCRs can be achieved with mutant stabilized protein targets.

Structure of the complement C5a receptor bound to the extra-helical antagonist NDT9513727.

The complement system is a crucial component of the host response to infection and tissue damage. Activation of the complement cascade generates anaphylatoxins including C5a and C3a. C5a exerts a pro-inflammatory effect via the G-protein-coupled receptor C5a anaphylatoxin chemotactic receptor 1 (C5aR1, also known as CD88) that is expressed on cells of myeloid origin. Inhibitors of the complement system have long been of interest as potential drugs for the treatment of diseases such as sepsis, rheumatoid arthritis, Crohn's disease and ischaemia-reperfusion injuries. More recently, a role of C5a in neurodegenerative conditions such as Alzheimer's disease has been identified. Peptide antagonists based on the C5a ligand have progressed to phase 2 trials in psoriasis and rheumatoid arthritis; however, these compounds exhibited problems with off-target activity, production costs, potential immunogenicity and poor oral bioavailability. Several small-molecule competitive antagonists for C5aR1, such as W-54011 and NDT9513727, have been identified by C5a radioligand-binding assays. NDT9513727 is a non-peptide inverse agonist of C5aR1, and is highly selective for the primate and gerbil receptors over those of other species. Here, to study the mechanism of action of C5a antagonists, we determine the structure of a thermostabilized C5aR1 (known as C5aR1 StaR) in complex with NDT9513727. We found that the small molecule bound between transmembrane helices 3, 4 and 5, outside the helical bundle. One key interaction between the small molecule and residue Trp2135.49 seems to determine the species selectivity of the compound. The structure demonstrates that NDT9513727 exerts its inverse-agonist activity through an extra-helical mode of action.

Towards high throughput GPCR crystallography: In Meso soaking of Adenosine A2A Receptor crystals.

Here we report an efficient method to generate multiple co-structures of the A2A G protein-coupled receptor (GPCR) with small-molecules from a single preparation of a thermostabilised receptor crystallised in Lipidic Cubic Phase (LCP). Receptor crystallisation is achieved following purification using a low affinity "carrier" ligand (theophylline) and crystals are then soaked in solutions containing the desired (higher affinity) compounds. Complete datasets to high resolution can then be collected from single crystals and seven structures are reported here of which three are novel. The method significantly improves structural throughput for ligand screening using stabilised GPCRs, thereby actively driving Structure-Based Drug Discovery (SBDD).

Structural Mapping of Adenosine Receptor Mutations: Ligand Binding and Signaling Mechanisms.

The four adenosine receptors (ARs), A1, A2A, A2B, and A3, constitute a subfamily of G protein-coupled receptors (GPCRs) with exceptional foundations for structure-based ligand design. The vast amount of mutagenesis data, accumulated in the literature since the 1990s, has been recently supplemented with structural information, currently consisting of several inactive and active structures of the A2A and inactive conformations of the A1 ARs. We provide the first integrated view of the pharmacological, biochemical, and structural data available for this receptor family, by mapping onto the relevant crystal structures all site-directed mutagenesis data, curated and deposited at the GPCR database (available through http://www.gpcrdb.org). This analysis provides novel insights into ligand binding, allosteric modulation, and signaling of the AR family.

Serial millisecond crystallography for routine room-temperature structure determination at synchrotrons.

Historically, room-temperature structure determination was succeeded by cryo-crystallography to mitigate radiation damage. Here, we demonstrate that serial millisecond crystallography at a synchrotron beamline equipped with high-viscosity injector and high frame-rate detector allows typical crystallographic experiments to be performed at room-temperature. Using a crystal scanning approach, we determine the high-resolution structure of the radiation sensitive molybdenum storage protein, demonstrate soaking of the drug colchicine into tubulin and native sulfur phasing of the human G protein-coupled adenosine receptor. Serial crystallographic data for molecular replacement already converges in 1,000-10,000 diffraction patterns, which we collected in 3 to maximally 82 minutes. Compared with serial data we collected at a free-electron laser, the synchrotron data are of slightly lower resolution, however fewer diffraction patterns are needed for de novo phasing. Overall, the data we collected by room-temperature serial crystallography are of comparable quality to cryo-crystallographic data and can be routinely collected at synchrotrons.Serial crystallography was developed for protein crystal data collection with X-ray free-electron lasers. Here the authors present several examples which show that serial crystallography using high-viscosity injectors can also be routinely employed for room-temperature data collection at synchrotrons.

Structures of Human A1 and A2A Adenosine Receptors with Xanthines Reveal Determinants of Selectivity.

The adenosine A1 and A2A receptors belong to the purinergic family of G protein-coupled receptors, and regulate diverse functions of the cardiovascular, respiratory, renal, inflammation, and CNS. Xanthines such as caffeine and theophylline are weak, non-selective antagonists of adenosine receptors. Here we report the structure of a thermostabilized human A1 receptor at 3.3 Å resolution with PSB36, an A1-selective xanthine-based antagonist. This is compared with structures of the A2A receptor with PSB36 (2.8 Å resolution), caffeine (2.1 Å), and theophylline (2.0 Å) to highlight features of ligand recognition which are common across xanthines. The structures of A1R and A2AR were analyzed to identify the differences that are important selectivity determinants for xanthine ligands, and the role of T2707.35 in A1R (M2707.35 in A2AR) in conferring selectivity was confirmed by mutagenesis. The structural differences confirmed to lead to selectivity can be utilized in the design of new subtype-selective A1R or A2AR antagonists.

Corrigendum: Crystal structure of the GLP-1 receptor bound to a peptide agonist.

This corrects the article DOI: 10.1038/nature22800.

Crystal structure of the GLP-1 receptor bound to a peptide agonist.

Glucagon-like peptide 1 (GLP-1) regulates glucose homeostasis through the control of insulin release from the pancreas. GLP-1 peptide agonists are efficacious drugs for the treatment of diabetes. To gain insight into the molecular mechanism of action of GLP-1 peptides, here we report the crystal structure of the full-length GLP-1 receptor bound to a truncated peptide agonist. The peptide agonist retains an α-helical conformation as it sits deep within the receptor-binding pocket. The arrangement of the transmembrane helices reveals hallmarks of an active conformation similar to that observed in class A receptors. Guided by this structural information, we design peptide agonists with potent in vivo activity in a mouse model of diabetes.

Structural insight into allosteric modulation of protease-activated receptor 2.

Protease-activated receptors (PARs) are a family of G-protein-coupled receptors (GPCRs) that are irreversibly activated by proteolytic cleavage of the N terminus, which unmasks a tethered peptide ligand that binds and activates the transmembrane receptor domain, eliciting a cellular cascade in response to inflammatory signals and other stimuli. PARs are implicated in a wide range of diseases, such as cancer and inflammation. PARs have been the subject of major pharmaceutical research efforts but the discovery of small-molecule antagonists that effectively bind them has proved challenging. The only marketed drug targeting a PAR is vorapaxar, a selective antagonist of PAR1 used to prevent thrombosis. The structure of PAR1 in complex with vorapaxar has been reported previously. Despite sequence homology across the PAR isoforms, discovery of PAR2 antagonists has been less successful, although GB88 has been described as a weak antagonist. Here we report crystal structures of PAR2 in complex with two distinct antagonists and a blocking antibody. The antagonist AZ8838 binds in a fully occluded pocket near the extracellular surface. Functional and binding studies reveal that AZ8838 exhibits slow binding kinetics, which is an attractive feature for a PAR2 antagonist competing against a tethered ligand. Antagonist AZ3451 binds to a remote allosteric site outside the helical bundle. We propose that antagonist binding prevents structural rearrangements required for receptor activation and signalling. We also show that a blocking antibody antigen-binding fragment binds to the extracellular surface of PAR2, preventing access of the tethered ligand to the peptide-binding site. These structures provide a basis for the development of selective PAR2 antagonists for a range of therapeutic uses.

Kinetics for Drug Discovery: an industry-driven effort to target drug residence time.

A considerable number of approved drugs show non-equilibrium binding characteristics, emphasizing the potential role of drug residence times for in vivo efficacy. Therefore, a detailed understanding of the kinetics of association and dissociation of a target-ligand complex might provide crucial insight into the molecular mechanism-of-action of a compound. This deeper understanding will help to improve decision making in drug discovery, thus leading to a better selection of interesting compounds to be profiled further. In this review, we highlight the contributions of the Kinetics for Drug Discovery (K4DD) Consortium, which targets major open questions related to binding kinetics in an industry-driven public-private partnership.

Impact of protein-ligand solvation and desolvation on transition state thermodynamic properties of adenosine A2A ligand binding kinetics.

Ligand-protein binding kinetic rates are growing in importance as parameters to consider in drug discovery and lead optimization. In this study we analysed using surface plasmon resonance (SPR) the transition state (TS) properties of a set of six adenosine A2A receptor inhibitors, belonging to both the xanthine and the triazolo-triazine scaffolds. SPR highlighted interesting differences among the ligands in the enthalpic and entropic components of the TS energy barriers for the binding and unbinding events. To better understand at a molecular level these differences, we developed suMetaD, a novel molecular dynamics (MD)-based approach combining supervised MD and metadynamics. This method allows simulation of the ligand unbinding and binding events. It also provides the system conformation corresponding to the highest energy barrier the ligand is required to overcome to reach the final state. For the six ligands evaluated in this study their TS thermodynamic properties were linked in particular to the role of water molecules in solvating/desolvating the pocket and the small molecules. suMetaD identified kinetic bottleneck conformations near the bound state position or in the vestibule area. In the first case the barrier is mainly enthalpic, requiring the breaking of strong interactions with the protein. In the vestibule TS location the kinetic bottleneck is instead mainly of entropic nature, linked to the solvent behaviour.

Intracellular allosteric antagonism of the CCR9 receptor.

Chemokines and their G-protein-coupled receptors play a diverse role in immune defence by controlling the migration, activation and survival of immune cells. They are also involved in viral entry, tumour growth and metastasis and hence are important drug targets in a wide range of diseases. Despite very significant efforts by the pharmaceutical industry to develop drugs, with over 50 small-molecule drugs directed at the family entering clinical development, only two compounds have reached the market: maraviroc (CCR5) for HIV infection and plerixafor (CXCR4) for stem-cell mobilization. The high failure rate may in part be due to limited understanding of the mechanism of action of chemokine antagonists and an inability to optimize compounds in the absence of structural information. CC chemokine receptor type 9 (CCR9) activation by CCL25 plays a key role in leukocyte recruitment to the gut and represents a therapeutic target in inflammatory bowel disease. The selective CCR9 antagonist vercirnon progressed to phase 3 clinical trials in Crohn's disease but efficacy was limited, with the need for very high doses to block receptor activation. Here we report the crystal structure of the CCR9 receptor in complex with vercirnon at 2.8 Å resolution. Remarkably, vercirnon binds to the intracellular side of the receptor, exerting allosteric antagonism and preventing G-protein coupling. This binding site explains the need for relatively lipophilic ligands and describes another example of an allosteric site on G-protein-coupled receptors that can be targeted for drug design, not only at CCR9, but potentially extending to other chemokine receptors.

Controlling the Dissociation of Ligands from the Adenosine A2A Receptor through Modulation of Salt Bridge Strength.

The association and dissociation kinetics of ligands binding to proteins vary considerably, but the mechanisms behind this variability are poorly understood, limiting their utilization for drug discovery. This is particularly so for G protein-coupled receptors (GPCRs) where high resolution structural information is only beginning to emerge. Engineering the human A2A adenosine receptor has allowed structures to be solved in complex with the reference compound ZM241385 and four related ligands at high resolution. Differences between the structures are limited, with the most pronounced being the interaction of each ligand with a salt bridge on the extracellular side of the receptor. Mutagenesis experiments confirm the role of this salt bridge in controlling the dissociation kinetics of the ligands from the receptor, while molecular dynamics simulations demonstrate the ability of ligands to modulate salt bridge stability. These results shed light on a structural determinant of ligand dissociation kinetics and identify a means by which this property may be optimized.

Extra-helical binding site of a glucagon receptor antagonist.

Glucagon is a 29-amino-acid peptide released from the α-cells of the islet of Langerhans, which has a key role in glucose homeostasis. Glucagon action is transduced by the class B G-protein-coupled glucagon receptor (GCGR), which is located on liver, kidney, intestinal smooth muscle, brain, adipose tissue, heart and pancreas cells, and this receptor has been considered an important drug target in the treatment of diabetes. Administration of recently identified small-molecule GCGR antagonists in patients with type 2 diabetes results in a substantial reduction of fasting and postprandial glucose concentrations. Although an X-ray structure of the transmembrane domain of the GCGR has previously been solved, the ligand (NNC0640) was not resolved. Here we report the 2.5 Å structure of human GCGR in complex with the antagonist MK-0893 (ref. 4), which is found to bind to an allosteric site outside the seven transmembrane (7TM) helical bundle in a position between TM6 and TM7 extending into the lipid bilayer. Mutagenesis of key residues identified in the X-ray structure confirms their role in the binding of MK-0893 to the receptor. The unexpected position of the binding site for MK-0893, which is structurally similar to other GCGR antagonists, suggests that glucagon activation of the receptor is prevented by restriction of the outward helical movement of TM6 required for G-protein coupling. Structural knowledge of class B receptors is limited, with only one other ligand-binding site defined--for the corticotropin-releasing hormone receptor 1 (CRF1R)--which was located deep within the 7TM bundle. We describe a completely novel allosteric binding site for class B receptors, providing an opportunity for structure-based drug design for this receptor class and furthering our understanding of the mechanisms of activation of these receptors.

Purification of Stabilized GPCRs for Structural and Biophysical Analyses.

G protein-coupled receptors (GPCRs) are of particular importance for drug discovery, being the targets of many existing drugs, and being linked to many diseases where new therapies are required. However, as integral membrane proteins, they are generally unstable when removed from their membrane environment, precluding them from the wide range of structural and biophysical techniques which can be applied to soluble proteins such as kinases. Through the use of protein engineering methods, mutations can be identified which both increase the thermostability of GPCRs when purified in detergent, as well as biasing the receptor toward a specific physiologically relevant conformational state. The resultant stabilized receptor (known as a StaR) can be purified in multiple-milligram quantities, whilst retaining correct folding, thus enabling the generation of reagents suitable for a broad range of structural and biophysical studies. Example protocols for the purification of StaR proteins for analysis, ligand screening with the thiol-specific fluorochrome N-[4-(7-diethylamino-4-methyl-3-coumarinyl)phenyl]maleimide (CPM), surface plasmon resonance (SPR), and crystallization for structural studies are presented.