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

G Protein-Coupled Receptor-Ligand Pose and Functional Class Prediction

G protein-coupled receptor (GPCR) transmembrane protein family members play essential roles in physiology. Numerous pharmaceuticals target GPCRs, and many drug discovery programs utilize virtual screening (VS) against GPCR targets. Improvements in the accuracy of predicting new molecules that bind to and either activate or inhibit GPCR function would accelerate such drug discovery programs. This work addresses two significant research questions. First, do ligand interaction fingerprints provide a substantial advantage over automated methods of binding site selection for classical docking? Second, can the functional status of prospective screening candidates be predicted from ligand interaction fingerprints using a random forest classifier? Ligand interaction fingerprints were found to offer modest advantages in sampling accurate poses, but no substantial advantage in the final set of top-ranked poses after scoring, and, thus, were not used in the generation of the ligand-receptor complexes used to train and test the random forest classifier. A binary classifier which treated agonists, antagonists, and inverse agonists as active and all other ligands as inactive proved highly effective in ligand function prediction in an external test set of GPR31 and TAAR2 candidate ligands with a hit rate of 82.6% actual actives within the set of predicted actives.

Rational Design of Drugs Targeting G-Protein-Coupled Receptors: A Structural Biology Perspective

G protein-coupled receptors (GPCRs) play a key role in the transduction of extracellular signals to cells and regulation of many biological processes, which makes these membrane proteins one of the most important targets for pharmacological agents. A significant increase in the number of resolved atomic structures of GPCRs has opened the possibility of developing pharmaceuticals targeting these receptors via structure-based drug design (SBDD). SBDD employs information on the structure of receptor-ligand complexes to search for selective ligands without the need for an extensive high-throughput experimental ligand screening and can significantly expand the chemical space for ligand search. In this review, we describe the process of deciphering GPCR structures using X-ray diffraction analysis and cryoelectron microscopy as an important stage in the rational design of drugs targeting this receptor class. Our main goal was to present modern developments and key features of experimental methods used in SBDD of GPCR-targeting agents to a wide range of specialists.

Experimental phasing opportunities for macromolecular crystallography at very long wavelengths

Despite recent advances in cryo-electron microscopy and artificial intelligence-based model predictions, a significant fraction of structure determinations by macromolecular crystallography still requires experimental phasing, usually by means of single-wavelength anomalous diffraction (SAD) techniques. Most synchrotron beamlines provide highly brilliant beams of X-rays of between 0.7 and 2 Å wavelength. Use of longer wavelengths to access the absorption edges of biologically important lighter atoms such as calcium, potassium, chlorine, sulfur and phosphorus for native-SAD phasing is attractive but technically highly challenging. The long-wavelength beamline I23 at Diamond Light Source overcomes these limitations and extends the accessible wavelength range to λ = 5.9 Å. Here we report 22 macromolecular structures solved in this extended wavelength range, using anomalous scattering from a range of elements which demonstrate the routine feasibility of lighter atom phasing. We suggest that, in light of its advantages, long-wavelength crystallography is a compelling option for experimental phasing.

A versatile approach to high-density microcrystals in lipidic cubic phase for room-temperature serial crystallography

Serial crystallography has emerged as an important tool for structural studies of integral membrane proteins. The ability to collect data from micrometre-sized weakly diffracting crystals at room temperature with minimal radiation damage has opened many new opportunities in time-resolved studies and drug discovery. However, the production of integral membrane protein microcrystals in lipidic cubic phase at the desired crystal density and quantity is challenging. This paper introduces VIALS (versatile approach to high-density microcrystals in lipidic cubic phase for serial crystallography), a simple, fast and efficient method for preparing hundreds of microlitres of high-density microcrystals suitable for serial X-ray diffraction experiments at both synchrotron and free-electron laser sources. The method is also of great benefit for rational structure-based drug design as it facilitates in situ crystal soaking and rapid determination of many co-crystal structures. Using the VIALS approach, room-temperature structures are reported of (i) the archaerhodopsin-3 protein in its dark-adapted state and 110 ns photocycle intermediate, determined to 2.2 and 1.7 Å, respectively, and (ii) the human A2A adenosine receptor in complex with two different ligands determined to a resolution of 3.5 Å.

Molecular Biophysics of Class A G Protein Coupled Receptors-Lipids Interactome at a Glance-Highlights from the A2A Adenosine Receptor

G protein-coupled receptors (GPCRs) are embedded in phospholipid membrane bilayers with cholesterol representing 34% of the total lipid content in mammalian plasma membranes. Membrane lipids interact with GPCRs structures and modulate their function and drug-stimulated signaling through conformational selection. It has been shown that anionic phospholipids form strong interactions between positively charged residues in the G protein and the TM5-TM6-TM 7 cytoplasmic interface of class A GPCRs stabilizing the signaling GPCR-G complex. Cholesterol with a high content in plasma membranes can be identified in more specific sites in the transmembrane region of GPCRs, such as the Cholesterol Consensus Motif (CCM) and Cholesterol Recognition Amino Acid Consensus (CRAC) motifs and other receptor dependent and receptor state dependent sites. Experimental biophysical methods, atomistic (AA) MD simulations and coarse-grained (CG) molecular dynamics simulations have been applied to investigate these interactions. We emphasized here the impact of phosphatidyl inositol-4,5-bisphosphate (PtdIns(4,5)P₂ or PIP₂), a minor phospholipid component and of cholesterol on the function-related conformational equilibria of the human A_2A adenosine receptor (A_2AR), a representative receptor in class A GPCR. Several GPCRs of class A interacted with PIP₂ and cholesterol and in many cases the mechanism of the modulation of their function remains unknown. This review provides a helpful comprehensive overview for biophysics that enter the field of GPCRs-lipid systems.

