Ligand binding and micro-switches in 7TM receptor structures

Primary cilia, A-kinase anchoring proteins and constitutive activity at the orphan G protein-coupled receptor GPR161: A tale about a tail

Primary cilia are non-motile antennae-like structures responsible for sensing environmental changes in most mammalian cells. Ciliary signalling is largely mediated by the Sonic Hedgehog (Shh) pathway, which acts as a master regulator of ciliary protein transit and is essential for normal embryonic development. One particularly important player in primary cilia is the orphan G protein-coupled receptor, GPR161. In this review, we introduce GPR161 in the context of Shh signalling and describe the unique features on its C-terminus such as PKA phosphorylation sites and an A-kinase anchoring protein motif, which may influence the function of the receptor, cAMP compartmentalisation and/or trafficking within primary cilia. We discuss the recent putative pairing of GPR161 and spexin-1, highlighting the additional steps needed before GPR161 could be considered 'deorphanised'. Finally, we speculate that the marked constitutive activity and unconventional regulation of GPR161 may indicate that the receptor may not require an endogenous ligand. LINKED ARTICLES: This article is part of a themed issue Therapeutic Targeting of G Protein-Coupled Receptors: hot topics from the Australasian Society of Clinical and Experimental Pharmacologists and Toxicologists 2021 Virtual Annual Scientific Meeting. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.14/issuetoc.

Molecular Basis of MC1R Activation: Mutation-Induced Alterations in Structural Dynamics

The MC1R protein is a receptor found in melanocytes that plays a role in melanin synthesis. Mutations in this protein can impact hair color, skin tone, tanning ability, and increase the risk of skin cancer. The MC1R protein is activated by the alpha-melanocyte-stimulating hormone (α-MSH). Previous studies have shown that mutations affect the interaction between MC1R and α-MSH; however, the mechanism behind this process is poorly understood. Our study aims to shed light on this mechanism using molecular dynamics (MD) simulations to analyze the Asp84Glu and Asp294His variants. We simulated both the wild-type (WT) protein and the mutants with and without ligand. Our results reveal that mutations induce unique conformations during state transitions, hindering the switch between active and inactive states and decreasing cellular levels of cAMP. Interestingly, Asp294His showed increased ligand affinity but decreased protein activity, highlighting that tighter binding does not always lead to increased activation. Our study provides insights into the molecular mechanisms underlying the impact of MC1R mutations on protein activity.

Structural Insights into Partial Activation of the Prototypic G Protein-Coupled Adenosine A2A Receptor

The adenosine A2A receptor (A2AAR) belongs to the rhodopsin-like G protein-coupled receptor (GPCR) family, which constitutes the largest class of GPCRs. Partial agonists show reduced efficacy as compared to physiological agonists and can even act as antagonists in the presence of a full agonist. Here, we determined an X-ray crystal structure of the partial A2AAR agonist 2-amino-6-[(1H-imidazol-2-ylmethyl)sulfanyl]-4-p-hydroxyphenyl-3,5-pyridinedicarbonitrile (LUF5834) in complex with the A2AAR construct A2A-PSB2-bRIL, stabilized in its inactive conformation and being devoid of any mutations in the ligand binding pocket. The determined high-resolution structure (2.43 Å) resolved water networks and crucial binding pocket interactions. A direct hydrogen bond of the p-hydroxy group of LUF5834 with T883.36 was observed, an amino acid that was mutated to alanine in the most frequently used A2AAR crystallization constructs thus preventing the discovery of its interactions in most of the previous A2AAR co-crystal structures. G protein dissociation studies confirmed partial agonistic activity of LUF5834 as compared to that of the full agonist N-ethylcarboxamidoadenosine (NECA). In contrast to NECA, the partial agonist was still able to bind to the receptor construct locked in its inactive conformation by an S913.39K mutation, although with an affinity lower than that at the native receptor. This could explain the compound's partial agonistic activity: while full A2AAR agonists bind exclusively to the active conformation, likely following conformational selection, partial agonists bind to active as well as inactive conformations, showing higher affinity for the active conformation. This might be a general mechanism of partial agonism also applicable to other GPCRs.

GPER binding site detection and description: A flavonoid-based docking and molecular dynamics simulations study

Flavonoids, a phenolic compounds class widely distributed in the plant kingdom, have attracted much interest for their implications on several health and disease processes. Usually, the consumption of this type of compounds is approximately 1 g/d, primarily obtained from cereals, chocolate, and dry legumes ensuring its beneficial role in maintaining the homeostasis of the human body. In this context, in cancer disease prominent data points to the role of flavonoids as adjuvant treatment aimed at the regression of the disease. GPER, an estrogen receptor on the cell surface, has been postulated as a probable orchestrator of the beneficial effects of several flavonoids through modulation/inhibition of various mechanisms that lead to cancer progression. Therefore, applying pocket and cavity protein detection and docking and molecular dynamics simulations (MD), we generate, from a cluster composed of 39 flavonoids, crucial insights into the potential role as GPER ligands, of Puerarin, Isoquercetin, Kaempferol 3-O-glucoside and Petunidin 3-O-glucoside, aglycones whose sugar moiety delimits a new described sugar-acceptor sub-cavity into the cavity binding site on the receptor, as well as of the probable activation mechanism of the receptor and the pivotal residues involved in it. Altogether, our results shed light on the potential use of the aforementioned flavonoids as GPER ligands and for further evaluations in in vitro and in vivo assays to elucidate their probable anti-cancer activity.

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.

Molecular dynamics-based identification of binding pathways and two distinct high-affinity sites for succinate in succinate receptor 1/GPR91

SUCNR1 is an auto- and paracrine sensor of the metabolic stress signal succinate. Using unsupervised molecular dynamics (MD) simulations (170.400 ns) and mutagenesis across human, mouse, and rat SUCNR1, we characterize how a five-arginine motif around the extracellular pole of TM-VI determines the initial capture of succinate in the extracellular vestibule (ECV) to either stay or move down to the orthosteric site. Metadynamics demonstrate low-energy succinate binding in both sites, with an energy barrier corresponding to an intermediate stage during which succinate, with an associated water cluster, unlocks the hydrogen-bond-stabilized conformationally constrained extracellular loop (ECL)-2b. Importantly, simultaneous binding of two succinate molecules through either a "sequential" or "bypassing" mode is a frequent endpoint. The mono-carboxylate NF-56-EJ40 antagonist enters SUCNR1 between TM-I and -II and does not unlock ECL-2b. It is proposed that occupancy of both high-affinity sites is required for selective activation of SUCNR1 by high local succinate concentrations.

Characterizing conformational states in GPCR structures using machine learning

G protein-coupled receptors (GPCRs) play a pivotal role in signal transduction and represent attractive targets for drug development. Recent advances in structural biology have provided insights into GPCR conformational states, which are critical for understanding their signaling pathways and facilitating structure-based drug discovery. In this study, we introduce a machine learning approach for conformational state annotation of GPCRs. We represent GPCR conformations as high-dimensional feature vectors, incorporating information about amino acid residue pairs involved in the activation pathway. Using a dataset of GPCR conformations in inactive and active states obtained through molecular dynamics simulations, we trained machine learning models to distinguish between inactive-like and active-like conformations. The developed model provides interpretable predictions and can be used for the large-scale analysis of molecular dynamics trajectories of GPCRs.

