Biophysical mapping of the adenosine A2A receptor

A chromatographic method for determining the interaction between a drug and two target proteins by fabricating a dual-heterogeneous surface

Characterization of the drug-target interactions is pivotal throughout the whole procedure of drug development. Most of the current assays, particularly, chromatographic methods lack the capacity to reveal drug adsorption on the muti-target surface. To this end, we derived a reliable and workable mathematical equation for revealing drug bindings to dual targets on the heterogeneous surface starting from the mass balance equation. The derivatization relied on the correlation of drug injection amounts with their retention factors. Experimental validation was performed by determining the binding parameters of three canonical drugs on a heterogeneous surface, which was fabricated by fusing angiotensin receptor type I and type II receptors (AT1R and AT2R) at the terminuses of circularly permuted HaloTag (cpHaloTag) and immobilizing the whole fusion protein onto 6-bromohexanoic acid modified silica gel. We proved that immobilized AT1R-cpHalo-AT2R maintained the original ligand- and antibody-binding activities of the two receptors in three weeks. The association constants of valsartan, candesartan, and telmisartan to AT1R were (6.26±0.14) × 105, (9.66±0.71) × 105, and (3.17±0.03) × 105 L/mol. In the same column, their association constants to AT2R were (1.25±0.04) × 104, (2.30±0.08) × 104, and (8.51±0.06) × 103 L/mol. The patterns of the association constants to AT1R/AT2R (candesartan>valsartan>telmisartan) were in good line with the data by performing nonlinear chromatography on control columns containing immobilized AT1R or AT2R alone. This provided proof of the fact that the derivatization allowed the determination of drug bindings on the heterogeneous surface with the utilization of a single series of injections and linear regression. We reasoned that is simple enough to model the bindings of drug adsorption on commercially available adsorbents in fundamental or industrial fields, thus having the potential to become a universal method for analyzing the bindings of a drug to the heterogeneous surface containing multiple targets.

Direct Binding Methods to Measure Receptor-Ligand Interactions

G-protein-coupled receptors (GPCRs) contribute to numerous physiological processes via complex network mechanisms. While indirect signaling assays (Ca2+ mobilization, cAMP production, and GTPγS binding) have been useful in identifying and characterizing downstream signaling mechanisms of GPCRs, these methods lack measurements of direct binding affinities, kinetics, binding specificity, and selectivity that are important parameters in GPCR drug discovery. In comparison to existing direct methods that use radio- or fluorescent labels, label-free techniques can closely emulate the native interactions around binding partners. Surface plasmon resonance (SPR) is a label-free technique that utilizes the refractive index (RI) property and is applied widely in quantitative GPCR-ligand binding kinetics measurement including small molecules screening. However, purified GPCRs are further embedded in a synthetic lipid environment which is immobilized through different tags to the SPR sensor surface, resulting in a non-native environment. Here, we introduced a methodology that also uses the RI property to measure binding interactions in a label-free, immobilization-free arrangement. The free-solution technique is successfully applied in quantifying the interaction of bioactive lipids to cognate lipid GPCRs, which is not purified but rather present in near-native conditions, i.e., in milieu of other cytoplasmic lipids and proteins. To further consider the wide applicability of these free-solution approaches in biomolecular interaction research, additional applications on a variety of receptor-ligand pairs are imperative.

Anti-nanodisc antibodies specifically capture nanodiscs and facilitate molecular interaction kinetics studies for membrane protein

Nanodisc technology has dramatically advanced the analysis of molecular interactions for membrane proteins. A nanodisc is designed as a vehicle for membrane proteins that provide a native-like phospholipid environment and better thermostability in a detergent-free buffer. This enables the determination of the thermodynamic and kinetic parameters of small molecule binding by surface plasmon resonance. In this study, we generated a nanodisc specific anti-MSP (membrane scaffold protein) monoclonal antibody biND5 for molecular interaction analysis of nanodiscs. The antibody, biND5 bound to various types of nanodiscs with sub-nanomolar to nanomolar affinity. Epitope mapping analysis revealed specific recognition of 8 amino acid residues in the exposed helix-4 structure of MSP. Further, we performed kinetics binding analysis between adenosine A_2a receptor reconstituted nanodiscs and small molecule antagonist ZM241385 using biND5 immobilized sensor chips. These results show that biND5 facilitates the molecular interaction kinetics analysis of membrane proteins substituted in nanodiscs.

Molecular docking/dynamics simulations, MEP analysis, bioisosteric replacement and ADME/T prediction for identification of dual targets inhibitors of Parkinson's disease with novel scaffold

