Structure-based drug design for G protein-coupled receptors

Investigating the Unbinding of Muscarinic Antagonists from the Muscarinic 3 Receptor

Patient symptom relief is often heavily influenced by the residence time of the inhibitor-target complex. For the human muscarinic receptor 3 (hMR3), tiotropium is a long-acting bronchodilator used in conditions such as asthma or chronic obstructive pulmonary disease (COPD). The mechanistic insights into this inhibitor remain unclear; specifically, the elucidation of the main factors determining the unbinding rates could help develop the next generation of antimuscarinic agents. Using our novel unbinding algorithm, we were able to investigate ligand dissociation from hMR3. The unbinding paths of tiotropium and two of its analogues, N-methylscopolamin and homatropine methylbromide, show a consistent qualitative mechanism and allow us to identify the structural bottleneck of the process. Furthermore, our machine learning-based analysis identified key roles of the ECL2/TM5 junction involved in the transition state. Additionally, our results point to relevant changes at the intracellular end of the TM6 helix leading to the ICL3 kinase domain, highlighting the closest residue L482. This residue is located right between two main protein binding sites involved in signal transduction for hMR3's activation and regulation. We also highlight key pharmacophores of tiotropium that play determining roles in the unbinding kinetics and could aid toward drug design and lead optimization.

Effect of Ions and Sequence Variants on the Antagonist Binding Properties of the Histamine H1 Receptor

The histamine H₁ receptor (H₁R) is a G protein-coupled receptor (GPCR) and represents a main target in the treatment of allergic reactions as well as inflammatory reactions and depressions. Although the overall effect of antagonists on H₁ function has been extensively investigated, rather little is known about the potential modulatory effect of ions or sequence variants on antagonist binding. We investigated the dynamics of a phosphate ion present in the crystal structure and of a sodium ion, for which we determined the position in the allosteric pocket by metadynamics simulations. Both types of ions exhibit significant dynamics within their binding site; however, some key contacts remain stable over the simulation time, which might be exploited to develop more potent drugs targeting these sites. The dynamics of the ions is almost unaffected by the presence or absence of doxepin, as also reflected in their small effect (less than 1 kcal·mol^-1) on doxepin binding affinity. We also examined the effect of four H₁R sequence variants observed in the human population on doxepin binding. These variants cause a reduction in doxepin affinity of up to 2.5 kcal·mol^-1, indicating that personalized medical treatments that take into account individual mutation patterns could increase precision in the dosage of GPCR-targeting drugs.

Molecular Regulation of Lysophosphatidic Acid Receptor 1 Maturation and Desensitization

Lysophosphatidic acid receptor 1 (LPA1) belongs to the G protein-coupled receptor family. The ligand for LPA1 is LPA, the simplest lysophospholipid. LPA is considered a growth factor and induces cell proliferation, anti-apoptosis, and cell migration. The pro-inflammatory and pro-fibrotic roles of LPA have also been well-demonstrated. Most of the biological functions of LPA are mostly executed through LPA1. The mature form of LPA1 is glycosylated and localized on the plasma membrane. LPA1 is bound to heterotrimetric G proteins and transduces intracellular signaling in response to ligation to LPA. Desensitization of LPA1 negatively regulates LPA1-mediated signaling and the resulting biological functions. Phosphorylation and ubiquitination are well-demonstrated posttranslational modifications of GPCR. In this review, we will discuss our knowledge of LPA1 glycosylation, maturation, and trafficking from the endoplasmic reticulum (ER)/Golgi to the plasma membrane. Moreover, in light of recent findings, we will also discuss molecular regulation of LPA1 internalization and stability.

Molecular Modeling of Histamine Receptors-Recent Advances in Drug Discovery

The recent developments of fast reliable docking, virtual screening and other algorithms gave rise to discovery of many novel ligands of histamine receptors that could be used for treatment of allergic inflammatory disorders, central nervous system pathologies, pain, cancer and obesity. Furthermore, the pharmacological profiles of ligands clearly indicate that these receptors may be considered as targets not only for selective but also for multi-target drugs that could be used for treatment of complex disorders such as Alzheimer's disease. Therefore, analysis of protein-ligand recognition in the binding site of histamine receptors and also other molecular targets has become a valuable tool in drug design toolkit. This review covers the period 2014-2020 in the field of theoretical investigations of histamine receptors mostly based on molecular modeling as well as the experimental characterization of novel ligands of these receptors.

