The importance of interactions with helix 5 in determining the efficacy of β-adrenoceptor ligands

Determinants of Neutral Antagonism and Inverse Agonism in the β2-Adrenergic Receptor

Free-energy profiles for the activation/deactivation of the β2-adrenergic receptor (ADRB2) with neutral antagonist and inverse agonist ligands have been determined with well-tempered multiple-walker (MW) metadynamics simulations. The inverse agonists carazolol and ICI118551 clearly favor single inactive conformational minima in both the binary and ternary ligand-receptor-G-protein complexes, in accord with the inverse-agonist activity of the ligands. The behavior of neutral antagonists is more complex, as they seem also to affect the recruitment of the G-protein. The results are analyzed in terms of the conformational states of the well-known microswitches that have been proposed as indicators of receptor activity.

Conformational selection guides β-arrestin recruitment at a biased G protein-coupled receptor

G protein-coupled receptors (GPCRs) recruit β-arrestins to coordinate diverse cellular processes, but the structural dynamics driving this process are poorly understood. Atypical chemokine receptors (ACKRs) are intrinsically biased GPCRs that engage β-arrestins but not G proteins, making them a model system for investigating the structural basis of β-arrestin recruitment. Here, we performed nuclear magnetic resonance (NMR) experiments on ¹³CH₃-ε-methionine-labeled ACKR3, revealing that β-arrestin recruitment is associated with conformational exchange at key regions of the extracellular ligand-binding pocket and intracellular β-arrestin-coupling region. NMR studies of ACKR3 mutants defective in β-arrestin recruitment identified an allosteric hub in the receptor core that coordinates transitions among heterogeneously populated and selected conformational states. Our data suggest that conformational selection guides β-arrestin recruitment by tuning receptor dynamics at intracellular and extracellular regions.

An Interpretable Convolutional Neural Network Framework for Analyzing Molecular Dynamics Trajectories: a Case Study on Functional States for G-Protein-Coupled Receptors

Molecular dynamics (MD) simulations have made great contribution to revealing structural and functional mechanisms for many biomolecular systems. However, how to identify functional states and important residues from vast conformation space generated by MD remains challenging; thus an intelligent navigation is highly desired. Despite intelligent advantages of deep learning exhibited in analyzing MD trajectory, its black-box nature limits its application. To address this problem, we explore an interpretable convolutional neural network (CNN)-based deep learning framework to automatically identify diverse active states from the MD trajectory for G-protein-coupled receptors (GPCRs), named the ICNNMD model. To avoid the information loss in representing the conformation structure, the pixel representation is introduced, and then the CNN module is constructed to efficiently extract features followed by a fully connected neural network to realize the classification task. More importantly, we design a local interpretable model-agnostic explanation interpreter for the classification result by local approximation with a linear model, through which important residues underlying distinct active states can be quickly identified. Our model showcases higher than 99% classification accuracy for three important GPCR systems with diverse active states. Notably, some important residues in regulating different biased activities are successfully identified, which are beneficial to elucidating diverse activation mechanisms for GPCRs. Our model can also serve as a general tool to analyze MD trajectory for other biomolecular systems. All source codes are freely available at https://github.com/Jane-Liu97/ICNNMD for aiding MD studies.

N-Palmitoylglycine and other N-acylamides activate the lipid receptor G2A/GPR132

The G-protein-coupled receptor GPR132, also known as G2A, is activated by 9-hydroxyoctadecadienoic acid (9-HODE) and other oxidized fatty acids. Other suggested GPR132 agonists including lysophosphatidylcholine (LPC) have not been readily reproduced. Here, we identify N-acylamides in particular N-acylglycines, as lipid activators of GPR132 with comparable activity to 9-HODE. The order-of-potency is N-palmitoylglycine > 9-HODE ≈ N-linoleoylglycine > linoleamide > N-oleoylglycine ≈ N-stereoylglycine > N-arachidonoylglycine > N-docosehexanoylglycine. Physiological concentrations of N-acylglycines in tissue are sufficient to activate GPR132. N-linoleoylglycine and 9-HODE also activate rat and mouse GPR132, despite limited sequence conservation to human. We describe pharmacological tools for GPR132, identified through drug screening. SKF-95667 is a novel GPR132 agonist. SB-583831 and SB-583355 are peptidomimetic molecules containing core amino acids (glycine and phenylalanine, respectively), and structurally related to previously described ligands. A telmisartan analog, GSK1820795A, antagonizes the actions of N-acylamides at GPR132. The synthetic cannabinoid CP-55 940 also activates GPR132. Molecular docking to a homology model suggested a site for lipid binding, predicting the acyl side-chain to extend into the membrane bilayer between TM4 and TM5 of GPR132. Small-molecule ligands are envisaged to occupy a "classical" site encapsulated in the 7TM bundle. Structure-directed mutagenesis indicates a critical role for arginine at position 203 in transmembrane domain 5 to mediate GPR132 activation by N-acylamides. Our data suggest distinct modes of binding for small-molecule and lipid agonists to the GPR132 receptor. Antagonists, such as those described here, will be vital to understand the physiological role of this long-studied target.

