Molecular basis for high-affinity agonist binding in GPCRs

Nanobodies against the myelin enzyme CNPase as tools for structural and functional studies

2',3'-cyclic nucleotide 3'-phosphodiesterase (CNPase) is an abundant constituent of central nervous system non-compact myelin, frequently used as a marker antigen for myelinating cells. The catalytic activity of CNPase, the 3'-hydrolysis of 2',3'-cyclic nucleotides, is well characterised in vitro, but the in vivo function of CNPase remains unclear. CNPase interacts with the actin cytoskeleton to counteract the developmental closure of cytoplasmic channels that travel through compact myelin; its enzymatic activity may be involved in adenosine metabolism and RNA degradation. We developed a set of high-affinity nanobodies recognizing the phosphodiesterase domain of CNPase, and the crystal structures of each complex show that the five nanobodies have distinct epitopes. One of the nanobodies bound deep into the CNPase active site and acted as an inhibitor. Moreover, the nanobodies were characterised in imaging applications and as intrabodies, expressed in mammalian cells, such as primary oligodendrocytes. Fluorescently labelled nanobodies functioned in imaging of teased nerve fibers and whole brain tissue sections, as well as super-resolution microscopy. These anti-CNPase nanobodies provide new tools for structural and functional biology of myelination, including high-resolution imaging of nerve tissue.

Binding kinetics drive G protein subtype selectivity at the β1-adrenergic receptor

G protein-coupled receptors (GPCRs) bind to different G protein α-subtypes with varying degrees of selectivity. The mechanism by which GPCRs achieve this selectivity is still unclear. Using 13C methyl methionine and 19F NMR, we investigate the agonist-bound active state of β1AR and its ternary complexes with different G proteins in solution. We find the receptor in the ternary complexes adopts very similar conformations. In contrast, the full agonist-bound receptor active state assumes a conformation differing from previously characterised activation intermediates or from β1AR in ternary complexes. Assessing the kinetics of binding for the agonist-bound receptor with different G proteins, we find the increased affinity of β1AR for Gs results from its much faster association with the receptor. Consequently, we suggest a kinetic-driven selectivity gate between canonical and secondary coupling which arises from differential favourability of G protein binding to the agonist-bound receptor active state.

Identification of a Novel Subtype-Selective α1B-Adrenoceptor Antagonist

α1A-, α1B-, and α1D-adrenoceptors (α1-ARs) are members of the adrenoceptor G protein-coupled receptor family that are activated by adrenaline (epinephrine) and noradrenaline. α1-ARs are clinically targeted using antagonists that have minimal subtype selectivity, such as prazosin and tamsulosin, to treat hypertension and benign prostatic hyperplasia, respectively. Abundant expression of α1-ARs in the heart and central nervous system (CNS) makes these receptors potential targets for the treatment of cardiovascular and CNS disorders, such as heart failure, epilepsy, and Alzheimer's disease. Our understanding of the precise physiological roles of α1-ARs, however, and their involvement in disease has been hindered by the lack of sufficiently subtype-selective tool compounds, especially for α1B-AR. Here, we report the discovery of 4-[(2-hydroxyethyl)amino]-6-methyl-2H-chromen-2-one (Cpd1), as an α1B-AR antagonist that has 10-15-fold selectivity over α1A-AR and α1D-AR. Through computational and site-directed mutagenesis studies, we have identified the binding site of Cpd1 in α1B-AR and propose the molecular basis of α1B-AR selectivity, where the nonconserved V19745.52 residue plays a major role, with contributions from L3146.55 within the α1B-AR pocket. By exploring the structure-activity relationships of Cpd1 at α1B-AR, we have also identified 3-[(cyclohexylamino)methyl]-6-methylquinolin-2(1H)-one (Cpd24), which has a stronger binding affinity than Cpd1, albeit with reduced selectivity for α1B-AR. Cpd1 and Cpd24 represent potential leads for α1B-AR-selective drug discovery and novel tool molecules to further study the physiology of α1-ARs.

Molecular determinants of ligand efficacy and potency in GPCR signaling

Heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs) bind to extracellular ligands and drugs and modulate intracellular responses through conformational changes. Despite their importance as drug targets, the molecular origins of pharmacological properties such as efficacy (maximum signaling response) and potency (the ligand concentration at half-maximal response) remain poorly understood for any ligand-receptor-signaling system. We used the prototypical adrenaline-β2 adrenergic receptor-G protein system to reveal how specific receptor residues decode and translate the information encoded in a ligand to mediate a signaling response. We present a data science framework to integrate pharmacological and structural data to uncover structural changes and allosteric networks relevant for ligand pharmacology. These methods can be tailored to study any ligand-receptor-signaling system, and the principles open possibilities for designing orthosteric and allosteric compounds with defined signaling properties.

Unravelling the mechanism of neurotensin recognition by neurotensin receptor 1

The conformational ensembles of G protein-coupled receptors (GPCRs) include inactive and active states. Spectroscopy techniques, including NMR, show that agonists, antagonists and other ligands shift the ensemble toward specific states depending on the pharmacological efficacy of the ligand. How receptors recognize ligands and the kinetic mechanism underlying this population shift is poorly understood. Here, we investigate the kinetic mechanism of neurotensin recognition by neurotensin receptor 1 (NTS1) using 19F-NMR, hydrogen-deuterium exchange mass spectrometry and stopped-flow fluorescence spectroscopy. Our results indicate slow-exchanging conformational heterogeneity on the extracellular surface of ligand-bound NTS1. Numerical analysis of the kinetic data of neurotensin binding to NTS1 shows that ligand recognition follows an induced-fit mechanism, in which conformational changes occur after neurotensin binding. This approach is applicable to other GPCRs to provide insight into the kinetic regulation of ligand recognition by GPCRs.

Quantitative flow cytometric selection of tau conformational nanobodies specific for pathological aggregates

Single-domain antibodies, also known as nanobodies, are broadly important for studying the structure and conformational states of several classes of proteins, including membrane proteins, enzymes, and amyloidogenic proteins. Conformational nanobodies specific for aggregated conformations of amyloidogenic proteins are particularly needed to better target and study aggregates associated with a growing class of associated diseases, especially neurodegenerative disorders such as Alzheimer's and Parkinson's diseases. However, there are few reported nanobodies with both conformational and sequence specificity for amyloid aggregates, especially for large and complex proteins such as the tau protein associated with Alzheimer's disease, due to difficulties in selecting nanobodies that bind to complex aggregated proteins. Here, we report the selection of conformational nanobodies that selectively recognize aggregated (fibrillar) tau relative to soluble (monomeric) tau. Notably, we demonstrate that these nanobodies can be directly isolated from immune libraries using quantitative flow cytometric sorting of yeast-displayed libraries against tau aggregates conjugated to quantum dots, and this process eliminates the need for secondary nanobody screening. The isolated nanobodies demonstrate conformational specificity for tau aggregates in brain samples from both a transgenic mouse model and human tauopathies. We expect that our facile approach will be broadly useful for isolating conformational nanobodies against diverse amyloid aggregates and other complex antigens.