Exploring the structural basis and atomistic binding mechanistic of the selective antagonist blockade at D3 dopamine receptor over D2 dopamine receptor

More recently, there has been a paradigm shift toward selective drug targeting in the treatment of neurological disorders, including drug addiction, schizophrenia, and Parkinson's disease mediated by the different dopamine receptor subtypes. Antagonists with higher selectivity for D₃ dopamine receptor (D3DR) over D₂ dopamine receptor (D2DR) have been shown to attenuate drug-seeking behavior and associated side effects compared to non-subtype selective antagonists. However, high conservations among constituent residues of both proteins, particularly at the ligand-binding pockets, remain a challenge to therapeutic drug design. Recent studies have reported the discovery of two small-molecules R-VK4-40 and Y-QA31 which substantially inhibited D3DR with >180-fold selectivity over D2DR. Therefore, in this study, we seek to provide molecular and structural insights into these differential binding mechanistic using meta-analytic computational simulation methods. Findings revealed that R-VK4-40 and Y-QA31 adopted shallow binding modes and were more surface-exposed at D3DR while on the contrary, they exhibited deep hydrophobic pocket binding at D2DR. Also, two non-conserved residues; Tyr36^1.39 and Ser182^45.51 were identified in D3DR, based on their crucial roles and contributions to the selective binding of R-VK4-40 and Y-QA31. Importantly, both antagonists exhibited high affinities in complex with D3DR compared to D2DR, while van der Waals energies contributed majorly to their binding and stability. Structural analyses also revealed the distinct stabilizing effects of both compounds on D3DR secondary architecture relative to D2DR. Therefore, findings herein pinpointed the origin and mechanistic of selectivity of the compounds, which may assist in the rational design of potential small molecules of the D₂ -like dopamine family receptor subtype with improved potency and selectivity.

Homology Modeling of Class A G-Protein-Coupled Receptors in the Age of the Structure Boom

With 700 members, G protein-coupled receptors (GPCRs) of the rhodopsin family (class A) form the largest membrane receptor family in humans and are the target of about 30% of presently available pharmaceutical drugs. The recent boom in GPCR structures led to the structural resolution of 57 unique receptors in different states (39 receptors in inactive state only, 2 receptors in active state only and 16 receptors in different activation states). In spite of these tremendous advances, most computational studies on GPCRs, including molecular dynamics simulations, virtual screening and drug design, rely on GPCR models obtained by homology modeling. In this protocol, we detail the different steps of homology modeling with the MODELLER software, from template selection to model evaluation. The present structure boom provides closely related templates for most receptors. If, in these templates, some of the loops are not resolved, in most cases, the numerous available structures enable to find loop templates with similar length for equivalent loops. However, simultaneously, the large number of putative templates leads to model ambiguities that may require additional information based on multiple sequence alignments or molecular dynamics simulations to be resolved. Using the modeling of the human bradykinin receptor B1 as a case study, we show how several templates are managed by MODELLER, and how the choice of template(s) and of template fragments can improve the quality of the models. We also give examples of how additional information and tools help the user to resolve ambiguities in GPCR modeling.

Structure-based discovery and development of metabotropic glutamate receptor 5 negative allosteric modulators

The metabotropic glutamate (mGlu) receptors are a family of eight class C G protein-coupled receptors (GPCRs) which modulate cell signaling and synaptic transmission to the major excitatory neurotransmitter l-glutamate (l-glutamic acid). Due to their role in modulating glutamate response, their widespread distribution in the central nervous system (CNS) and some evidence of dysregulation in disease, the mGlu receptors have become attractive pharmacological targets. As the orthosteric (glutamate) binding site is highly conserved across the eight mGlu receptors, it is difficult not only to generate ligands with subtype selectivity but, due to the nature of the binding site, with suitable drug-like properties to allow oral bioavailability and CNS penetration. Selective pharmacological targeting of a single receptor subtype can be achieved by targeting alternative (allosteric) binding sites. The nature of the allosteric binding pockets allows ligands to be developed that have good physical chemical properties as evidenced by several allosteric modulators of mGlu receptors entering clinical trials. The first negative allosteric modulators of the metabotropic glutamate 5 (mGlu₅) receptor were discovered from high throughput screening activities. An alternative approach to drug discovery is to use structural knowledge to enable structure-based drug design (SBDD), which allows the design of molecules in a more rational, rather than empirical, fashion. Here we will describe the process of SBDD in the discovery of the mGlu₅ negative allosteric modulator HTL0014242 and describe how knowledge of receptor structure can also be used to gain insights into the receptor activation mechanisms.