Monoamine oxidase B and Adenosine A2A receptors are used as key targets for Parkinson's disease. Recently, hMAO-B and hA2AR Dual-targets inhibitory potential of a novel series of Phenylxanthine derivatives has been established in experimental findings. Hence, the current study examines the interactions between 38 compounds of this series with hMAO-B and hA2AR targets using different molecular modeling techniques to investigate the binding mode and stability of the formed complexes. A molecular docking study revealed that the compounds L24 ((E)-3-(3-Chlorophenyl)-N-(4-(1,3-dimethyl-2,6-dioxo-2,3,6,7-tetrahydro-1H-purin-8-yl) phenyl) acrylamide and L32 ((E)-3-(3-Chlorophenyl)-N-(3-(1,3-dimethyl-2,6-dioxo-2,3,6,7-tetrahydro-1H-purin-8-yl)phenyl)acrylamide) had a high affinity (S-score: -10.160 and -7.344 kcal/mol) with the pocket of hMAO-B and hA2AR targets respectively, and the stability of the studied complexes was confirmed during MD simulations. Also, the MEP maps of compounds 24 and 32 were used to identify the nucleophilic and electrophilic attack regions. Moreover, the bioisosteric replacement approach was successfully applied to design two new analogs of each compound with similar biological activities and low energy scores. Furthermore, ADME-T and Drug-likeness results revealed the promising pharmacokinetic properties and oral bioavailability of these compounds. Thus, compounds L24, L32, and their analogs can undergo further analysis and optimization in order to design new lead compounds with higher efficacy toward Parkinson's disease.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-023-00139-3.

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.

Cancer-Associated Mutations of the Adenosine A2A Receptor Have Diverse Influences on Ligand Binding and Receptor Functions

The adenosine A_2A receptor (A_2AAR) is a class A G-protein-coupled receptor (GPCR). It is an immune checkpoint in the tumor micro-environment and has become an emerging target for cancer treatment. In this study, we aimed to explore the effects of cancer-patient-derived A_2AAR mutations on ligand binding and receptor functions. The wild-type A_2AAR and 15 mutants identified by Genomic Data Commons (GDC) in human cancers were expressed in HEK293T cells. Firstly, we found that the binding affinity for agonist NECA was decreased in six mutants but increased for the V275A mutant. Mutations A165V and A265V decreased the binding affinity for antagonist ZM241385. Secondly, we found that the potency of NECA (EC₅₀) in an impedance-based cell-morphology assay was mostly correlated with the binding affinity for the different mutants. Moreover, S132L and H278N were found to shift the A_2AAR towards the inactive state. Importantly, we found that ZM241385 could not inhibit the activation of V275A and P285L stimulated by NECA. Taken together, the cancer-associated mutations of A_2AAR modulated ligand binding and receptor functions. This study provides fundamental insights into the structure–activity relationship of the A_2AAR and provides insights for A_2AAR-related personalized treatment in cancer.

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-Based Virtual Screening for Ligands of G Protein-Coupled Receptors: What Can Molecular Docking Do for You?

G protein-coupled receptors (GPCRs) constitute the largest family of membrane proteins in the human genome and are important therapeutic targets. During the last decade, the number of atomic-resolution structures of GPCRs has increased rapidly, providing insights into drug binding at the molecular level. These breakthroughs have created excitement regarding the potential of using structural information in ligand design and initiated a new era of rational drug discovery for GPCRs. The molecular docking method is now widely applied to model the three-dimensional structures of GPCR-ligand complexes and screen for chemical probes in large compound libraries. In this review article, we first summarize the current structural coverage of the GPCR superfamily and the understanding of receptor-ligand interactions at atomic resolution. We then present the general workflow of structure-based virtual screening and strategies to discover GPCR ligands in chemical libraries. We assess the state of the art of this research field by summarizing prospective applications of virtual screening based on experimental structures. Strategies to identify compounds with specific efficacy and selectivity profiles are discussed, illustrating the opportunities and limitations of the molecular docking method. Our overview shows that structure-based virtual screening can discover novel leads and will be essential in pursuing the next generation of GPCR drugs. SIGNIFICANCE STATEMENT: Extraordinary advances in the structural biology of G protein-coupled receptors have revealed the molecular details of ligand recognition by this large family of therapeutic targets, providing novel avenues for rational drug design. Structure-based docking is an efficient computational approach to identify novel chemical probes from large compound libraries, which has the potential to accelerate the development of drug candidates.

Nonantimicrobial Actions of Macrolides: Overview and Perspectives for Future Development