Molecular recognition of an aversive odorant by the murine trace amine-associated receptor TAAR7f

There are two main families of G protein-coupled receptors that detect odours in humans, the odorant receptors (ORs) and the trace amine-associated receptors (TAARs). Their amino acid sequences are distinct, with the TAARs being most similar to the aminergic receptors such as those activated by adrenaline, serotonin and histamine. To elucidate the structural determinants of ligand recognition by TAARs, we have determined the cryo-EM structure of a murine receptor, mTAAR7f, coupled to the heterotrimeric G protein Gs and bound to the odorant N,N-dimethylcyclohexylamine (DMCH) to an overall resolution of 2.9 Å. DMCH is bound in a hydrophobic orthosteric binding site primarily through van der Waals interactions and a strong charge-charge interaction between the tertiary amine of the ligand and an aspartic acid residue. This site is distinct and non-overlapping with the binding site for the odorant propionate in the odorant receptor OR51E2. The structure, in combination with mutagenesis data and molecular dynamics simulations suggests that the activation of the receptor follows a similar pathway to that of the β-adrenoceptors, with the significant difference that DMCH interacts directly with one of the main activation microswitch residues.

Comparative Saturation Binding Analysis of 64Cu-Labeled Somatostatin Analogues Using Cell Homogenates and Intact Cells

The development of novel ligands for G-protein-coupled receptors (GPCRs) typically entails the characterization of their binding affinity, which is often performed with radioligands in a competition or saturation binding assay format. Since GPCRs are transmembrane proteins, receptor samples for binding assays are prepared from tissue sections, cell membranes, cell homogenates, or intact cells. As part of our investigations on modulating the pharmacokinetics of radiolabeled peptides for improved theranostic targeting of neuroendocrine tumors with a high abundance of the somatostatin receptor sub-type 2 (SST₂), we characterized a series of ⁶⁴Cu-labeled [Tyr³]octreotate (TATE) derivatives in vitro in saturation binding assays. Herein, we report on the SST₂ binding parameters measured toward intact mouse pheochromocytoma cells and corresponding cell homogenates and discuss the observed differences taking the physiology of SST₂ and GPCRs in general into account. Furthermore, we point out method-specific advantages and limitations.

Clinical pharmacology of β-3 adrenergic receptor agonists for cardiovascular diseases

INTRODUCTION

Few agonists of the third isotype of beta-adrenergic receptors, the β3-adrenoreceptor, are currently used clinically, and new agonists are under development for the treatment of overactive bladder disease. As the receptor is expressed in human cardiac and vascular tissues, it is important to understand their beneficial (or adverse) effect(s) on these targets.

AREAS COVERED

We discuss the most recent results of clinical trials testing the benefit and safety of β3-adrenoreceptor activation on cardiovascular outcomes in light of current knowledge on the receptor biology, genetic polymorphisms, and agonist pharmacology.

EXPERT OPINION

While evidence from small clinical trials is limited so far, the β3-agonist, mirabegron seems to be safe in patients at high cardiovascular risk but produces benefits on selected cardiovascular outcomes only at higher than standard doses. Activation of cardiovascular β3-adrenoreceptors deserves to be tested with more potent agonists, such as vibegron.

Non-traditional approaches for control of antibiotic resistance

INTRODUCTION

The drying up of antibiotic pipeline has necessitated the development of alternative therapeutic strategies to control the problem of antimicrobial resistance (AMR) that is expected to kill 10-million people annually by 2050. Newer therapeutic approaches address the shortcomings of traditional small-molecule antibiotics - the lack of specificity, evolvability, and susceptibility to mutation-based resistance. These 'non-traditional' molecules are biologicals having a complex structure and mode(s) of action that makes them resilient to resistance.

AREAS COVERED

This review aims to provide information about the non-traditional drug development approaches to tackle the problem of antimicrobial resistance, from the pre-antibiotic era to the latest developments. We have covered the molecules under development in the clinic with literature sourced from reviewed scholarly articles, official company websites involved in innovation of concerned therapeutics, press releases from the regulatory bodies, and clinical trial databases.

EXPERT OPINION

Formal introduction of non-traditional therapies in general practice can be quick and feasible only if supported with companion diagnostics and used in conjunction with established therapies. Owing to relatively higher development costs, non-traditional therapeutics require more funding as well as well as clarity in regulatory and clinical path. We are hopeful these issues are adequately addressed before AMR develops into a pandemic.

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.

Computational Analysis of Histamine Protonation Effects on H1R Binding

Despite numerous studies investigating histamine and its receptors, the impact of histamine protonation states on binding to the histamine H1-receptor (H1R) has remained elusive. Therefore, we assessed the influence of different histamine tautomers (τ-tautomer, π-tautomer) and charge states (mono- vs. dicationic) on the interaction with the ternary histamine-H1R-G_q complex. In atomistic molecular dynamics simulations, the τ-tautomer formed stable interactions with the receptor, while the π-tautomer induced a rotation of the histamine ring by 180° and formed only weaker hydrogen bonding interactions. This suggests that the τ-tautomer is more relevant for stabilization of the active ternary histamine-H1R-G_q complex. In addition to the two monocationic tautomers, the binding of dicationic histamine was investigated, whose interaction with the H1R had been observed in a previous experimental study. Our simulations showed that the dication is less compatible with the ternary histamine-H1R-G_q complex and rather induces an inactive conformation in the absence of the G_q protein. Our data thus indicate that the charge state of histamine critically affects its interactions with the H1R. Ultimately these findings might have implications for the future development of new ligands that stabilize distinct H1R activation states.

Improving the Modeling of Extracellular Ligand Binding Pockets in RosettaGPCR for Conformational Selection