Therapeutic Targets for Treatment of Heart Failure: Focus on GRKs and β-Arrestins Affecting βAR Signaling

Heart failure (HF) is a heart disease that is classified into two main types: HF with reduced ejection fraction (HFrEF) and HF with preserved ejection fraction (HFpEF). Both types of HF lead to significant risk of mortality and morbidity. Pharmacological treatment with β-adrenergic receptor (βAR) antagonists (also called β-blockers) has been shown to reduce the overall hospitalization and mortality rates and improve the clinical outcomes in HF patients with HFrEF but not HFpEF. Although, the survival rate of patients suffering from HF continues to drop, the management of HF still faces several limitations and discrepancies highlighting the need to develop new treatment strategies. Overstimulation of the sympathetic nervous system is an adaptive neurohormonal response to acute myocardial injury and heart damage, whereas prolonged exposure to catecholamines causes defects in βAR regulation, including a reduction in the amount of βARs and an increase in βAR desensitization due to the upregulation of G protein-coupled receptor kinases (GRKs) in the heart, contributing in turn to the progression of HF. Several studies show that myocardial GRK2 activity and expression are raised in the failing heart. Furthermore, β-arrestins play a pivotal role in βAR desensitization and, interestingly, can mediate their own signal transduction without any G protein-dependent pathway involved. In this review, we provide new insight into the role of GRKs and β-arrestins on how they affect βAR signaling regarding the molecular and cellular pathophysiology of HF. Additionally, we discuss the therapeutic potential of targeting GRKs and β-arrestins for the treatment of HF.

Structure-inspired design of β-arrestin-biased ligands for aminergic GPCRs

Development of biased ligands targeting G protein-coupled receptors (GPCRs) is a promising approach for current drug discovery. Although structure-based drug design of biased agonists remains challenging even with an abundance of GPCR crystal structures, we present an approach for translating GPCR structural data into β-arrestin-biased ligands for aminergic GPCRs. We identified specific amino acid-ligand contacts at transmembrane helix 5 (TM5) and extracellular loop 2 (EL2) responsible for Gi/o and β-arrestin signaling, respectively, and targeted those residues to develop biased ligands. For these ligands, we found that bias is conserved at other aminergic GPCRs that retain similar residues at TM5 and EL2. Our approach provides a template for generating arrestin-biased ligands by modifying predicted ligand interactions that block TM5 interactions and promote EL2 interactions. This strategy could facilitate the structure-guided design of arrestin-biased ligands at other GPCRs, including polypharmacological biased ligands.

Evidence for a shape-based recognition of odorants in vivo in the human nose from an analysis of the molecular mechanism of lily-of-the-valley odorants detection in the Lilial and Bourgeonal family using the C/Si/Ge/Sn switch strategy

We performed an analysis of possible mechanisms of ligand recognition in the human nose. The analysis is based on in vivo odor threshold determination and in vitro Ca2+ imaging assays with a C/Si/Ge/Sn switch strategy applied to the compounds Lilial and Bourgeonal, to differentiate between different molecular mechanisms of odorant detection. Our results suggest that odorant detection under threshold conditions is mainly based on the molecular shape, i.e. the van der Waals surface, and electrostatics of the odorants. Furthermore, we show that a single olfactory receptor type is responsible for odor detection of Bourgeonal at the threshold level in humans in vivo. Carrying out a QM analysis of vibrational energies contained in the odorants, there is no evidence for a vibration-based recognition.

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