Macrolides are among the most widely prescribed broad spectrum antibacterials, particularly for respiratory infections. It is now recognized that these drugs, in particular azithromycin, also exert time-dependent immunomodulatory actions that contribute to their therapeutic benefit in both infectious and other chronic inflammatory diseases. Their increased chronic use in airway inflammation and, more recently, of azithromycin in COVID-19, however, has led to a rise in bacterial resistance. An additional crucial aspect of chronic airway inflammation, such as chronic obstructive pulmonary disease, as well as other inflammatory disorders, is the loss of epithelial barrier protection against pathogens and pollutants. In recent years, azithromycin has been shown with time to enhance the barrier properties of airway epithelial cells, an action that makes an important contribution to its therapeutic efficacy. In this article, we review the background and evidence for various immunomodulatory and time-dependent actions of macrolides on inflammatory processes and on the epithelium and highlight novel nonantibacterial macrolides that are being studied for immunomodulatory and barrier-strengthening properties to circumvent the risk of bacterial resistance that occurs with macrolide antibacterials. We also briefly review the clinical effects of macrolides in respiratory and other inflammatory diseases associated with epithelial injury and propose that the beneficial epithelial effects of nonantibacterial azithromycin derivatives in chronic inflammation, even given prophylactically, are likely to gain increasing attention in the future. SIGNIFICANCE STATEMENT: Based on its immunomodulatory properties and ability to enhance the protective role of the lung epithelium against pathogens, azithromycin has proven superior to other macrolides in treating chronic respiratory inflammation. A nonantibiotic azithromycin derivative is likely to offer prophylactic benefits against inflammation and epithelial damage of differing causes while preserving the use of macrolides as antibiotics.

Binding Characterization of Agonists and Antagonists by MCCS: A Case Study from Adenosine A2A Receptor

Characterizing the structural basis of ligand recognition of adenosine A_2A receptor (AA_2AR) will facilitate its rational design and development of small molecules with high affinity and selectivity, as well as optimal therapeutic effects for pain, cancers, drug abuse disorders, etc. In the present work, we applied our reported algorithm, molecular complex characterizing system (MCCS), to characterize the binding features of AA_2AR based on its reported 3D structures of protein-ligand complexes. First, we compared the binding score to the reported experimental binding affinities of each compound. Then, we calculated an output example of residue energy contribution using MCCS and compared the results with data obtained from MM/GBSA. The consistency in results indicated that MCCS is a powerful, fast, and accurate method. Sequentially, using a receptor-ligand data set of 57 crystallized structures of AA_2ARs, we characterized the binding features of the binding pockets in AA_2AR, summarized the key residues that distinguish antagonist from agonist, produced heatmaps of residue energy contribution for clustering various statuses of AA_2ARs, explored the selectivity between AA_2AR and AA_1AR, etc. All the information provided new insights into the protein features of AA_2AR and will facilitate its rational drug design.

Imine bond transformation of a dynamic Sm(III) macrocycle-based chemosensor: The indirect approach for detecting cyanuric chloride

Herein, we report our strategy to develop the efficient chemosensor and real-time monitoring technique for cyanuric chloride (TCT) detection. A luminescent macrocyclic mononuclear Sm(III) complex Sm-2_k bearing with two dynamic imine bonds has been constructed via the template synthesis between dialdehyde H₂Q_k and matched diamine 1,2-bis(2-aminoethoxy)ethane. Sensing experiments reveal that complex Sm-2_k exhibits the turn-off fluorescent and colorimetric response for TCT in CH₃OH. It is especially encouraging that this optical sensing process is not only rapid within 60 s but also high-efficient in the presence of TCT analogues as well as sensitive with the low limit of detection (LOD, 1.74 μM) and wide linear sensing range. Mechanism studies demonstrate that TCT sensing is mainly based on the imine bond transformation of probe Sm-2_k, which is due to the increased acidity induced by TCT. Meanwhile, a smartphone-based analytical method was developed to make complex Sm-2_k accessible for the real-time TCT detection by RGB value outputs. It is believed that this work can shed some constructive lights on design of chemosensors and convenient detection technique for highly reactive analytes.

Free Energy Calculations for Protein-Ligand Binding Prediction

Computational prediction of protein-ligand binding involves initial determination of the binding mode and subsequent evaluation of the strength of the protein-ligand interactions, which directly correlates with ligand binding affinities. As a consequence of increasing computer power, rigorous approaches to calculate protein-ligand binding affinities, such as free energy perturbation (FEP) methods, are becoming an essential part of the toolbox of computer-aided drug design. In this chapter, we provide a general overview of these methods and introduce the QFEP modules, which are open-source API workflows based on our molecular dynamics (MD) package Q. The module QligFEP allows estimation of relative binding affinities along ligand series, while QresFEP is a module to estimate binding affinity shifts caused by single-point mutations of the protein. We herein provide guidelines for the use of each of these modules based on data extracted from ligand-design projects. While these modules are stand-alone, the combined use of the two workflows in a drug-design project yields complementary perspectives of the ligand binding problem, providing two sides of the same coin. The selected case studies illustrate how to use QFEP to approach the two key questions associated with ligand binding prediction: identifying the most favorable binding mode from different alternatives and establishing structure-affinity relationships that allow the rational optimization of hit compounds.

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.