Advances in therapeutic peptides targeting G protein-coupled receptors

Dysregulation of peptide-activated pathways causes a range of diseases, fostering the discovery and clinical development of peptide drugs. Many endogenous peptides activate G protein-coupled receptors (GPCRs) - nearly 50 GPCR peptide drugs have been approved to date, most of them for metabolic disease or oncology, and more than 10 potentially first-in-class peptide therapeutics are in the pipeline. The majority of existing peptide therapeutics are agonists, which reflects the currently dominant strategy of modifying the endogenous peptide sequence of ligands for peptide-binding GPCRs. Increasingly, novel strategies are being employed to develop both agonists and antagonists, to both introduce chemical novelty and improve drug-like properties. Pharmacodynamic improvements are evolving to allow biasing ligands to activate specific downstream signalling pathways, in order to optimize efficacy and reduce side effects. In pharmacokinetics, modifications that increase plasma half-life have been revolutionary. Here, we discuss the current status of the peptide drugs targeting GPCRs, with a focus on evolving strategies to improve pharmacokinetic and pharmacodynamic properties.

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.

Distinct inactive conformations of the dopamine D2 and D3 receptors correspond to different extents of inverse agonism

By analyzing and simulating inactive conformations of the highly homologous dopamine D2 and D3 receptors (D2R and D3R), we find that eticlopride binds D2R in a pose very similar to that in the D3R/eticlopride structure but incompatible with the D2R/risperidone structure. In addition, risperidone occupies a sub-pocket near the Na+ binding site, whereas eticlopride does not. Based on these findings and our experimental results, we propose that the divergent receptor conformations stabilized by Na+-sensitive eticlopride and Na+-insensitive risperidone correspond to different degrees of inverse agonism. Moreover, our simulations reveal that the extracellular loops are highly dynamic, with spontaneous transitions of extracellular loop 2 from the helical conformation in the D2R/risperidone structure to an extended conformation similar to that in the D3R/eticlopride structure. Our results reveal previously unappreciated diversity and dynamics in the inactive conformations of D2R. These findings are critical for rational drug discovery, as limiting a virtual screen to a single conformation will miss relevant ligands.

Molecular Mechanism of S1P Binding and Activation of the S1P1 Receptor

Sphingosine-1-phosphate (S1P) is a lipidic mediator in mammals that functions either as a second messenger or as a ligand. In the latter case, it is transported by its HDL-associated apoM carrier and circulated in blood where it binds to specific S1P receptors on cell membranes and induces downstream reactions. Although S1P signaling pathways are essential for many biological processes, they are poorly understood at the molecular level. Here, the solved crystal structures of the S1P₁ receptor were used to evaluate molecular dynamics (MD) simulations to generate greater detailed molecular insights into the mechanism of S1P signaling. The MD simulations provided observations at the coarse-grained and atomic levels indicating that S1P may access the receptor binding pocket directly from solvents. Lifting of the bulky N-terminal cap region of the receptor precedes initial S1P binding. Glu121^3.29 guides S1P penetration, and together with Arg292^7.34 is responsible for the stabilization of S1P in the binding pocket, which is consistent with experimental predictions. The complete binding of S1P is followed by receptor activation, wherein Trp269^6.48 moves toward the transmembrane helix (TM) 7, resulting in the formation of an enhanced hydrogen bond network in the lower region of TM7. The distance between TM3 and TM6 is subsequently increased, resulting in the opening of the intracellular binding pocket that enables G protein binding. Further analysis of the force distribution network in the receptor yielded a detailed molecular understanding of the signal transmission network that is activated upon agonist binding.

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.

Interplay of cysteine exposure and global protein dynamics in small-molecule recognition by a regulator of G-protein signaling protein