G protein-coupled receptors (GPCRs) are the largest class of drug targets and undergo substantial conformational changes in response to ligand binding. Despite recent progress in GPCR structure determination, static snapshots fail to reflect the conformational space of putative binding pocket geometries to which small molecule ligands can bind. In comparative modeling of GPCRs in the absence of a ligand, often a shrinking of the orthosteric binding pocket is observed. However, the exact prediction of the flexible orthosteric binding site is crucial for adequate structure-based drug discovery. In order to improve ligand docking and guide virtual screening experiments in computer-aided drug discovery, we developed RosettaGPCRPocketSize. The algorithm creates a conformational ensemble of biophysically realistic conformations of the GPCR binding pocket between the TM bundle, which is consistent with a knowledge base of expected pocket geometries. Specifically, tetrahedral volume restraints are defined based on information about critical residues in the orthosteric binding site and their experimentally observed range of Cα-Cα-distances. The output of RosettaGPCRPocketSize is an ensemble of binding pocket geometries that are filtered by energy to ensure biophysically probable arrangements, which can be used for docking simulations. In a benchmark set, pocket shrinkage observed in the default RosettaGPCR was reduced by up to 80% and the binding pocket volume range and geometric diversity were increased. Compared to models from four different GPCR homology model databases (RosettaGPCR, GPCR-Tasser, GPCR-SSFE, and GPCRdb), the here-created models showed more accurate volumes of the orthosteric pocket when evaluated with respect to the crystallographic reference structure. Furthermore, RosettaGPCRPocketSize was able to generate an improved realistic pocket distribution. However, while being superior to other homology models, the accuracy of generated model pockets was comparable to AlphaFold2 models. Furthermore, in a docking benchmark using small-molecule ligands with a higher molecular weight between 400 and 700 Da, a higher success rate in creating native-like binding poses was observed. In summary, RosettaGPCRPocketSize can generate GPCR models with realistic orthosteric pocket volumes, which are useful for structure-based drug discovery applications.

Constrained catecholamines gain β2AR selectivity through allosteric effects on pocket dynamics

G protein-coupled receptors (GPCRs) within the same subfamily often share high homology in their orthosteric pocket and therefore pose challenges to drug development. The amino acids that form the orthosteric binding pocket for epinephrine and norepinephrine in the β1 and β2 adrenergic receptors (β1AR and β2AR) are identical. Here, to examine the effect of conformational restriction on ligand binding kinetics, we synthesized a constrained form of epinephrine. Surprisingly, the constrained epinephrine exhibits over 100-fold selectivity for the β2AR over the β1AR. We provide evidence that the selectivity may be due to reduced ligand flexibility that enhances the association rate for the β2AR, as well as a less stable binding pocket for constrained epinephrine in the β1AR. The differences in the amino acid sequence of the extracellular vestibule of the β1AR allosterically alter the shape and stability of the binding pocket, resulting in a marked difference in affinity compared to the β2AR. These studies suggest that for receptors containing identical binding pocket residues, the binding selectivity may be influenced in an allosteric manner by surrounding residues, like those of the extracellular loops (ECLs) that form the vestibule. Exploiting these allosteric influences may facilitate the development of more subtype-selective ligands for GPCRs.

Mandibulofacial dysostosis with alopecia results from ETAR gain-of-function mutations via allosteric effects on ligand binding

Mutations of G protein-coupled receptors (GPCRs) cause various human diseases, but the mechanistic details are limited. Here, we establish p.E303K in the gene encoding the endothelin receptor type A (ETAR/EDNRA) as a recurrent mutation causing mandibulofacial dysostosis with alopecia (MFDA), with craniofacial changes similar to those caused by p.Y129F. Mouse models carrying either of these missense mutations exhibited a partial maxillary-to-mandibular transformation, which was rescued by deleting the ligand endothelin 3 (ET3/EDN3). Pharmacological experiments confirmed the causative ETAR mutations as gain of function, dependent on ET3. To elucidate how an amino acid substitution far from the ligand binding site can increase ligand affinity, we used molecular dynamics (MD) simulations. E303 is located at the intracellular end of transmembrane domain 6, and its replacement by a lysine increased flexibility of this portion of the helix, thus favoring G protein binding and leading to G protein-mediated enhancement of agonist affinity. The Y129F mutation located under the ligand binding pocket reduced the sodium-water network, thereby affecting the extracellular portion of helices in favor of ET3 binding. These findings provide insight into the pathogenesis of MFDA and into allosteric mechanisms regulating GPCR function, which may provide the basis for drug design targeting GPCRs.

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.

Rational design, synthesis, and pharmacological evaluation of a cohort of novel beta-adrenergic receptors ligands enables an assessment of structure-activity relationships

Biomedical applications of molecules that are able to modulate β-adrenergic signaling have become increasingly attractive over the last decade, revealing that β-adrenergic receptors (β-ARs) are key targets for a plethora of therapeutic interventions, including cancer. Despite successes in β-AR drug discovery, identification of β-AR ligands that are useful as selective chemical tools in pharmacological studies of the three β-AR subtypes, or lead compounds for drug development is still a highly challenging task. This is mainly due to the intrinsic plasticity of β-ARs as G protein-coupled receptors in conjunction with the requirement for functional receptor subtype selectivity, tissue specificity and minimal off-target effects. With the aim to provide insight into structure-activity relationships for the three β-AR subtypes, we have synthesized and obtained the pharmacological profile of a series of structurally diverse compounds (named MC) that were designed based on the aryloxy-propanolamine scaffold of SR59230A. Comparative analysis of their predicted binding mode within the active and inactive states of the receptors in combination with their pharmacological profile revealed key structural elements that control their activity as agonists or antagonists, in addition to clues about substituents that mediate selectivity for one receptor subtype over the others. We anticipate that these results will facilitate selective β-AR drug development efforts.

The β2-adrenergic receptor in the apical membrane of intestinal enterocytes senses sugars to stimulate glucose uptake from the gut

The presence of sugar in the gut causes induction of SGLT1, the sodium/glucose cotransporter in intestinal epithelial cells (enterocytes), and this is accompanied by stimulation of sugar absorption. Sugar sensing was suggested to involve a G-protein coupled receptor and cAMP - protein kinase A signalling, but the sugar receptor has remained unknown. We show strong expression and co-localization with SGLT1 of the β2-adrenergic receptor (β ₂-AR) at the enterocyte apical membrane and reveal its role in stimulating glucose uptake from the gut by the sodium/glucose-linked transporter, SGLT1. Upon heterologous expression in different reporter systems, the β ₂-AR responds to multiple sugars in the mM range, consistent with estimated gut sugar levels after a meal. Most adrenergic receptor antagonists inhibit sugar signaling, while some differentially inhibit epinephrine and sugar responses. However, sugars did not inhibit binding of I¹²⁵-cyanopindolol, a β ₂-AR antagonist, to the ligand-binding site in cell-free membrane preparations. This suggests different but interdependent binding sites. Glucose uptake into everted sacs from rat intestine was stimulated by epinephrine and sugars in a β ₂-AR-dependent manner. STD-NMR confirmed direct physical binding of glucose to the β ₂-AR. Oral administration of glucose with a non-bioavailable β ₂-AR antagonist lowered the subsequent increase in blood glucose levels, confirming a role for enterocyte apical β ₂-ARs in stimulating gut glucose uptake, and suggesting enterocyte β ₂-AR as novel drug target in diabetic and obese patients. Future work will have to reveal how glucose sensing by enterocytes and neuroendocrine cells is connected, and whether β ₂-ARs mediate glucose sensing also in other tissues.