Allosteric modulation of G protein-coupled receptors as a novel therapeutic strategy in neuropathic pain

Neuropathic pain is a debilitating pathological condition that presents significant therapeutic challenges in clinical practice. Unfortunately, current pharmacological treatments for neuropathic pain lack clinical efficacy and often lead to harmful adverse reactions. As G protein-coupled receptors (GPCRs) are widely distributed throughout the body, including the pain transmission pathway and descending inhibition pathway, the development of novel neuropathic pain treatments based on GPCRs allosteric modulation theory is gaining momentum. Extensive research has shown that allosteric modulators targeting GPCRs on the pain pathway can effectively alleviate symptoms of neuropathic pain while reducing or eliminating adverse effects. This review aims to provide a comprehensive summary of the progress made in GPCRs allosteric modulators in the treatment of neuropathic pain, and discuss the potential benefits and adverse factors of this treatment. We will also concentrate on the development of biased agonists of GPCRs, and based on important examples of biased agonist development in recent years, we will describe universal strategies for designing structure-based biased agonists. It is foreseeable that, with the continuous improvement of GPCRs allosteric modulation and biased agonist theory, effective GPCRs allosteric drugs will eventually be available for the treatment of neuropathic pain with acceptable safety.

Structure, function and drug discovery of GPCR signaling

G protein-coupled receptors (GPCRs) are versatile and vital proteins involved in a wide array of physiological processes and responses, such as sensory perception (e.g., vision, taste, and smell), immune response, hormone regulation, and neurotransmission. Their diverse and essential roles in the body make them a significant focus for pharmaceutical research and drug development. Currently, approximately 35% of marketed drugs directly target GPCRs, underscoring their prominence as therapeutic targets. Recent advances in structural biology have substantially deepened our understanding of GPCR activation mechanisms and interactions with G-protein and arrestin signaling pathways. This review offers an in-depth exploration of both traditional and recent methods in GPCR structure analysis. It presents structure-based insights into ligand recognition and receptor activation mechanisms and delves deeper into the mechanisms of canonical and noncanonical signaling pathways downstream of GPCRs. Furthermore, it highlights recent advancements in GPCR-related drug discovery and development. Particular emphasis is placed on GPCR selective drugs, allosteric and biased signaling, polyphamarcology, and antibody drugs. Our goal is to provide researchers with a thorough and updated understanding of GPCR structure determination, signaling pathway investigation, and drug development. This foundation aims to propel forward-thinking therapeutic approaches that target GPCRs, drawing upon the latest insights into GPCR ligand selectivity, activation, and biased signaling mechanisms.

Structural basis of antibody inhibition and chemokine activation of the human CC chemokine receptor 8

The C-C motif chemokine receptor 8 (CCR8) is a class A G-protein coupled receptor that has emerged as a promising therapeutic target in cancer. Targeting CCR8 with an antibody has appeared to be an attractive therapeutic approach, but the molecular basis for chemokine-mediated activation and antibody-mediated inhibition of CCR8 are not fully elucidated. Here, we obtain an antagonist antibody against human CCR8 and determine structures of CCR8 in complex with either the antibody or the endogenous agonist ligand CCL1. Our studies reveal characteristic antibody features allowing recognition of the CCR8 extracellular loops and CCL1-CCR8 interaction modes that are distinct from other chemokine receptor - ligand pairs. Informed by these structural insights, we demonstrate that CCL1 follows a two-step, two-site binding sequence to CCR8 and that antibody-mediated inhibition of CCL1 signaling can occur by preventing the second binding event. Together, our results provide a detailed structural and mechanistic framework of CCR8 activation and inhibition that expands our molecular understanding of chemokine - receptor interactions and offers insight into the development of therapeutic antibodies targeting chemokine GPCRs.

Structure-Based Discovery of Mouse Trace Amine-Associated Receptor 5 Antagonists

Trace amine-associated receptors (TAARs) were discovered in 2001 as new members of class A G protein-coupled receptors (GPCRs). With the only exception of TAAR1, TAAR members (TAAR2-9, also known as noncanonical olfactory receptors) were originally described exclusively in the olfactory epithelium and believed to mediate the innate perception of volatile amines. However, most noncanonical olfactory receptors are still orphan receptors. Given its recently discovered nonolfactory expression and therapeutic potential, TAAR5 has been the focus of deorphanization campaigns that led to the discovery of a few druglike antagonists. Here, we report four novel TAAR5 antagonists identified through high-throughput screening, which, along with the four ligands published in the literature, constituted our starting point to design a computational strategy for the identification of TAAR5 ligands. We developed a structure-based virtual screening protocol that allowed us to identify three new TAAR5 antagonists with a hit rate of 10%. Despite lacking an experimental structure, we accurately modeled the TAAR5 binding site by integrating comparative sequence- and structure-based analyses of serotonin receptors with homology modeling and side-chain optimization. In summary, we have identified seven new TAAR5 antagonists that could serve as lead candidates for the development of new treatments for depression, anxiety, and neurodegenerative diseases.

Activation/Deactivation Free-Energy Profiles for the β2-Adrenergic Receptor: Ligand Modes of Action

We use enhanced-sampling simulations with an effective collective variable to study the activation of the β2-adrenergic receptor in the presence of ligands with different efficacy. The free-energy profiles are computed for the ligand-free (apo) receptor and binary (apo-receptor + G-protein α-subunit and receptor + ligand) and ternary complexes. The results are not only compatible with available experiments but also allow unprecedented structural insight into the nature of GPCR conformations along the activation pathway and their role in the activation mechanism. In particular, the simulations reveal an unexpected mode of action of partial agonists such as salmeterol and salbutamol that arises already in the binary complex without the G-protein. Specific differences in the polar interactions with residues in TM5, which are required to stabilize an optimal TM6 conformation that facilitates G-protein binding and receptor activation, play a major role in differentiating them from full agonists.

Pharmacological hallmarks of allostery at the M4 muscarinic receptor elucidated through structure and dynamics

Allosteric modulation of G protein-coupled receptors (GPCRs) is a major paradigm in drug discovery. Despite decades of research, a molecular-level understanding of the general principles that govern the myriad pharmacological effects exerted by GPCR allosteric modulators remains limited. The M₄ muscarinic acetylcholine receptor (M₄ mAChR) is a validated and clinically relevant allosteric drug target for several major psychiatric and cognitive disorders. In this study, we rigorously quantified the affinity, efficacy, and magnitude of modulation of two different positive allosteric modulators, LY2033298 (LY298) and VU0467154 (VU154), combined with the endogenous agonist acetylcholine (ACh) or the high-affinity agonist iperoxo (Ipx), at the human M₄ mAChR. By determining the cryo-electron microscopy structures of the M₄ mAChR, bound to a cognate G_i1 protein and in complex with ACh, Ipx, LY298-Ipx, and VU154-Ipx, and applying molecular dynamics simulations, we determine key molecular mechanisms underlying allosteric pharmacology. In addition to delineating the contribution of spatially distinct binding sites on observed pharmacology, our findings also revealed a vital role for orthosteric and allosteric ligand-receptor-transducer complex stability, mediated by conformational dynamics between these sites, in the ultimate determination of affinity, efficacy, cooperativity, probe dependence, and species variability. There results provide a holistic framework for further GPCR mechanistic studies and can aid in the discovery and design of future allosteric drugs.

