Pharmacological Analysis and Structure Determination of 7-Methylcyanopindolol-Bound β1-Adrenergic Receptor

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

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

Structures of Adrenoceptors

The first structure of an adrenoceptor (AR), the human β2-adrenoceptor (hβ2AR) was published in 2007 and since then a total of 78 structures (up to June 2022) have been determined by X-ray crystallography and electron cryo-microscopy (cryo-EM) of all three βARs (β1, β2 and β3) and four out of six αARs (α1B, α2A, α2B, α2C). The structures are in a number of different conformational states, including the inactive state bound to an antagonist, an intermediate state bound to agonist and active states bound to agonist and an intracellular transducer (G protein or arrestin) or transducer mimetic (nanobody). The structures identify molecular details of how ligands bind in the orthosteric binding pocket (OBP; 19 antagonists, 18 agonists) and also how three different small molecule allosteric modulators bind. The structures have been used to define the molecular details of receptor activation and also the molecular determinants for transducer coupling. This chapter will give a brief overview of the structures, receptor activation, a comparison across the different subfamilies and commonalities of ligand-receptor interactions.

Adrenoceptors: Receptors, Ligands and Their Clinical Uses, Molecular Pharmacology and Assays

The nine G protein-coupled adrenoceptor subtypes are where the endogenous catecholamines adrenaline and noradrenaline interact with cells. Since they are important therapeutic targets, over a century of effort has been put into developing drugs that modify their activity. This chapter provides an outline of how we have arrived at current knowledge of the receptors, their physiological roles and the methods used to develop ligands. Initial studies in vivo and in vitro with isolated organs and tissues progressed to cell-based techniques and the use of cloned adrenoceptor subtypes together with high-throughput assays that allow close examination of receptors and their signalling pathways. The crystal structures of many of the adrenoceptor subtypes have now been determined opening up new possibilities for drug development.

IUPHAR themed review: Opioid efficacy, bias, and selectivity

Drugs acting at the opioid receptor family are clinically used to treat chronic and acute pain, though they represent the second line of treatment behind GABA analogs, antidepressants and SSRI's. Within the opioid family mu and kappa opioid receptor are commonly targeted. However, activation of the mu opioid receptor has side effects of constipation, tolerance, dependence, euphoria, and respiratory depression; activation of the kappa opioid receptor leads to dysphoria and sedation. The side effects of mu opioid receptor activation have led to mu receptor drugs being widely abused with great overdose risk. For these reasons, newer safer opioid analgesics are in high demand. For many years a focus within the opioid field was finding drugs that activated the G protein pathway at mu opioid receptor, without activating the β-arrestin pathway, known as biased agonism. Recent advances have shown that this may not be the way forward to develop safer analgesics at mu opioid receptor, though there is still some promise at the kappa opioid receptor. Here we discuss recent novel approaches to develop safer opioid drugs including efficacy vs bias and fine-tuning receptor activation by targeting sub-pockets in the orthosteric site, we explore recent works on the structural basis of bias, and we put forward the suggestion that Gα subtype selectivity may be an exciting new area of interest.

A Novel "Activation Switch" Motif Common to All Aminergic Receptors

Aminergic receptors are G protein-coupled receptors (GPCRs) that transduce signals from small endogenous biogenic amines to regulate intracellular signaling pathways. Agonist binding in the ligand binding pocket on the extracellular side opens and prepares a cavity on the intracellular face of the receptors to interact with and activate G proteins and β-arrestins. Here, by reviewing and analyzing all available aminergic receptor structures, we seek to identify activation-related conformational changes that are independent of the specific scaffold of the bound agonist, which we define as "activation conformational changes" (ACCs). While some common intracellular ACCs have been well-documented, identifying common extracellular ACCs, including those in the ligand binding pocket, is complicated by local adjustments to different ligand scaffolds. Our analysis shows no common ACCs at the extracellular ends of the transmembrane helices. Furthermore, the restricted access to the ligand binding pocket identified previously in some receptors is not universal. Notably, the Trp6.48 toggle switch and the Pro5.50-Ile3.40-Phe6.44 (PIF) motif at the bottom of the ligand binding pocket have previously been proposed to mediate the conformational consequences of ligand binding to the intracellular side of the receptors. Our analysis shows that common ACCs in the ligand binding pocket are associated with the PIF motif and nearby residues, including Trp6.48, but fails to support a shared rotamer toggle associated with activation. However, we identify two common rearrangements between the extracellular and middle subsegments, and propose a novel "activation switch" motif common to all aminergic receptors. This motif includes the middle subsegments of transmembrane helices 3, 5, and 6 and integrates both the PIF motif and Trp6.48.