Mass spectrometry captures biased signalling and allosteric modulation of a G-protein-coupled receptor

G-protein-coupled receptors signal through cognate G proteins. Despite the widespread importance of these receptors, their regulatory mechanisms for G-protein selectivity are not fully understood. Here we present a native mass spectrometry-based approach to interrogate both biased signalling and allosteric modulation of the β1-adrenergic receptor in response to various ligands. By simultaneously capturing the effects of ligand binding and receptor coupling to different G proteins, we probed the relative importance of specific interactions with the receptor through systematic changes in 14 ligands, including isoprenaline derivatives, full and partial agonists, and antagonists. We observed enhanced dynamics of the intracellular loop 3 in the presence of isoprenaline, which is capable of acting as a biased agonist. We also show here that endogenous zinc ions augment the binding in receptor-Gs complexes and propose a zinc ion-binding hotspot at the TM5/TM6 intracellular interface of the receptor-Gs complex. Further interrogation led us to propose a mechanism in which zinc ions facilitate a structural transition of the intermediate complex towards the stable state.

Assessment of the Generalization Abilities of Machine-Learning Scoring Functions for Structure-Based Virtual Screening

In structure-based virtual screening (SBVS), it is critical that scoring functions capture protein-ligand atomic interactions. By focusing on the local domains of ligand binding pockets, a standardized pocket Pfam-based clustering (Pfam-cluster) approach was developed to assess the cross-target generalization ability of machine-learning scoring functions (MLSFs). Subsequently, 12 typical MLSFs were evaluated using random cross-validation (Random-CV), protein sequence similarity-based cross-validation (Seq-CV), and pocket Pfam-based cross-validation (Pfam-CV) methods. Surprisingly, all of the tested models showed decreased performances from Random-CV to Seq-CV to Pfam-CV experiments, not showing satisfactory generalization capacity. Our interpretable analysis suggested that the predictions on novel targets by MLSFs were dependent on buried solvent-accessible surface area (SASA)-related features of complex structures, with greater predicted binding affinities on complexes owning larger protein-ligand interfaces. By combining buried SASA-related features with target-specific patterns that were only shared among structurally similar compounds in the same cluster, the random forest (RF)-Score attained a good performance in the Random-CV test. Based on these findings, we strongly advise assessing the generalization ability of MLSFs with the Pfam-cluster approach and being cautious with the features learned by MLSFs.

Intracellular VHHs to monitor and modulate GPCR signaling

Single-domain antibody fragments, also known as VHHs or nanobodies, have opened promising avenues in therapeutics and in exploration of intracellular processes. Because of their unique structural properties, they can reach cryptic regions in their cognate antigen. Intracellular VHHs/antibodies primarily directed against cytosolic proteins or transcription factors have been described. In contrast, few of them target membrane proteins and even less recognize G protein-coupled receptors. These receptors are major therapeutic targets, which reflects their involvement in a plethora of physiological responses. Hence, they elicit a tremendous interest in the scientific community and in the industry. Comprehension of their pharmacology has been obscured by their conformational complexity, that has precluded deciphering their structural properties until the early 2010's. To that respect, intracellular VHHs have been instrumental in stabilizing G protein-coupled receptors in active conformations in order to solve their structure, possibly bound to their primary transducers, G proteins or β-arrestins. In contrast, the modulatory properties of VHHs recognizing the intracellular regions of G protein-coupled receptors on the induced signaling network have been poorly studied. In this review, we will present the advances that the intracellular VHHs have permitted in the field of GPCR signaling and trafficking. We will also discuss the methodological hurdles that linger the discovery of modulatory intracellular VHHs directed against GPCRs, as well as the opportunities they open in drug discovery.

The full activation mechanism of the adenosine A1 receptor revealed by GaMD and Su-GaMD simulations

The full activation process of G protein-coupled receptor (GPCR) plays an important role in cellular signal transduction. However, it remains challenging to simulate the whole process in which the GPCR is recognized and activated by a ligand and then couples to the G protein on a reasonable simulation timescale. Here, we developed a molecular dynamics (MD) approach named supervised (Su) Gaussian accelerated MD (GaMD) by incorporating a tabu-like supervision algorithm into a standard GaMD simulation. By using this Su-GaMD method, from the active and inactive structure of adenosine A₁ receptor (A₁R), we successfully revealed the full activation mechanism of A₁R, including adenosine (Ado)-A₁R recognition, preactivation of A₁R, and A₁R-G protein recognition, in hundreds of nanoseconds of simulations. The binding of Ado to the extracellular side of A₁R initiates conformational changes and the preactivation of A₁R. In turn, the binding of G_i2 to the intracellular side of A₁R causes a decrease in the volume of the extracellular orthosteric site and stabilizes the binding of Ado to A₁R. Su-GaMD could be a useful tool to reconstruct or even predict ligand-protein and protein-protein recognition pathways on a short timescale. The intermediate states revealed in this study could provide more detailed complementary structural characterizations to facilitate the drug design of A₁R in the future.

Single Domain Antibody application in bacterial infection diagnosis and neutralization

Increasing antibiotic resistance to bacterial infections causes a serious threat to human health. Efficient detection and treatment strategies are the keys to preventing and reducing bacterial infections. Due to the high affinity and antigen specificity, antibodies have become an important tool for diagnosis and treatment of various human diseases. In addition to conventional antibodies, a unique class of “heavy-chain-only” antibodies (HCAbs) were found in the serum of camelids and sharks. HCAbs binds to the antigen through only one variable domain Referred to as VHH (variable domain of the heavy chain of HCAbs). The recombinant format of the VHH is also called single domain antibody (sdAb) or nanobody (Nb). Sharks might also have an ancestor HCAb from where SdAbs or V-NAR might be engineered. Compared with traditional Abs, Nbs have several outstanding properties such as small size, high stability, strong antigen-binding affinity, high solubility and low immunogenicity. Furthermore, they are expressed at low cost in microorganisms and amenable to engineering. These superior properties make Nbs a highly desired alternative to conventional antibodies, which are extensively employed in structural biology, unravelling biochemical mechanisms, molecular imaging, diagnosis and treatment of diseases. In this review, we summarized recent progress of nanobody-based approaches in diagnosis and neutralization of bacterial infection and further discussed the challenges of Nbs in these fields.