Small molecule-mediated targeting of microRNAs for drug discovery: Experiments, computational techniques, and disease implications

Small molecules have been providing medical breakthroughs for human diseases for more than a century. Recently, identifying small molecule inhibitors that target microRNAs (miRNAs) has gained importance, despite the challenges posed by labour-intensive screening experiments and the significant efforts required for medicinal chemistry optimization. Numerous experimentally-verified cases have demonstrated the potential of miRNA-targeted small molecule inhibitors for disease treatment. This new approach is grounded in their posttranscriptional regulation of the expression of disease-associated genes. Reversing dysregulated gene expression using this mechanism may help control dysfunctional pathways. Furthermore, the ongoing improvement of algorithms has allowed for the integration of computational strategies built on top of laboratory-based data, facilitating a more precise and rational design and discovery of lead compounds. To complement the use of extensive pharmacogenomics data in prioritising potential drugs, our previous work introduced a computational approach based on only molecular sequences. Moreover, various computational tools for predicting molecular interactions in biological networks using similarity-based inference techniques have been accumulated in established studies. However, there are a limited number of comprehensive reviews covering both computational and experimental drug discovery processes. In this review, we outline a cohesive overview of both biological and computational applications in miRNA-targeted drug discovery, along with their disease implications and clinical significance. Finally, utilizing drug-target interaction (DTIs) data from DrugBank, we showcase the effectiveness of deep learning for obtaining the physicochemical characterization of DTIs.

General Metadynamics Protocol To Simulate Activation/Deactivation of Class A GPCRs: Proof of Principle for the Serotonin Receptor

We present a generally applicable metadynamics protocol for characterizing the activation free-energy profiles of class A G-protein coupled receptors and a proof-of-principle study for the 5HT_1A-receptor. The almost universal A¹⁰⁰ activation index, which depends on five inter-helix distances, is used as the single collective variable in well-tempered multiple-walker metadynamics simulations. Here, we show free-energy profiles for the serotonin receptor as binary (apo-receptor + G-protein-α-subunit and receptor + ligand) and ternary complexes with two prototypical orthosteric ligands: the full agonist serotonin and the partial agonist aripiprazole. Our results are not only compatible with previously reported experimental and computational data, but they also allow differences between active and inactive conformations to be determined in unprecedented atomic detail, and with respect to the so-called microswitches that have been suggested as determinants of activation, giving insight into their role in the activation mechanism.

Prospects for targeting ACKR1 in cancer and other diseases

The chemokine network is comprised of a family of signal proteins that encode messages for cells displaying chemokine G-protein coupled receptors (GPCRs). The diversity of effects on cellular functions, particularly directed migration of different cell types to sites of inflammation, is enabled by different combinations of chemokines activating signal transduction cascades on cells displaying a combination of receptors. These signals can contribute to autoimmune disease or be hijacked in cancer to stimulate cancer progression and metastatic migration. Thus far, three chemokine receptor-targeting drugs have been approved for clinical use: Maraviroc for HIV, Plerixafor for hematopoietic stem cell mobilization, and Mogalizumab for cutaneous T-cell lymphoma. Numerous compounds have been developed to inhibit specific chemokine GPCRs, but the complexity of the chemokine network has precluded more widespread clinical implementation, particularly as anti-neoplastic and anti-metastatic agents. Drugs that block a single signaling axis may be rendered ineffective or cause adverse reactions because each chemokine and receptor often have multiple context-specific functions. The chemokine network is tightly regulated at multiple levels, including by atypical chemokine receptors (ACKRs) that control chemokine gradients independently of G-proteins. ACKRs have numerous functions linked to chemokine immobilization, movement through and within cells, and recruitment of alternate effectors like β-arrestins. Atypical chemokine receptor 1 (ACKR1), previously known as the Duffy antigen receptor for chemokines (DARC), is a key regulator that binds chemokines involved in inflammatory responses and cancer proliferation, angiogenesis, and metastasis. Understanding more about ACKR1 in different diseases and populations may contribute to the development of therapeutic strategies targeting the chemokine network.

Structural Basis for Agonistic Activity and Selectivity toward Melatonin Receptors hMT1 and hMT2

Glaucoma, a major ocular neuropathy originating from a progressive degeneration of retinal ganglion cells, is often associated with increased intraocular pressure (IOP). Daily IOP fluctuations are physiologically influenced by the antioxidant and signaling activities of melatonin. This endogenous modulator has limited employment in treating altered IOP disorders due to its low stability and bioavailability. The search for low-toxic compounds as potential melatonin agonists with higher stability and bioavailability than melatonin itself could start only from knowing the molecular basis of melatonergic activity. Thus, using a computational approach, we studied the melatonin binding toward its natural macromolecular targets, namely melatonin receptors 1 (MT1) and 2 (MT2), both involved in IOP signaling regulation. Besides, agomelatine, a melatonin-derivative agonist and, at the same time, an atypical antidepressant, was also included in the study due to its powerful IOP-lowering effects. For both ligands, we evaluated both stability and ligand positioning inside the orthosteric site of MTs, mapping the main molecular interactions responsible for receptor activation. Affinity values in terms of free binding energy (ΔG_bind) were calculated for the selected poses of the chosen compounds after stabilization through a dynamic molecular docking protocol. The results were compared with experimental in vivo effects, showing a higher potency and more durable effect for agomelatine with respect to melatonin, which could be ascribed both to its higher affinity for hMT2 and to its additional activity as an antagonist for the serotonin receptor 5-HT2c, in agreement with the in silico results.

Structure-function relationship and physiological role of apelin and its G protein coupled receptor

Apelin receptor (APJR) is a class A peptide (apelin) binding G protein-coupled receptor (GPCR) that plays a significant role in regulating blood pressure, cardiac output, and maintenance of fluid homeostasis. It is activated by a wide range of endogenous peptide isoforms of apelin and elabela. The apelin peptide isoforms contain distinct structural features that aid in ligand recognition and activation of the receptor. Site-directed mutagenesis and structure-based studies have revealed the involvement of extracellular and transmembrane regions of the receptor in binding to the peptide isoforms. The structural features of APJR activation of the receptor as well as mediating G-protein and β-arrestin-mediated signaling are delineated by multiple mutagenesis studies. There is increasing evidence that the structural requirements of APJR to activate G-proteins and β-arrestins are different, leading to biased signaling. APJR also responds to mechanical stimuli in a ligand-independent manner. A multitude of studies has focused on developing both peptide and non-peptide agonists and antagonists specific to APJR. Apelin/elabela-activated APJR orchestrates major signaling pathways such as extracellular signal-regulated kinase (ERKs), protein kinase B (PKB/Akt), and p70S. This review focuses on the structural and functional characteristics of apelin, elabela, APJR, and their interactions involved in the binding and activation of the downstream signaling cascade. We also focus on the diverse signaling profile of APJR and its ligands and their involvement in various physiological systems.