Comparisons between structures of the β1-adrenergic receptor (AR) bound to either agonists, partial agonists, or weak partial agonists led to the proposal that rotamer changes of Ser(5.46), coupled to a contraction of the binding pocket, are sufficient to increase the probability of receptor activation. (RS)-4-[3-(tert-butylamino)-2-hydroxypropoxy]-1H-indole-2-carbonitrile (cyanopindolol) is a weak partial agonist of β1AR and, based on the hypothesis above, we predicted that the addition of a methyl group to form 4-[(2S)-3-(tert-butylamino)-2-hydroxypropoxy]-7-methyl-1H-indole-2-carbonitrile (7-methylcyanopindolol) would dramatically reduce its efficacy. An eight-step synthesis of 7-methylcyanopindolol was developed and its pharmacology was analyzed. 7-Methylcyanopindolol bound with similar affinity to cyanopindolol to both β1AR and β2AR. As predicted, the efficacy of 7-methylcyanopindolol was reduced significantly compared with cyanopindolol, acting as a very weak partial agonist of turkey β1AR and an inverse agonist of human β2AR. The structure of 7-methylcyanopindolol-bound β1AR was determined to 2.4-Å resolution and found to be virtually identical to the structure of cyanopindolol-bound β1AR. The major differences in the orthosteric binding pocket are that it has expanded by 0.3 Å in 7-methylcyanopindolol-bound β1AR and the hydroxyl group of Ser(5.46) is positioned 0.8 Å further from the ligand, with respect to the position of the Ser(5.46) side chain in cyanopindolol-bound β1AR. Thus, the molecular basis for the reduction in efficacy of 7-methylcyanopindolol compared with cyanopindolol may be regarded as the opposite of the mechanism proposed for the increase in efficacy of agonists compared with antagonists.

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.

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.

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.

Modeling G protein-coupled receptors in complex with biased agonists

The biological response to the activation of G protein-coupled receptors (GPCRs) typically originates from the simultaneous modulation of various signaling pathways that lead to distinct biological consequences. Hence, 'biased agonists' (i.e., compounds that selectively activate one of the pathways while blocking the others) are highly sought-after molecules to provide fine-tuned pharmacological interventions. This review describes strategies that can be deployed to model the conformation of GPCRs in complex with ligands endowed with specific signaling profiles useful for the generation of hypotheses on the structural requirements for the activation of different signaling pathways or for rational computer-aided ligand discovery campaigns. In particular, it focuses on strategies potentially applicable to model the global or local conformational states of GPCRs stabilized by specific ligands.

The 2.1 Å resolution structure of cyanopindolol-bound β1-adrenoceptor identifies an intramembrane Na+ ion that stabilises the ligand-free receptor

The β₁-adrenoceptor (β₁AR) is a G protein-coupled receptor (GPCR) that is activated by the endogenous agonists adrenaline and noradrenaline. We have determined the structure of an ultra-thermostable β₁AR mutant bound to the weak partial agonist cyanopindolol to 2.1 Å resolution. High-quality crystals (100 μm plates) were grown in lipidic cubic phase without the assistance of a T4 lysozyme or BRIL fusion in cytoplasmic loop 3, which is commonly employed for GPCR crystallisation. An intramembrane Na⁺ ion was identified co-ordinated to Asp87^2.50, Ser128^3.39 and 3 water molecules, which is part of a more extensive network of water molecules in a cavity formed between transmembrane helices 1, 2, 3, 6 and 7. Remarkably, this water network and Na⁺ ion is highly conserved between β₁AR and the adenosine A_2A receptor (rmsd of 0.3 Å), despite an overall rmsd of 2.4 Å for all Cα atoms and only 23% amino acid identity in the transmembrane regions. The affinity of agonist binding and nanobody Nb80 binding to β₁AR is unaffected by Na⁺ ions, but the stability of the receptor is decreased by 7.5°C in the absence of Na⁺. Mutation of amino acid side chains that are involved in the co-ordination of either Na⁺ or water molecules in the network decreases the stability of β₁AR by 5–10°C. The data suggest that the intramembrane Na⁺ and associated water network stabilise the ligand-free state of β₁AR, but still permits the receptor to form the activated state which involves the collapse of the Na⁺ binding pocket on agonist binding.

Unifying family A GPCR theories of activation

Several new pairs of active and inactive GPCR structures have recently been solved enabling detailed structural insight into the activation process, not only of rhodopsin but now also of the β2 adrenergic, M2 muscarinic and adenosine A2A receptors. Combined with structural analyses they have enabled us to examine the different recent theories proposed for GPCR activation and show that they are all indeed parts of the same process, and are intrinsically related through their effect on the central hydrophobic core of GPCRs. This new unifying general process of activation is consistent with the identification of known constitutively active mutants and an in-depth conservational analysis of significant residues implicated in the process.

Molecular pharmacology and ligand docking studies reveal a single amino acid difference between mouse and human serotonin 5-HT2A receptors that impacts behavioral translation of novel 4-phenyl-2-dimethylaminotetralin ligands