Finding the Perfect Fit: Conformational Biosensors to Determine the Efficacy of GPCR Ligands

G protein-coupled receptors (GPCRs) are highly druggable targets that adopt numerous conformations. A ligand's ability to stabilize specific conformation(s) of its cognate receptor determines its efficacy or ability to produce a biological response. Identifying ligands that produce different receptor conformations and potentially discrete pharmacological effects (e.g., biased agonists, partial agonists, antagonists, allosteric modulators) is a major goal in drug discovery and necessary to develop drugs with better effectiveness and fewer side effects. Fortunately, direct measurements of ligand efficacy, via receptor conformational changes are possible with the recent development of conformational biosensors. In this review, we discuss classical efficacy models, including the two-state model, the ternary-complex model, and multistate models. We describe how nanobody-, transducer-, and receptor-based conformational biosensors detect and/or stabilize specific GPCR conformations to identify ligands with different levels of efficacy. In particular, conformational biosensors provide the potential to identify and/or characterize therapeutically desirable but often difficult to measure conformations of receptors faster and better than current methods. For drug discovery/development, several recent proof-of-principle studies have optimized conformational biosensors for high-throughput screening (HTS) platforms. However, their widespread use is limited by the fact that few sensors are reliably capable of detecting low-frequency conformations and technically demanding assay conditions. Nonetheless, conformational biosensors do help identify desirable ligands such as allosteric modulators, biased ligands, or partial agonists in a single assay, representing a distinct advantage over classical methods.

Nanobodies as Probes and Modulators of Cardiovascular G Protein-Coupled Receptors

ABSTRACT

Understanding the activation of G protein-coupled receptors (GPCRs) is of paramount importance to the field of cardiovascular medicine due to the critical physiological roles of these receptors and their prominence as drug targets. Although many cardiovascular GPCRs have been extensively studied as model receptors for decades, new complexities in their regulation continue to emerge. As a result, there is an ongoing need to develop novel approaches to monitor and to modulate GPCR activation. In less than a decade, nanobodies, or recombinant single-domain antibody fragments from camelids, have become indispensable tools for interrogating GPCRs both in purified systems and in living cells. Nanobodies have gained traction rapidly due to their biochemical tractability and their ability to recognize defined states of native proteins. Here, we review how nanobodies have been adopted to elucidate the structure, pharmacology, and signaling of cardiovascular GPCRs, resolving long-standing mysteries and revealing unexpected mechanisms. We also discuss how advancing technologies to discover nanobodies with tailored specificities may expand the impact of these tools for both basic science and therapeutic applications.

Filling of a water-free void explains the allosteric regulation of the β1-adrenergic receptor by cholesterol

Recent high-pressure NMR results indicate that the preactive conformation of the β₁-adrenergic receptor (β₁AR) harbours completely empty cavities of ~100 ų volume, which disappear in the active conformation of the receptor. Here we have localized these cavities using X-ray crystallography of xenon-derivatized β₁AR crystals. One of the cavities is in direct contact with the cholesterol-binding pocket. Solution NMR shows that addition of the cholesterol analogue cholesteryl hemisuccinate impedes the formation of the active conformation of detergent-solubilized β₁AR by blocking conserved G protein-coupled receptor microswitches, concomitant with an affinity reduction of both isoprenaline and G protein-mimicking nanobody Nb80 for β₁AR detected by isothermal titration calorimetry. This wedge-like action explains the function of cholesterol as a negative allosteric modulator of β₁AR. A detailed understanding of G protein-coupled receptor regulation by cholesterol by filling of a dry void and the easy scouting for such voids by xenon may provide new routes for the development of allosteric drugs.

Structures of β1-adrenergic receptor in complex with Gs and ligands of different efficacies

G-protein-coupled receptors (GPCRs) receive signals from ligands with different efficacies, and transduce to heterotrimeric G-proteins to generate different degrees of physiological responses. Previous studies revealed how ligands with different efficacies activate GPCRs. Here, we investigate how a GPCR activates G-proteins upon binding ligands with different efficacies. We report the cryo-EM structures of β₁-adrenergic receptor (β₁-AR) in complex with Gs (Gα_sGβ₁Gγ₂) and a partial agonist or a very weak partial agonist, and compare them to the β₁-AR-Gs structure in complex with a full agonist. Analyses reveal similar overall complex architecture, with local conformational differences. Cellular functional studies with mutations of β₁-AR residues show effects on the cellular signaling from β₁-AR to the cAMP response initiated by the three different ligands, with residue-specific functional differences. Biochemical investigations uncover that the intermediate state complex comprising β₁-AR and nucleotide-free Gs is more stable when binding a full agonist than a partial agonist. Molecular dynamics simulations support the local conformational flexibilities and different stabilities among the three complexes. These data provide insights into the ligand efficacy in the activation of GPCRs and G-proteins.

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

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

Molecular insights into the allosteric coupling mechanism between an agonist and two different transducers for μ-opioid receptors

G protein-coupled receptors (GPCRs) as the most important class of pharmacological targets regulate G-protein and β-arrestin-mediated signaling through allosteric interplay, which are responsible for different biochemical and physiological actions like therapeutic efficacy and side effects. However, the allosteric mechanism underlying preferentially recruiting one transducer versus the other has been poorly understood, limiting drug design. Motivated by this issue, we utilize accelerated molecular dynamics simulation coupled with potential of mean force (PMF), molecular mechanics Poisson Boltzmann surface area (MM/PBSA) and protein structure network (PSN) to study two ternary complex systems of a representative class A GPCR (μ-opioid receptor (μOR)) bound by an agonist and one specific transducer (G-protein and β-arrestin). The results show that no significant difference exists in the whole structure of μOR between two transducer couplings, but displays transducer-dependent changes in the intracellular binding region of μOR, where the β-arrestin coupling results in a narrower crevice with TM7 inward movement compared with the G-protein. In addition, both the G-protein and β-arrestin coupling can increase the binding affinity of the agonist to the receptor. However, the interactions between the agonist and μOR also exhibit transducer-specific changes, in particular for the interaction with ECL2 that plays an important role in recruiting β-arrestin. The allosteric network analysis further indicates that Y148^3.33, F152^3.37, F156^3.41, N191^4.49, T160^3.45, Y106^2.42, W293^6.48, F289^6.44, I248^5.54 and Y252^5.58 play important roles in equally activating G-protein and β-arrestin. In contrast, M161^3.46 and R165^3.50 devote important contributions to preferentially recruit G-protein while D164^3.49 and R179^ICL2 are revealed to be important for selectively activating β-arrestin. The observations provide useful information for understanding the biased activation mechanism.