Ligand Pose Predictions for Human G Protein-Coupled Receptors: Insights from the Amber-Based Hybrid Molecular Mechanics/Coarse-Grained Approach

Human G protein-coupled receptors (hGPCRs) are the most frequent targets of Food and Drug Administration (FDA)-approved drugs. Structural bioinformatics, along with molecular simulation, can support structure-based drug design targeting hGPCRs. In this context, several years ago, we developed a hybrid molecular mechanics (MM)/coarse-grained (CG) approach to predict ligand poses in low-resolution hGPCR models. The approach was based on the GROMOS96 43A1 and PRODRG united-atom force fields for the MM part. Here, we present a new MM/CG implementation using, instead, the Amber 14SB and GAFF all-atom potentials for proteins and ligands, respectively. The new implementation outperforms the previous one, as shown by a variety of applications on models of hGPCR/ligand complexes at different resolutions, and it is also more user-friendly. Thus, it emerges as a useful tool to predict poses in low-resolution models and provides insights into ligand binding similarly to all-atom molecular dynamics, albeit at a lower computational cost.

Perspective: Implications of Ligand-Receptor Binding Kinetics for Therapeutic Targeting of G Protein-Coupled Receptors

The concept of ligand-receptor binding kinetics has been broadly applied in drug development pipelines focusing on G protein-coupled receptors (GPCRs). The ligand residence time (RT) for a receptor describes how long a ligand-receptor complex exists, and is defined as the reciprocal of the dissociation rate constant (k _off). RT has turned out to be a valuable parameter for GPCR researchers focusing on drug development as a good predictor of in vivo efficacy. The positive correlation between RT and in vivo efficacy has been established for several drugs targeting class A GPCRs (e.g., the neurokinin-1 receptor (NK₁R), the β₂ adrenergic receptor (β₂AR), and the muscarinic 3 receptor (M₃R)) and for drugs targeting class B1 (e.g., the glucagon-like peptide 1 receptor (GLP-1R)). Recently, the association rate constant (k _on) has gained similar attention as another parameter affecting in vivo efficacy. In the current perspective, we address the importance of studying ligand-receptor binding kinetics for therapeutic targeting of GPCRs, with an emphasis on how binding kinetics can be altered by subtle molecular changes in the ligands and/or the receptors and how such changes affect treatment outcome. Moreover, we speculate on the impact of binding kinetic parameters for functional selectivity and sustained receptor signaling from endosomal compartments; phenomena that have gained increasing interest in attempts to improve therapeutic targeting of GPCRs.

Drug-Target Association Kinetics in Drug Discovery

The important role of ligand-receptor binding kinetics in drug design and discovery is increasingly recognized by the drug research community. Over the past decade, accumulating evidence has shown that optimizing the ligand's dissociation rate constant can lead to desirable duration of in vivo target occupancy and, hence, improved pharmacodynamic properties. However, the association rate constant as a pharmacological principle remains less investigated, whereas it can play an equally important role in the selection of drug candidates. This review provides a compilation and discussion of otherwise scarce and dispersed information on this topic, bringing to light the importance of drug-target association in kinetics-directed drug design and discovery.

Exploring Molecular-Biomembrane Interactions with Surface Plasmon Resonance and Dual Polarization Interferometry Technology: Expanding the Spotlight onto Biomembrane Structure

The molecular analysis of biomolecular-membrane interactions is central to understanding most cellular systems but has emerged as a complex technical challenge given the complexities of membrane structure and composition across all living cells. We present a review of the application of surface plasmon resonance and dual polarization interferometry-based biosensors to the study of biomembrane-based systems using both planar mono- or bilayers or liposomes. We first describe the optical principals and instrumentation of surface plasmon resonance, including both linear and extraordinary transmission modes and dual polarization interferometry. We then describe the wide range of model membrane systems that have been developed for deposition on the chips surfaces that include planar, polymer cushioned, tethered bilayers, and liposomes. This is followed by a description of the different chemical immobilization or physisorption techniques. The application of this broad range of engineered membrane surfaces to biomolecular-membrane interactions is then overviewed and how the information obtained using these techniques enhance our molecular understanding of membrane-mediated peptide and protein function. We first discuss experiments where SPR alone has been used to characterize membrane binding and describe how these studies yielded novel insight into the molecular events associated with membrane interactions and how they provided a significant impetus to more recent studies that focus on coincident membrane structure changes during binding of peptides and proteins. We then discuss the emerging limitations of not monitoring the effects on membrane structure and how SPR data can be combined with DPI to provide significant new information on how a membrane responds to the binding of peptides and proteins.