Crystal structure of adenosine A2A receptor in complex with clinical candidate Etrumadenant reveals unprecedented antagonist interaction

The G_s protein-coupled adenosine A_2A receptor (A_2AAR) represents an emerging drug target for cancer immunotherapy. The clinical candidate Etrumadenant was developed as an A_2AAR antagonist with ancillary blockade of the A_2BAR subtype. It constitutes a unique chemotype featuring a poly-substituted 2-amino-4-phenyl-6-triazolylpyrimidine core structure. Herein, we report two crystal structures of the A_2AAR in complex with Etrumadenant, obtained with differently thermostabilized A_2AAR constructs. This led to the discovery of an unprecedented interaction, a hydrogen bond of T88^3.36 with the cyano group of Etrumadenant. T88^3.36 is mutated in most A_2AAR constructs used for crystallization, which has prevented the discovery of its interactions. In-vitro characterization of Etrumadenant indicated low selectivity versus the A₁AR subtype, which can be rationalized by the structural data. These results will facilitate the future design of AR antagonists with desired selectivity. Moreover, they highlight the advantages of the employed A_2AAR crystallization construct that is devoid of ligand binding site mutations.

Comparative Study of Receptor-, Receptor State-, and Membrane-Dependent Cholesterol Binding Sites in A2A and A1 Adenosine Receptors Using Coarse-Grained Molecular Dynamics Simulations

We used coarse-grained molecular dynamics (CG MD) simulations to study protein-cholesterol interactions for different activation states of the A_2A adenosine receptor (A_2AR) and the A₁ adenosine receptor (A₁R) and predict new cholesterol binding sites indicating amino acid residues with a high residence time in three biologically relevant membranes. Compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-cholesterol and POPC-phosphatidylinositol-bisphosphate (PIP₂)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A_2AR and revealed the binding sites with long-lasting amino acid residues. We observed that using the plasma mimetic membrane and plotting residues with cholesterol residence time ≥2 μs, our CG MD simulations captured most obviously the cholesterol-protein interactions. For the inactive A_2AR, we identified one more binding site in which cholesterol is bound to residues with a long residence time compared to the previously detected, for the active A₁R, three binding sites, and for the inactive A₁R, two binding sites. We calculated that for the active states, cholesterol binds to residues with a much longer residence time compared to the inactive state for both A_2AR and A₁R. The stability of the identified binding sites to A₁R or A_2AR with CG MD simulations was additionally investigated with potential of mean force calculations using umbrella sampling. We observed that the binding sites with residues to which cholesterol has a long residence time in A_2AR have shallow binding free energy minima compared to the related binding sites in A₁R, suggesting a stronger binding for cholesterol to A₁R. The differences in binding sites in which cholesterol is stabilized and interacts with residues with a long residence time between active and inactive states of A₁R and A_2AR can be important for differences in functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators, which can bind through lipid pathways. We observed a stronger binding for cholesterol to A₁R (i.e., generally higher association rates) compared to A_2AR, which remains to be demonstrated. For the active states, cholesterol binds to residues with much longer residence times compared to the inactive state for both A_2AR and A₁R. Taken together, binding sites of active A₁R may be considered as promising allosteric targets.

Systematic analyses of the sequence conservation and ligand interaction patterns of purinergic P1 and P2Y receptors provide a structural basis for receptor selectivity

Purinergic receptors are membrane proteins that regulate numerous cellular functions by catalyzing reactions involving purine nucleotides or nucleosides. Among the three receptor families, i.e., P1, P2X, and P2Y, the P1 and P2Y receptors share common structural features of class A GPCR. Comprehensive sequence and structural analysis revealed that the P1 and P2Y receptors belong to two distinct groups. They exhibit different ligand-binding site features that can distinguish between specific activators. These specific amino acid residues in the binding cavity may be involved in the selectivity and unique pharmacological behavior of each subtype. In this study, we conducted a structure-based analysis of purinergic P1 and P2Y receptors to identify their evolutionary signature and obtain structural insights into ligand recognition and selectivity. The structural features of the P1 and P2Y receptor classes were compared based on sequence conservation and ligand interaction patterns. Orthologous protein sequences were collected for the P1 and P2Y receptors, and sequence conservation was calculated based on Shannon entropy to identify highly conserved residues. To analyze the ligand interaction patterns, we performed docking studies on the P1 and P2Y receptors using known ligand information extracted from the ChEMBL database. We analyzed how the conserved residues are related to ligand-binding sites and how the key interacting residues differ between P1 and P2Y receptors, or between agonists and antagonists. We extracted new similarities and differences between the receptor subtypes, and the results can be used for designing new ligands by predicting hotspot residues that are important for functional selectivity.

GPCR Agonist-to-Antagonist Conversion: Enabling the Design of Nucleoside Functional Switches for the A2A Adenosine Receptor

Modulators of the G protein-coupled A_2A adenosine receptor (A_2AAR) have been considered promising agents to treat Parkinson’s disease, inflammation, cancer, and central nervous system disorders. Herein, we demonstrate that a thiophene modification at the C8 position in the common adenine scaffold converted an A_2AAR agonist into an antagonist. We synthesized and characterized a novel A_2AAR antagonist, 2 (LJ-4517), with K_i = 18.3 nM. X-ray crystallographic structures of 2 in complex with two thermostabilized A_2AAR constructs were solved at 2.05 and 2.80 Å resolutions. In contrast to A_2AAR agonists, which simultaneously interact with both Ser277^7.42 and His278^7.43, 2 only transiently contacts His278^7.43, which can be direct or water-mediated. The n-hexynyl group of 2 extends into an A_2AAR exosite. Structural analysis revealed that the introduced thiophene modification restricted receptor conformational rearrangements required for subsequent activation. This approach can expand the repertoire of adenosine receptor antagonists that can be designed based on available agonist scaffolds.