Activation mechanism of the class D fungal GPCR dimer Ste2

The fungal class D1 G-protein-coupled receptor (GPCR) Ste2 has a different arrangement of transmembrane helices compared with mammalian GPCRs and a distinct mode of coupling to the heterotrimeric G protein Gpa1-Ste2-Ste181. In addition, Ste2 lacks conserved sequence motifs such as DRY, PIF and NPXXY, which are associated with the activation of class A GPCRs2. This suggested that the activation mechanism of Ste2 may also differ. Here we determined structures of Saccharomyces cerevisiae Ste2 in the absence of G protein in two different conformations bound to the native agonist α-factor, bound to an antagonist and without ligand. These structures revealed that Ste2 is indeed activated differently from other GPCRs. In the inactive state, the cytoplasmic end of transmembrane helix H7 is unstructured and packs between helices H1-H6, blocking the G protein coupling site. Agonist binding results in the outward movement of the extracellular ends of H6 and H7 by 6 Å. On the intracellular surface, the G protein coupling site is formed by a 20 Å outward movement of the unstructured region in H7 that unblocks the site, and a 12 Å inward movement of H6. This is a distinct mechanism in GPCRs, in which the movement of H6 and H7 upon agonist binding facilitates G protein coupling.

Ligand-directed bias of G protein signaling at the dopamine D2 receptor

G-protein-coupled receptors (GPCRs) represent the largest family of drug targets. Upon activation, GPCRs signal primarily via a diverse set of heterotrimeric G proteins. Most GPCRs can couple to several different G protein subtypes. However, how drugs act at GPCRs contributing to the selectivity of G protein recognition is poorly understood. Here, we examined the G protein selectivity profile of the dopamine D₂ receptor (D₂), a GPCR targeted by antipsychotic drugs. We show that D₂ discriminates between six individual members of the G_i/o family, and its profile of functional selectivity is remarkably different across its ligands, which all engaged D₂ with a distinct G protein coupling pattern. Using structural modeling, receptor mutagenesis, and pharmacological evaluation, we identified residues in the D₂ binding pocket that shape these ligand-directed biases. We further provide pharmacogenomic evidence that natural variants in D₂ differentially affect its G protein biases in response to different ligands.

Easily Established and Multifunctional Synthetic Nanobody Libraries as Research Tools

Nanobodies, or VHHs, refer to the antigen-binding domain of heavy-chain antibodies (HCAbs) from camelids. They have been widely used as research tools for protein purification and structure determination due to their small size, high specificity, and high stability, overcoming limitations with conventional antibody fragments. However, animal immunization and subsequent retrieval of antigen-specific nanobodies are expensive and complicated. Construction of synthetic nanobody libraries using DNA oligonucleotides is a cost-effective alternative for immunization libraries and shows great potential in identifying antigen-specific or even conformation-specific nanobodies. This review summarizes and analyses synthetic nanobody libraries in the current literature, including library design and biopanning methods, and further discusses applications of antigen-specific nanobodies obtained from synthetic libraries to research.

Crystal structure of the α1B-adrenergic receptor reveals molecular determinants of selective ligand recognition

α-adrenergic receptors (αARs) are G protein-coupled receptors that regulate vital functions of the cardiovascular and nervous systems. The therapeutic potential of αARs, however, is largely unexploited and hampered by the scarcity of subtype-selective ligands. Moreover, several aminergic drugs either show off-target binding to αARs or fail to interact with the desired subtype. Here, we report the crystal structure of human α_1BAR bound to the inverse agonist (+)-cyclazosin, enabled by the fusion to a DARPin crystallization chaperone. The α_1BAR structure allows the identification of two unique secondary binding pockets. By structural comparison of α_1BAR with α₂ARs, and by constructing α_1BAR-α_2CAR chimeras, we identify residues 3.29 and 6.55 as key determinants of ligand selectivity. Our findings provide a basis for discovery of α_1BAR-selective ligands and may guide the optimization of aminergic drugs to prevent off-target binding to αARs, or to elicit a selective interaction with the desired subtype.

This study reports the X-ray structure of the α_1B-adrenergic G protein-coupled receptor bound to an inverse agonist, and unveils key determinants of subtype-selective ligand binding that may help the design of aminergic drugs with fewer side-effects.

Free-Energy Simulations Support a Lipophilic Binding Route for Melatonin Receptors

The effects of the neurohormone melatonin are mediated by the activation of the GPCRs MT₁ and MT₂ in a variety of tissues. Crystal structures suggest ligand access to the orthosteric binding site of MT₁ and MT₂ receptors through a lateral channel between transmembrane (TM) helices IV and V. We investigated the feasibility of this lipophilic entry route for 2-iodomelatonin, a nonselective agonist with a slower dissociation rate from the MT₂ receptor, applying enhanced sampling simulations and free-energy calculations. 2-Iodomelatonin unbinding was investigated with steered molecular dynamics simulations which revealed different trajectories passing through the gap between TM helices IV and V for both receptors. For one of these unbinding trajectories from the MT₁ receptor, an umbrella-sampling protocol with path-collective variables provided a calculated energy barrier consistent with the experimental dissociation rate. The side-chain flexibility of Tyr5.38 was significantly different in the two receptor subtypes, as assessed by metadynamics simulations, and during ligand unbinding it frequently assumes an open conformation in the MT₁ but not in the MT₂ receptor, favoring 2-iodomelatonin egress. Taken together, our simulations are consistent with the possibility that the gap between TM IV and V is a way of connecting the orthosteric binding site and the membrane core for lipophilic melatonin receptor ligands. Our simulations also suggest that the open state of Tyr5.38 generates a small pocket on the surface of MT₁ receptor, which could participate in the recognition of MT₁-selective ligands and may be exploited in the design of new selective compounds.

Vibrational spectroscopy analysis of ligand efficacy in human M2 muscarinic acetylcholine receptor (M2R)

The intrinsic efficacy of ligand binding to G protein-coupled receptors (GPCRs) reflects the ability of the ligand to differentially activate its receptor to cause a physiological effect. Here we use attenuated total reflection-Fourier transform infrared (ATR-FTIR) spectroscopy to examine the ligand-dependent conformational changes in the human M2 muscarinic acetylcholine receptor (M2R). We show that different ligands affect conformational alteration appearing at the C=O stretch of amide-I band in M2R. Notably, ATR-FTIR signals strongly correlated with G-protein activation levels in cells. Together, we propose that amide-I band serves as an infrared probe to distinguish the ligand efficacy in M2R and paves the path to rationally design ligands with varied efficacy towards the target GPCR.

Structure determination of GPCRs: cryo-EM compared with X-ray crystallography