Regulator of G protein signaling (RGS) proteins play a pivotal role in regulation of G protein-coupled receptor (GPCR) signaling and are therefore becoming an increasingly important therapeutic target. Recently discovered thiadiazolidinone (TDZD) compounds that target cysteine residues have shown different levels of specificities and potencies for the RGS4 protein, thereby suggesting intrinsic differences in dynamics of this protein upon binding of these compounds. In this work, we investigated using atomistic molecular dynamics (MD) simulations the effect of binding of several small-molecule inhibitors on perturbations and dynamical motions in RGS4. Specifically, we studied two conformational models of RGS4 in which a buried cysteine residue is solvent-exposed due to side-chain motions or due to flexibility in neighboring helices. We found that TDZD compounds with aromatic functional groups perturb the RGS4 structure more than compounds with aliphatic functional groups. Moreover, small-molecules with aromatic functional groups but lacking sulfur atoms only transiently reside within the protein and spontaneously dissociate to the solvent. We further measured inhibitory effects of TDZD compounds using a protein-protein interaction assay on a single-cysteine RGS4 protein showing trends in potencies of compounds consistent with our simulation studies. Thermodynamic analyses of RGS4 conformations in the apo-state and on binding to TDZD compounds revealed links between both conformational models of RGS4. The exposure of cysteine side-chains appears to facilitate initial binding of TDZD compounds followed by migration of the compound into a bundle of four helices, thereby causing allosteric perturbations in the RGS/Gα protein-protein interface.

Structure-guided development of selective M3 muscarinic acetylcholine receptor antagonists

The development of selective antagonists for muscarinic acetylcholine receptors is challenging due to high homology in orthosteric binding sites among subtypes. Starting from a single amino acid difference in the orthosteric pockets in M2 muscarinic acetylcholine receptor (M2R) and M3R, we developed an M3R-selective antagonist using molecular docking and structure-based design. The resulting M3R antagonist showed up to 100-fold selectivity over the M2R in affinity and 1,000-fold selectivity in vivo. The docking-predicted geometry was further confirmed by a 3.1 Å crystal structure of M3R in complex with the selective antagonist. The potential of structure-based design to develop selective drugs with reduced off-target effects is supported by this study.

Drugs that treat chronic obstructive pulmonary disease by antagonizing the M3 muscarinic acetylcholine receptor (M3R) have had a significant effect on health, but can suffer from their lack of selectivity against the M2R subtype, which modulates heart rate. Beginning with the crystal structures of M2R and M3R, we exploited a single amino acid difference in their orthosteric binding pockets using molecular docking and structure-based design. The resulting M3R antagonists had up to 100-fold selectivity over M2R in affinity and over 1,000-fold selectivity in vivo. The crystal structure of the M3R-selective antagonist in complex with M3R corresponded closely to the docking-predicted geometry, providing a template for further optimization.

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

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

Computational Methods Used in Hit-to-Lead and Lead Optimization Stages of Structure-Based Drug Discovery

GPCR modeling approaches are widely used in the hit-to-lead (H2L) and lead optimization (LO) stages of drug discovery. The aims of these modeling approaches are to predict the 3D structures of the receptor-ligand complexes, to explore the key interactions between the receptor and the ligand and to utilize these insights in the design of new molecules with improved binding, selectivity or other pharmacological properties. In this book chapter, we present a brief survey of key computational approaches integrated with hierarchical GPCR modeling protocol (HGMP) used in hit-to-lead (H2L) and in lead optimization (LO) stages of structure-based drug discovery (SBDD). We outline the differences in modeling strategies used in H2L and LO of SBDD and illustrate how these tools have been applied in three drug discovery projects.

GPCR Homology Model Generation for Lead Optimization

The vast increase of recently solved GPCR X-ray structures forms the basis for GPCR homology modeling to atomistic accuracy. Nowadays, homology models can be employed for GPCR-ligand optimization and have been reported as invaluable tools for drug design in the last few years. Elucidation of the complex GPCR pharmacology and the associated GPCR conformations made clear that different homology models have to be constructed for different activation states of the GPCRs. Therefore, templates have to be chosen accordingly to their sequence homology as well as to their activation state. The subsequent ligand placement is nontrivial, as some recent X-ray structures show very unusual ligand binding sites and solvent involvement, expanding the space of the putative ligand binding site from the generic retinal binding pocket to the whole receptor. In the present study, a workflow is presented starting from the selection of the target sequence, guiding through the GPCR modeling process, and finishing with ligand placement and pose validation.

Open Boundary Simulations of Proteins and Their Hydration Shells by Hamiltonian Adaptive Resolution Scheme