Adenosine A2A receptor antagonists: from caffeine to selective non-xanthines

A long evolution of knowledge of the psychostimulant caffeine led in the 1960s to another purine natural product, adenosine and its A2A receptor. Adenosine is a short-lived autocrine/paracrine mediator that acts pharmacologically at four different adenosine receptors in a manner opposite to the pan-antagonist caffeine and serves as an endogenous allostatic regulator. Although detrimental in the developing brain, caffeine appears to be cerebroprotective in aging. Moderate caffeine consumption in adults, except in pregnancy, may also provide benefit in pain, diabetes, and kidney and liver disorders. Inhibition of A2A receptors is one of caffeine's principal effects and we now understand this interaction at the atomic level. The A2A receptor has become a prototypical example of utilizing high-resolution structures of GPCRs for the rational design of chemically diverse drug molecules. The previous focus on discovery of selective A2A receptor antagonists for neurodegenerative diseases has expanded to include immunotherapy for cancer, and clinical trials have ensued. LINKED ARTICLES: This article is part of a themed issue on Structure Guided Pharmacology of Membrane Proteins (BJP 75th Anniversary). To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v179.14/issuetoc.

Translocation Mechanism of Allosteric Sodium Ions in β2-Adrenoceptor

The allosteric modulation of G-protein-coupled receptors (GPCRs) by sodium ions has received significant attention as the crystal structures of several receptors show the binding of sodium ions (Na⁺) at the conserved D^2.50. Theoretical studies have shown that extracellular Na⁺ would enter the allosteric D^2.50 via the orthosteric site. However, it remains unclear how the bound allosteric Na⁺ would leave the GPCRs. In this study, we performed molecular dynamics (MD) simulations to illustrate the energy barriers of Na⁺ transfer through the transmembrane helix bundle of β₂AR. In contrast to the postulations from other GPCRs, the translocation of this allosteric Na⁺ into the intracellular side is found to be significantly difficult. Hence, the translocation direction could be receptor-specific.

Single Stabilizing Point Mutation Enables High-Resolution Co-Crystal Structures of the Adenosine A2A Receptor with Preladenant Conjugates

The G protein-coupled adenosine A_2A receptor (A_2A AR) is an important new (potential) drug target in immuno-oncology, and for neurodegenerative diseases. Preladenant and its derivatives belong to the most potent A_2A AR antagonists displaying exceptional selectivity. While crystal structures of the human A_2A AR have been solved, mostly using the A_2A -StaR2 protein that bears 9 point mutations, co-crystallization with Preladenant derivatives has so far been elusive. We developed a new A_2A AR construct harboring a single point mutation (S91^3.39 K) which renders it extremely thermostable. This allowed the co-crystallization of two novel Preladenant derivatives, the polyethylene glycol-conjugated (PEGylated) PSB-2113, and the fluorophore-labeled PSB-2115. The obtained crystal structures (2.25 Å and 2.6 Å resolution) provide explanations for the high potency and selectivity of Preladenant derivatives. They represent the first crystal structures of a GPCR in complex with PEG- and fluorophore-conjugated ligands. The applied strategy is predicted to be applicable to further class A GPCRs.

Yeast-based directed-evolution for high-throughput structural stabilization of G protein-coupled receptors (GPCRs)

The immense potential of G protein-coupled receptors (GPCRs) as targets for drug discovery is not fully realized due to the enormous difficulties associated with structure elucidation of these profoundly unstable membrane proteins. The existing methods of GPCR stability-engineering are cumbersome and low-throughput; in addition, the scope of GPCRs that could benefit from these techniques is limited. Here, we present a yeast-based screening platform for a single-step isolation of GRCR variants stable in the presence of short-chain detergents, a feature essential for their successful crystallization using vapor diffusion method. The yeast detergent-resistant cell wall presents a unique opportunity for compartmentalization, to physically link the receptor's phenotype to its encoding DNA, and thus enable discovery of stable GPCR variants with unprecedent efficiency. The scope of mutations identified by the method reveals a surprising amenability of the GPCR scaffold to stabilization, and suggests an intriguing possibility of amending the stability properties of GPCR by varying the structural status of the C-terminus.

Purinergic receptors modulators: An emerging pharmacological tool for disease management