Optimization and 13CH3 methionine labeling of a signaling competent neurotensin receptor 1 variant for NMR studies

Neurotensin is a 13-residue peptide that acts as a neuromodulator of classical neurotransmitters such as dopamine and glutamate in the mammalian central nervous system, mainly by activating the G protein-coupled receptor (GPCR), neurotensin receptor 1 (NTS₁). Agonist binding to GPCRs shifts the conformational equilibrium of the transmembrane helices towards distinct, thermodynamically favorable conformations that favor effector protein interactions and promotes cell signaling. The introduction of site specific labels for NMR spectroscopy has proven useful for investigating this dynamic process, but the low expression levels and poor stability of GPCRs is a hindrance to solution NMR experiments. Several thermostabilized mutants of NTS₁ have been engineered to circumvent this, with the crystal structures of four of these published. The conformational dynamics of NTS₁ however, has not been thoroughly investigated with NMR. It is generally accepted that stabilized GPCRs exhibit attenuated signaling, thus we thoroughly characterized the signaling characteristics of several thermostabilized NTS₁ variants to identify an optimal variant for protein NMR studies. A variant termed enNTS₁ exhibited the best combination of signaling capability and stability upon solubilization with detergents. enNTS₁ was subsequently labeled with ¹³CH₃-methionine in E. coli and purified to homogeneity in the absence of bound ligands. Using solution NMR spectroscopy we observed several well dispersed ¹³CH₃-methionine resonances, many of which exhibited chemical shift changes upon the addition of the high affinity agonist peptide, NT8-13. Thus, enNTS₁ represents a novel tool for investigating ligand induced conformational changes in NTS₁ to gain insights into the molecular mechanisms underlying neurotensin signaling.

Elucidation of the Binding Mode of the Carboxyterminal Region of Peptide YY to the Human Y2 Receptor

Understanding the agonist-receptor interactions in the neuropeptide Y (NPY)/peptide YY (PYY) signaling system is fundamental for the design of novel modulators of appetite regulation. We report here the results of a multidisciplinary approach to elucidate the binding mode of the native peptide agonist PYY to the human Y₂ receptor, based on computational modeling, peptide chemistry and in vitro pharmacological analyses. The preserved binding orientation proposed for full-length PYY and five analogs, truncated at the amino terminus, explains our pharmacological results where truncations of the N-terminal proline helix showed little effect on peptide affinity. This was followed by receptor mutagenesis to investigate the roles of several receptor positions suggested by the modeling. As a complement, PYY-(3-36) analogs were synthesized with modifications at different positions in the common PYY/NPY C-terminal fragment (³²TRQRY³⁶-amide). The results were assessed and interpreted by molecular dynamics and Free Energy Perturbation (FEP) simulations of selected mutants, providing a detailed map of the interactions of the PYY/NPY C-terminal fragment with the transmembrane cavity of the Y₂ receptor. The amidated C-terminus would be stabilized by polar interactions with Gln288^6.55 and Tyr219^5.39, while Gln130^3.32 contributes to interactions with Q³⁴ in the peptide and T³² is close to the tip of TM7 in the receptor. This leaves the core, α-helix of the peptide exposed to make potential interactions with the extracellular loops. This model agrees with most experimental data available for the Y₂ system and can be used as a basis for optimization of Y₂ receptor agonists.

Characterization of Ligand Binding to GPCRs Through Computational Methods

The recent increase in available G protein-coupled receptor structures now contributes decisively to the structure-based ligand design. In this context, computational approaches in combination with medicinal chemistry and pharmacology are extremely helpful. Here, we provide an update on our structure-based computational protocols, used to answer key questions related to GPCR-ligand binding. All combined, these techniques can shed light on ligand binding modes, determine the molecular basis of conformational selection, for agonists and antagonists, as well as of subtype selectivity. To illustrate each of these questions, we will consider examples from existing projects on three families of class A (rhodopsin-like) GPCRs: one small-molecule (nucleotide-like) family, i.e., the adenosine receptors, and two peptide-binding receptors: neuropeptide-Y and angiotensin II receptors. The successful application of the same computational protocols to investigate this diverse group of receptor families gives an idea of the general applicability of our methodology in the characterization of GPCR-ligand binding.

Analysis of Protein Interactions by Surface Plasmon Resonance