G protein-coupled receptors (GPCRs) are the largest single family of cell surface receptors encoded by the human genome and they play pivotal roles in co-ordinating cellular systems throughout the human body, making them ideal drug targets. Structural biology has played a key role in defining how receptors are activated and signal through G proteins and β-arrestins. The application of structure-based drug design (SBDD) is now yielding novel compounds targeting GPCRs. There is thus significant interest from both academia and the pharmaceutical industry in the structural biology of GPCRs as currently only about one quarter of human non-odorant receptors have had their structure determined. Initially, all the structures were determined by X-ray crystallography, but recent advances in electron cryo-microscopy (cryo-EM) now make GPCRs tractable targets for single-particle cryo-EM with comparable resolution to X-ray crystallography. So far this year, 78% of the 99 GPCR structures deposited in the PDB (Jan-Jul 2021) were determined by cryo-EM. Cryo-EM has also opened up new possibilities in GPCR structural biology, such as determining structures of GPCRs embedded in a lipid nanodisc and multiple GPCR conformations from a single preparation. However, X-ray crystallography still has a number of advantages, particularly in the speed of determining many structures of the same receptor bound to different ligands, an essential prerequisite for effective SBDD. We will discuss the relative merits of cryo-EM and X-ray crystallography for the structure determination of GPCRs and the future potential of both techniques.

Partial agonists of dopamine receptors: receptor theory and the dopamine hypothesis of psychosis

This article discusses dopamine partial agonism, which is the main mechanism of action of the psychiatric drugs aripiprazole, brexpiprazole and cariprazine. It outlines the principles of receptor theory and the structure of dopamine receptors; characterises agonists, antagonists and partial agonists; and summarises the dopamine hypothesis of psychosis and the role of dopamine and serotonin in depression.

The pesticide chlorpyrifos promotes obesity by inhibiting diet-induced thermogenesis in brown adipose tissue

Obesity results from a caloric imbalance between energy intake, absorption and expenditure. In both rodents and humans, diet-induced thermogenesis contributes to energy expenditure and involves the activation of brown adipose tissue (BAT). We hypothesize that environmental toxicants commonly used as food additives or pesticides might reduce BAT thermogenesis through suppression of uncoupling protein 1 (UCP1) and this may contribute to the development of obesity. Using a step-wise screening approach, we discover that the organophosphate insecticide chlorpyrifos suppresses UCP1 and mitochondrial respiration in BAT at concentrations as low as 1 pM. In mice housed at thermoneutrality and fed a high-fat diet, chlorpyrifos impairs BAT mitochondrial function and diet-induced thermogenesis, promoting greater obesity, non-alcoholic fatty liver disease (NAFLD) and insulin resistance. This is associated with reductions in cAMP; activation of p38MAPK and AMPK; protein kinases critical for maintaining UCP1 and mitophagy, respectively in BAT. These data indicate that the commonly used pesticide chlorpyrifos, suppresses diet-induced thermogenesis and the activation of BAT, suggesting its use may contribute to the obesity epidemic.

High-mass MALDI-MS unravels ligand-mediated G protein-coupling selectivity to GPCRs

G protein–coupled receptors (GPCRs) are important pharmaceutical targets for the treatment of a broad spectrum of diseases. Upon ligand binding, GPCRs initiate intracellular signaling pathways by interacting with partner proteins. Assays that quantify the interplay between ligand binding and initiation of downstream signaling cascades are critical in the early stages of drug development. We have developed a high-throughput mass spectrometry method to unravel GPCR–protein complex interplay and demonstrated its use with three GPCRs to provide quantitative information about ligand-modulated coupling selectivity. This method provides insights into the molecular details of GPCR interactions and could serve as an approach for discovery of drugs that initiate specific cell-signaling pathways.

G protein–coupled receptors (GPCRs) are important pharmaceutical targets for the treatment of a broad spectrum of diseases. Although there are structures of GPCRs in their active conformation with bound ligands and G proteins, the detailed molecular interplay between the receptors and their signaling partners remains challenging to decipher. To address this, we developed a high-sensitivity, high-throughput matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) method to interrogate the first stage of signal transduction. GPCR–G protein complex formation is detected as a proxy for the effect of ligands on GPCR conformation and on coupling selectivity. Over 70 ligand–GPCR–partner protein combinations were studied using as little as 1.25 pmol protein per sample. We determined the selectivity profile and binding affinities of three GPCRs (rhodopsin, beta-1 adrenergic receptor [β1AR], and angiotensin II type 1 receptor) to engineered Gα-proteins (mGs, mGo, mGi, and mGq) and nanobody 80 (Nb80). We found that GPCRs in the absence of ligand can bind mGo, and that the role of the G protein C terminus in GPCR recognition is receptor-specific. We exemplified our quantification method using β1AR and demonstrated the allosteric effect of Nb80 binding in assisting displacement of nadolol to isoprenaline. We also quantified complex formation with wild-type heterotrimeric Gα_iβγ and β-arrestin-1 and showed that carvedilol induces an increase in coupling of β-arrestin-1 and Gα_iβγ to β1AR. A normalization strategy allows us to quantitatively measure the binding affinities of GPCRs to partner proteins. We anticipate that this methodology will find broad use in screening and characterization of GPCR-targeting drugs.

Cryo-EM structure of the β3-adrenergic receptor reveals the molecular basis of subtype selectivity

The β₃-adrenergic receptor (β₃AR) is predominantly expressed in adipose tissue and urinary bladder and has emerged as an attractive drug target for the treatment of type 2 diabetes, obesity, and overactive bladder (OAB). Here, we report the cryogenic electron microscopy structure of the β₃AR-G_s signaling complex with the selective agonist mirabegron, a first-in-class drug for OAB. Comparison of this structure with the previously reported β₁AR and β₂AR structures reveals a receptor activation mechanism upon mirabegron binding to the orthosteric site. Notably, the narrower exosite in β₃AR creates a perpendicular pocket for mirabegron. Mutational analyses suggest that a combination of both the exosite shape and the amino-acid-residue substitutions defines the drug selectivity of the βAR agonists. Our findings provide a molecular basis for βAR subtype selectivity, allowing the design of more-selective agents with fewer adverse effects.

Purinergic GPCR transmembrane residues involved in ligand recognition and dimerization