The recently proposed Hamiltonian adaptive resolution scheme (H-AdResS) allows the performance of molecular simulations in an open boundary framework. It allows changing, on the fly, the resolution of specific subsets of molecules (usually the solvent), which are free to diffuse between the atomistic region and the coarse-grained reservoir. So far, the method has been successfully applied to pure liquids. Coupling the H-AdResS methodology to hybrid models of proteins, such as the molecular mechanics/coarse-grained (MM/CG) scheme, is a promising approach for rigorous calculations of ligand binding free energies in low-resolution protein models. Toward this goal, here we apply for the first time H-AdResS to two atomistic proteins in dual-resolution solvent, proving its ability to reproduce structural and dynamic properties of both the proteins and the solvent, as obtained from atomistic simulations.

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

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

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

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

A staged screening of registered drugs highlights remyelinating drug candidates for clinical trials

There is no treatment for the myelin loss in multiple sclerosis, ultimately resulting in the axonal degeneration that leads to the progressive phase of the disease. We established a multi-tiered platform for the sequential screening of drugs that could be repurposed as remyelinating agents. We screened a library of 2,000 compounds (mainly Food and Drug Administration (FDA)-approved compounds and natural products) for cellular metabolic activity on mouse oligodendrocyte precursors (OPC), identifying 42 molecules with significant stimulating effects. We then characterized the effects of these compounds on OPC proliferation and differentiation in mouse glial cultures, and on myelination and remyelination in organotypic cultures. Three molecules, edaravone, 5-methyl-7-methoxyisoflavone and lovastatin, gave positive results in all screening tiers. We validated the results by retesting independent stocks of the compounds, analyzing their purity, and performing dose-response curves. To identify the chemical features that may be modified to enhance the compounds' activity, we tested chemical analogs and identified, for edaravone, the functional groups that may be essential for its activity. Among the selected remyelinating candidates, edaravone appears to be of strong interest, also considering that this drug has been approved as a neuroprotective agent for acute ischemic stroke and amyotrophic lateral sclerosis in Japan.

Visualizing GPCR 'Megaplexes' Which Enable Sustained Intracellular Signaling

A range of cutting-edge techniques have been employed to visualize 'megaplexes' consisting of a G protein-coupled receptor (GPCR) bound to β-arrestin in intracellular endosomes following agonist-induced internalization. Surprisingly, the complex includes simultaneous binding of the heterotrimeric G protein, which retains full functional activity and supports sustained signaling from within the cell.

Discovery of nanomolar ligands with novel scaffolds for the histamine H4 receptor by virtual screening

The involvement of histamine H4 receptor (H4R) in immune cells chemotaxis and mediator release makes it an attractive target for the treatment of inflammation disorders. A decade of medicinal chemistry efforts has led to several promising ligands, although the chemical structures described so far possesses a singular limited diversity. We report here the discovery of novel structures, belonging to completely different scaffolds. The virtual screening was planed as a two-steps process. First, using a "scout screening" methodology, we have experimentally probed the H4R ligand binding site using a small size chemical library with very diverse structures, and identified a hit that further assist us in refining a raw 3D homology model. Second, the refined 3D model was used to conduct a widened virtual screening. This two-steps strategy proved to be very successful, both in terms of structural diversity and hit rate (23%). Moreover, the hits have high affinity for the H4R, with most potent ligands in the nanomolar range.

Integrating structural and mutagenesis data to elucidate GPCR ligand binding

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins, as well as drug targets. A recent boom in GPCR structural biology has provided detailed images of receptor ligand binding sites and interactions on the molecular level. An ever-increasing number of ligands is reported that exhibit activity through multiple receptors, binding in allosteric sites, and bias towards different intracellular signalling pathways. Furthermore, a wealth of single point mutants has accumulated in literature and public databases. Integrating these structural and mutagenesis data will help elucidate new GPCR ligand binding sites, and ultimately design drugs with tailored pharmacological activity.

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.

Function-specific virtual screening for GPCR ligands using a combined scoring method

The ability of scoring functions to correctly select and rank docking poses of small molecules in protein binding sites is highly target dependent, which presents a challenge for structure-based drug discovery. Here we describe a virtual screening method that combines an energy-based docking scoring function with a molecular interaction fingerprint (IFP) to identify new ligands based on G protein-coupled receptor (GPCR) crystal structures. The consensus scoring method is prospectively evaluated by: 1) the discovery of chemically novel, fragment-like, high affinity histamine H₁ receptor (H₁R) antagonists/inverse agonists, 2) the selective structure-based identification of ß₂-adrenoceptor (ß₂R) agonists, and 3) the experimental validation and comparison of the combined and individual scoring approaches. Systematic retrospective virtual screening simulations allowed the definition of scoring cut-offs for the identification of H₁R and ß₂R ligands and the selection of an optimal ß-adrenoceptor crystal structure for the discrimination between ß₂R agonists and antagonists. The consensus approach resulted in the experimental validation of 53% of the ß₂R and 73% of the H₁R virtual screening hits with up to nanomolar affinities and potencies. The selective identification of ß₂R agonists shows the possibilities of structure-based prediction of GPCR ligand function by integrating protein-ligand binding mode information.

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.