Surface plasmon resonance is an optical technique that is utilized for detecting molecular interactions, such as interactions that occur between proteins or other classes of molecules. Binding of a mobile molecule (analyte) to a molecule immobilized on a thin metal film (ligand) changes the refractive index of the film. The angle of extinction of light that is completely reflected after polarized light impinges upon the film, is altered and monitored as a change in detector position for a dip in reflected intensity (the surface plasmon resonance phenomenon). Because the method strictly detects mass, there is no need to label the interacting components, thus eliminating possible changes of their molecular properties. In this chapter, we review essential SPR methodology and present applications to basic science and human disease.

Discovery of Novel and Selective Adenosine A2A Receptor Antagonists for Treating Parkinson's Disease through Comparative Structure-Based Virtual Screening

Among non-dopaminergic strategies for combating Parkinson's disease (PD), antagonism of the A_2A adenosine receptor (AR) has emerged to show great potential. In this study, on the basis of two crystal structures of the A_2A AR with the best capability to distinguish known antagonists from decoys, docking-based virtual screening (VS) was conducted to identify novel A_2A AR antagonists. A total of 63 structurally diverse compounds identified by VS were submitted to experimental testing, and 11 of them exhibited substantial activity against the A_2A AR (K_i < 10 μM), including two compounds with K_i below 1 μM (compound 43, 0.42 μM; compound 51, 0.27 μM) and good A_2A/A₁ selectivity (fold < 0.1). Compounds 43 and 51 demonstrated antagonistic activity according to the results of cAMP measurements (cAMP IC₅₀ = 1.67 and 1.80 μM, respectively) and showed good efficacy in the haloperidol-induced catalepsy (HIC) rat model for PD at doses of up to 30 mg/kg. Further lead optimization based on a substructure searching strategy led to the discovery of compound 84 as an excellent A_2A AR antagonist (A_2A K_i = 54 nM, A_2A/A₁ fold < 0.1, cAMP IC₅₀ = 0.3 μM) that exhibited significant improvement in anti-PD efficacy in the HIC rat model.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Design of Next-Generation G Protein-Coupled Receptor Drugs: Linking Novel Pharmacology and In Vivo Animal Models

Despite the fact that G protein-coupled receptors (GPCRs) are the most successful drug targets in history, this supergene family of cell surface receptors has yet to be fully exploited as targets in the treatment of human disease. Here, we present optimism that this may change in the future by reviewing the substantial progress made in the understanding of GPCR molecular pharmacology that has generated an extensive toolbox of ligand types that include orthosteric, allosteric, and bitopic ligands, many of which show signaling bias. We discuss how combining these advances with recently described transgenic, chemical genetic, and optogenetic animal models will provide the framework to allow for the rational design of next-generation GPCR drugs that possess increased therapeutic efficacy and decreased adverse/toxic responses.

A2AR Binding Kinetics in the Ligand Depletion Regime

Ligand binding plays a fundamental role in stimulating the downstream signaling of membrane receptors. Here, ligand-binding kinetics of the full-length human adenosine A_2A receptor (A_2AR) reconstituted in detergent micelles were measured using a fluorescently labeled ligand via fluorescence anisotropy. Importantly, to optimize the signal-to-noise ratio, these experiments were conducted in the ligand depletion regime. In the ligand depletion regime, the assumptions used to determine analytical solutions for one-site binding models for either one or two ligands in competition are no longer valid. We therefore implemented a numerical solution approach to analyze kinetic binding data as experimental conditions approach the ligand depletion regime. By comparing the results from the numerical and the analytical solutions, we highlight the ligand-receptor ratios at which the analytical solution begins to lose predictive accuracy. Using the numerical solution approach, we determined the kinetic rate constants of the fluorescent ligand, FITC-APEC, and those for three unlabeled ligands using competitive association experiments. The association and dissociation rate constants of the unlabeled ligands determined from the competitive association experiments were then independently validated using competitive dissociation data. Based on this study, a numerical solution is recommended to determine kinetic ligand-binding parameters for experiments conducted in the ligand-depletion regime.

In search of novel ligands using a structure-based approach: a case study on the adenosine A2A receptor

In this work, we present a case study to explore the challenges associated with finding novel molecules for a receptor that has been studied in depth and has a wealth of chemical information available. Specifically, we apply a previously described protocol that incorporates explicit water molecules in the ligand binding site to prospectively screen over 2.5 million drug-like and lead-like compounds from the commercially available eMolecules database in search of novel binders to the adenosine A_2A receptor (A_2AAR). A total of seventy-one compounds were selected for purchase and biochemical assaying based on high ligand efficiency and high novelty (Tanimoto coefficient ≤0.25 to any A_2AAR tested compound). These molecules were then tested for their affinity to the adenosine A_2A receptor in a radioligand binding assay. We identified two hits that fulfilled the criterion of ~50 % radioligand displacement at a concentration of 10 μM. Next we selected an additional eight novel molecules that were predicted to make a bidentate interaction with Asn253^6.55, a key interacting residue in the binding pocket of the A_2AAR. None of these eight molecules were found to be active. Based on these results we discuss the advantages of structure-based methods and the challenges associated with finding chemically novel molecules for well-explored targets.