Protein Design Strategies for the Structural–Functional Studies of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are an important family of membrane proteins responsible for many physiological functions in human body. High resolution GPCR structures are required to understand their molecular mechanisms and perform rational drug design, as GPCRs play a crucial role in a variety of diseases. That is difficult to obtain for the wild-type proteins because of their low stability. In this review, we discuss how this problem can be solved by using protein design strategies developed to obtain homogeneous stabilized GPCR samples for crystallization and cryoelectron microscopy.

Agonists in the Extended Conformation Stabilize the Active State of β-Adrenoceptors

In this study, we conducted a comparative analysis of the structure of agonists and antagonists of transmembrane (TM) β-adrenoceptors (β-ARs) and their interactions with the β-ARs and proposed the mechanism of receptor activation. A characteristic feature of agonist and antagonist molecules is the presence of a hydrophobic head (most often, one or two aromatic rings) and a tail with a positively charged amino group. All β-adrenergic agonists have two carbon atoms between the aromatic ring of the head and the nitrogen atom of the amino group. In antagonist molecules, this fragment can be either reduced or increased to four atoms due to the additional carbon and oxygen atoms. The agonist head, as a rule, has two H-bond donors or acceptors in the para- and meta-positions of the aromatic rings, while in the antagonist heads, these donors/acceptors are absent or located in other positions. Analysis of known three-dimensional structures of β-AR complexes with agonists showed that the agonist head forms two H-bonds with the TM5 helix, and the tail forms an ionic bond with the D3.32 residue of the TM3 helix and one or two H-bonds with the TM7 helix. The tail of the antagonist can form similar bonds, but the interaction between the head and the TM5 helix is much weaker. As a result of these interactions, the agonist molecule acquires an extended "strained string" conformation, in contrast to the antagonist molecule, which has a longer, bended, and flexible tail. The "strained string" of the agonist interacts with the TM6 helix (primarily with the W6.48 residue) and turns it, which leads to the opening of the G protein-binding site on the intracellular side of the receptor, while flexible and larger antagonist molecules do not have the same effect on the receptor.

Ligands of Adrenergic Receptors: A Structural Point of View

Adrenergic receptors are G protein-coupled receptors for epinephrine and norepinephrine. They are targets of many drugs for various conditions, including treatment of hypertension, hypotension, and asthma. Adrenergic receptors are intensively studied in structural biology, displayed for binding poses of different types of ligands. Here, we summarized molecular mechanisms of ligand recognition and receptor activation exhibited by structure. We also reviewed recent advances in structure-based ligand discovery against adrenergic receptors.

Molecular Mechanism for Ligand Recognition and Subtype Selectivity of α2C Adrenergic Receptor

Adrenergic G-protein-coupled receptors (GPCRs) mediate different cellular signaling pathways in the presence of endogenous catecholamines and play important roles in both physiological and pathological conditions. Extensive studies have been carried out to investigate the structure and function of β adrenergic receptors (βARs). However, the structure of α adrenergic receptors (αARs) remains to be determined. Here, we report the structure of the human α_2C adrenergic receptor (α_2CAR) with the non-selective antagonist, RS79948, at 2.8 Å. Our structure, mutations, modeling, and functional experiments indicate that a α_2CAR-specific D206^ECL2-R409^ECL3-Y405^6.58 network plays a role in determining α₂ adrenergic subtype selectivity. Furthermore, our results show that a specific loosened helix at the top of TM4 in α_2CAR is involved in receptor activation. Together, our structure of human α_2CAR-RS79948 provides key insight into the mechanism underlying the α₂ adrenergic receptor activation and subtype selectivity.

A Systematic Review of Inverse Agonism at Adrenoceptor Subtypes

As many, if not most, ligands at G protein-coupled receptor antagonists are inverse agonists, we systematically reviewed inverse agonism at the nine adrenoceptor subtypes. Except for β3-adrenoceptors, inverse agonism has been reported for each of the adrenoceptor subtypes, most often for β2-adrenoceptors, including endogenously expressed receptors in human tissues. As with other receptors, the detection and degree of inverse agonism depend on the cells and tissues under investigation, i.e., they are greatest when the model has a high intrinsic tone/constitutive activity for the response being studied. Accordingly, they may differ between parts of a tissue, for instance, atria vs. ventricles of the heart, and within a cell type, between cellular responses. The basal tone of endogenously expressed receptors is often low, leading to less consistent detection and a lesser extent of observed inverse agonism. Extent inverse agonism depends on specific molecular properties of a compound, but inverse agonism appears to be more common in certain chemical classes. While inverse agonism is a fascinating facet in attempts to mechanistically understand observed drug effects, we are skeptical whether an a priori definition of the extent of inverse agonism in the target product profile of a developmental candidate is a meaningful option in drug discovery and development.