Development of Pharmacophore Models for the Important Off-Target 5-HT2B Receptor

Toxicity is a major cause of attrition in the development of pharmaceuticals, and the off-target effects are a frequent contributor. The 5-HT_2B receptor agonism is known to be responsible for a variety of safety concerns including valvular heart disease which was the cause for the withdrawal of several compounds from the market. An early detection of potential binding to this receptor is thus desirable. Herein, we present the identification of key amino acid residues in the active site of 5-HT_2B by molecular dynamics simulations, the development of pharmacophore models and their performance on in-house data, and a structurally highly diverse subset of Enamine REAL labeled for 5-HT_2B activity by a machine learning model. These models may be used as filters employed on screening compound sets for the early filtration of compounds with potential 5-HT_2B off-target liabilities.

MUG: A mutation overview of GPCR subfamily A17 receptors

G protein-coupled receptors (GPCRs) mediate several signaling pathways through a general mechanism that involves their activation, upholding a chain of events that lead to the release of molecules responsible for cytoplasmic action and further regulation. These physiological functions can be severely altered by mutations in GPCR genes. GPCRs subfamily A17 (dopamine, serotonin, adrenergic and trace amine receptors) are directly related with neurodegenerative diseases, and as such it is crucial to explore known mutations on these systems and their impact in structure and function. A comprehensive and detailed computational framework - MUG (Mutations Understanding GPCRs) - was constructed, illustrating key reported mutations and their effect on receptors of the subfamily A17 of GPCRs. We explored the type of mutations occurring overall and in the different families of subfamily A17, as well their localization within the receptor and potential effects on receptor functionality. The mutated residues were further analyzed considering their pathogenicity. The results reveal a high diversity of mutations in the GPCR subfamily A17 structures, drawing attention to the considerable number of mutations in conserved residues and domains. Mutated residues were typically hydrophobic residues enriched at the ligand binding pocket and known activating microdomains, which may lead to disruption of receptor function. MUG as an interactive web application is available for the management and visualization of this dataset. We expect that this interactive database helps the exploration of GPCR mutations, their influence, and their familywise and receptor-specific effects, constituting the first step in elucidating their structures and molecules at the atomic level.

Superconserved receptors expressed in the brain: Expression, function, motifs and evolution of an orphan receptor family

GPR27, GPR85 and GPR173 constitute a small family of G protein-coupled receptors (GPCR) that share the distinctive characteristics of being highly conserved throughout vertebrate evolution and predominantly expressed in the brain. Accordingly, they have been coined as "Superconserved Receptors Expressed in the Brain" (SREB), although their expression profile is more complex than what was originally thought. SREBs have no known validated endogenous ligands and are thus labeled as "orphan" receptors. The investigation of this particular category of uncharacterized receptors holds great promise both in terms of physiology and drug development. In the largest GPCR family, the Rhodopsin-like or Class A, around 100 receptors are considered orphans. Because GPCRs are the most successful source of drug targets, the discovery of a novel function or ligand most likely will lead to significant breakthroughs for the discovery of innovative therapies. The high level of conservation is one of the characteristic features of the SREBs. We propose herein a detailed analysis of the putative evolutionary origin of this family. We highlight the properties that distinguish SREBs from other rhodopsin-like GPCRs. We present the current evidence for these receptors downstream signaling pathways and functions. We discuss the pharmacological challenge for the identification of natural or synthetic ligands of orphan receptors like SREBs. The different SREB-related scientific questions are presented with a highlight on what should be addressed in the near future, including the confirmation of published evidence and their validation as drug targets. In particular, we discuss in which pathological conditions these receptors may be of great relevance to solve unmet medical needs.

Class A and C GPCR Dimers in Neurodegenerative Diseases

Neurodegenerative diseases affect over 30 million people worldwide with an ascending trend. Most individuals suffering from these irreversible brain damages belong to the elderly population, with onset between 50 and 60 years. Although the pathophysiology of such diseases is partially known, it remains unclear upon which point a disease turns degenerative. Moreover, current therapeutics can treat some of the symptoms but often have severe side effects and become less effective in long-term treatment. For many neurodegenerative diseases, the involvement of G proteincoupled receptors (GPCRs), which are key players of neuronal transmission and plasticity, has become clearer and holds great promise in elucidating their biological mechanism. With this review, we introduce and summarize class A and class C GPCRs, known to form heterodimers or oligomers to increase their signalling repertoire. Additionally, the examples discussed here were shown to display relevant alterations in brain signalling and had already been associated with the pathophysiology of certain neurodegenerative diseases. Lastly, we classified the heterodimers into two categories of crosstalk, positive or negative, for which there is known evidence.

Structural basis for receptor selectivity and inverse agonism in S1P5 receptors

The bioactive lysophospholipid sphingosine-1-phosphate (S1P) acts via five different subtypes of S1P receptors (S1PRs) - S1P_1-5. S1P₅ is predominantly expressed in nervous and immune systems, regulating the egress of natural killer cells from lymph nodes and playing a role in immune and neurodegenerative disorders, as well as carcinogenesis. Several S1PR therapeutic drugs have been developed to treat these diseases; however, they lack receptor subtype selectivity, which leads to side effects. In this article, we describe a 2.2 Å resolution room temperature crystal structure of the human S1P₅ receptor in complex with a selective inverse agonist determined by serial femtosecond crystallography (SFX) at the Pohang Accelerator Laboratory X-Ray Free Electron Laser (PAL-XFEL) and analyze its structure-activity relationship data. The structure demonstrates a unique ligand-binding mode, involving an allosteric sub-pocket, which clarifies the receptor subtype selectivity and provides a template for structure-based drug design. Together with previously published S1PR structures in complex with antagonists and agonists, our structure with S1P₅-inverse agonist sheds light on the activation mechanism and reveals structural determinants of the inverse agonism in the S1PR family.

S1P₅ is a sphingosine-1-phosphate (S1P) receptor implicated in immune and neurodegenerative disorders. Here, authors report a crystal structure of the S1P₅ receptor in complex with a selective inverse agonist, revealing an allosteric subpocket and shedding light on inverse agonism in S1P receptors.

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

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

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.

Recent Molecular Insights into Agonist-specific Binding to the Mu-Opioid Receptor

Opioid agonists produce their analgesic effects primarily by acting at the µ-opioid receptor (µOR). µOR agonists with different efficacies exert diverse molecular changes in the µOR which dictate the faith of the receptor’s signaling pathway and possibly it’s the degree of desensitization. Since the development of the active conformations of the µOR, growing data have been published in relation to ligand-specific changes in µOR activation. In this regard, this review summarizes recent data regarding the most studied opioid agonists in in silico µOR activation, including how these ligands are recognized by the µOR, how their binding signal is transmitted toward the intracellular parts of the µOR, and finally, what type of large-scale movements do these changes trigger in the µOR’s domains.

Angiotensin and Endothelin Receptor Structures With Implications for Signaling Regulation and Pharmacological Targeting