During translational studies to develop 4-phenyl-2-dimethylaminotetralin (PAT) compounds for neuropsychiatric disorders, the (2R,4S)-trans-(+)- and (2S,4R)-trans-(-)-enantiomers of the analog 6-hydroxy-7-chloro-PAT (6-OH-7-Cl-PAT) demonstrated unusual pharmacology at serotonin (5-HT) 5-HT2 G protein-coupled receptors (GPCRs). The enantiomers had similar affinities (Ki) at human (h) 5-HT2A receptors (≈ 70 nM). In an in vivo mouse model of 5-HT2A receptor activation [(±)-(2,5)-dimethoxy-4-iodoamphetamine (DOI)-elicited head twitch], however, (-)-6-OH-7-Cl-PAT was about 5-fold more potent than the (+)-enantiomer at attenuating the DOI-elicited response. It was discovered that (+)-6-OH-7-Cl-PAT (only) had ≈ 40-fold-lower affinity at mouse (m) compared with h5-HT2A receptors. Molecular modeling and computational ligand docking studies indicated that the 6-OH moiety of (+)- but not (-)-6-OH-7-Cl-PAT could form a hydrogen bond with serine residue 5.46 of the h5-HT2A receptor. The m5-HT2A as well as m5-HT2B, h5-HT2B, m5-HT2C, and h5-HT2C receptors have alanine at position 5.46, obviating this interaction; (+)-6-OH-7-Cl-PAT also showed ≈ 50-fold lower affinity than (-)-6-OH-7-Cl-PAT at m5-HT2C and h5-HT2C receptors. Mutagenesis studies confirmed that 5-HT2A S5.46 is critical for (+)- but not (-)-6-OH-7-Cl-PAT binding, as well as function. The (+)-6-OH-7-Cl-PAT enantiomer showed partial agonist effects at h5-HT2A wild-type (WT) and m5-HT2A A5.46S point-mutated receptors but did not activate m5-HT2A WT and h5-HT2A S5.46A point-mutated receptors, or h5-HT2B, h5-HT2C, and m5-HT2C receptors; (-)-6-OH-7-Cl-PAT did not activate any of the 5-HT2 receptors. Experiments also included the (2R,4S)-trans-(+)- and (2S,4R)-trans-(-)-enantiomers of 6-methoxy-7-chloro-PAT to validate hydrogen bonding interactions proposed for the corresponding 6-OH analogs. Results indicate that PAT ligand three-dimensional structure impacts target receptor binding and translational outcomes, supporting the hypothesis that GPCR ligand structure governs orthosteric binding pocket molecular determinants and resulting pharmacology.

The role of ECL2 in CGRP receptor activation: a combined modelling and experimental approach

The calcitonin gene-related peptide (CGRP) receptor is a complex of a calcitonin receptor-like receptor (CLR), which is a family B G-protein-coupled receptor (GPCR) and receptor activity modifying protein 1. The role of the second extracellular loop (ECL2) of CLR in binding CGRP and coupling to Gs was investigated using a combination of mutagenesis and modelling. An alanine scan of residues 271–294 of CLR showed that the ability of CGRP to produce cAMP was impaired by point mutations at 13 residues; most of these also impaired the response to adrenomedullin (AM). These data were used to select probable ECL2-modelled conformations that are involved in agonist binding, allowing the identification of the likely contacts between the peptide and receptor. The implications of the most likely structures for receptor activation are discussed.

Biophysical fragment screening of the β1-adrenergic receptor: identification of high affinity arylpiperazine leads using structure-based drug design

Biophysical fragment screening of a thermostabilized β1-adrenergic receptor (β1AR) using surface plasmon resonance (SPR) enabled the identification of moderate affinity, high ligand efficiency (LE) arylpiperazine hits 7 and 8. Subsequent hit to lead follow-up confirmed the activity of the chemotype, and a structure-based design approach using protein-ligand crystal structures of the β1AR resulted in the identification of several fragments that bound with higher affinity, including indole 19 and quinoline 20. In the first example of GPCR crystallography with ligands derived from fragment screening, structures of the stabilized β1AR complexed with 19 and 20 were determined at resolutions of 2.8 and 2.7 Å, respectively.

G-Protein-Coupled Receptors: from Structural Insights to Functional Mechanisms

The papers resulting from the recent Biochemical Society Focused Meeting ‘G-Protein-Coupled Receptors: from Structural Insights to Functional Mechanisms’ held in Prato in September 2012 are introduced in the present overview. A number of future goals for GPCR (G-protein-coupled receptor) research are considered, including the need to develop biophysical and computational methods to explore the full range of GPCR conformations and their dynamics, the need to develop methods to take this into account for drug discovery and the importance of relating observations on isolated receptors or receptors expressed in model systems to receptor function in vivo.