A benchmark study of loop modeling methods applied to G protein-coupled receptors

G protein-coupled receptors (GPCR) are important drug discovery targets. Despite progress, many GPCR structures have not yet been solved. For these targets, comparative modeling is used in virtual ligand screening to prioritize experimental efforts. However, the structure of extracellular loop 2 (ECL2) is often poorly predicted. This is significant due to involvement of ECL2 in ligand binding for many Class A GPCR. Here we examine the performance of loop modeling protocols available in the Rosetta (cyclic coordinate descent [CCD], KIC with fragments [KICF] and next generation KIC [NGK]) and Molecular Operating Environment (MOE) software suites (de novo search). ECL2 from GPCR crystal structures served as the structure prediction targets and were divided into four sets depending on loop length. Results suggest that KICF and NGK sampled and scored more loop models with sub-angstrom and near-atomic accuracy than CCD or de novo search for loops of 24 or fewer residues. None of the methods were able to sample loop conformations with near-atomic accuracy for the longest targets ranging from 25 to 32 residues based on 1000 models generated. For these long loop targets, increased conformational sampling is necessary. The strongly conserved disulfide bond between Cys3.25 and Cys45.50 in ECL2 proved an effective filter. Setting an upper limit of 5.1 Å on the S-S distance improved the lowest RMSD model included in the top 10 scored structures in Groups 1-4 on average between 0.33 and 1.27 Å. Disulfide bond formation and geometry optimization of ECL2 provided an additional incremental benefit in structure quality.

Development of covalent antagonists for β1- and β2-adrenergic receptors

The selective covalent tethering of ligands to a specific GPCR binding site has attracted considerable interest in structural biology, molecular pharmacology and drug design. We recently reported on a covalently binding noradrenaline analog (FAUC37) facilitating crystallization of the β₂-adrenergic receptor (β₂AR^H2.64C) in an active state. We herein present the stereospecific synthesis of covalently binding disulfide ligands based on the pharmacophores of adrenergic β₁- and β₂ receptor antagonists. Radioligand depletion experiments revealed that the disulfide-functionalized ligands were able to rapidly form a covalent bond with a specific cysteine residue of the receptor mutants β₁AR^I2.64C and β₂AR^H2.64C. The propranolol derivative (S)-1a induced nearly complete irreversible blockage of the β₂AR^H2.64C within 30 min incubation. The CGP20712A-based ligand (S)-4 showed efficient covalent blocking of the β₂AR^H2.64C at very low concentrations. The analog (S)-5a revealed extraordinary covalent cross-linking at the β₁AR^I2.64C and β₂AR^H2.64C mutant while retaining a 41-fold selectivity for the β₁AR wild type over β₂AR. These compounds may serve as valuable molecular tools for studying β₁/β₂ subtype selectivity or investigations on GPCR trafficking and dimerization.

Dynamic Role of the G Protein in Stabilizing the Active State of the Adenosine A2A Receptor

Agonist binding in the extracellular region of the G protein-coupled adenosine A2A receptor increases its affinity to the G proteins in the intracellular region, and vice versa. The structural basis for this effect is not evident from the crystal structures of A_2AR in various conformational states since it stems from the receptor dynamics. Using atomistic molecular dynamics simulations on four different conformational states of the adenosine A_2A receptor, we observed that the agonists show decreased ligand mobility, lower entropy of the extracellular loops in the active-intermediate state compared with the inactive state. In contrast, the entropy of the intracellular region increases to prime the receptor for coupling the G protein. Coupling of the G protein to A_2AR shrinks the agonist binding site, making tighter receptor agonist contacts with an increase in the strength of allosteric communication compared with the active-intermediate state. These insights provide a strong basis for structure-based ligand design studies.

•

GPCR conformation dynamics reveals the forward and backward allosteric mechanism

•

Agonist binding increases the entropy in the intracellular region of the GPCR

•

G protein binding shrinks the receptor-ligand contacts in the extracellular region