A mutagenesis and screening strategy to generate optimally thermostabilized membrane proteins for structural studies

The thermostability of an integral membrane protein (MP) in detergent solution is a key parameter that dictates the likelihood of obtaining well-diffracting crystals that are suitable for structure determination. However, many mammalian MPs are too unstable for crystallization. We developed a thermostabilization strategy based on systematic mutagenesis coupled to a radioligand-binding thermostability assay that can be applied to receptors, ion channels and transporters. It takes ∼6-12 months to thermostabilize a G-protein-coupled receptor (GPCR) containing 300 amino acid (aa) residues. The resulting thermostabilized MPs are more easily crystallized and result in high-quality structures. This methodology has facilitated structure-based drug design applied to GPCRs because it is possible to determine multiple structures of the thermostabilized receptors bound to low-affinity ligands. Protocols and advice are given on how to develop thermostability assays for MPs and how to combine mutations to make an optimally stable mutant suitable for structural studies. The steps in the procedure include the generation of ∼300 site-directed mutants by Ala/Leu scanning mutagenesis, the expression of each mutant in mammalian cells by transient transfection and the identification of thermostable mutants using a thermostability assay that is based on binding of an (125)I-labeled radioligand to the unpurified, detergent-solubilized MP. Individual thermostabilizing point mutations are then combined to make an optimally stable MP that is suitable for structural biology and other biophysical studies.

Impact, determination and prediction of drug-receptor residence times for GPCRs

The residence time of a ligand on a GPCR of interest has become an optimization parameter in many examples in drug design. Long residence times can counterbalance unfavorable pharmacokinetic parameters, contributing to compound safety, and short residence times can be a tool to avoid target related side effects. Unlike the prediction and interpretation of the structure-activity relationship (SAR) of a ligand class on a receptor, the understanding and prediction of the structure-kinetics relationship (SKR) is much more demanding. Experimental and computational approaches are described, which serve to either rationalize SKR or to predict the kinetic parameters such as on-rates and off-rates.

The use of conformationally thermostabilised GPCRs in drug discovery: application to fragment, structure and biophysical techniques

Recent developments in receptor stabilisation have facilitated major advances in G protein-coupled receptor (GPCR) research, notably structural biology, over the past eight years. Here we review the application of fragment, structure and biophysical techniques using stabilised GPCRs (StaR proteins), and their impact in the drug discovery process. These techniques have, most recently, been utilised in the discovery of the non-alkyne mGlu₅ negative allosteric modulator HTL14242, in addition to the dual orexin receptor antagonist HTL6641, with differentiated residence time kinetics.

Structurally Enabled Discovery of Adenosine A2A Receptor Antagonists

Over the past decade there has been a revolution in the field of G protein-coupled receptor (GPCR) structural biology. Many years of innovative research from different areas have come together to fuel this significant change in the fortunes of this field, which for many years was characterized by the paucity of high-resolution structures. The determination to succeed has been in part due to the recognized importance of these proteins as drug targets, and although the pharmaceutical industry has been focusing on these receptors, it can be justifiably argued and demonstrated that many of the approved and commercially successful GPCR drugs can be significantly improved to increase efficacy and/or reduce undesired side effects. In addition, many validated targets in this class remain to be drugged. It is widely recognized that application of structure-based drug design approaches can help medicinal chemists a long way toward discovering better drugs. The achievement of structural biologists in providing high-resolution insight is beginning to transform drug discovery efforts, and there are a number of GPCR drugs that have been discovered by use of structural information that are in clinical development. This review aims to highlight the key developments that have brought success to GPCR structure resolution efforts and exemplify the practical application of structural information for the discovery of adenosine A_2A receptor antagonists that have potential to treat multiple conditions.

New paradigms in GPCR drug discovery