Application of advanced X-ray methods in life sciences

BACKGROUND

Synchrotron radiation (SR) sources provide diverse X-ray methods for the investigation of structure-function relationships in biological macromolecules.

SCOPE OF REVIEW

Recent developments in SR sources and in the X-ray tools they offer for life sciences are reviewed. Specifically, advances in macromolecular crystallography, small angle X-ray solution scattering, X-ray absorption and fluorescence spectroscopy, and imaging are discussed with examples.

MAJOR CONCLUSIONS

SR sources offer a range of X-ray techniques that can be used in a complementary fashion in studies of biological systems at a wide range of resolutions from atomic to cellular scale. Emerging applications of X-ray techniques include the characterization of disordered proteins, noncrystalline and nonequilibrium systems, elemental imaging of tissues, cells and organs, and detection of time-resolved changes in molecular structures.

GENERAL SIGNIFICANCE

X-ray techniques are in the center of hybrid approaches that are used to gain insight into complex problems relating to biomolecular mechanisms, disease and possible therapeutic solutions. This article is part of a Special Issue entitled "Science for Life". Guest Editors: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Small Molecule CXCR3 Antagonists

Chemokines and their receptors are known to play important roles in disease. More than 40 chemokine ligands and 20 chemokine receptors have been identified, but, to date, only two small molecule chemokine receptor antagonists have been approved by the FDA. The chemokine receptor CXCR3 was identified in 1996, and nearly 20 years later, new areas of CXCR3 disease biology continue to emerge. Several classes of small molecule CXCR3 antagonists have been developed, and two have shown efficacy in preclinical models of inflammatory disease. However, only one CXCR3 antagonist has been evaluated in clinical trials, and there remain many opportunities to further investigate known classes of CXCR3 antagonists and to identify new chemotypes. This Perspective reviews the known CXCR3 antagonists and considers future opportunities for the development of small molecules for clinical evaluation.

A Molecular Pharmacologist's Guide to G Protein-Coupled Receptor Crystallography

G protein-coupled receptor (GPCR) structural biology has progressed dramatically in the last decade. There are now over 120 GPCR crystal structures deposited in the Protein Data Bank of 32 different receptors from families scattered across the phylogenetic tree, including class B, C, and Frizzled GPCRs. These structures have been obtained in combination with a wide variety of ligands and captured in a range of conformational states. This surge in structural knowledge has enlightened research into the molecular recognition of biologically active molecules, the mechanisms of receptor activation, the dynamics of functional selectivity, and fueled structure-based drug design efforts for GPCRs. Here we summarize the innovations in both protein engineering/molecular biology and crystallography techniques that have led to these advances in GPCR structural biology and discuss how they may influence the resulting structural models. We also provide a brief molecular pharmacologist's guide to GPCR X-ray crystallography, outlining some key aspects in the process of structure determination, with the goal to encourage noncrystallographers to interrogate structures at the molecular level. Finally, we show how chemogenomics approaches can be used to marry the wealth of existing receptor pharmacology data with the expanding repertoire of structures, providing a deeper understanding of the mechanistic details of GPCR function.

Optimized Proteomic Mass Spectrometry Characterization of Recombinant Human μ-Opioid Receptor Functionally Expressed in Pichia pastoris Cell Lines

Human μ-opioid receptor (hMOR) is a class-A G-protein-coupled receptor (GPCR), a prime therapeutic target for the management of moderate and severe pain. A chimeric form of the receptor has been cocrystallized with an opioid antagonist and resolved by X-ray diffraction; however, further direct structural analysis is still required to identify the active form of the receptor to facilitate the rational design of hMOR-selective agonist and antagonists with therapeutic potential. Toward this goal and in spite of the intrinsic difficulties posed by the highly hydrophobic transmembrane motives of hMOR, we have comprehensively characterized by mass spectrometry (MS) analysis the primary sequence of the functional hMOR. Recombinant hMOR was overexpressed as a C-terminal c-myc and 6-his tagged protein using an optimized expression procedure in Pichia pastoris cells. After membrane solubilization and metal-affinity chromatography purification, a procedure was devised to enhance the concentration of the receptor. Subsequent combinations of in-solution and in-gel digestions using either trypsin, chymotrypsin, or proteinase K, followed by matrix-assisted laser desorption ionization time-of-flight MS or nanoliquid chromatography coupled with tandem MS analyses afforded an overall sequence coverage of up to >80%, a level of description first attained for an opioid receptor and one of the six such high-coverage MS-based analyses of any GPCR.