•

Increased allostery between G protein and agonist in the GPCR-G protein complex

Three-dimensional descriptors for aminergic GPCRs: dependence on docking conformation and crystal structure

Three-dimensional descriptors are often used to search for new biologically active compounds, in both ligand- and structure-based approaches, capturing the spatial orientation of molecules. They frequently constitute an input for machine learning-based predictions of compound activity or quantitative structure-activity relationship modeling; however, the distribution of their values and the accuracy of depicting compound orientations might have an impact on the power of the obtained predictive models. In this study, we analyzed the distribution of three-dimensional descriptors calculated for docking poses of active and inactive compounds for all aminergic G protein-coupled receptors with available crystal structures, focusing on the variation in conformations for different receptors and crystals. We demonstrated that the consistency in compound orientation in the binding site is rather not correlated with the affinity itself, but is more influenced by other factors, such as the number of rotatable bonds and crystal structure used for docking studies. The visualizations of the descriptors distributions were prepared and made available online at http://chem.gmum.net/vischem_stability , which enables the investigation of chemical structures referring to particular data points depicted in the figures. Moreover, the performed analysis can assist in choosing crystal structure for docking studies, helping in selection of conditions providing the best discrimination between active and inactive compounds in machine learning-based experiments.

Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

The aminergic family of G protein-coupled receptors (GPCRs) plays an important role in various diseases and represents a major drug discovery target class. Structure determination of all major aminergic subfamilies has enabled structure-based ligand design for these receptors. Site-directed mutagenesis data provides an invaluable complementary source of information for elucidating the structural determinants of binding of different ligand chemotypes. The current study provides a comparative analysis of 6692 mutation data points on 34 aminergic GPCR subtypes, covering the chemical space of 540 unique ligands from mutagenesis experiments and information from experimentally determined structures of 52 distinct aminergic receptor-ligand complexes. The integrated analysis enables detailed investigation of structural receptor-ligand interactions and assessment of the transferability of combined binding mode and mutation data across ligand chemotypes and receptor subtypes. An overview is provided of the possibilities and limitations of using mutation data to guide the design of novel aminergic receptor ligands.

Structural determinants of 5-HT2B receptor activation and biased agonism

Serotonin (5-hydroxytryptamine; 5-HT) receptors modulate a variety of physiological processes ranging from perception, cognition and emotion to vascular and smooth muscle contraction, platelet aggregation, gastrointestinal function and reproduction. Drugs that interact with 5-HT receptors effectively treat diseases as diverse as migraine headaches, depression and obesity. Here we present four structures of a prototypical serotonin receptor-the human 5-HT_2B receptor-in complex with chemically and pharmacologically diverse drugs, including methysergide, methylergonovine, lisuride and LY266097. A detailed analysis of these structures complemented by comprehensive interrogation of signaling illuminated key structural determinants essential for activation. Additional structure-guided mutagenesis experiments revealed binding pocket residues that were essential for agonist-mediated biased signaling and β-arrestin2 translocation. Given the importance of 5-HT receptors for a large number of therapeutic indications, insights derived from these studies should accelerate the design of safer and more effective medications.

A Structural Framework for GPCR Chemogenomics: What's In a Residue Number?

The recent surge of crystal structures of G protein-coupled receptors (GPCRs), as well as comprehensive collections of sequence, structural, ligand bioactivity, and mutation data, has enabled the development of integrated chemogenomics workflows for this important target family. This chapter will focus on cross-family and cross-class studies of GPCRs that have pinpointed the need for, and the implementation of, a generic numbering scheme for referring to specific structural elements of GPCRs. Sequence- and structure-based numbering schemes for different receptor classes will be introduced and the remaining caveats will be discussed. The use of these numbering schemes has facilitated many chemogenomics studies such as consensus binding site definition, binding site comparison, ligand repurposing (e.g. for orphan receptors), sequence-based pharmacophore generation for homology modeling or virtual screening, and class-wide chemogenomics studies of GPCRs.

In Silico Studies Targeting G-protein Coupled Receptors for Drug Research Against Parkinson's Disease