Purinergic signaling is mediated through extracellular nucleotides (adenosine 5'-triphosphate, uridine-5'-triphosphate, adenosine diphosphate, uridine-5'-diphosphate, and adenosine) that serve as signaling molecules. In the early 1990s, purines and pyrimidine receptors were cloned and characterized drawing the attention of scientists toward this aspect of cellular signaling. This signaling pathway is comprised of four subtypes of adenosine receptors (P1), eight subtypes of G-coupled protein receptors (P2YRs), and seven subtypes of ligand-gated ionotropic receptors (P2XRs). In current studies, the pathophysiology and therapeutic potentials of these receptors have been focused on. Various ligands, modulating the functions of purinergic receptors, are in current clinical practices for the treatment of various neurodegenerative disorders and cardiovascular diseases. Moreover, several purinergic receptors ligands are in advanced phases of clinical trials as a remedy for depression, epilepsy, autism, osteoporosis, atherosclerosis, myocardial infarction, diabetes, irritable bowel syndrome, and cancers. In the present study, agonists and antagonists of purinergic receptors have been summarized that may serve as pharmacological tools for drug design and development.

International Union of Basic and Clinical Pharmacology. CXII: Adenosine Receptors: A Further Update

Our previous International Union of Basic and Clinical Pharmacology report on the nomenclature and classification of adenosine receptors (2011) contained a number of emerging developments with respect to this G protein-coupled receptor subfamily, including protein structure, protein oligomerization, protein diversity, and allosteric modulation by small molecules. Since then, a wealth of new data and results has been added, allowing us to explore novel concepts such as target binding kinetics and biased signaling of adenosine receptors, to examine a multitude of receptor structures and novel ligands, to gauge new pharmacology, and to evaluate clinical trials with adenosine receptor ligands. This review should therefore be considered a further update of our previous reports from 2001 and 2011. SIGNIFICANCE STATEMENT: Adenosine receptors (ARs) are of continuing interest for future treatment of chronic and acute disease conditions, including inflammatory diseases, neurodegenerative afflictions, and cancer. The design of AR agonists ("biased" or not) and antagonists is largely structure based now, thanks to the tremendous progress in AR structural biology. The A2A- and A2BAR appear to modulate the immune response in tumor biology. Many clinical trials for this indication are ongoing, whereas an A2AAR antagonist (istradefylline) has been approved as an anti-Parkinson agent.

Acoustic levitation and rotation of thin films and their application for room temperature protein crystallography

Acoustic levitation has attracted attention in terms of chemical and biochemical analysis in combination with various analytical methods because of its unique container-less environment for samples that is not reliant on specific material characteristics. However, loading samples with very high viscosity is difficult. To expand the scope, we propose the use of polymer thin films as sample holders, whereby the sample is dispensed on a film that is subsequently loaded onto an acoustic levitator. When applied for protein crystallography experiments, rotation controllability and positional stability are important prerequisites. We therefore study the acoustic levitation and rotation of thin films with an aspect ratio (the diameter-to-thickness ratio) of 80-240, which is an order of magnitude larger than those reported previously. For films with empirically optimized shapes, we find that it is possible to control the rotation speed in the range of 1-4 rotations per second while maintaining a positional stability of 12 ± 5 µm. The acoustic radiation force acting on the films is found to be a factor of 26-30 higher than that for same-volume water droplets. We propose use cases of the developed films for protein crystallography experiments and demonstrate data collections for large single crystal samples at room temperature.

Dual-acting antitumor agents targeting the A2A adenosine receptor and histone deacetylases: Design and synthesis of 4-(furan-2-yl)-1H-pyrazolo[3,4-d]pyrimidin-6-amine derivatives

Based on its inhibition by antagonists, the A_2A adenosine receptor (A_2AAR) has attracted attention as an anti-tumor drug target; however, in preclinical models and clinical trials, A_2AAR antagonists have so far shown only limited efficacy as standalone therapies. The design of dual-acting compounds, targeting the A_2AAR and histone deacetylases (HDACs), is used here as an approach to the discovery of novel and more potent antitumor agents. Based on the core structures of the A_2AAR antagonists V-2006 and CPI-444, novel 4-(furan-2-yl)-1H-pyrazolo[3,4-d]pyrimidin-6-amine derivatives were designed as such dual-acting compounds. The binding affinities for A_2AAR of all the new compounds were tested, and their HDAC inhibitory activity was evaluated. Compounds with balanced A_2AAR antagonism and HDAC inhibition were tested for their in vitro anti-proliferative activity and pharmacokinetic properties. One of the compounds, 14c (4-(2-(6-Amino-4-(furan-2-yl)-1H-pyrazolo[3,4-d]pyrimidin-1-yl)ethyl)-N-(2-amino-phenyl)benzamide) showed an overall favorable pharmacokinetic profile; in the mouse MC38 xenograft model, it showed potent anti-tumor effects with inhibition rates of 44% (90 mg/kg, po, bid) and 85% (60 mg/kg, ip, bid), respectively.

Molecular Simulations and Drug Discovery of Adenosine Receptors

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.

Structure determination of GPCRs: cryo-EM compared with X-ray crystallography

G protein-coupled receptors (GPCRs) are the largest single family of cell surface receptors encoded by the human genome and they play pivotal roles in co-ordinating cellular systems throughout the human body, making them ideal drug targets. Structural biology has played a key role in defining how receptors are activated and signal through G proteins and β-arrestins. The application of structure-based drug design (SBDD) is now yielding novel compounds targeting GPCRs. There is thus significant interest from both academia and the pharmaceutical industry in the structural biology of GPCRs as currently only about one quarter of human non-odorant receptors have had their structure determined. Initially, all the structures were determined by X-ray crystallography, but recent advances in electron cryo-microscopy (cryo-EM) now make GPCRs tractable targets for single-particle cryo-EM with comparable resolution to X-ray crystallography. So far this year, 78% of the 99 GPCR structures deposited in the PDB (Jan-Jul 2021) were determined by cryo-EM. Cryo-EM has also opened up new possibilities in GPCR structural biology, such as determining structures of GPCRs embedded in a lipid nanodisc and multiple GPCR conformations from a single preparation. However, X-ray crystallography still has a number of advantages, particularly in the speed of determining many structures of the same receptor bound to different ligands, an essential prerequisite for effective SBDD. We will discuss the relative merits of cryo-EM and X-ray crystallography for the structure determination of GPCRs and the future potential of both techniques.

An automated platform for structural analysis of membrane proteins through serial crystallography

Membrane proteins are central to many pathophysiological processes, yet remain very difficult to analyze structurally. Moreover, high-throughput structure-based drug discovery has not yet been exploited for membrane proteins because of lack of automation. Here, we present a facile and versatile platform for in meso membrane protein crystallization, enabling rapid atomic structure determination at both cryogenic and room temperatures. We apply this approach to human integral membrane proteins, which allowed us to identify different conformational states of intramembrane enzyme-product complexes and analyze by molecular dynamics simulations the structural dynamics of the ADIPOR2 integral membrane protein. Finally, we demonstrate an automated pipeline combining high-throughput microcrystal soaking, automated laser-based harvesting, and serial crystallography, enabling screening of small-molecule libraries with membrane protein crystals grown in meso. This approach brings needed automation to this important class of drug targets and enables high-throughput structure-based ligand discovery with membrane proteins.

Structure-Guided Design of G-Protein-Coupled Receptor Polypharmacology

Many diseases are polygenic and can only be treated efficiently with drugs that modulate multiple targets. However, rational design of compounds with multi-target profiles is rarely pursued because it is considered too difficult, in particular if the drug must enter the central nervous system. Here, a structure-based strategy to identify dual-target ligands of G-protein-coupled receptors is presented. We use this approach to design compounds that both antagonize the A_2A adenosine receptor and activate the D₂ dopamine receptor, which have excellent potential as antiparkinson drugs. Atomic resolution models of the receptors guided generation of a chemical library with compounds designed to occupy orthosteric and secondary binding pockets in both targets. Structure-based virtual screens identified ten compounds, of which three had affinity for both targets. One of these scaffolds was optimized to nanomolar dual-target activity and showed the predicted pharmacodynamic effect in a rat model of Parkinsonism.

Purinergic GPCR transmembrane residues involved in ligand recognition and dimerization

We compare the GPCR-ligand interactions and highlight important residues for recognition in purinergic receptors-from both X-ray crystallographic and cryo-EM structures. These include A1 and A2A adenosine receptors, and P2Y1 and P2Y12 receptors that respond to ADP and other nucleotides. These receptors are important drug discovery targets for immune, metabolic and nervous system disorders. In most cases, orthosteric ligands are represented, except for one allosteric P2Y1 antagonist. This review catalogs the residues and regions that engage in contacts with ligands or with other GPCR protomers in dimeric forms. Residues that are in proximity to bound ligands within purinergic GPCR families are correlated. There is extensive conservation of recognition motifs between adenosine receptors, but the P2Y1 and P2Y12 receptors are each structurally distinct in their ligand recognition. Identifying common interaction features for ligand recognition within a receptor class that has multiple structures available can aid in the drug discovery process.

Pathways and Mechanism of Caffeine Binding to Human Adenosine A2A Receptor

Caffeine (CFF) is a common antagonist to the four subtypes of adenosine G-protein-coupled receptors (GPCRs), which are critical drug targets for treating heart failure, cancer, and neurological diseases. However, the pathways and mechanism of CFF binding to the target receptors remain unclear. In this study, we have performed all-atom-enhanced sampling simulations using a robust Gaussian-accelerated molecular dynamics (GaMD) method to elucidate the binding mechanism of CFF to human adenosine A_2A receptor (A_2AAR). Multiple 500–1,000 ns GaMD simulations captured both binding and dissociation of CFF in the A_2AAR. The GaMD-predicted binding poses of CFF were highly consistent with the x-ray crystal conformations with a characteristic hydrogen bond formed between CFF and residue N6.55 in the receptor. In addition, a low-energy intermediate binding conformation was revealed for CFF at the receptor extracellular mouth between ECL2 and TM1. While the ligand-binding pathways of the A_2AAR were found similar to those of other class A GPCRs identified from previous studies, the ECL2 with high sequence divergence serves as an attractive target site for designing allosteric modulators as selective drugs of the A_2AAR.

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 A_2A 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.

Engineering of Challenging G Protein-Coupled Receptors for Structure Determination and Biophysical Studies

Membrane proteins such as G protein-coupled receptors (GPCRs) exert fundamental biological functions and are involved in a multitude of physiological responses, making these receptors ideal drug targets. Drug discovery programs targeting GPCRs have been greatly facilitated by the emergence of high-resolution structures and the resulting opportunities to identify new chemical entities through structure-based drug design. To enable the determination of high-resolution structures of GPCRs, most receptors have to be engineered to overcome intrinsic hurdles such as their poor stability and low expression levels. In recent years, multiple engineering approaches have been developed to specifically address the technical difficulties of working with GPCRs, which are now beginning to make more challenging receptors accessible to detailed studies. Importantly, successfully engineered GPCRs are not only valuable in X-ray crystallography, but further enable biophysical studies with nuclear magnetic resonance spectroscopy, surface plasmon resonance, native mass spectrometry, and fluorescence anisotropy measurements, all of which are important for the detailed mechanistic understanding, which is the prerequisite for successful drug design. Here, we summarize engineering strategies based on directed evolution to reduce workload and enable biophysical experiments of particularly challenging GPCRs.

Structural determinants of cholesterol recognition in helical integral membrane proteins

Cholesterol is an integral component of mammalian membranes. It has been shown to modulate membrane fluidity and dynamics and alter integral membrane protein function. However, understanding the molecular mechanisms of how cholesterol impacts protein function is complicated by limited and conflicting structural data. Because of the nature of the crystallization and cryo-EM structure determination, it is difficult to distinguish between specific and biologically relevant interactions and a nonspecific association. The only widely recognized search algorithm for cholesterol-integral-membrane-protein interaction sites is sequence based, i.e., searching for the so-called "Cholesterol Recognition/interaction Amino acid Consensus" motif. Although these motifs are present in numerous integral membrane proteins, there is inconclusive evidence to support their necessity or sufficiency for cholesterol binding. Here, we leverage the increasing number of experimental cholesterol-integral-membrane-protein structures to systematically analyze putative interaction sites based on their spatial arrangement and evolutionary conservation. This analysis creates three-dimensional representations of general cholesterol interaction sites that form clusters across multiple integral membrane protein classes. We also classify cholesterol-integral-membrane-protein interaction sites as either likely-specific or nonspecific. Information gleaned from our characterization will eventually enable a structure-based approach to predict and design cholesterol-integral-membrane-protein interaction sites.

The Specificity of Downstream Signaling for A1 and A2AR Does Not Depend on the C-Terminus, Despite the Importance of This Domain in Downstream Signaling Strength

Recent efforts to determine the high-resolution crystal structures for the adenosine receptors (A₁R and A_2AR) have utilized modifications to the native receptors in order to facilitate receptor crystallization and structure determination. One common modification is a truncation of the unstructured C-terminus, which has been utilized for all the adenosine receptor crystal structures obtained to date. Ligand binding for this truncated receptor has been shown to be similar to full-length receptor for A_2AR. However, the C-terminus has been identified as a location for protein-protein interactions that may be critical for the physiological function of these important drug targets. We show that variants with A_2AR C-terminal truncations lacked cAMP-linked signaling compared to the full-length receptor constructs transfected into mammalian cells (HEK-293). In addition, we show that in a humanized yeast system, the absence of the full-length C-terminus affected downstream signaling using a yeast MAPK response-based fluorescence assay, though full-length receptors showed native-like G-protein coupling. To further study the G protein coupling, we used this humanized yeast platform to explore coupling to human-yeast G-protein chimeras in a cellular context. Although the C-terminus was essential for Gα protein-associated signaling, chimeras of A₁R with a C-terminus of A_2AR coupled to the A₁R-specific Gα (i.e., Gαi1 versus Gαs). This surprising result suggests that the C-terminus is important in the signaling strength, but not specificity, of the Gα protein interaction. This result has further implications in drug discovery, both in enabling the experimental use of chimeras for ligand design, and in the cautious interpretation of structure-based drug design using truncated receptors.

A comprehensive approach to X-ray crystallography for drug discovery at a synchrotron facility - The example of Diamond Light Source

A detailed understanding of the interactions between drugs and their targets is crucial to develop the best possible therapeutic agents. Structure-based drug design relies on the availability of high-resolution structures obtained primarily through X-ray crystallography. Collecting and analysing quickly a large quantity of structural data is crucial to accelerate drug discovery pipelines. Researchers from academia and industry can access the highly automated macromolecular crystallography (MX) beamlines of Diamond Light Source, the UK national synchrotron, to rapidly collect diffraction data from large numbers of crystals. With seven beamlines dedicated to MX, Diamond offers bespoke solutions for a wide variety of user requirements. Working in synergy with state-of-the-art laboratories and other life science instruments to provide an integrated offering, the MX beamlines provide innovative and multidisciplinary approaches to the determination of structures of new pharmacological targets as well as the efficient study of protein-ligand complexes.

Probing Membrane Protein Assembly into Nanodiscs by In Situ Dynamic Light Scattering: A2A Receptor as a Case Study

Simple Summary

Located within the biological cell membranes, integral membrane proteins are responsible for a large variety of vital cellular processes. In humans, nearly a quarter of the genome codes integral membrane proteins, therefore malfunction of these proteins is associated with a variety of symptoms and diseases such as obesity, cancer and Parkinson’s disease. Clearly, knowledge of membrane proteins behaviour, in both structural and functional terms, is important not only in medicine but also in the design of better drugs with improved pharmaceutical properties. Nevertheless, much still remains unknown about these proteins, mainly because of the technical challenges associated with their production and stability in vitro once removed from their native lipidic environment. Recently, several membrane mimetic systems have been developed including nanodisc lipid particles. Nanodiscs are self-assembled lipidic structures that “trap” membrane proteins into a disc shaped phospholipid bilayer that is stabilised by a belt made of a protein know as membrane scaffold protein (MSP). Membrane proteins assembled into lipidic nanodiscs can maintain their structural and functional integrity and are compatible with most biophysical methods. Here we demonstrate the use of in situ dynamic light scattering as a high-throughput screening tool to assess the best conditions for nanodisc assembly and protein incorporation.

Abstract

Membrane proteins play a crucial role in cell physiology by participating in a variety of essential processes such as transport, signal transduction and cell communication. Hence, understanding their structure–function relationship is vital for the improvement of therapeutic treatments. Over the last decade, based on the development of detergents, amphipoles and styrene maleic-acid lipid particles (SMALPs), remarkable accomplishments have been made in the field of membrane protein structural biology. Nevertheless, there are still many drawbacks associated with protein–detergent complexes, depending on the protein in study or experimental application. Recently, newly developed membrane mimetic systems have become very popular for allowing a structural and functional characterisation of membrane proteins in vitro. The nanodisc technology is one such valuable tool, which provides a more native-like membrane environment than detergent micelles or liposomes. In addition, it is also compatible with many biophysical and biochemical methods. Here we describe the use of in situ dynamic light scattering to accurately and rapidly probe membrane proteins’ reconstitution into nanodiscs. The adenosine type 2A receptor (A_2AR) was used as a case study.

Isoprenoid-chained lipid EROCOC17+4: a new matrix for membrane protein crystallization and a crystal delivery medium in serial femtosecond crystallography

In meso crystallization of membrane proteins relies on the use of lipids capable of forming a lipidic cubic phase (LCP). However, almost all previous crystallization trials have used monoacylglycerols, with 1-(cis-9-octadecanoyl)-rac-glycerol (MO) being the most widely used lipid. We now report that EROCOC₁₇₊₄ mixed with 10% (w/w) cholesterol (Fig. 1) serves as a new matrix for crystallization and a crystal delivery medium in the serial femtosecond crystallography of Adenosine A_2A receptor (A_2AR). The structures of EROCOC₁₇₊₄-matrix grown A_2AR crystals were determined at 2.0 Å resolution by serial synchrotron rotation crystallography at a cryogenic temperature, and at 1.8 Å by LCP-serial femtosecond crystallography, using an X-ray free-electron laser at 4 and 20 °C sample temperatures, and are comparable to the structure of the MO-matrix grown A_2AR crystal (PDB ID: 4EIY). Moreover, X-ray scattering measurements indicated that the EROCOC₁₇₊₄/water system did not form the crystalline L_C phase at least down to - 20 °C, in marked contrast to the equilibrium MO/water system, which transforms into the crystalline L_C phase below about 17 °C. As the L_C phase formation within the LCP-matrix causes difficulties in protein crystallography experiments in meso, this feature of EROCOC₁₇₊₄ will expand the utility of the in meso method.

Harnessing the power of an X-ray laser for serial crystallography of membrane proteins crystallized in lipidic cubic phase

Serial femtosecond crystallography (SFX) with X-ray free-electron lasers (XFELs) has proven highly successful for structure determination of challenging membrane proteins crystallized in lipidic cubic phase; however, like most techniques, it has limitations. Here we attempt to address some of these limitations related to the use of a vacuum chamber and the need for attenuation of the XFEL beam, in order to further improve the efficiency of this method. Using an optimized SFX experimental setup in a helium atmosphere, the room-temperature structure of the adenosine A2A receptor (A2AAR) at 2.0 Å resolution is determined and compared with previous A2AAR structures determined in vacuum and/or at cryogenic temperatures. Specifically, the capability of utilizing high XFEL beam transmissions is demonstrated, in conjunction with a high dynamic range detector, to collect high-resolution SFX data while reducing crystalline material consumption and shortening the collection time required for a complete dataset. The experimental setup presented herein can be applied to future SFX applications for protein nanocrystal samples to aid in structure-based discovery efforts of therapeutic targets that are difficult to crystallize.

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.

Tumor Immunotherapy Using A2A Adenosine Receptor Antagonists

The A_2A adenosine receptor (A_2AAR) plays critical roles in human physiology and pathophysiology, which makes it an important drug target. Previous drug-discovery efforts targeting the A_2AAR have been focused on the use of A_2AAR antagonists for the treatment of Parkinson's disease. More recently, the A_2AAR has attracted additional attention for its roles in immuno-oncology, and a number of A_2AAR antagonists are currently used as lead compounds for antitumor drugs in both preclinical models and clinical trials. This review surveys recent advances in the development of A_2AAR antagonists for cancer immunotherapy. The therapeutic potential of representative A_2AAR antagonists is discussed based on both animal efficacy studies and clinical data.

Adenosine A2A Receptor Antagonists for Cancer Immunotherapy

Currently, the most promising therapeutic modality for cancer treatment is the blockade of immune checkpoint pathways, which has revolutionized cancer therapy in the past 15 years. Strategies targeting and modulating adenosine A_2A receptor (A_2AR), an emerging alternative immune checkpoint, have shown the potential to produce significant therapeutic effects. In this review, we describe the immunosuppressive activities of A_2AR and A_2BR in the tumor microenvironment (TME), followed by a summary and discussion of the structure-activity relationship (SAR) of the A_2AR (and dual A_2AR/A_2BR) antagonists that have been experimentally confirmed to exert oncoimmunological effects. This review also provides an update on the compounds under clinical evaluation and insights into the ligand binding modes of the receptor.

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation

The computational prediction of relative binding free energies is a crucial goal for drug discovery, and G protein-coupled receptors (GPCRs) are arguably the most important drug target class. However, they present increased complexity to model compared to soluble globular proteins. Despite breakthroughs, experimental X-ray crystal and cryo-EM structures are challenging to attain, meaning computational models of the receptor and ligand binding mode are sometimes necessary. This leads to uncertainty in understanding ligand-protein binding induced changes such as, water positioning and displacement, side chain positioning, hydrogen bond networks, and the overall structure of the hydration shell around the ligand and protein. In other words, the very elements that define structure activity relationships (SARs) and are crucial for accurate binding free energy calculations are typically more uncertain for GPCRs. In this work we use free energy perturbation (FEP) to predict the relative binding free energies for ligands of two different GPCRs. We pinpoint the key aspects for success such as the important role of key water molecules, amino acid ionization states, and the benefit of equilibration with specific ligands. Initial calculations following typical FEP setup and execution protocols delivered no correlation with experiment, but we show how results are improved in a logical and systematic way. This approach gave, in the best cases, a coefficient of determination (R²) compared with experiment in the range of 0.6-0.9 and mean unsigned errors compared to experiment of 0.6-0.7 kcal/mol. We anticipate that our findings will be applicable to other difficult-to-model protein ligand data sets and be of wide interest for the community to continue improving FE binding energy predictions.

3D-printed holders for in meso in situ fixed-target serial X-ray crystallography

The in meso in situ serial X-ray crystallography method was developed to ease the handling of small fragile crystals of membrane proteins and for rapid data collection on hundreds of microcrystals directly in the growth medium without the need for crystal harvesting. To facilitate mounting of these in situ samples on a goniometer at cryogenic or at room temperatures, two new 3D-printed holders have been developed. They provide for cubic and sponge phase sample stability in the X-ray beam and are compatible with sample-changing robots. The holders can accommodate a variety of window material types, as well as bespoke samples for diffraction screening and data collection at conventional macromolecular crystallography beamlines. They can be used for convenient post-crystallization treatments such as ligand and heavy-atom soaking. The design, assembly and application of the holders for in situ serial crystallography are described. Files for making the holders using a 3D printer are included as supporting information.

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.

Quantitative prediction of selectivity between the A1 and A2A adenosine receptors

The development of drugs is often hampered due to off-target interactions leading to adverse effects. Therefore, computational methods to assess the selectivity of ligands are of high interest. Currently, selectivity is often deduced from bioactivity predictions of a ligand for multiple targets (individual machine learning models). Here we show that modeling selectivity directly, by using the affinity difference between two drug targets as output value, leads to more accurate selectivity predictions. We test multiple approaches on a dataset consisting of ligands for the A₁ and A_2A adenosine receptors (among others classification, regression, and we define different selectivity classes). Finally, we present a regression model that predicts selectivity between these two drug targets by directly training on the difference in bioactivity, modeling the selectivity-window. The quality of this model was good as shown by the performances for fivefold cross-validation: ROC A₁AR-selective 0.88 ± 0.04 and ROC A_2AAR-selective 0.80 ± 0.07. To increase the accuracy of this selectivity model even further, inactive compounds were identified and removed prior to selectivity prediction by a combination of statistical models and structure-based docking. As a result, selectivity between the A₁ and A_2A adenosine receptors was predicted effectively using the selectivity-window model. The approach presented here can be readily applied to other selectivity cases.

Insights into adenosine A2A receptor activation through cooperative modulation of agonist and allosteric lipid interactions

The activation process of G protein-coupled receptors (GPCRs) has been extensively studied, both experimentally and computationally. In particular, Molecular Dynamics (MD) simulations have proven useful in exploring GPCR conformational space. The typical behaviour of class A GPCRs, when subjected to unbiased MD simulations from their crystallized inactive state, is to fluctuate between inactive and intermediate(s) conformations, even with bound agonist. Fully active conformation(s) are rarely stabilized unless a G protein is also bound. Despite several crystal structures of the adenosine A2a receptor (A2aR) having been resolved in complex with co-crystallized agonists and Gs protein, its agonist-mediated activation process is still not completely understood. In order to thoroughly examine the conformational landscape of A2aR activation, we performed unbiased microsecond-length MD simulations in quadruplicate, starting from the inactive conformation either in apo or with bound agonists: endogenous adenosine or synthetic NECA, embedded in two homogeneous phospholipid membranes: 1,2-dioleoyl-sn-glycerol-3-phosphoglycerol (DOPG) or 1,2-dioleoyl-sn-glycerol-3-phosphocholine (DOPC). In DOPC with bound adenosine or NECA, we observe transition to an intermediate receptor conformation consistent with the known adenosine-bound crystal state. In apo state in DOPG, two different intermediate conformations are obtained. One is similar to that observed with bound adenosine in DOPC, while the other is closer to the active state but not yet fully active. Exclusively, in DOPG with bound adenosine or NECA, we reproducibly identify receptor conformations with fully active features, which are able to dock Gs protein. These different receptor conformations can be attributed to the action/absence of agonist and phospholipid-mediated allosteric effects on the intracellular side of the receptor.

Towards an Optimal Sample Delivery Method for Serial Crystallography at XFEL

The advent of the X-ray free electron laser (XFEL) in the last decade created the discipline of serial crystallography but also the challenge of how crystal samples are delivered to X-ray. Early sample delivery methods demonstrated the proof-of-concept for serial crystallography and XFEL but were beset with challenges of high sample consumption, jet clogging and low data collection efficiency. The potential of XFEL and serial crystallography as the next frontier of structural solution by X-ray for small and weakly diffracting crystals and provision of ultra-fast time-resolved structural data spawned a huge amount of scientific interest and innovation. To utilize the full potential of XFEL and broaden its applicability to a larger variety of biological samples, researchers are challenged to develop better sample delivery methods. Thus, sample delivery is one of the key areas of research and development in the serial crystallography scientific community. Sample delivery currently falls into three main systems: jet-based methods, fixed-target chips, and drop-on-demand. Huge strides have since been made in reducing sample consumption and improving data collection efficiency, thus enabling the use of XFEL for many biological systems to provide high-resolution, radiation damage-free structural data as well as time-resolved dynamics studies. This review summarizes the current main strategies in sample delivery and their respective pros and cons, as well as some future direction.

Toward G protein-coupled receptor structure-based drug design using X-ray lasers

Rational structure-based drug design (SBDD) relies on the availability of a large number of co-crystal structures to map the ligand-binding pocket of the target protein and use this information for lead-compound optimization via an iterative process. While SBDD has proven successful for many drug-discovery projects, its application to G protein-coupled receptors (GPCRs) has been limited owing to extreme difficulties with their crystallization. Here, a method is presented for the rapid determination of multiple co-crystal structures for a target GPCR in complex with various ligands, taking advantage of the serial femtosecond crystallography approach, which obviates the need for large crystals and requires only submilligram quantities of purified protein. The method was applied to the human β2-adrenergic receptor, resulting in eight room-temperature co-crystal structures with six different ligands, including previously unreported structures with carvedilol and propranolol. The generality of the proposed method was tested with three other receptors. This approach has the potential to enable SBDD for GPCRs and other difficult-to-crystallize membrane proteins.

Cryo-Electron Microscopy: Moving Beyond X-Ray Crystal Structures for Drug Receptors and Drug Development

Electron cryo-microscopy (cryo-EM) has revolutionized structure determination of membrane proteins and holds great potential for structure-based drug discovery. Here we discuss the potential of cryo-EM in the rational design of therapeutics for membrane proteins compared to X-ray crystallography. We also detail recent progress in the field of drug receptors, focusing on cryo-EM of two protein families with established therapeutic value, the γ-aminobutyric acid A receptors (GABA_ARs) and G protein-coupled receptors (GPCRs). GABA_ARs are pentameric ion channels, and cryo-EM structures of physiological heteromeric receptors in a lipid environment have uncovered the molecular basis of receptor modulation by drugs such as diazepam. The structures of ten GPCR-G protein complexes from three different classes of GPCRs have now been determined by cryo-EM. These structures give detailed insights into molecular interactions with drugs, GPCR-G protein selectivity, how accessory membrane proteins alter receptor-ligand pharmacology, and the mechanism by which HIV uses GPCRs to enter host cells.