In conjunction with the endothelin (ET) type A (ETAR) and type B (ETBR) receptors, angiotensin (AT) type 1 (AT1R) and type 2 (AT2R) receptors, are peptide-binding class A G-protein-coupled receptors (GPCRs) acting in a physiologically overlapping context. Angiotensin receptors (ATRs) are involved in regulating cell proliferation, as well as cardiovascular, renal, neurological, and endothelial functions. They are important therapeutic targets for several diseases or pathological conditions, such as hypertrophy, vascular inflammation, atherosclerosis, angiogenesis, and cancer. Endothelin receptors (ETRs) are expressed primarily in blood vessels, but also in the central nervous system or epithelial cells. They regulate blood pressure and cardiovascular homeostasis. Pathogenic conditions associated with ETR dysfunctions include cancer and pulmonary hypertension. While both receptor groups are activated by their respective peptide agonists, pathogenic autoantibodies (auto-Abs) can also activate the AT1R and ETAR accompanied by respective clinical conditions. To date, the exact mechanisms and differences in binding and receptor-activation mediated by auto-Abs as opposed to endogenous ligands are not well understood. Further, several questions regarding signaling regulation in these receptors remain open. In the last decade, several receptor structures in the apo- and ligand-bound states were determined with protein X-ray crystallography using conventional synchrotrons or X-ray Free-Electron Lasers (XFEL). These inactive and active complexes provide detailed information on ligand binding, signal induction or inhibition, as well as signal transduction, which is fundamental for understanding properties of different activity states. They are also supportive in the development of pharmacological strategies against dysfunctions at the receptors or in the associated signaling axis. Here, we summarize current structural information for the AT1R, AT2R, and ETBR to provide an improved molecular understanding.

The Many Faces of G Protein-Coupled Receptor 143, an Atypical Intracellular Receptor

GPCRs transform extracellular stimuli into a physiological response by activating an intracellular signaling cascade initiated via binding to G proteins. Orphan G protein-coupled receptors (GPCRs) hold the potential to pave the way for development of new, innovative therapeutic strategies. In this review we will introduce G protein-coupled receptor 143 (GPR143), an enigmatic receptor in terms of classification within the GPCR superfamily and localization. GPR143 has not been assigned to any of the GPCR families due to the lack of common structural motifs. Hence we will describe the most important motifs of classes A and B and compare them to the protein sequence of GPR143. While a precise function for the receptor has yet to be determined, the protein is expressed abundantly in pigment producing cells. Many GPR143 mutations cause X-linked Ocular Albinism Type 1 (OA1, Nettleship-Falls OA), which results in hypopigmentation of the eyes and loss of visual acuity due to disrupted visual system development and function. In pigment cells of the skin, loss of functional GPR143 results in abnormally large melanosomes (organelles in which pigment is produced). Studies have shown that the receptor is localized internally, including at the melanosomal membrane, where it may function to regulate melanosome size and/or facilitate protein trafficking to the melanosome through the endolysosomal system. Numerous additional roles have been proposed for GPR143 in determining cancer predisposition, regulation of blood pressure, development of macular degeneration and signaling in the brain, which we will briefly describe as well as potential ligands that have been identified. Furthermore, GPR143 is a promiscuous receptor that has been shown to interact with multiple other melanosomal proteins and GPCRs, which strongly suggests that this orphan receptor is likely involved in many different physiological actions.

Identification of a conserved chemokine receptor motif that enables ligand discrimination

Extensive ligand-receptor promiscuity in the chemokine signaling system balances beneficial redundancy and specificity. However, this feature poses a major challenge to selectively modulate the system pharmacologically. Here, we identified a conserved cluster of three aromatic receptor residues that anchors the second extracellular loop (ECL2) to the top of receptor transmembrane helices (TM) 4 and 5 and enables recognition of both shared and specific characteristics of interacting chemokines. This cluster was essential for the activation of several chemokine receptors. Furthermore, characteristic motifs of the ß₁ strand and 30s loop make the two main CC-chemokine subgroups-the macrophage inflammatory proteins (MIPs) and monocyte chemoattractant proteins (MCPs)-differentially dependent on this cluster in the promiscuous receptors CCR1, CCR2, and CCR5. The cluster additionally enabled CCR1 and CCR5 to discriminate between closely related MIPs based on the N terminus of the chemokine. G protein signaling and β-arrestin2 recruitment assays confirmed the importance of the conserved cluster in receptor discrimination of chemokine ligands. This extracellular site may facilitate the development of chemokine-related therapeutics.

Into the Tissues: Extracellular Matrix and Its Artificial Substitutes: Cell Signalling Mechanisms

The existence of orderly structures, such as tissues and organs is made possible by cell adhesion, i.e., the process by which cells attach to neighbouring cells and a supporting substance in the form of the extracellular matrix. The extracellular matrix is a three-dimensional structure composed of collagens, elastin, and various proteoglycans and glycoproteins. It is a storehouse for multiple signalling factors. Cells are informed of their correct connection to the matrix via receptors. Tissue disruption often prevents the natural reconstitution of the matrix. The use of appropriate implants is then required. This review is a compilation of crucial information on the structural and functional features of the extracellular matrix and the complex mechanisms of cell–cell connectivity. The possibilities of regenerating damaged tissues using an artificial matrix substitute are described, detailing the host response to the implant. An important issue is the surface properties of such an implant and the possibilities of their modification.

Endomembrane-Based Signaling by GPCRs and G-Proteins

G-protein-coupled receptors (GPCRs) and G-proteins have a range of roles in many physiological and pathological processes and are among the most studied signaling proteins. A plethora of extracellular stimuli can activate the GPCR and can elicit distinct intracellular responses through the activation of specific transduction pathways. For many years, biologists thought that GPCR signaling occurred entirely on the plasma membrane. However, in recent decades, many lines of evidence have proved that the GPCRs and G-proteins may reside on endomembranes and can start or propagate signaling pathways through the organelles that form the secretory route. How these alternative intracellular signaling pathways of the GPCR and G-proteins influence the physiological and pathological function of the endomembranes is still under investigation. Here, we review the general role and classification of GPCRs and G-proteins with a focus on their signaling pathways in the membrane transport apparatus.

Unique features of different classes of G-protein-coupled receptors revealed from sequence coevolutionary and structural analysis

G-protein-coupled receptors (GPCRs) are the largest family of human membrane proteins and represent the primary targets of about one third of currently marketed drugs. Despite the critical importance, experimental structures have been determined for only a limited portion of GPCRs and functional mechanisms of GPCRs remain poorly understood. Here, we have constructed novel sequence coevolutionary models of the A and B classes of GPCRs and compared them with residue contact frequency maps generated with available experimental structures. Significant portions of structural residue contacts were successfully detected in the sequence-based covariational models. "Exception" residue contacts predicted from sequence coevolutionary models but not available structures added missing links that were important for GPCR activation and allosteric modulation. Moreover, we identified distinct residue contacts involving different sets of functional motifs for GPCR activation, such as the Na⁺ pocket, CWxP, DRY, PIF, and NPxxY motifs in the class A and the HETx and PxxG motifs in the class B. Finally, we systematically uncovered critical residue contacts tuned by allosteric modulation in the two classes of GPCRs, including those from the activation motifs and particularly the extracellular and intracellular loops in class A GPCRs. These findings provide a promising framework for rational design of ligands to regulate GPCR activation and allosteric modulation.

Understanding the role of the CB1 toggle switch in interaction networks using molecular dynamics simulation