We compare the GPCR-ligand interactions and highlight important residues for recognition in purinergic receptors-from both X-ray crystallographic and cryo-EM structures. These include A1 and A2A adenosine receptors, and P2Y1 and P2Y12 receptors that respond to ADP and other nucleotides. These receptors are important drug discovery targets for immune, metabolic and nervous system disorders. In most cases, orthosteric ligands are represented, except for one allosteric P2Y1 antagonist. This review catalogs the residues and regions that engage in contacts with ligands or with other GPCR protomers in dimeric forms. Residues that are in proximity to bound ligands within purinergic GPCR families are correlated. There is extensive conservation of recognition motifs between adenosine receptors, but the P2Y1 and P2Y12 receptors are each structurally distinct in their ligand recognition. Identifying common interaction features for ligand recognition within a receptor class that has multiple structures available can aid in the drug discovery process.

Ligands of Adrenergic Receptors: A Structural Point of View

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

Optimized damping parameters for empirical dispersion corrections to symmetry-adapted perturbation theory

Symmetry-adapted perturbation theory (SAPT) has become an invaluable tool for studying the fundamental nature of non-covalent interactions by directly computing the electrostatics, exchange (steric) repulsion, induction (polarization), and London dispersion contributions to the interaction energy using quantum mechanics. Further application of SAPT is primarily limited by its computational expense, where even its most affordable variant (SAPT0) scales as the fifth power of system size [O(N⁵)] due to the dispersion terms. The algorithmic scaling of SAPT0 is reduced from O(N⁵)→O(N⁴) by replacing these terms with the empirical D3 dispersion correction of Grimme and co-workers, forming a method that may be termed SAPT0-D3. Here, we optimize the damping parameters for the -D3 terms in SAPT0-D3 using a much larger training set than has previously been considered, namely, 8299 interaction energies computed at the complete-basis-set limit of coupled cluster through perturbative triples [CCSD(T)/CBS]. Perhaps surprisingly, with only three fitted parameters, SAPT0-D3 improves on the accuracy of SAPT0, reducing mean absolute errors from 0.61 to 0.49 kcal mol^-1 over the full set of complexes. Additionally, SAPT0-D3 exhibits a nearly 2.5× speedup over conventional SAPT0 for systems with ∼300 atoms and is applied here to systems with up to 459 atoms. Finally, we have also implemented a functional group partitioning of the approach (F-SAPT0-D3) and applied it to determine important contacts in the binding of salbutamol to G-protein coupled β₁-adrenergic receptor in both active and inactive forms. SAPT0-D3 capabilities have been added to the open-source Psi4 software.

Binding pathway determines norepinephrine selectivity for the human β1AR over β2AR

Beta adrenergic receptors (βARs) mediate physiologic responses to the catecholamines epinephrine and norepinephrine released by the sympathetic nervous system. While the hormone epinephrine binds β₁AR and β₂AR with similar affinity, the smaller neurotransmitter norepinephrine is approximately tenfold selective for the β₁AR. To understand the structural basis for this physiologically important selectivity, we solved the crystal structures of the human β₁AR bound to an antagonist carazolol and different agonists including norepinephrine, epinephrine and BI-167107. Structural comparison revealed that the catecholamine-binding pockets are identical between β₁AR and β₂AR, but the extracellular vestibules have different shapes and electrostatic properties. Metadynamics simulations and mutagenesis studies revealed that these differences influence the path norepinephrine takes to the orthosteric pocket and contribute to the different association rates and thus different affinities.

Allosteric communication regulates ligand-specific GPCR activity

G protein-coupled receptors (GPCRs) are membrane-bound proteins that are ubiquitously expressed in many cell types and take part in mediating multiple signaling pathways. GPCRs are dynamic proteins and exist in an equilibrium between an ensemble of conformational states such as inactive and fully active states. This dynamic nature of GPCRs is one of the factors that confers their basal activity even in the absence of any ligand-mediated activation. Ligands selectively bind and stabilize a subset of the conformations from the ensemble leading to a shift in the equilibrium toward the inactive or the active state depending on the nature of the ligand. This ligand-selective effect is achieved through allosteric communication between the ligand binding site and G protein or β-arrestin coupling site. Similarly, the G protein coupling to the receptor exerts the allosteric effect on the ligand binding region leading to increased binding affinity for agonists and decreased affinity for antagonists or inverse agonists. In this review, we enumerate the current state of our understanding of the mechanism of allosteric communication in GPCRs with a specific focus on the critical role of computational methods in delineating the residues involved in allosteric communication. Analyzing allosteric communication mechanism using molecular dynamics simulations has revealed (a) a structurally conserved mechanism of allosteric communication that regulates the G protein coupling, (b) a rational structure-based approach to designing selective ligands, and (c) an approach to designing allosteric GPCR mutants that are either ligand and G protein or β-arrestin selective.

Structural Insights into Ligand-Receptor Interactions Involved in Biased Agonism of G-Protein Coupled Receptors

G protein-coupled receptors (GPCRs) are versatile signaling proteins that mediate complex cellular responses to hormones and neurotransmitters. Ligand directed signaling is observed when agonists, upon binding to the same receptor, trigger significantly different configuration of intracellular events. The current work reviews the structurally defined ligand – receptor interactions that can be related to specific molecular mechanisms of ligand directed signaling across different receptors belonging to class A of GPCRs. Recent advances in GPCR structural biology allow for mapping receptors’ binding sites with residues particularly important in recognition of ligands’ structural features that are responsible for biased signaling. Various studies show particular role of specific residues lining the extended ligand binding domains, biased agonists may alternatively affect their interhelical interactions and flexibility what can be translated into intracellular loop rearrangements. Studies on opioid and angiotensin receptors indicate importance of residues located deeper within the binding cavity and direct interactions with receptor residues linking the ortosteric ligand binding site with the intracellular transducer binding domain. Collection of results across different receptors may suggest elements of common molecular mechanisms which are responsible for passing alternative signals from biased agonists.

Complexes of the neurotensin receptor 1 with small-molecule ligands reveal structural determinants of full, partial, and inverse agonism

Neurotensin receptor 1 (NTSR1) and related G protein-coupled receptors of the ghrelin family are clinically unexploited, and several mechanistic aspects of their activation and inactivation have remained unclear. Enabled by a new crystallization design, we present five new structures: apo-state NTSR1 as well as complexes with nonpeptide inverse agonists SR48692 and SR142948A, partial agonist RTI-3a, and the novel full agonist SRI-9829, providing structural rationales on how ligands modulate NTSR1. The inverse agonists favor a large extracellular opening of helices VI and VII, undescribed so far for NTSR1, causing a constriction of the intracellular portion. In contrast, the full and partial agonists induce a binding site contraction, and their efficacy correlates with the ability to mimic the binding mode of the endogenous agonist neurotensin. Providing evidence of helical and side-chain rearrangements modulating receptor activation, our structural and functional data expand the mechanistic understanding of NTSR1 and potentially other peptidergic receptors.