Adrenaline-activated structure of β2-adrenoceptor stabilized by an engineered nanobody

Distinct binding conformations of epinephrine with α- and β-adrenergic receptors

Agonists targeting α2-adrenergic receptors (ARs) are used to treat diverse conditions, including hypertension, attention-deficit/hyperactivity disorder, pain, panic disorders, opioid and alcohol withdrawal symptoms, and cigarette cravings. These receptors transduce signals through heterotrimeric Gi proteins. Here, we elucidated cryo-EM structures that depict α2A-AR in complex with Gi proteins, along with the endogenous agonist epinephrine or the synthetic agonist dexmedetomidine. Molecular dynamics simulations and functional studies reinforce the results of the structural revelations. Our investigation revealed that epinephrine exhibits different conformations when engaging with α-ARs and β-ARs. Furthermore, α2A-AR and β1-AR (primarily coupled to Gs, with secondary associations to Gi) were compared and found to exhibit different interactions with Gi proteins. Notably, the stability of the epinephrine-α2A-AR-Gi complex is greater than that of the dexmedetomidine-α2A-AR-Gi complex. These findings substantiate and improve our knowledge on the intricate signaling mechanisms orchestrated by ARs and concurrently shed light on the regulation of α-ARs and β-ARs by epinephrine.

SARS-CoV-2 spike protein acts as a β-adrenergic receptor agonist: A potential mechanism for cardiac sequelae of long COVID

BACKGROUND

Currently, pathophysiological mechanisms of post-acute sequelae of coronavirus disease-19-cardiovascular syndrome (PASC-CVS) remain unknown.

METHODS AND RESULTS

Patients with PASC-CVS exhibited significantly higher circulating levels of severe acute respiratory syndrome-coronavirus-2 spike protein S1 than the non-PASC-CVS patients and healthy controls. Moreover, individuals with high plasma spike protein S1 concentrations exhibited elevated heart rates and normalized low frequency, suggesting cardiac β-adrenergic receptor (β-AR) hyperactivity. Microscale thermophoresis (MST) assay revealed that the spike protein bound to β1- and β2-AR, but not to D1-dopamine receptor. These interactions were blocked by β1- and β2-AR blockers. Molecular docking and MST assay of β-AR mutants revealed that the spike protein interacted with the extracellular loop 2 of both β-ARs. In cardiomyocytes, spike protein dose-dependently increased the cyclic adenosine monophosphate production with or without epinephrine, indicating its allosteric effects on β-ARs.

CONCLUSION

Severe acute respiratory syndrome-coronavirus-2 spike proteins act as an allosteric β-AR agonist, leading to cardiac β-AR hyperactivity, thus contributing to PASC-CVS.

Structural basis of psychedelic LSD recognition at dopamine D1 receptor

Understanding the kinetics of LSD in receptors and subsequent induced signaling is crucial for comprehending both the psychoactive and therapeutic effects of LSD. Despite extensive research on LSD's interactions with serotonin 2A and 2B receptors, its behavior on other targets, including dopamine receptors, has remained elusive. Here, we present cryo-EM structures of LSD/PF6142-bound dopamine D1 receptor (DRD1)-legobody complexes, accompanied by a β-arrestin-mimicking nanobody, NBA3, shedding light on the determinants of G protein coupling versus β-arrestin coupling. Structural analysis unveils a distinctive binding mode of LSD in DRD1, particularly with the ergoline moiety oriented toward TM4. Kinetic investigations uncover an exceptionally rapid dissociation rate of LSD in DRD1, attributed to the flexibility of extracellular loop 2 (ECL2). Moreover, G protein can stabilize ECL2 conformation, leading to a significant slowdown in ligand's dissociation rate. These findings establish a solid foundation for further exploration of G protein-coupled receptor (GPCR) dynamics and their relevance to signal transduction.

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

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

Intracellular Allosteric Antagonist of the Olfactory Receptor OR51E2

Olfactory receptors are members of class A (rhodopsin-like) family of G protein-coupled receptors (GPCRs). Their expression and function have been increasingly studied in nonolfactory tissues, and many have been identified as potential therapeutic targets. In this manuscript, we focus on the discovery of novel ligands for the olfactory receptor family 51 subfamily E2 (OR51E2). We performed an artificial intelligence-based virtual drug screen of a ∼2.2 million small molecule library. Cell-based functional assay identified compound 80 (C80) as an antagonist and inverse agonist, and detailed pharmacological analysis revealed C80 acts as a negative allosteric modulator by significantly decreasing the agonist efficacy, while having a minimal effect on receptor affinity for agonist. C80 binds to an allosteric binding site formed by a network of nine residues localized in the intracellular parts of transmembrane domains 3, 5, 6, 7, and H8, which also partially overlaps with a G protein binding site. Mutational experiments of residues involved in C80 binding uncovered the significance of the C2406.37 position in blocking the activation-related conformational change and keeping the receptor in the inactive form. Our study provides a mechanistic understanding of the negative allosteric action of C80 on agonist-ctivated OR51E2. We believe the identification of the antagonist of OR51E2 will enable a multitude of studies aiming to determine the functional role of this receptor in specific biologic processes. SIGNIFICANCE STATEMENT: OR51E2 has been implicated in various biological processes, and its antagonists that can effectively modulate its activity have therapeutic potential. Here we report the discovery of a negative allosteric modulator of OR51E2 and provide a mechanistic understanding of its action. We demonstrate that this modulator has an inhibitory effect on the efficacy of the agonist for the receptor and reveal a network of nine residues that constitute its binding pocket, which also partially overlaps with the G protein binding site.

Fast, accurate ranking of engineered proteins by target-binding propensity using structure modeling