The cannabinoid receptor 1 (CB1) is a class A G-protein coupled receptor (GPCR) that can exert various effects on the human body through the endocannabinoid system. Understanding CB1 activation has many benefits for the medical use of cannabinoids. A previous study reported that CB1 has two notable residues referred to as the toggle switch, F3.36 and W6.48, which are important for its activation mechanism. We performed a molecular dynamics simulation with a mutation in the toggle switch to examine its role in active and inactive states. We also examined structural changes, the residue-residue interaction network, and the interaction network among helices and loops of wildtype and mutant CB1 for both activation states. As a result, we found that the energetic changes in the hydrogen-bond network of the Na⁺ pocket, extracellular N-terminus-TM2-ECL1-TM3 interface including D2.63-K3.28 salt-bridge, and extracellular ECL2-TM5-ECL3-TM6 interface directly linked to the toggle switch contribute to the stability of CB1 by the broken aromatic interaction of the toggle switch. It makes the conformation of inactive CB1 receptor to be unstable. Our study explained the role of the toggle switch regarding the energetic interactions related to the Na⁺ pocket and extracellular loop interfaces, which could contribute to a better understanding of the activation mechanism of CB1.

Probing Allosteric Regulation Mechanism of W7.35 on Agonist-Induced Activity for μOR by Mutation Simulation

The residue located at 15 positions before the most conserved residue in TM7 (7.35 of Ballesteros-Weinstein number) plays important roles in ligand binding and the receptor activity for class A GPCRs. Nevertheless, its regulation mechanism has not been clearly clarified in experiments, and some controversies also exist for its impact on μ-opioid receptors (μOR) bound by agonists. Thus, we chose the μ-opioid receptor (μOR) of class A GPCRs as a representative and utilized a microsecond accelerated molecular dynamics simulation (aMD) coupled with a protein structure network (PSN) to explore the effect of W318^7.35 on its functional activity induced by the agonist endomorphin2 mainly by a comparison of the wild system and its W7.35A mutant. When endomorphin2 binds to the wild-type μOR, TM6 in μOR moves outward to form an open intracellular conformation that is beneficial to accommodating the β-arrestin transducer, rather than the G-protein transducer due to the clash with the α5 helix of G-protein, thus acting as a β-arrestin biased agonist. However, the W318A mutation induces the intracellular part of μOR to form a closed state, which disfavors coupling with either G-protein or β-arrestin. The allosteric pathway analysis further reveals that the binding of endomorphin2 to the wild-type μOR transmits more activation signals to the β-arrestin binding site while the W318A mutation induces more structural signals to transmit to common binding residues of the G protein and β-arrestin. More interestingly, the residue at the 7.35 position regulates the shortest allosteric pathway in indirect ways by influencing the interactions between other ligand-binding residues and endomorphin2. W293^6.48 and F289^6.44 are important for regulating the different activities of μOR induced either by the agonist or by the mutation. Y336^7.53, F343^8.50, and D340^8.47 play crucial roles in activating the β-arrestin biased signal induced by the agonist endomorphin2, while L158^3.43 and V286^6.41 devote important contributions to the change in the activity of endomorphin2 from the β-arrestin biased agonist to the antagonist upon the W318A mutation.

Characterization of a new type of neuronal 5-HT G- protein coupled receptor in the cestode nervous system

Cestodes are platyhelminth parasites with a wide range of hosts that cause neglected diseases. Neurotransmitter signaling is of critical importance for these parasites which lack circulatory, respiratory and digestive systems. For example, serotonin (5-HT) and serotonergic G-protein coupled receptors (5-HT GPCRs) play major roles in cestode motility, development and reproduction. In previous work, we deorphanized a group of 5-HT7 type GPCRs from cestodes. However, little is known about another type of 5-HT GPCR, the 5-HT1 clade, which has been studied in several invertebrate phyla but not in platyhelminthes. Three putative 5-HT GPCRs from Echinococcus canadensis, Mesocestoides vogae (syn. M. corti) and Hymenolepis microstoma were cloned, sequenced and bioinformatically analyzed. Evidence grouped these new sequences within the 5-HT1 clade of GPCRs but differences in highly conserved GPCR motifs were observed. Transcriptomic analysis, heterologous expression and immunolocalization studies were performed to characterize the E. canadensis receptor, called Eca-5-HT_1a. Functional heterologous expression studies showed that Eca-5-HT_1a is highly specific for serotonin. 5-Methoxytryptamine and α-methylserotonin, both known 5-HT GPCR agonists, give stimulatory responses whereas methysergide, a known 5-HT GPCR ligand, give an antagonist response in Eca-5-HT_1a. Mutants obtained by the substitution of key predicted residues resulted in severe impairment of receptor activity, confirming that indeed, these residues have important roles in receptor function. Immunolocalization studies on the protoscolex stage from E. canadensis, showed that Eca-5-HT_1a is localized in branched fibers which correspond to the nervous system of the parasite. The patterns of immunoreactive fibers for Eca-5-HT_1a and for serotonin were intimately intertwined but not identical, suggesting that they are two separate groups of fibers. These data provide the first functional, pharmacological and localization report of a serotonergic receptor that putatively belongs to the 5-HT1 type of GPCRs in cestodes. The serotonergic GPCR characterized here may represent a new target for antiparasitic intervention.

GPCR activation mechanisms across classes and macro/microscales

Two-thirds of human hormones and one-third of clinical drugs activate ~350 G-protein-coupled receptors (GPCR) belonging to four classes: A, B1, C and F. Whereas a model of activation has been described for class A, very little is known about the activation of the other classes, which differ by being activated by endogenous ligands bound mainly or entirely extracellularly. Here we show that, although they use the same structural scaffold and share several 'helix macroswitches', the GPCR classes differ in their 'residue microswitch' positions and contacts. We present molecular mechanistic maps of activation for each GPCR class and methods for contact analysis applicable for any functional determinants. This provides a superfamily residue-level rationale for conformational selection and allosteric communication by ligands and G proteins, laying the foundation for receptor-function studies and drugs with the desired modality.

Subtle Chemical Changes Cross the Boundary between Agonist and Antagonist: New A3 Adenosine Receptor Homology Models and Structural Network Analysis Can Predict This Boundary

Distinguishing compounds' agonistic or antagonistic behavior would be of great utility for the rational discovery of selective modulators. We synthesized truncated nucleoside derivatives and discovered 6c (K_i = 2.40 nM) as a potent human A₃ adenosine receptor (hA₃AR) agonist, and subtle chemical modification induced a shift from antagonist to agonist. We elucidated this shift by developing new hA₃AR homology models that consider the pharmacological profiles of the ligands. Taken together with molecular dynamics (MD) simulation and three-dimensional (3D) structural network analysis of the receptor-ligand complex, the results indicated that the hydrogen bonding with Thr94^3.36 and His272^7.43 could make a stable interaction between the 3'-amino group with TM3 and TM7, and the corresponding induced-fit effects may play important roles in rendering the agonistic effect. Our results provide a more precise understanding of the compounds' actions at the atomic level and a rationale for the design of new drugs with specific pharmacological profiles.