Selectivity and anti-Parkinson's potential of thiadiazolidinone RGS4 inhibitors

Many current therapies target G protein coupled receptors (GPCR), transporters, or ion channels. In addition to directly targeting these proteins, disrupting the protein-protein interactions that localize or regulate their function could enhance selectivity and provide unique pharmacologic actions. Regulators of G protein signaling (RGS) proteins, especially RGS4, play significant roles in epilepsy and Parkinson's disease. Thiadiazolidinone (TDZD) inhibitors of RGS4 are nanomolar potency blockers of the biochemical actions of RGS4 in vitro. Here, we demonstrate the substantial selectivity (8- to >5000-fold) of CCG-203769 for RGS4 over other RGS proteins. It is also 300-fold selective for RGS4 over GSK-3β, another target of this class of chemical scaffolds. It does not inhibit the cysteine protease papain at 100 μM. CCG-203769 enhances Gαq-dependent cellular Ca(2+) signaling in an RGS4-dependent manner. TDZD inhibitors also enhance Gαi-dependent δ-OR inhibition of cAMP production in SH-SY-5Y cells, which express endogenous receptors and RGS4. Importantly, CCG-203769 potentiates the known RGS4 mechanism of Gαi-dependent muscarinic bradycardia in vivo. Furthermore, it reverses raclopride-induced akinesia and bradykinesia in mice, a model of some aspects of the movement disorder in Parkinson's disease. A broad assessment of compound effects revealed minimal off-target effects at concentrations necessary for cellular RGS4 inhibition. These results expand our understanding of the mechanism and specificity of TDZD RGS inhibitors and support the potential for therapeutic targeting of RGS proteins in Parkinson's disease and other neural disorders.

Structural insights into ligand recognition and selectivity for classes A, B, and C GPCRs

The G protein-coupled receptor (GPCR) superfamily constitutes the largest collection of cell surface signaling proteins with approximately 800 members in the human genome. GPCRs regulate virtually all aspects of physiology and they are an important class of drug targets with ~30% of drugs on the market targeting a GPCR. Breakthroughs in GPCR structural biology in recent years have significantly expanded our understanding of GPCR structure and function and ushered in a new era of structure-based drug design for GPCRs. Crystal structures for nearly thirty distinct GPCRs are now available including receptors from each of the major classes, A, B, C, and F. These structures provide a foundation for understanding the molecular basis of GPCR pharmacology. Here, we review structural mechanisms of ligand recognition and selectivity of GPCRs with a focus on selected examples from classes A, B, and C, and we highlight major unresolved questions for future structural studies.

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.

Large-scale production and protein engineering of G protein-coupled receptors for structural studies

Structural studies of G protein-coupled receptors (GPCRs) gave insights into molecular mechanisms of their action and contributed significantly to molecular pharmacology. This is primarily due to technical advances in protein engineering, production and crystallization of these important receptor targets. On the other hand, NMR spectroscopy of GPCRs, which can provide information about their dynamics, still remains challenging due to difficulties in preparation of isotopically labeled receptors and their low long-term stabilities. In this review, we discuss methods used for expression and purification of GPCRs for crystallographic and NMR studies. We also summarize protein engineering methods that played a crucial role in obtaining GPCR crystal structures.

GPCR structure, function, drug discovery and crystallography: report from Academia-Industry International Conference (UK Royal Society) Chicheley Hall, 1-2 September 2014