Parkinson's Disease (PD) is a long-term neurodegenerative brain disorder that mainly affects the motor system. The causes are still unknown, and even though currently there is no cure, several therapeutic options are available to manage its symptoms. The development of novel antiparkinsonian agents and an understanding of their proper and optimal use are, indeed, highly demanding. For the last decades, L-3,4-DihydrOxyPhenylAlanine or levodopa (L-DOPA) has been the gold-standard therapy for the symptomatic treatment of motor dysfunctions associated to PD. However, the development of dyskinesias and motor fluctuations (wearing-off and on-off phenomena) associated with long-term L-DOPA replacement therapy have limited its antiparkinsonian efficacy. The investigation for non-dopaminergic therapies has been largely explored as an attempt to counteract the motor side effects associated with dopamine replacement therapy. Being one of the largest cell membrane protein families, G-Protein-Coupled Receptors (GPCRs) have become a relevant target for drug discovery focused on a wide range of therapeutic areas, including Central Nervous System (CNS) diseases. The modulation of specific GPCRs potentially implicated in PD, excluding dopamine receptors, may provide promising non-dopaminergic therapeutic alternatives for symptomatic treatment of PD. In this review, we focused on the impact of specific GPCR subclasses, including dopamine receptors, adenosine receptors, muscarinic acetylcholine receptors, metabotropic glutamate receptors, and 5-hydroxytryptamine receptors, on the pathophysiology of PD and the importance of structure- and ligand-based in silico approaches for the development of small molecules to target these receptors.

Insight into partial agonism by observing multiple equilibria for ligand-bound and Gs-mimetic nanobody-bound β1-adrenergic receptor

A complex conformational energy landscape determines G-protein-coupled receptor (GPCR) signalling via intracellular binding partners (IBPs), e.g., G_s and β-arrestin. Using ¹³C methyl methionine NMR for the β₁-adrenergic receptor, we identify ligand efficacy-dependent equilibria between an inactive and pre-active state and, in complex with G_s-mimetic nanobody, between more and less active ternary complexes. Formation of a basal activity complex through ligand-free nanobody–receptor interaction reveals structural differences on the cytoplasmic receptor side compared to the full agonist-bound nanobody-coupled form, suggesting that ligand-induced variations in G-protein interaction underpin partial agonism. Significant differences in receptor dynamics are observed ranging from rigid nanobody-coupled states to extensive μs-to-ms timescale dynamics when bound to a full agonist. We suggest that the mobility of the full agonist-bound form primes the GPCR to couple to IBPs. On formation of the ternary complex, ligand efficacy determines the quality of the interaction between the rigidified receptor and an IBP and consequently the signalling level.

β₁-adrenergic receptors are expressed in cardiac tissue and stimulated by the sympathetic nervous system. Here, the authors use NMR spectroscopy to unravel the conformational diversity upon β₁-adrenergic receptor activation and provide structural insights into partial agonism and basal activity.

G protein-coupled receptors: the evolution of structural insight

G protein-coupled receptors (GPCR) comprise a diverse superfamily of over 800 proteins that have gained relevance as biological targets for pharmaceutical drug design. Although these receptors have been investigated for decades, three-dimensional structures of GPCR have only recently become available. In this review, we focus on the technological advancements that have facilitated efforts to gain insights into GPCR structure. Progress in these efforts began with the initial crystal structure determination of rhodopsin (PDB: 1F88) in 2000 and has continued to the most recently published structure of the A_1AR (PDB: 5UEN) in 2017. Numerous experimental developments over the past two decades have opened the door for widespread GPCR structural characterization. These efforts have resulted in the determination of three-dimensional structures for over 40 individual GPCR family members. Herein we present a comprehensive list and comparative analysis of over 180 individual GPCR structures. This includes a summary of different GPCR functional states crystallized with agonists, dual agonists, partial agonists, inverse agonists, antagonists, and allosteric modulators.

Recent Advances in Structure-Based Drug Design Targeting Class A G Protein-Coupled Receptors Utilizing Crystal Structures and Computational Simulations

G protein-coupled receptors (GPCRs) represent the largest and most physiologically important integral membrane protein family, and these receptors respond to a wide variety of physiological and environmental stimuli. GPCRs are among the most critical therapeutic targets for numerous human diseases, and approximately one-third of the currently marketed drugs target this receptor family. The recent breakthroughs in GPCR structural biology have significantly contributed to our understanding of GPCR function, ligand binding, and pharmacological action as well as to the design of new drugs. This perspective highlights the latest advances in GPCR structures with a focus on the receptor-ligand interactions of each receptor family in class A nonrhodopsin GPCRs as well as the structural features for their activation, biased signaling, and allosteric mechanisms. The current state-of-the-art methodologies of structure-based drug design (SBDD) approaches in the GPCR research field are also discussed.

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.