Activation-Induced Reorganization of Energy Transport Networks in the β2 Adrenergic Receptor

We compute energy exchange networks (EENs) through the β₂ adrenergic receptor (β₂AR), a G-protein coupled receptor (GPCR), in inactive and active states, based on the results of molecular dynamics simulations of this membrane bound protein. We introduce a new definition for the reorganization of EENs upon activation that depends on the relative change in rates of energy transfer across noncovalent contacts throughout the protein. On the basis of the reorganized network that we obtain for β₂AR upon activation, we identify a branched pathway between the agonist binding site and the cytoplasmic region, where a G-protein binds to the receptor when activated. The pathway includes all of the motifs containing molecular switches previously identified as contributing to the allosteric transition of β₂AR upon agonist binding. EENs and their reorganization upon activation are compared with structure-based contact networks computed for the inactive and active states of β₂AR.

Characterization of the G protein-coupled receptor family SREB across fish evolution

The SREB (Super-conserved Receptors Expressed in Brain) family of G protein-coupled receptors is highly conserved across vertebrates and consists of three members: SREB1 (orphan receptor GPR27), SREB2 (GPR85), and SREB3 (GPR173). Ligands for these receptors are largely unknown or only recently identified, and functions for all three are still beginning to be understood, including roles in glucose homeostasis, neurogenesis, and hypothalamic control of reproduction. In addition to the brain, all three are expressed in gonads, but relatively few studies have focused on this, especially in non-mammalian models or in an integrated approach across the entire receptor family. The purpose of this study was to more fully characterize sreb genes in fish, using comparative genomics and gonadal expression analyses in five diverse ray-finned (Actinopterygii) species across evolution. Several unique characteristics were identified in fish, including: (1) a novel, fourth euteleost-specific gene (sreb3b or gpr173b) that likely emerged from a copy of sreb3 in a separate event after the teleost whole genome duplication, (2) sreb3a gene loss in Order Cyprinodontiformes, and (3) expression differences between a gar species and teleosts. Overall, gonadal patterns suggested an important role for all sreb genes in teleost testicular development, while gar were characterized by greater ovarian expression that may reflect similar roles to mammals. The novel sreb3b gene was also characterized by several unique features, including divergent but highly conserved amino acid positions, and elevated brain expression in puffer (Dichotomyctere nigroviridis) that more closely matched sreb2, not sreb3a. These results demonstrate that SREBs may differ among vertebrates in genomic structure and function, and more research is needed to better understand these roles in fish.

Structural insights on ligand recognition at the human leukotriene B4 receptor 1

The leukotriene B4 receptor 1 (BLT1) regulates the recruitment and chemotaxis of different cell types and plays a role in the pathophysiology of infectious, allergic, metabolic, and tumorigenic human diseases. Here we present a crystal structure of human BLT1 (hBLT1) in complex with a selective antagonist MK-D-046, developed for the treatment of type 2 diabetes and other inflammatory conditions. Comprehensive analysis of the structure and structure-activity relationship data, reinforced by site-directed mutagenesis and docking studies, reveals molecular determinants of ligand binding and selectivity toward different BLT receptor subtypes and across species. The structure helps to identify a putative membrane-buried ligand access channel as well as potential receptor binding modes of endogenous agonists. These structural insights of hBLT1 enrich our understanding of its ligand recognition and open up future avenues in structure-based drug design.

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.

G protein-coupled receptor-G protein interactions: a single-molecule perspective

G protein-coupled receptors (GPCRs) regulate many cellular and physiological processes, responding to a diverse range of extracellular stimuli including hormones, neurotransmitters, odorants, and light. Decades of biochemical and pharmacological studies have provided fundamental insights into the mechanisms of GPCR signaling. Thanks to recent advances in structural biology, we now possess an atomistic understanding of receptor activation and G protein coupling. However, how GPCRs and G proteins interact in living cells to confer signaling efficiency and specificity remains insufficiently understood. The development of advanced optical methods, including single-molecule microscopy, has provided the means to study receptors and G proteins in living cells with unprecedented spatio-temporal resolution. The results of these studies reveal an unexpected level of complexity, whereby GPCRs undergo transient interactions among themselves as well as with G proteins and structural elements of the plasma membrane to form short-lived signaling nanodomains that likely confer both rapidity and specificity to GPCR signaling. These findings may provide new strategies to pharmaceutically modulate GPCR function, which might eventually pave the way to innovative drugs for common diseases such as diabetes or heart failure.

The Dynamics of the Neuropeptide Y Receptor Type 1 Investigated by Solid-State NMR and Molecular Dynamics Simulation

We report data on the structural dynamics of the neuropeptide Y (NPY) G-protein-coupled receptor (GPCR) type 1 (Y1R), a typical representative of class A peptide ligand GPCRs, using a combination of solid-state NMR and molecular dynamics (MD) simulation. First, the equilibrium dynamics of Y1R were studied using 15N-NMR and quantitative determination of 1H-13C order parameters through the measurement of dipolar couplings in separated-local-field NMR experiments. Order parameters reporting the amplitudes of the molecular motions of the C-H bond vectors of Y1R in DMPC membranes are 0.57 for the Cα sites and lower in the side chains (0.37 for the CH2 and 0.18 for the CH3 groups). Different NMR excitation schemes identify relatively rigid and also dynamic segments of the molecule. In monounsaturated membranes composed of longer lipid chains, Y1R is more rigid, attributed to a higher hydrophobic thickness of the lipid membrane. The presence of an antagonist or NPY has little influence on the amplitude of motions, whereas the addition of agonist and arrestin led to a pronounced rigidization. To investigate Y1R dynamics with site resolution, we conducted extensive all-atom MD simulations of the apo and antagonist-bound state. In each state, three replicas with a length of 20 μs (with one exception, where the trajectory length was 10 μs) were conducted. In these simulations, order parameters of each residue were determined and showed high values in the transmembrane helices, whereas the loops and termini exhibit much lower order. The extracellular helix segments undergo larger amplitude motions than their intracellular counterparts, whereas the opposite is observed for the loops, Helix 8, and termini. Only minor differences in order were observed between the apo and antagonist-bound state, whereas the time scale of the motions is shorter for the apo state. Although these relatively fast motions occurring with correlation times of ns up to a few µs have no direct relevance for receptor activation, it is believed that they represent the prerequisite for larger conformational transitions in proteins.

Deciphering collaborative sidechain motions in proteins during molecular dynamics simulations

The dynamic structure of proteins is essential for their functions and may include large conformational transitions which can be studied by molecular dynamics (MD) simulations. However, details of these transitions are difficult to automatically track. To facilitate their analysis, we developed two scores of correlation between sidechain dihedral angles. The CIRCULAR and OMES scores are computed from, respectively, dihedral angle values and rotamer distributions. As a case study, we applied our methods to an activation-like transition of the chemokine receptor CXCR4, observed during accelerated MD simulations. The principal component analysis of the correlation matrices was consistent with the networking structure of the top ranking pairs. Both scores identify a set of residues whose “collaborative” sidechain rotamerization immediately preceded or accompanied the conformational transition of CXCR4. Detailed analysis of the sequential order of these rotamerizations suggests that an allosteric mechanism, involving the outward motion of an asparagine residue in transmembrane helix 3, might be a prerequisite to the large scale conformational transition of CXCR4. This case study provides the proof-of-concept that the correlation methods developed here are valuable exploratory techniques to help decipher complex reactional pathways.