G protein-coupled receptors (GPCRs) remain a major domain of pharmaceutical discovery. The identification of GPCR lead compounds and their optimization are now structure-based, thanks to advances in X-ray crystallography, molecular modeling, protein engineering and biophysical techniques. In silico screening provides useful hit molecules. New pharmacological approaches to tuning the pleotropic action of GPCRs include: allosteric modulators, biased ligands, GPCR heterodimer-targeted compounds, manipulation of polypharmacology, receptor antibodies and tailoring of drug molecules to fit GPCR pharmacogenomics. Measurements of kinetics and drug efficacy are factors influencing clinical success. With the exception of inhibitors of GPCR kinases, targeting of intracellular GPCR signaling or receptor cycling for therapeutic purposes remains a futuristic concept. New assay approaches are more efficient and multidimensional: cell-based, label-free, fluorescence-based assays, and biosensors. Tailoring GPCR drugs to a patient's genetic background is now being considered. Chemoinformatic tools can predict ADME-tox properties. New imaging technology visualizes drug action in vivo. Thus, there is reason to be optimistic that new technology for GPCR ligand discovery will help reverse the current narrowing of the pharmaceutical pipeline.

Analytical methods for kinetic studies of biological interactions: A review

The rates at which biological interactions occur can provide important information concerning the mechanism and behavior of these processes in living systems. This review discusses several analytical methods that can be used to examine the kinetics of biological interactions. These techniques include common or traditional methods such as stopped-flow analysis and surface plasmon resonance spectroscopy, as well as alternative methods based on affinity chromatography and capillary electrophoresis. The general principles and theory behind these approaches are examined, and it is shown how each technique can be utilized to provide information on the kinetics of biological interactions. Examples of applications are also given for each method. In addition, a discussion is provided on the relative advantages or potential limitations of each technique regarding its use in kinetic studies.

Biosensor-based affinities and binding kinetics of small molecule antagonists to the adenosine A(2A) receptor reconstituted in HDL like particles

The options for investigating solubilised G protein-coupled receptors (GPCRs) by biophysical techniques have long been hampered by their instability. A thermostabilised adenosine A2A receptor expressed in insect cells, purified in detergent and reconstituted into high-density lipoprotein (HDL) particles was immobilised onto a Surface Plasmon Resonance sensor chip. This allowed measurement of affinities and kinetics for A2A antagonists with affinities ranging from 50 pM to almost 2 μM. Compared with other formats, reproduction of affinities, and dissociation and association rate constants are good, reasonable and poor respectively, indicating stabilised receptors in HDL particles are useful for investigating specific aspects of GPCR-ligand interactions.

Stabilization of G protein-coupled receptors by point mutations

G protein-coupled receptors (GPCRs) are flexible integral membrane proteins involved in transmembrane signaling. Their involvement in many physiological processes makes them interesting targets for drug development. Determination of the structure of these receptors will help to design more specific drugs, however, their structural characterization has so far been hampered by the low expression and their inherent instability in detergents which made protein engineering indispensable for structural and biophysical characterization. Several approaches to stabilize the receptors in a particular conformation have led to breakthroughs in GPCR structure determination. These include truncations of the flexible regions, stabilization by antibodies and nanobodies, fusion partners, high affinity and covalently bound ligands as well as conformational stabilization by mutagenesis. In this review we focus on stabilization of GPCRs by insertion of point mutations, which lead to increased conformational and thermal stability as well as improved expression levels. We summarize existing mutagenesis strategies with different coverage of GPCR sequence space and depth of information, design and transferability of mutations and the molecular basis for stabilization. We also discuss whether mutations alter the structure and pharmacological properties of GPCRs.

Real-time monitoring of binding events on a thermostabilized human A2A receptor embedded in a lipid bilayer by surface plasmon resonance

Membrane proteins (MPs) are prevalent drug discovery targets involved in many cell processes. Despite their high potential as drug targets, the study of MPs has been hindered by limitations in expression, purification and stabilization in order to acquire thermodynamic and kinetic parameters of small molecules binding. These bottlenecks are grounded on the mandatory use of detergents to isolate and extract MPs from the cell plasma membrane and the coexistence of multiple conformations, which reflects biochemical versatility and intrinsic instability of MPs. In this work ,we set out to define a new strategy to enable surface plasmon resonance (SPR) measurements on a thermostabilized and truncated version of the human adenosine (A2A) G-protein-coupled receptor (GPCR) inserted in a lipid bilayer nanodisc in a label- and detergent-free manner by using a combination of affinity tags and GFP-based fluorescence techniques. We were able to detect and characterize small molecules binding kinetics on a GPCR fully embedded in a lipid environment. By providing a comparison between different binding assays in membranes, nanodiscs and detergent micelles, we show that nanodiscs can be used for small molecule binding studies by SPR to enhance the MP stability and to trigger a more native-like behaviour when compared to kinetics on A2A receptors isolated in detergent. This work provides thus a new methodology in drug discovery to characterize the binding kinetics of small molecule ligands for MPs targets in a lipid environment.