G-protein coupled receptors (GPCRs) are the targets of over half of all prescribed drugs today. The UniProt database has records for about 800 proteins classified as GPCRs, but drugs have only been developed against 50 of these. Thus, there is huge potential in terms of the number of targets for new therapies to be designed. Several breakthroughs in GPCRs biased pharmacology, structural biology, modelling and scoring have resulted in a resurgence of interest in GPCRs as drug targets. Therefore, an international conference, sponsored by the Royal Society, with world-renowned researchers from industry and academia was recently held to discuss recent progress and highlight key areas of future research needed to accelerate GPCR drug discovery. Several key points emerged. Firstly, structures for all three major classes of GPCRs have now been solved and there is increasing coverage across the GPCR phylogenetic tree. This is likely to be substantially enhanced with data from x-ray free electron sources as they move beyond proof of concept. Secondly, the concept of biased signalling or functional selectivity is likely to be prevalent in many GPCRs, and this presents exciting new opportunities for selectivity and the control of side effects, especially when combined with increasing data regarding allosteric modulation. Thirdly, there will almost certainly be some GPCRs that will remain difficult targets because they exhibit complex ligand dependencies and have many metastable states rendering them difficult to resolve by crystallographic methods. Subtle effects within the packing of the transmembrane helices are likely to mask and contribute to this aspect, which may play a role in species dependent behaviour. This is particularly important because it has ramifications for how we interpret pre-clinical data. In summary, collaborative efforts between industry and academia have delivered significant progress in terms of structure and understanding of GPCRs and will be essential for resolving problems associated with the more difficult targets in the future.

Fluorescent knock-in mice to decipher the physiopathological role of G protein-coupled receptors

G protein-coupled receptors (GPCRs) modulate most physiological functions but are also critically involved in numerous pathological states. Approximately a third of marketed drugs target GPCRs, which places this family of receptors in the main arena of pharmacological pre-clinical and clinical research. The complexity of GPCR function demands comprehensive appraisal in native environment to collect in-depth knowledge of receptor physiopathological roles and assess the potential of therapeutic molecules. Identifying neurons expressing endogenous GPCRs is therefore essential to locate them within functional circuits whereas GPCR visualization with subcellular resolution is required to get insight into agonist-induced trafficking. Both remain frequently poorly investigated because direct visualization of endogenous receptors is often hampered by the lack of appropriate tools. Also, monitoring intracellular trafficking requires real-time visualization to gather in-depth knowledge. In this context, knock-in mice expressing a fluorescent protein or a fluorescent version of a GPCR under the control of the endogenous promoter not only help to decipher neuroanatomical circuits but also enable real-time monitoring with subcellular resolution thus providing invaluable information on their trafficking in response to a physiological or a pharmacological challenge. This review will present the animal models and discuss their contribution to the understanding of the physiopathological role of GPCRs. We will also address the drawbacks associated with this methodological approach and browse future directions.

Purification and Crystallization of a Thermostabilized Agonist-Bound Conformation of the Human Adenosine A(2A) Receptor

Crystallization of G protein-coupled receptors (GPCRs) is successful due to the development of generic protein engineering strategies, which has resulted in the structure determination of more than 25 GPCRs, including representatives from class A, B, C, and F. Most of the X-ray structures available correspond to an inactive conformation of the receptor bound to an antagonist. Only a few high-resolution structures of agonist-bound conformations of GPCRs have been determined over the last 6 years. Here, we describe the purification and crystallization protocols of a thermostabilized agonist-bound conformation of the human adenosine A2A receptor.

Insights into a defined secondary binding region on β-adrenoceptors and putative roles in ligand binding and drug design

Putative roles of a secondary binding region shared among beta-adrenoceptors.

Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models

We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.

Structure-based virtual screening for fragment-like ligands of the G protein-coupled histamine H4 receptor

Structure-based virtual screening using H₁R- and β₂R-based histamine H₄R homology models identified 9 fragments with an affinity ranging from 0.14 to 6.3 μm for H₄R.

The receptor concept in 3D: from hypothesis and metaphor to GPCR-ligand structures

The first mentioning of the word "receptor" for the structure with which a bioactive compound should react for obtaining its specific influence on a physiological system goes back to the years around 1900. The receptor concept was adapted from the lock and key theory for the enzyme substrate and blockers interactions. Through the years the concept, in the beginning rather being a metaphor, not a model, was refined and became reality in recent years. Not only the structures of receptors were elucidated, also the receptor machineries were unraveled. Following a brief historical review we will describe how the recent breakthroughs in the experimental determination of G protein-coupled receptor (GPCR) crystal structures can be complemented by computational modeling, medicinal chemistry, biochemical, and molecular pharmacological studies to obtain new insights into the molecular determinants of GPCR-ligand binding and activation. We will furthermore discuss how this information can be used for structure-based discovery of novel GPCR ligands that bind specific (allosteric) binding sites with desired effects on GPCR functional activity.