The precious fluorine on the ring: fluorine NMR for biological systems

The fluorine-19 nucleus was recognized early to harbor exceptional properties for NMR spectroscopy. With 100% natural abundance, a high gyromagnetic ratio (83% sensitivity compared to 1H), a chemical shift that is extremely sensitive to its surroundings and near total absence in biological systems, it was destined to become a favored NMR probe, decorating small and large molecules. However, after early excitement, where uptake of fluorinated aromatic amino acids was explored in a series of animal studies, 19F-NMR lost popularity, especially in large molecular weight systems, due to chemical shift anisotropy (CSA) induced line broadening at high magnetic fields. Recently, two orthogonal approaches, (i) CF3 labeling and (ii) aromatic 19F-13C labeling leveraging the TROSY (Transverse Relaxation Optimized Spectroscopy) effect have been successfully applied to study large biomolecular systems. In this perspective, we will discuss the fascinating early work with fluorinated aromatic amino acids, which reveals the enormous potential of these non-natural amino acids in biological NMR and the potential of 19F-NMR to characterize protein and nucleic acid structure, function and dynamics in the light of recent developments. Finally, we explore how fluorine NMR might be exploited to implement small molecule or fragment screens that resemble physiological conditions and discuss the opportunity to follow the fate of small molecules in living cells.

Molecular evolution of a collage of cholesterol interaction motifs in transmembrane helix V of the serotonin1A receptor

The human serotonin_1A receptor is a representative member of the superfamily of G protein-coupled receptors (GPCRs) and an important drug target for neurological disorders. Using a combination of biochemical, biophysical and molecular dynamics simulation approaches, we and others have shown that membrane cholesterol modulates the organization, dynamics and function of vertebrate serotonin_1A receptors. Previous studies have shown that the cytoplasmic portion of transmembrane helix V (TM V) and the extramembraneous intracellular loop 3 are critical for G-protein coupling, phosphorylation and desensitization of the receptor. We have recently resolved a collage of putative cholesterol interaction motifs from the amino acid sequence overlapping this region. In this paper, we explore the sequence plasticity of this fragment that may have adapted to altered membrane lipidome, after vertebrates evolved from primordial invertebrates. Since invertebrates have lower levels of membrane cholesterol relative to vertebrates, we compared TM V sequence fragments from invertebrate serotonin₁ receptors with vertebrate orthologs to infer the sequence plasticity in TM V. We report that the average number of cholesterol interaction motifs in TM V for diverse phyla represents an increasing trend that could mirror vertebrate evolution from primordial invertebrates. By statistical modeling, we propose that the collage of cholesterol interaction motifs in TM V of the human serotonin_1A receptor may have evolved from rudimentary collages, reminiscent of primordial invertebrate orthologs. Taken together, we propose that a repertoire of cholesterol-philic nonsynonymous substitutions may have enhanced collage complexity in TM V during vertebrate evolution.

Structural Complexity and Plasticity of Signaling Regulation at the Melanocortin-4 Receptor

The melanocortin-4 receptor (MC4R) is a class A G protein-coupled receptor (GPCR), essential for regulation of appetite and metabolism. Pathogenic inactivating MC4R mutations are the most frequent cause of monogenic obesity, a growing medical and socioeconomic problem worldwide. The MC4R mediates either ligand-independent or ligand-dependent signaling. Agonists such as α-melanocyte-stimulating hormone (α-MSH) induce anorexigenic effects, in contrast to the endogenous inverse agonist agouti-related peptide (AgRP), which causes orexigenic effects by suppressing high basal signaling activity. Agonist action triggers the binding of different subtypes of G proteins and arrestins, leading to concomitant induction of diverse intracellular signaling cascades. An increasing number of experimental studies have unraveled molecular properties and mechanisms of MC4R signal transduction related to physiological and pathophysiological aspects. In addition, the MC4R crystal structure was recently determined at 2.75 Å resolution in an inactive state bound with a peptide antagonist. Underpinned by structural homology models of MC4R complexes simulating a presumably active-state conformation compared to the structure of the inactive state, we here briefly summarize the current understanding and key players involved in the MC4R switching process between different activity states. Finally, these perspectives highlight the complexity and plasticity in MC4R signaling regulation and identify gaps in our current knowledge.

An investigation into the allosteric mechanism of GPCR A2A adenosine receptor with trajectory-based information theory and complex network model

G protein-coupled receptors (GPCRs), a large superfamily of transmembrane (TM) proteins, allosterically transduce the signal of ligand binding in the extracellular (EC) domain to couple to effector proteins in the intracellular (IC) domain, therefore forming the largest class of drug targets. The A_2A adenosine receptor (A_2AAR), a class-A GPCR, has been extensively studied as it offers numerous possibilities for therapeutic applications. However, the mechanism of allosteric communication between EC and IC domains is not completely clear. In this work, we utilize torsional mutual information to quantify the correlated motions of residue pairs from its molecular dynamics (MD) simulation trajectories, and further use the complex network model to obtain allosteric pipelines and hubs. The identified allosteric communication pipelines mainly transmit the signal from EC domain to the cytoplasmic ends of TM helix 5 (TM5), TM6 and TM7. The allosteric hubs, mostly located at TM5, TM6 and TM7, play an important role in mediating allosteric signal transmission to keep the receptor rigid and prevent G protein from binding to IC domain, which can explain the reason why their mutations distant from ligand-binding site do not affect the ligand binding affinity but affect the ligand efficacy. Additionally, we identify the key residues located in antagonist ZM241385 binding pocket which mediate multiple allosteric pathways and have been experimentally proven to play a critical role in affecting the ligand potency. This study is helpful for understanding the allosteric communication mechanism of A_2AAR, and can provide valuable information for the structure-based drug design of GPCRs.Communicated by Ramaswamy H. Sarma.

Structural insights into emergent signaling modes of G protein-coupled receptors

G protein-coupled receptors (GPCRs) represent the largest family of cell membrane proteins, with >800 GPCRs in humans alone, and recognize highly diverse ligands, ranging from photons to large protein molecules. Very important to human medicine, GPCRs are targeted by about 35% of prescription drugs. GPCRs are characterized by a seven-transmembrane α-helical structure, transmitting extracellular signals into cells to regulate major physiological processes via heterotrimeric G proteins and β-arrestins. Initially viewed as receptors whose signaling via G proteins is delimited to the plasma membrane, it is now recognized that GPCRs signal also at various intracellular locations, and the mechanisms and (patho)physiological relevance of such signaling modes are actively investigated. The propensity of GPCRs to adopt different signaling modes is largely encoded in the structural plasticity of the receptors themselves and of their signaling complexes. Here, we review emerging modes of GPCR signaling via endosomal membranes and the physiological implications of such signaling modes. We further summarize recent structural insights into mechanisms of GPCR activation and signaling. We particularly emphasize the structural mechanisms governing the continued GPCR signaling from endosomes and the structural aspects of the GPCR resensitization mechanism and discuss the recently uncovered and important roles of lipids in these processes.