Deep-learning-based methods for protein structure prediction have achieved unprecedented accuracy, yet their utility in the engineering of protein-based binders remains constrained due to a gap between the ability to predict the structures of candidate proteins and the ability toprioritize proteins by their potential to bind to a target. To bridge this gap, we introduce Automated Pairwise Peptide-Receptor Analysis for Screening Engineered proteins (APPRAISE), a method for predicting the target-binding propensity of engineered proteins. After generating structural models of engineered proteins competing for binding to a target using an established structure prediction tool such as AlphaFold-Multimer or ESMFold, APPRAISE performs a rapid (under 1 CPU second per model) scoring analysis that takes into account biophysical and geometrical constraints. As proof-of-concept cases, we demonstrate that APPRAISE can accurately classify receptor-dependent vs. receptor-independent adeno-associated viral vectors and diverse classes of engineered proteins such as miniproteins targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike, nanobodies targeting a G-protein-coupled receptor, and peptides that specifically bind to transferrin receptor or programmed death-ligand 1 (PD-L1). APPRAISE is accessible through a web-based notebook interface using Google Colaboratory (https://tiny.cc/APPRAISE). With its accuracy, interpretability, and generalizability, APPRAISE promises to expand the utility of protein structure prediction and accelerate protein engineering for biomedical applications.

Can the heptapeptide ASSIVSF of the β2-adrenoceptor recognize ephedrine and pseudoephedrine epimers in a complex system?

Epimer separation is crucial in the field of analytical chemistry, separation science, and the pharmaceutical industry. No reported methods could separate simultaneously epimers or even isomers and remove other unwanted, co-existing, interfering substances from complex systems like herbal extracts. Herein, we prepared a heptapeptide-modified stationary phase for the separation of 1R,2S-(-)-ephedrine [(-)-Ephe] and 1S,2S-(+)-pseudoephedrine [(+)-Pse] epimers from Ephedra sinica Stapf extract and blood samples. The heptapeptide stationary phase was comprehensively characterized by scanning electron microscopy, X-ray photoelectron spectroscopy, and Fourier transform infrared spectroscopy. The separation efficiency of the heptapeptide column was compared with an affinity column packed with full-length β2-AR functionalized silica gel (β2-AR column). The binding affinity of the heptapeptide with (+)-Pse was 3-fold greater than that with (-)-Ephe. Their binding mechanisms were extensively characterized by chromatographic analysis, ultraviolet spectra, circular dichroism analysis, isothermal titration calorimetry, and molecule docking. An enhanced hydrogen bonding was clearly observed in the heptapeptide-(+)-Pse complex. Such results demonstrated that the heptapeptide can recognize (+)-Pse and (-)-Ephe epimers in a complex system. This work, we believe, was the first report to simultaneously separate epimers and remove non-specific interfering substances from complex samples. The method was potentially applicable to more challenging sample separation, such as chiral separation from complex systems.

Molecular basis of proton-sensing by G protein-coupled receptors

Three proton-sensing G protein-coupled receptors (GPCRs), GPR4, GPR65, and GPR68, respond to changes in extracellular pH to regulate diverse physiology and are implicated in a wide range of diseases. A central challenge in determining how protons activate these receptors is identifying the set of residues that bind protons. Here, we determine structures of each receptor to understand the spatial arrangement of putative proton sensing residues in the active state. With a newly developed deep mutational scanning approach, we determined the functional importance of every residue in proton activation for GPR68 by generating ~9,500 mutants and measuring effects on signaling and surface expression. This unbiased screen revealed that, unlike other proton-sensitive cell surface channels and receptors, no single site is critical for proton recognition in GPR68. Instead, a network of titratable residues extend from the extracellular surface to the transmembrane region and converge on canonical class A GPCR activation motifs to activate proton-sensing GPCRs. More broadly, our approach integrating structure and unbiased functional interrogation defines a new framework for understanding the rich complexity of GPCR signaling.

Structure and dynamics determine G protein coupling specificity at a class A GPCR

G protein coupled receptors (GPCRs) exhibit varying degrees of selectivity for different G protein isoforms. Despite the abundant structures of GPCR-G protein complexes, little is known about the mechanism of G protein coupling specificity. The β2-adrenergic receptor is an example of GPCR with high selectivity for Gαs, the stimulatory G protein for adenylyl cyclase, and much weaker for the Gαi family of G proteins inhibiting adenylyl cyclase. By developing a new Gαi-biased agonist (LM189), we provide structural and biophysical evidence supporting that distinct conformations at ICL2 and TM6 are required for coupling of the different G protein subtypes Gαs and Gαi. These results deepen our understanding of G protein specificity and bias and can accelerate the design of ligands that select for preferred signaling pathways.

Mechano-sensitivity of β2-adrenoceptors enhances constitutive activation of cAMP generation that is inhibited by inverse agonists

The concept of agonist-independent signalling that can be attenuated by inverse agonists is a fundamental element of the cubic ternary complex model of G protein-coupled receptor (GPCR) activation. This model shows how a GPCR can exist in two conformational states in the absence of ligands; an inactive R state and an active R* state that differ in their affinities for agonists, inverse agonists, and G-protein alpha subunits. The proportion of R* receptors that exist in the absence of agonists determines the level of constitutive receptor activity. In this study we demonstrate that mechanical stimulation can induce β2-adrenoceptor agonist-independent Gs-mediated cAMP signalling that is sensitive to inhibition by inverse agonists such as ICI-118551 and propranolol. The size of the mechano-sensitive response is dependent on the cell surface receptor expression level in HEK293G cells, is still observed in a ligand-binding deficient D113A mutant β2-adrenoceptor and can be attenuated by site-directed mutagenesis of the extracellular N-glycosylation sites on the N-terminus and second extracellular loop of the β2-adrenoceptor. Similar mechano-sensitive agonist-independent responses are observed in HEK293G cells overexpressing the A2A-adenosine receptor. These data provide new insights into how agonist-independent constitutive receptor activity can be enhanced by mechanical stimulation and regulated by inverse agonists.

GPR161 structure uncovers the redundant role of sterol-regulated ciliary cAMP signaling in the Hedgehog pathway

The orphan G protein-coupled receptor (GPCR) GPR161 plays a central role in development by suppressing Hedgehog signaling. The fundamental basis of how GPR161 is activated remains unclear. Here, we determined a cryogenic-electron microscopy structure of active human GPR161 bound to heterotrimeric Gs. This structure revealed an extracellular loop 2 that occupies the canonical GPCR orthosteric ligand pocket. Furthermore, a sterol that binds adjacent to transmembrane helices 6 and 7 stabilizes a GPR161 conformation required for Gs coupling. Mutations that prevent sterol binding to GPR161 suppress Gs-mediated signaling. These mutants retain the ability to suppress GLI2 transcription factor accumulation in primary cilia, a key function of ciliary GPR161. By contrast, a protein kinase A-binding site in the GPR161 C terminus is critical in suppressing GLI2 ciliary accumulation. Our work highlights how structural features of GPR161 interface with the Hedgehog pathway and sets a foundation to understand the role of GPR161 function in other signaling pathways.

Computational Modeling and Characterization of Peptides Derived from Nanobody Complementary-Determining Region 2 (CDR2) Targeting Active-State Conformation of the β2-Adrenergic Receptor (β2AR)

This study assessed the suitability of the complementarity-determining region 2 (CDR2) of the nanobody (Nb) as a template for the derivation of nanobody-derived peptides (NDPs) targeting active-state β2-adrenergic receptor (β2AR) conformation. Sequences of conformationally selective Nbs favoring the agonist-occupied β2AR were initially analyzed by the informational spectrum method (ISM). The derived NDPs in complex with β2AR were subjected to protein-peptide docking, molecular dynamics (MD) simulations, and metadynamics-based free-energy binding calculations. Computational analyses identified a 25-amino-acid-long CDR2-NDP of Nb71, designated P4, which exhibited the following binding free-energy for the formation of the β2AR:P4 complex (ΔG = -6.8 ± 0.8 kcal/mol or a Ki = 16.5 μM at 310 K) and mapped the β2AR:P4 amino acid interaction network. In vitro characterization showed that P4 (i) can cross the plasma membrane, (ii) reduces the maximum isoproterenol-induced cAMP level by approximately 40% and the isoproterenol potency by up to 20-fold at micromolar concentration, (iii) has a very low affinity to interact with unstimulated β2AR in the cAMP assay, and (iv) cannot reduce the efficacy and potency of the isoproterenol-mediated β2AR/β-arrestin-2 interaction in the BRET2-based recruitment assay. In summary, the CDR2-NDP, P4, binds preferentially to agonist-activated β2AR and disrupts Gαs-mediated signaling.

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.

GPR101: Modeling a constitutively active receptor linked to X-linked acrogigantism

GPR101 is a G protein-coupled receptor (GPCR) implicated in a rare form of genetic gigantism known as X-linked acrogigantism, or X-LAG. In particular, X-LAG patients harbor microduplications in the long arm of the X-chromosome that invariably include the GPR101 gene. Duplications of the GPR101 gene lead to the formation of a new chromatin domain that causes over-expression of the receptor in the pituitary tumors of the patients. Notably, GPR101 is a constitutively active receptor, which stimulates cells to produce the second messenger cyclic AMP (cAMP) in the absence of ligands. Moreover, GPR101 was recently reported to constitutively activate not only the cAMP pathway via Gs, but also other G protein subunits (Gq/11 and G12/13). Hence, chemicals that block the constitutive activity of GPR101, known as inverse agonists, have the potential to be useful for the development of pharmacological tools for the treatment of X-LAG. In this study, we provide structural insights into the putative structure of GPR101 based on in-house built homology models, as well as third party models based on the machine learning methods AlphaFold and AlphaFold-Multistate. Moreover, we report a molecular dynamics study, meant to further probe the constitutive activity of GPR101. Finally, we provide a structural comparison with the closest GPCRs, which suggests that GPR101 does not share their natural ligands. While this manuscript was under review, cryo-electron microscopy structures of GPR101 were reported. These structures are expected to enable computer-aided ligand discovery efforts targeting GPR101.

CMV-encoded GPCRs in infection, disease, and pathogenesis

G protein coupled receptors (GPCRs) are seven-transmembrane domain proteins that modulate cellular processes in response to external stimuli. These receptors represent the largest family of membrane proteins, and in mammals, their signaling regulates important physiological functions, such as vision, taste, and olfaction. Many organisms, including yeast, slime molds, and viruses encode GPCRs. Cytomegaloviruses (CMVs) are large, betaherpesviruses, that encode viral GPCRs (vGPCRs). Human CMV (HCMV) encodes four vGPCRs, including UL33, UL78, US27, and US28. Each of these vGPCRs, as well as their rodent and primate orthologues, have been investigated for their contributions to viral infection and disease. Herein, we discuss how the CMV vGPCRs function during lytic and latent infection, as well as our understanding of how they impact viral pathogenesis.

Mechanisms of biased agonism by Gαi/o-biased stapled peptide agonists of the relaxin-3 receptor

The neuropeptide relaxin-3 is composed of an A chain and a B chain held together by disulfide bonds, and it modulates functions such as anxiety and food intake by binding to and activating its cognate receptor RXFP3, mainly through the B chain. Biased ligands of RXFP3 would help to determine the molecular mechanisms underlying the activation of G proteins and β-arrestins downstream of RXFP3 that lead to such diverse functions. We showed that the i, i+4 stapled relaxin-3 B chains, 14s18 and d(1-7)14s18, were Gαi/o-biased agonists of RXFP3. These peptides did not induce recruitment of β-arrestin1/2 to RXFP3 by GPCR kinases (GRKs), in contrast to relaxin-3, which enabled the GRK2/3-mediated recruitment of β-arrestin1/2 to RXFP3. Relaxin-3 and the previously reported peptide 4 (an i, i+4 stapled relaxin-3 B chain) did not exhibit biased signaling. The staple linker of peptide 4 and parts of both the A chain and B chain of relaxin-3 interacted with extracellular loop 3 (ECL3) of RXFP3, moving it away from the binding pocket, suggesting that unbiased ligands promote a more open conformation of RXFP3. These findings highlight roles for the A chain and the N-terminal residues of the B chain of relaxin-3 in inducing conformational changes in RXFP3, which will help in designing selective biased ligands with improved therapeutic efficacy.

The chemokine receptor CCR5: multi-faceted hook for HIV-1

Chemokines are cytokines whose primary role is cellular activation and stimulation of leukocyte migration. They perform their various functions by interacting with G protein-coupled cell surface receptors (GPCRs) and are involved in the regulation of many biological processes such as apoptosis, proliferation, angiogenesis, hematopoiesis or organogenesis. They contribute to the maintenance of the homeostasis of lymphocytes and coordinate the function of the immune system. However, chemokines and their receptors are sometimes hijacked by some pathogens to infect the host organism. For a given chemokine receptor, there is a wide structural, organizational and conformational diversity. In this review, we describe the evidence for structural variety reported for the chemokine receptor CCR5, how this variability can be exploited by HIV-1 to infect its target cells and what therapeutic solutions are currently being developed to overcome this problem.

Applications and challenges in designing VHH-based bispecific antibodies: leveraging machine learning solutions

The development of bispecific antibodies that bind at least two different targets relies on bringing together multiple binding domains with different binding properties and biophysical characteristics to produce a drug-like therapeutic. These building blocks play an important role in the overall quality of the molecule and can influence many important aspects from potency and specificity to stability and half-life. Single-domain antibodies, particularly camelid-derived variable heavy domain of heavy chain (VHH) antibodies, are becoming an increasingly popular choice for bispecific construction due to their single-domain modularity, favorable biophysical properties, and potential to work in multiple antibody formats. Here, we review the use of VHH domains as building blocks in the construction of multispecific antibodies and the challenges in creating optimized molecules. In addition to exploring traditional approaches to VHH development, we review the integration of machine learning techniques at various stages of the process. Specifically, the utilization of machine learning for structural prediction, lead identification, lead optimization, and humanization of VHH antibodies.

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.

Covalent crosslinking in gas-phase biomolecular ions. An account and perspective

Photochemical crosslinking in gas-phase ion complexes has been introduced as a method to study biomolecular structures and dynamics. Emphasis has been on carbene-based crosslinking induced by photodissociation of diazirine-tagged ions. The features that characterize gas-phase crosslinking include (1) complex formation in electrospray droplets that allows for library-type screening; (2) well defined stoichiometry of the complexes due to mass-selective isolation; (3) facile reaction monitoring and yield determination, and (4) post-crosslinking structure analysis by tandem mass spectrometry that has been combined with hydrogen-deuterium exchange, UV-vis action spectroscopy, and ion mobility measurements. In this account, examples are given of peptide-peptide, peptide-nucleotide, and peptide-ligand crosslinking that chiefly used carbene-based reactions. The pros and cons of gas-phase crosslinking are discussed. Nitrile-imine based crosslinking in gas-phase ions is introduced as a promising new approach to ion structure analysis that offers high efficiency and has the potential for wide ranging applications.

Coevolution-Driven Method for Efficiently Simulating Conformational Changes in Proteins Reveals Molecular Details of Ligand Effects in the β2AR Receptor

With the advent of AI-powered structure prediction, the scientific community is inching closer to solving protein folding. An unresolved enigma, however, is to accurately, reliably, and deterministically predict alternative conformational states that are crucial for the function of, e.g., transporters, receptors, or ion channels where conformational cycling is innately coupled to protein function. Accurately discovering and exploring all conformational states of membrane proteins has been challenging due to the need to retain atomistic detail while enhancing the sampling along interesting degrees of freedom. The challenges include but are not limited to finding which degrees of freedom are relevant, how to accelerate the sampling along them, and then quantifying the populations of each micro- and macrostate. In this work, we present a methodology that finds relevant degrees of freedom by combining evolution and physics through machine learning and apply it to the conformational sampling of the β2 adrenergic receptor. In addition to predicting new conformations that are beyond the training set, we have computed free energy surfaces associated with the protein's conformational landscape. We then show that the methodology is able to quantitatively predict the effect of an array of ligands on the β2 adrenergic receptor activation through the discovery of new metastable states not present in the training set. Lastly, we also stake out the structural determinants of activation and inactivation pathway signaling through different ligands and compare them to functional experiments to validate our methodology and potentially gain further insights into the activation mechanism of the β2 adrenergic receptor.

Selective activation of intracellular β1AR using a spatially restricted antagonist

G-protein-coupled receptors (GPCRs) regulate several physiological and pathological processes and represent the target of approximately 30% of FDA-approved drugs. GPCR-mediated signaling was thought to occur exclusively at the plasma membrane. However, recent studies have unveiled their presence and function at subcellular membrane compartments. There is a growing interest in studying compartmentalized signaling of GPCRs. This requires development of novel tools to separate GPCRs signaling at the plasma membrane from the ones initiated at intracellular compartments. We took advantage of the structural and pharmacological information available for β1-adrenergic receptor (β1AR), an exemplary GPCR that functions at subcellular compartments, and rationally designed spatially restricted antagonists. We generated a cell impermeable β1AR antagonist by conjugating a suitable pharmacophore to a sulfonate-containing fluorophore. This cell-impermeable antagonist only inhibited β1AR on the plasma membrane. In contrast, a cell permeable β1AR agonist containing a non-sulfonated fluorophore, efficiently inhibited both the plasma membrane and Golgi pools of β1ARs. Furthermore, the cell impermeable antagonist selectively inhibited the phosphorylation of downstream effectors of PKA proximal to the plasma membrane in adult cardiomyocytes while β1AR intracellular pool remained active. Our tools offer promising avenues for investigating compartmentalized β1AR signaling in various context, potentially advancing our understanding of β1AR-mediated cellular responses in health and disease. They also offer a general strategy to study compartmentalized signaling for other GPCRs in various biological systems.

Molecular dynamics and free energy calculations of clozapine bound to D2 and H1 receptors reveal a cardiometabolic mitigated derivative

Most atypical antipsychotics derive from a high dropout of drug treatments due to adverse cardiometabolic side effects. These side effects are caused, in part, by the H1 receptor blockade. The current work sought a clozapine derivative with a reduced affinity for the H1 receptor while maintaining its therapeutic effect linked to D2 receptor binding. Explicit molecular dynamics simulations and end-point free energy calculations of clozapine in complex with the D2 and H1 receptors embedded in cholesterol-rich lipid bilayers were performed to analyze the intermolecular interactions and address the relevance of clozapine-functional groups. Based on that, free energy perturbation calculations were performed to measure the change in free energy of clozapine structural modifications. Our results indicate the best clozapine derivative is the iodine atom substitution for chlorine. The latter is mainly due to electrostatic interaction loss for the H1 receptor, while the halogen orientation out of the D2 active site reduces the impact on the affinity.Communicated by Ramaswamy H. Sarma.

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.

Accessible and Generalizable in Vitro Luminescence Assay for Detecting GPCR Activation

G protein-coupled receptors (GPCRs) serve critical physiological roles as the most abundant family of receptors. Here, we describe the design of a generalizable and cell lysate-based method that leverages the interaction between an agonist-activated GPCR and a conformation-specific binder to reconstitute split nanoluciferase (NanoLuc) in vitro. This tool, In vitro GPCR split NanoLuc ligand Triggered Reporter (IGNiTR), has broad applications. We have demonstrated IGNiTR's use with three Gs-coupled GPCRs, two Gi-coupled GPCRs and three classes of conformation-specific binders: nanobodies, miniG proteins, and G protein peptidomimetics. As an in vitro method, IGNiTR enables the use of synthetic G protein peptidomimetics and provides easily scalable and portable reagents for characterizing GPCRs and ligands. We tested three diverse applications of IGNiTR: (1) proof-of-concept GPCR ligand screening using dopamine receptor D1 IGNiTR; (2) detection of opioids for point-of-care testing; and (3) characterizing GPCR functionality during Nanodisc-based reconstitution processes. Due to IGNiTR's unique advantages and the convenience of its cell lysate-based format, this tool will find extensive applications in GPCR ligand detection, screening, and GPCR characterization.

Interaction of Synthetic Cannabinoid Receptor Agonists with Cannabinoid Receptor I: Insights into Activation Molecular Mechanism

Synthetic cannabinoid receptor agonists (SCRAs) have become a wide group of new psychoactive substances since the 2010s. For the last few years, the X-ray structures of the complexes of cannabinoid receptor I (CB1) with SCRAs as well as the complexes of CB1 with its antagonist have been published. Based on those data, SCRA-CB1 interactions are analyzed in detail, using molecular modeling and molecular dynamics simulations. The molecular mechanism of the conformational transformation of the transmembrane domain of CB1 caused by its interaction with SCRA is studied. These conformational changes allosterically modulate the CB1-Gi complex, providing activation of the Gi protein. Based on the X-ray-determined structures of the CB1-ligand complexes, a stable apo conformation of inactive CB1 with a relatively low potential barrier of receptor activation was modeled. For that model, molecular dynamic simulations of SCRA binding to CB1 led to the active state of CB1, which allowed us to explore the key features of this activation and the molecular mechanism of the receptor's structural transformation. The simulated CB1 activation is in accordance with the previously published experimental data for the activation at protein mutations or structural changes of ligands. The key feature of the suggested activation mechanism is the determination of the stiff core of the CB1 transmembrane domain and the statement that the entire conformational transformation of the receptor to the active state is caused by a shift of alpha helix TM7 relative to this core. The shift itself is caused by protein-ligand interactions. It was verified via steered molecular dynamics simulations of the X-ray-determined structures of the inactive receptor, which resulted in the active conformation of CB1 irrespective of the placement of agonist ligand in the receptor's active site.

Molecular Basis of Human Trace Amine-Associated Receptor 1 Activation

The human trace amine-associated receptor 1 (hTAAR1, hTA1) is a key regulator of monoaminergic neurotransmission and the actions of psychostimulants. Despite preclinical research demonstrating its tractability as a drug target, its molecular mechanisms of activation remain unclear. Moreover, poorly understood pharmacological differences between rodent and human TA1 complicate the translation of findings from preclinical disease models into novel pharmacotherapies. To elucidate hTA1's mechanisms on the molecular scale and investigate the underpinnings of its divergent pharmacology from rodent orthologs, we herein report the structure of the human TA1 receptor in complex with a Gαs heterotrimer. Our structure reveals shared structural elements with other TAARs, as well as with its closest monoaminergic ortholog, the serotonin receptor 5-HT4R. We further find that a single mutation dramatically shifts the selectivity of hTA1 towards that of its rodent orthologs, and report on the effects of substituting residues to those found in serotonin and dopamine receptors. Strikingly, we also discover that the atypical antipsychotic medication and pan-monoaminergic antagonist asenapine potently and efficaciously activates hTA1. Together our studies provide detailed insight into hTA1 structure and function, contrast its molecular pharmacology with that of related receptors, and uncover off-target activities of monoaminergic drugs at hTA1.

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

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

Structural basis for the ligand recognition and signaling of free fatty acid receptors

Free fatty acid receptors 1-4 (FFA1-4) are class A G protein-coupled receptors (GPCRs). FFA1-3 share substantial sequence similarity whereas FFA4 is unrelated. Despite this FFA1 and FFA4 are activated by the same range of long chain fatty acids (LCFAs) whilst FFA2 and FFA3 are instead activated by short chain fatty acids (SCFAs) generated by the intestinal microbiota. Each of FFA1, 2 and 4 are promising targets for novel drug development in metabolic and inflammatory conditions. To gain insights into the basis of ligand interactions with, and molecular mechanisms underlying activation of, FFAs by LCFAs and SCFAs, we determined the active structures of FFA1 and FFA4 bound to the polyunsaturated LCFA docosahexaenoic acid (DHA), FFA4 bound to the synthetic agonist TUG-891, as well as SCFA butyrate-bound FFA2, each complexed with an engineered heterotrimeric Gq protein (miniGq), by cryo-electron microscopy. Together with computational simulations and mutagenesis studies, we elucidated the similarities and differences in the binding modes of fatty acid ligands with varying chain lengths to their respective GPCRs. Our findings unveil distinct mechanisms of receptor activation and G protein coupling. We anticipate that these outcomes will facilitate structure-based drug development and underpin future research to understand allosteric modulation and biased signaling of this group of GPCRs.

Identification of a Salt Bridge That Is Functionally Important for Chemokine Receptor CXCR1 but not CXCR2

CXC chemokine receptors 1 (CXCR1) and 2 (CXCR2) have high sequence similarity and overlapping chemokine ligand profiles. Residue positions 3.32 and 7.39 are critical for signal transduction in the related CXCR4, and in these positions CXCR1 and CXCR2 contain oppositely charged residues (Lys3.32 and Glu7.39). Experimental and computed receptor structures reveal the possible formation of a salt bridge between transmembrane (TM) helices 3 and 7 via these two residues. To investigate the functional importance of Lys1173.32 and Glu2917.39 in CXCR1, along with the flanking Glu1183.33, we performed a signaling study on 16 CXCR1 mutants using two different CXCL8 isoforms. While single Ala-mutation (K1173.32A, E2917.39A) and charge reversal (K1173.32E, E2917.39K) resulted in nonfunctional receptors, double (K1173.32E-E2917.39K) and triple (K1173.32E-E1183.33A-E2917.39K) mutants rescued CXCR1 function. In contrast, the corresponding mutations did not affect the CXCR2 function to the same extent. Our findings show that the Lys3.32-Glu7.39 salt bridge between TM3 and -7 is functionally important for CXCR1 but not for CXCR2, meaning that signal transduction for these highly homologous receptors is not conserved.

Structural basis of agonist specificity of α1A-adrenergic receptor

α1-adrenergic receptors (α1-ARs) play critical roles in the cardiovascular and nervous systems where they regulate blood pressure, cognition, and metabolism. However, the lack of specific agonists for all α1 subtypes has limited our understanding of the physiological roles of different α1-AR subtypes, and led to the stagnancy in agonist-based drug development for these receptors. Here we report cryo-EM structures of α1A-AR in complex with heterotrimeric G-proteins and either the endogenous common agonist epinephrine or the α1A-AR-specific synthetic agonist A61603. These structures provide molecular insights into the mechanisms underlying the discrimination between α1A-AR and α1B-AR by A61603. Guided by the structures and corresponding molecular dynamics simulations, we engineer α1A-AR mutants that are not responsive to A61603, and α1B-AR mutants that can be potently activated by A61603. Together, these findings advance our understanding of the agonist specificity for α1-ARs at the molecular level, opening the possibility of rational design of subtype-specific agonists.

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.

The relaxin receptor RXFP1 signals through a mechanism of autoinhibition

The relaxin family peptide receptor 1 (RXFP1) is the receptor for relaxin-2, an important regulator of reproductive and cardiovascular physiology. RXFP1 is a multi-domain G protein-coupled receptor (GPCR) with an ectodomain consisting of a low-density lipoprotein receptor class A (LDLa) module and leucine-rich repeats. The mechanism of RXFP1 signal transduction is clearly distinct from that of other GPCRs, but remains very poorly understood. In the present study, we determine the cryo-electron microscopy structure of active-state human RXFP1, bound to a single-chain version of the endogenous agonist relaxin-2 and the heterotrimeric Gs protein. Evolutionary coupling analysis and structure-guided functional experiments reveal that RXFP1 signals through a mechanism of autoinhibition. Our results explain how an unusual GPCR family functions, providing a path to rational drug development targeting the relaxin receptors.

The dynamics of agonist-β2-adrenergic receptor activation induced by binding of GDP-bound Gs protein

There is considerable uncertainty about the mechanism by which the β₂-adrenergic receptor (β₂AR) is activated. Here we use molecular metadynamics computations to predict the mechanism by which an agonist induces the activation of the β₂AR and its cognate Gs protein. We found that binding agonist alone to the inactive β₂AR does not break the ionic lock and hence does not drive the β₂AR towards the activated conformation. However, we found that attaching the inactive Gs protein to the agonist-bound inactive β₂AR (containing the ionic lock) leads to partial insertion of Gαs-α5 into the core of β₂AR, which breaks the ionic lock, leading to activation of the Gs protein coupled to β₂AR. Upon activation, the Gαs protein undergoes a remarkable opening of the GDP binding pocket, making the GDP available for exchange or release. Concomitantly, Gαs-α5 undergoes a remarkable expansion in the β₂AR cytoplasmic region after the ionic lock is broken, inducing TM6 to displace outward by ~5 Å from TM3.

Structural basis of α1A-adrenergic receptor activation and recognition by an extracellular nanobody

The α_1A-adrenergic receptor (α_1AAR) belongs to the family of G protein-coupled receptors that respond to adrenaline and noradrenaline. α_1AAR is involved in smooth muscle contraction and cognitive function. Here, we present three cryo-electron microscopy structures of human α_1AAR bound to the endogenous agonist noradrenaline, its selective agonist oxymetazoline, and the antagonist tamsulosin, with resolutions range from 2.9 Å to 3.5 Å. Our active and inactive α_1AAR structures reveal the activation mechanism and distinct ligand binding modes for noradrenaline compared with other adrenergic receptor subtypes. In addition, we identified a nanobody that preferentially binds to the extracellular vestibule of α_1AAR when bound to the selective agonist oxymetazoline. These results should facilitate the design of more selective therapeutic drugs targeting both orthosteric and allosteric sites in this receptor family.

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.

GPR161 structure uncovers the redundant role of sterol-regulated ciliary cAMP signaling in the Hedgehog pathway

The orphan G protein-coupled receptor (GPCR) GPR161 is enriched in primary cilia, where it plays a central role in suppressing Hedgehog signaling1. GPR161 mutations lead to developmental defects and cancers2,3,4. The fundamental basis of how GPR161 is activated, including potential endogenous activators and pathway-relevant signal transducers, remains unclear. To elucidate GPR161 function, we determined a cryogenic-electron microscopy structure of active GPR161 bound to the heterotrimeric G protein complex Gs. This structure revealed an extracellular loop 2 that occupies the canonical GPCR orthosteric ligand pocket. Furthermore, we identify a sterol that binds to a conserved extrahelical site adjacent to transmembrane helices 6 and 7 and stabilizes a GPR161 conformation required for Gs coupling. Mutations that prevent sterol binding to GPR161 suppress cAMP pathway activation. Surprisingly, these mutants retain the ability to suppress GLI2 transcription factor accumulation in cilia, a key function of ciliary GPR161 in Hedgehog pathway suppression. By contrast, a protein kinase A-binding site in the GPR161 C-terminus is critical in suppressing GLI2 ciliary accumulation. Our work highlights how unique structural features of GPR161 interface with the Hedgehog pathway and sets a foundation to understand the broader role of GPR161 function in other signaling pathways.

Computational structural-based GPCR optimization for user-defined ligand: Implications for the development of biosensors

Organisms have developed effective mechanisms to sense the external environment. Human-designed biosensors exploit this natural optimization, where different biological machinery have been adapted to detect the presence of user-defined molecules. Specifically, the pheromone pathway in the model organism Saccharomyces cerevisiae represents a suitable candidate as a synthetic signaling system. Indeed, it expresses just one G-Protein Coupled Receptor (GPCR), Ste2, able to recognize pheromone and initiate the expression of pheromone-dependent genes. To date, the standard procedure to engineer this system relies on the substitution of the yeast GPCR with another one and on the modification of the yeast G-protein to bind the inserted receptor. Here, we propose an innovative computational procedure, based on geometrical and chemical optimization of protein binding pockets, to select the amino acid substitutions required to make the native yeast GPCR able to recognize a user-defined ligand. This procedure would allow the yeast to recognize a wide range of ligands, without a-priori knowledge about a GPCR recognizing them or the corresponding G protein. We used Monte Carlo simulations to design on Ste2 a binding pocket able to recognize epinephrine, selected as a test ligand. We validated Ste2 mutants via molecular docking and molecular dynamics. We verified that the amino acid substitutions we identified make Ste2 able to accommodate and remain firmly bound to epinephrine. Our results indicate that we sampled efficiently the huge space of possible mutants, proposing such a strategy as a promising starting point for the development of a new kind of S.cerevisiae-based biosensors.

Cryo-EM structure of the endothelin-1-ETB-Gi complex

The endothelin ET_B receptor is a promiscuous G-protein coupled receptor that is activated by vasoactive peptide endothelins. ET_B signaling induces reactive astrocytes in the brain and vasorelaxation in vascular smooth muscle. Consequently, ET_B agonists are expected to be drugs for neuroprotection and improved anti-tumor drug delivery. Here, we report the cryo-electron microscopy structure of the endothelin-1-ET_B-G_i complex at 2.8 Å resolution, with complex assembly stabilized by a newly established method. Comparisons with the inactive ET_B receptor structures revealed how endothelin-1 activates the ET_B receptor. The NPxxY motif, essential for G-protein activation, is not conserved in ET_B, resulting in a unique structural change upon G-protein activation. Compared with other GPCR-G-protein complexes, ET_B binds G_i in the shallowest position, further expanding the diversity of G-protein binding modes. This structural information will facilitate the elucidation of G-protein activation and the rational design of ET_B agonists.

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.

Function and dynamics of the intrinsically disordered carboxyl terminus of β2 adrenergic receptor

Advances in structural biology have provided important mechanistic insights into signaling by the transmembrane core of G-protein coupled receptors (GPCRs); however, much less is known about intrinsically disordered regions such as the carboxyl terminus (CT), which is highly flexible and not visible in GPCR structures. The β2 adrenergic receptor's (β2AR) 71 amino acid CT is a substrate for GPCR kinases and binds β-arrestins to regulate signaling. Here we show that the β2AR CT directly inhibits basal and agonist-stimulated signaling in cell lines lacking β-arrestins. Combining single-molecule fluorescence resonance energy transfer (FRET), NMR spectroscopy, and molecular dynamics simulations, we reveal that the negatively charged β2AR-CT serves as an autoinhibitory factor via interacting with the positively charged cytoplasmic surface of the receptor to limit access to G-proteins. The stability of this interaction is influenced by agonists and allosteric modulators, emphasizing that the CT plays important role in allosterically regulating GPCR activation.

A comprehensive pharmacological analysis of fenoterol and its derivatives to unravel the role of β2-adrenergic receptor in zebrafish

β-adrenergic receptors (βARs) belong to a key molecular targets that regulate the most important processes occurring in the human organism. Although over the last decades a zebrafish model has been developed as a model complementary to rodents in biomedical research, the role of β₂AR in regulation of pathological and toxicological effects remains to elucidate. Therefore, the study aimed to clarify the role of β₂AR with a particular emphasis on the distinct role of subtypes A and B of zebrafish β₂AR. As model compounds selective β₂AR agonists - (R,R)-fenoterol ((R,R)-Fen) and its new derivatives: (R,R)-4'-methoxyfenoterol ((R,R)-MFen) and (R,R)-4'-methoxy-1-naphtylfenoterol ((R,R)-MNFen) - were tested. We described dose-dependent changes observed after fenoterols exposure in terms of general toxicity, cardiotoxicity and neurobehavioural responses. Subsequently, to better characterise the role of β₂-adrenergic stimulation in zebrafish, we have performed a series of molecular docking simulations. Our results indicate that (R,R)-Fen displays the highest affinity for subtype A of zebrafish β₂AR and β_2AAR might be involved in pigment depletion. (R,R)-MFen shows the lowest affinity for zebrafish β₂ARs out of the tested fenoterols and this might be associated with its cardiotoxic and anxiogenic effects. (R,R)-MNFen displays the highest affinity for subtype B of zebrafish β₂AR and modulation of this receptor might be associated with the development of malformations, increases locomotor activity and induces a negative chronotropic effect. Taken together, the presented data offer insights into the functional responses of the zebrafish β₂ARs confirming their intraspecies conservation, and support the translation of the zebrafish model in pharmacological and toxicological research.

Allosteric Regulation of G-Protein-Coupled Receptors: From Diversity of Molecular Mechanisms to Multiple Allosteric Sites and Their Ligands

Allosteric regulation is critical for the functioning of G protein-coupled receptors (GPCRs) and their signaling pathways. Endogenous allosteric regulators of GPCRs are simple ions, various biomolecules, and protein components of GPCR signaling (G proteins and β-arrestins). The stability and functional activity of GPCR complexes is also due to multicenter allosteric interactions between protomers. The complexity of allosteric effects caused by numerous regulators differing in structure, availability, and mechanisms of action predetermines the multiplicity and different topology of allosteric sites in GPCRs. These sites can be localized in extracellular loops; inside the transmembrane tunnel and in its upper and lower vestibules; in cytoplasmic loops; and on the outer, membrane-contacting surface of the transmembrane domain. They are involved in the regulation of basal and orthosteric agonist-stimulated receptor activity, biased agonism, GPCR-complex formation, and endocytosis. They are targets for a large number of synthetic allosteric regulators and modulators, including those constructed using molecular docking. The review is devoted to the principles and mechanisms of GPCRs allosteric regulation, the multiplicity of allosteric sites and their topology, and the endogenous and synthetic allosteric regulators, including autoantibodies and pepducins. The allosteric regulation of chemokine receptors, proteinase-activated receptors, thyroid-stimulating and luteinizing hormone receptors, and beta-adrenergic receptors are described in more detail.

Elucidation of a dynamic interplay between a beta-2 adrenergic receptor, its agonist, and stimulatory G protein

G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (β2AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (Gs). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to β2AR by Gs protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from β2AR and the conformational interconversions of Gs between closed and open conformations in the NE(+)-β2AR-Gs ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter Gs α subunit (Gsα) conformational transitions. Our simulations showed that the interdomain movement and the stacking of Gsα α1 and α5 helices are significant for increasing the distance between the Gsα and β2AR, which may indicate a partial dissociation of Gsα The distance increase commences when Gsα is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of β2AR interacting with Gsα, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.

Autoregulation of GPCR signalling through the third intracellular loop

The third intracellular loop (ICL3) of the G protein-coupled receptor (GPCR) fold is important for the signal transduction process downstream of receptor activation1-3. Despite this, the lack of a defined structure of ICL3, combined with its high sequence divergence among GPCRs, complicates characterization of its involvement in receptor signalling4. Previous studies focusing on the β2 adrenergic receptor (β2AR) suggest that ICL3 is involved in the structural process of receptor activation and signalling5-7. Here we derive mechanistic insights into the role of ICL3 in β2AR signalling, observing that ICL3 autoregulates receptor activity through a dynamic conformational equilibrium between states that block or expose the receptor's G protein-binding site. We demonstrate the importance of this equilibrium for receptor pharmacology, showing that G protein-mimetic effectors bias the exposed states of ICL3 to allosterically activate the receptor. Our findings additionally reveal that ICL3 tunes signalling specificity by inhibiting receptor coupling to G protein subtypes that weakly couple to the receptor. Despite the sequence diversity of ICL3, we demonstrate that this negative G protein-selection mechanism through ICL3 extends to GPCRs across the superfamily, expanding the range of known mechanisms by which receptors mediate G protein subtype selective signalling. Furthermore, our collective findings suggest ICL3 as an allosteric site for receptor- and signalling pathway-specific ligands.

Structural basis of odorant recognition by a human odorant receptor

Our sense of smell enables us to navigate a vast space of chemically diverse odour molecules. This task is accomplished by the combinatorial activation of approximately 400 odorant G protein-coupled receptors encoded in the human genome1-3. How odorants are recognized by odorant receptors remains unclear. Here we provide mechanistic insight into how an odorant binds to a human odorant receptor. Using cryo-electron microscopy, we determined the structure of the active human odorant receptor OR51E2 bound to the fatty acid propionate. Propionate is bound within an occluded pocket in OR51E2 and makes specific contacts critical to receptor activation. Mutation of the odorant-binding pocket in OR51E2 alters the recognition spectrum for fatty acids of varying chain length, suggesting that odorant selectivity is controlled by tight packing interactions between an odorant and an odorant receptor. Molecular dynamics simulations demonstrate that propionate-induced conformational changes in extracellular loop 3 activate OR51E2. Together, our studies provide a high-resolution view of chemical recognition of an odorant by a vertebrate odorant receptor, providing insight into how this large family of G protein-coupled receptors enables our olfactory sense.

Biophysical investigations of class A GPCRs

G protein-coupled receptors (GPCRs) play a central role in cellular communication, converting external stimuli into intracellular responses. GPCRs bind a very broad panel of ligands, such as hormones, neurotransmitters, peptides and lipids. Ligand binding triggers a series of receptor conformational rearrangements, enabling the coupling to intracellular partners and the activation of signaling cascades. The major breakthrough in GPCRs structural biology of the past decade has considerably advanced our understanding of GPCR activation. However, structural information cannot fully explain the molecular details of GPCRs pharmacology. Biophysical investigations reveal that GPCRs are very dynamic proteins, capable of exploring a wide range of conformational states. Binding to ligands of various pharmacological classes, as well as intracellular effectors and allosteric modulators, can shift the equilibrium between these states and the kinetic of interconversions among the different conformers. Investigation of GPCR dynamic interplay is therefore important to better understand the complex pharmacology and signaling profile of these receptors.

Molecular Modeling Study of a Receptor-Orthosteric Ligand-Allosteric Modulator Signaling Complex

Allosteric modulators (AMs) are considered as a perpetual hotspot in research for their higher selectivity and various effects on orthosteric ligands (OL). They are classified in terms of their functionalities as positive, negative, or silent allosteric modulators (PAM, NAM, or SAM, respectively). In the present work, 11 pairs of three-dimensional (3D) structures of receptor-orthosteric ligand and receptor-orthosteric ligand-allosteric modulator complexes have been collected for the studies, including three different systems: GPCR, enzyme, and ion channel. Molecular dynamics (MD) simulations are applied to quantify the dynamic interactions in both the orthosteric and allosteric binding pockets and the structural fluctuation of the involved proteins. Our results showed that MD simulations of moderately large molecules or peptides undergo insignificant changes compared to crystal structure results. Furthermore, we also studied the conformational changes of receptors that bound with PAM and NAM, as well as the different allosteric binding sites in a receptor. There should be no preference for the position of the allosteric binding pocket after comparing the allosteric binding pockets of these three systems. Finally, we aligned four distinct β2 adrenoceptor structures and three N-methyl-d-aspartate receptor (NMDAR) structures to investigate conformational changes. In the β2 adrenoceptor systems, the aligned results revealed that transmembrane (TM) helices 1, 5, and 6 gradually increased outward movement from an enhanced inactive state to an improved active state. TM6 endured the most significant conformational changes (around 11 Å). For NMDAR, the bottom section of NMDAR's ligand-binding domain (LBD) experienced an upward and outward shift during the gradually activating process. In conclusion, our research provides insight into receptor-orthosteric ligand-allosteric modulator studies and the design and development of allosteric modulator drugs using MD simulation.

Structural basis for Y2 receptor-mediated neuropeptide Y and peptide YY signaling

Neuropeptide Y (NPY) and its receptors are expressed in various human tissues including the brain where they regulate appetite and emotion. Upon NPY stimulation, the neuropeptide Y1 and Y2 receptors (Y₁R and Y₂R, respectively) activate G_I signaling, but their physiological responses to food intake are different. In addition, deletion of the two N-terminal amino acids of peptide YY (PYY(3-36)), the endogenous form found in circulation, can stimulate Y₂R but not Y₁R, suggesting that Y₁R and Y₂R may have distinct ligand-binding modes. Here, we report the cryo-electron microscopy structures of the PYY(3-36)‒Y₂R‒G_i and NPY‒Y₂R‒G_i complexes. Using cell-based assays, molecular dynamics simulations, and structural analysis, we revealed the molecular basis of the exclusive binding of PYY(3-36) to Y₂R. Furthermore, we demonstrated that Y₂R favors G protein signaling over β-arrestin signaling upon activation, whereas Y₁R does not show a preference between these two pathways.

Homology Modeling of the G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are therapeutically important family of membrane proteins. Despite growing number of experimental structures available for GPCRs, homology modeling remains a relevant method for studying these receptors and for discovering new ligands for them. Here we describe the state-of-the-art methods for modeling GPCRs, starting from template selection, through fine-tuning sequence alignment to model refinement.