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

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.

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

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

Modeling of Olfactory Receptors

Olfactory receptors (ORs) form the largest subfamily within class A G protein-coupled receptors (GPCRs). No experimental structural data of any OR is available to date. Homology modeling has become a popular strategy to propose plausible OR models, in order to study the structure-function relationships of the receptors and to aid the discovery and development of ligands capable of modulating receptor activity. In this chapter, we provide a general guideline for OR structure construction, including the collection of candidate templates, structure-based sequence alignment, 3D structure construction, ligand docking, and molecular dynamic simulation.

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.

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.

Energy Transport and Its Function in Heptahelical Transmembrane Proteins

Photoproteins such as bacteriorhodopsin (bR) and rhodopsin (Rho) need to effectively dissipate photoinduced excess energy to prevent themselves from damage. Another well-studied seven transmembrane (TM) helices protein is the β₂ adrenergic receptor (β₂AR), a G protein-coupled receptor for which energy dissipation paths have been linked with allosteric communication. To study the vibrational energy transport in the active and inactive states of these proteins, a master equation approach [J. Chem. Phys.2020, 152, 045103] is employed, which uses scaling rules that allow us to calculate energy transport rates solely based on the protein structure. Despite their overall structural similarity, the three 7TM proteins reveal quite different strategies to redistribute excess energy. While bR quickly removes the energy using the TM7 helix as a "lightning rod", Rho exhibits a rather poor energy dissipation, which might eventually require the hydrolysis of the Schiff base between the protein and the retinal chromophore to prevent overheating. Heating the ligand adrenaline of β₂AR, the resulting energy transport network of the protein is found to change significantly upon switching from the active state to the inactive state. While the energy flow may highlight aspects of the inter-residue couplings of β₂AR, it seems not particularly suited to explain allosteric phenomena.

One class classification for the detection of β2 adrenergic receptor agonists using single-ligand dynamic interaction data

G protein-coupled receptors are involved in many biological processes, relaying the extracellular signal inside the cell. Signaling is regulated by the interactions between receptors and their ligands, it can be stimulated by agonists, or inhibited by antagonists or inverse agonists. The development of a new drug targeting a member of this family requires to take into account the pharmacological profile of the designed ligands in order to elicit the desired response. The structure-based virtual screening of chemical libraries may prioritize a specific class of ligands by combining docking results and ligand binding information provided by crystallographic structures. The performance of the method depends on the relevance of the structural data, in particular the conformation of the targeted site, the binding mode of the reference ligand, and the approach used to compare the interactions formed by the docked ligand with those formed by the reference ligand in the crystallographic structure. Here, we propose a new method based on the conformational dynamics of a single protein–ligand reference complex to improve the biased selection of ligands with specific pharmacological properties in a structure-based virtual screening exercise. Interactions patterns between a reference agonist and the receptor, here exemplified on the β2 adrenergic receptor, were extracted from molecular dynamics simulations of the agonist/receptor complex and encoded in graphs used to train a one-class machine learning classifier. Different conditions were tested: low to high affinity agonists, varying simulation duration, considering or ignoring hydrophobic contacts, and tuning of the classifier parametrization. The best models applied to post-process raw data from retrospective virtual screening obtained by docking of test libraries effectively filtered out irrelevant poses, discarding inactive and non-agonist ligands while identifying agonists. Taken together, our results suggest that consistency of the binding mode during the simulation is a key to the success of the method.

Combined docking and machine learning identify key molecular determinants of ligand pharmacological activity on β2 adrenoceptor

G protein‐coupled receptors (GPCRs) are valuable therapeutic targets for many diseases. A central question of GPCR drug discovery is to understand what determines the agonism or antagonism of ligands that bind them. Ligands exert their action via the interactions in the ligand binding pocket. We hypothesized that there is a common set of receptor interactions made by ligands of diverse structures that mediate their action and that among a large dataset of different ligands, the functionally important interactions will be over‐represented. We computationally docked ~2700 known β2AR ligands to multiple β2AR structures, generating ca 75 000 docking poses and predicted all atomic interactions between the receptor and the ligand. We used machine learning (ML) techniques to identify specific interactions that correlate with the agonist or antagonist activity of these ligands. We demonstrate with the application of ML methods that it is possible to identify the key interactions associated with agonism or antagonism of ligands. The most representative interactions for agonist ligands involve K97^2.68×67, F194^ECL2, S203^5.42×43, S204^5.43×44, S207^5.46×641, H296^6.58×58, and K305^7.32×31. Meanwhile, the antagonist ligands made interactions with W286^6.48×48 and Y316^7.43×42, both residues considered to be important in GPCR activation. The interpretation of ML analysis in human understandable form allowed us to construct an exquisitely detailed structure‐activity relationship that identifies small changes to the ligands that invert their pharmacological activity and thus helps to guide the drug discovery process. This approach can be readily applied to any drug target.

L-DOPA and Droxidopa: From Force Field Development to Molecular Docking into Human β2-Adrenergic Receptor

The increasing interest in the molecular mechanism of the binding of different agonists and antagonists to β₂-adrenergic receptor (β₂AR) inactive and active states has led us to investigate protein–ligand interactions using molecular docking calculations. To perform this study, the 3.2 Å X-ray crystal structure of the active conformation of human β₂AR in the complex with the endogenous agonist adrenaline has been used as a template for investigating the binding of two exogenous catecholamines to this adrenergic receptor. Here, we show the derivation of L-DOPA and Droxidopa OPLS all atom (AA) force field (FF) parameters via quantum mechanical (QM) calculations, molecular dynamics (MD) simulations in aqueous solutions of the two catecholamines and the molecular docking of both ligands into rigid and flexible β₂AR models. We observe that both ligands share with adrenaline similar experimentally observed binding anchor sites, which are constituted by Asp113/Asn312 and Ser203/Ser204/Ser207 side chains. Moreover, both L-DOPA and Droxidopa molecules exhibit binding affinities comparable to that predicted for adrenaline, which is in good agreement with previous experimental and computational results. L-DOPA and Droxidopa OPLS AA FFs have also been tested by performing MD simulations of these ligands docked into β₂AR proteins embedded in lipid membranes. Both hydrogen bonds and hydrophobic interaction networks observed over the 1 μs MD simulation are comparable with those derived from molecular docking calculations and MD simulations performed with the CHARMM FF.

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.

Molecular Basis of Inhibitory Mechanism of Naltrexone and Its Metabolites through Structural and Energetic Analyses

Naltrexone is a potent opioid antagonist with good blood–brain barrier permeability, targeting different endogenous opioid receptors, particularly the mu-opioid receptor (MOR). Therefore, it represents a promising candidate for drug development against drug addiction. However, the details of the molecular interactions of naltrexone and its derivatives with MOR are not fully understood, hindering ligand-based drug discovery. In the present study, taking advantage of the high-resolution X-ray crystal structure of the murine MOR (mMOR), we constructed a homology model of the human MOR (hMOR). A solvated phospholipid bilayer was built around the hMOR and submitted to microsecond (µs) molecular dynamics (MD) simulations to obtain an optimized hMOR model. Naltrexone and its derivatives were docked into the optimized hMOR model and submitted to µs MD simulations in an aqueous membrane system. The MD simulation results were submitted to the molecular mechanics–generalized Born surface area (MMGBSA) binding free energy calculations and principal component analysis. Our results revealed that naltrexone and its derivatives showed differences in protein–ligand interactions; however, they shared contacts with residues at TM2, TM3, H6, and TM7. The binding free energy and principal component analysis revealed the structural and energetic effects responsible for the higher potency of naltrexone compared to its derivatives.

Exploring the kinetic selectivity of drugs targeting the β1 -adrenoceptor

In this study, we report the β1 -adrenoceptor binding kinetics of several clinically relevant β1/2 -adrenoceptor (β1/2 AR) agonists and antagonists. [3 H]-DHA was used to label CHO-β1 AR for binding studies. The kinetics of ligand binding was assessed using a competition association binding method. Ligand physicochemical properties, including logD7.4 and the immobilized artificial membrane partition coefficient (KIAM ), were assessed using column-based methods. Protein Data Bank (PDB) structures and hydrophobic and electrostatic surface maps were constructed in PyMOL. We demonstrate that the hydrophobic properties of a molecule directly affect its kinetic association rate (kon ) and affinity for the β1 AR. In contrast to our findings at the β2 -adrenoceptor, KIAM , reflecting both hydrophobic and electrostatic interactions of the drug with the charged surface of biological membranes, was no better predictor than simple hydrophobicity measurements such as clogP or logD7.4 , at predicting association rate. Bisoprolol proved kinetically selective for the β1 AR subtype, dissociating 50 times slower and partly explaining its higher measured affinity for the β1 AR. We speculate that the association of positively charged ligands at the β1 AR is curtailed somewhat by its predominantly neutral/positive charged extracellular surface. Consequently, hydrophobic interactions in the ligand-binding pocket dominate the kinetics of ligand binding. In comparison at the β2 AR, a combination of hydrophobicity and negative charge attracts basic, positively charged ligands to the receptor's surface promoting the kinetics of ligand binding. Additionally, we reveal the potential role kinetics plays in the on-target and off-target pharmacology of clinically used β-blockers.

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.

Development of an Allostery Responsive Chromatographic Method for Screening Potential Allosteric Modulator of Beta2-adrenoceptor from a Natural Product-Derived DNA-Encoded Chemical Library

Allosteric ligands are promising drugs owing to their remote regulations of the orthosteric ligand signaling pathway. There are few allosteric ligands due to the lack of handy and efficacious method for the screening. Herein, we developed an affinity chromatographic method for allosteric ligand screening by immobilizing purified beta2 adrenoceptor (β₂-AR) onto macroporous silica gel by a two-point tethering method. The method relies on the occupation of the orthosteric site by an antagonist and the chelation of N-terminal His-tag of the receptor and Ni²⁺ coated on the gel. The immobilized β₂-AR demonstrated the greatest allosteric responsive feature when Cmpd-15 (0.25 μM) was included in the mobile phase. Under the same conditions, the association constants of three agonists (salbutamol, terbutaline, and tulobuterol) reduced to 47%, 19%, and 27% compared with the data without the inclusion of Cmpd-15 in the mobile phase. APF was screened as a potential allosteric modulator of β₂-AR by applying the immobilized receptor in a natural product-derived DNA-encoded chemical library (DEL). Relying on these results, we reasoned that the current method has potential in screening allosteric ligands of the receptor. We expect that it is applicable for the discovery of new allosteric binding sites of a target protein and screening allosteric modulators of the other receptors from complex samples.

Accelerating GPCR Drug Discovery With Conformation-Stabilizing VHHs

The human genome encodes 850 G protein-coupled receptors (GPCRs), half of which are considered potential drug targets. GPCRs transduce extracellular stimuli into a plethora of vital physiological processes. Consequently, GPCRs are an attractive drug target class. This is underlined by the fact that approximately 40% of marketed drugs modulate GPCRs. Intriguingly 60% of non-olfactory GPCRs have no drugs or candidates in clinical development, highlighting the continued potential of GPCRs as drug targets. The discovery of small molecules targeting these GPCRs by conventional high throughput screening (HTS) campaigns is challenging. Although the definition of success varies per company, the success rate of HTS for GPCRs is low compared to other target families (Fujioka and Omori, 2012; Dragovich et al., 2022). Beyond this, GPCR structure determination can be difficult, which often precludes the application of structure-based drug design approaches to arising HTS hits. GPCR structural studies entail the resource-demanding purification of native receptors, which can be challenging as they are inherently unstable when extracted from the lipid matrix. Moreover, GPCRs are flexible molecules that adopt distinct conformations, some of which need to be stabilized if they are to be structurally resolved. The complexity of targeting distinct therapeutically relevant GPCR conformations during the early discovery stages contributes to the high attrition rates for GPCR drug discovery programs. Multiple strategies have been explored in an attempt to stabilize GPCRs in distinct conformations to better understand their pharmacology. This review will focus on the use of camelid-derived immunoglobulin single variable domains (VHHs) that stabilize disease-relevant pharmacological states (termed ConfoBodies by the authors) of GPCRs, as well as GPCR:signal transducer complexes, to accelerate drug discovery. These VHHs are powerful tools for supporting in vitro screening, deconvolution of complex GPCR pharmacology, and structural biology purposes. In order to demonstrate the potential impact of ConfoBodies on translational research, examples are presented of their role in active state screening campaigns and structure-informed rational design to identify de novo chemical space and, subsequently, how such matter can be elaborated into more potent and selective drug candidates with intended pharmacology.

The Val34Met, Thr164Ile and Ser220Cys Polymorphisms of the β2-Adrenergic Receptor and Their Consequences on the Receptor Conformational Features: A Molecular Dynamics Simulation Study

The gene encoding the β2-adrenergic receptor (β2-AR) is polymorphic, which results in possible differences in a primary structure of this protein. It has been shown that certain types of polymorphisms are correlated with some clinical features of asthma, including airways reactivity, whereas the influence of other is not yet understood. Among polymorphisms affecting amino acids at positions 16, 27, 34, 164 and 220, the latter three are present in the crystal structure of β2-AR, which facilitates studying them by means of molecular dynamics simulations. The current study was focused on investigating to what extent the three polymorphisms of β2-AR (i.e., Val34Met, Thr164Ile and Ser220Cys) affect the interaction of β2-AR with its natural molecular environment which includes: lipid bilayer (in the case of all three polymorphs) and Gs protein (which participates in β2-AR-mediated signaling; in the case of Ser220Cys). We have designed and carried out a series of molecular dynamics simulations at different level of resolution (i.e., either coarse-grained or atomistic simulations), accompanied by thermodynamic integration protocol, in order to identify potential polymorphism-induced alterations in structural, conformational or energetic features of β2-AR. The results indicate the lack of significant differences in the case of energies involved in the β2-AR-lipid bilayer interactions. Some differences have been observed when considering the polymorphism-induced alterations in β2-AR-Gs protein binding, but their magnitude is also negligible in relation to the absolute free energy difference correlated with the β2-AR-Gs affinity. The Val34Met and Thr164Ile polymorphisms are weakly correlated with alteration of the conformational features of the receptor around polymorphic sites. On the contrary, it has been concluded that the Ser220Cys polymorphism is correlated with several structural alterations located in the intracellular region of β2-AR, which can induce G-protein binding and, subsequently, the polymorphism-correlated therapeutic responses. More precisely, these alterations involve vicinity of intracellular loops and, in part, are the direct consequence of disturbed interactions of Ser/Cys220 sidechain within 5th transmembrane domain. Structurally, the dynamic structure exhibited by the β2-AR^Ser220 polymorph is closer to the Gs-compatible structure of β2-AR.

GPCR-mediated β-arrestin activation deconvoluted with single-molecule precision

β-arrestins bind G protein-coupled receptors to terminate G protein signaling and to facilitate other downstream signaling pathways. Using single-molecule fluorescence resonance energy transfer imaging, we show that β-arrestin is strongly autoinhibited in its basal state. Its engagement with a phosphopeptide mimicking phosphorylated receptor tail efficiently releases the β-arrestin tail from its N domain to assume distinct conformations. Unexpectedly, we find that β-arrestin binding to phosphorylated receptor, with a phosphorylation barcode identical to the isolated phosphopeptide, is highly inefficient and that agonist-promoted receptor activation is required for β-arrestin activation, consistent with the release of a sequestered receptor C tail. These findings, together with focused cellular investigations, reveal that agonism and receptor C-tail release are specific determinants of the rate and efficiency of β-arrestin activation by phosphorylated receptor. We infer that receptor phosphorylation patterns, in combination with receptor agonism, synergistically establish the strength and specificity with which diverse, downstream β-arrestin-mediated events are directed.

The organic cation Transporter 2 regulates dopamine D1 receptor signaling at the Golgi apparatus

Dopamine is a key catecholamine in the brain and kidney, where it is involved in a number of physiological functions such as locomotion, cognition, emotion, endocrine regulation, and renal function. As a membrane-impermeant hormone and neurotransmitter, dopamine is thought to signal by binding and activating dopamine receptors, members of the G protein coupled receptor (GPCR) family, only on the plasma membrane. Here, using novel nanobody-based biosensors, we demonstrate for the first time that the dopamine D1 receptor (D1DR), the primary mediator of dopaminergic signaling in the brain and kidney, not only functions on the plasma membrane but becomes activated at the Golgi apparatus in the presence of its ligand. We present evidence that activation of the Golgi pool of D1DR is dependent on organic cation transporter 2 (OCT2), a dopamine transporter, providing an explanation for how the membrane-impermeant dopamine accesses subcellular pools of D1DR. We further demonstrate that dopamine activates Golgi-D1DR in murine striatal medium spiny neurons, and this activity depends on OCT2 function. We also introduce a new approach to selectively interrogate compartmentalized D1DR signaling by inhibiting Gαs coupling using a nanobody-based chemical recruitment system. Using this strategy, we show that Golgi-localized D1DRs regulate cAMP production and mediate local protein kinase A activation. Together, our data suggest that spatially compartmentalized signaling hubs are previously unappreciated regulatory aspects of D1DR signaling. Our data provide further evidence for the role of transporters in regulating subcellular GPCR activity.

Cryo-EM structures of the β3 adrenergic receptor bound to solabegron and isoproterenol

The β₃-adrenergic receptor (β₃AR) is the most essential drug target for overactive bladder and has therapeutic potentials for the treatments of type 2 diabetes and obesity. Here, we report the cryo-electron microscopy structures of the β₃AR-G_s signaling complexes with the selective agonist, solabegron and the nonselective agonist, isoproterenol. Comparison of the isoproterenol-, mirabegron-, and solabegron-bound β₃AR structures revealed that the extracellular loop 2 changes its conformation depending on the bound agonist and plays an essential role in solabegron binding. Moreover, β₃AR has an intrinsically narrow exosite, regardless of the agonist type. This structural feature clearly explains why β₃AR prefers mirabegron and solabegron, as the narrow exosite is suitable for binding with agonists with elongated shapes. Our study deepens the understanding of the binding characteristics of β₃AR agonists and may pave the way for developing β₃AR-selective drugs.

Synthesis of 2,2-difluoro-2-arylethylamines as fluorinated analogs of octopamine and noradrenaline

A series of 2,2-difluoro-2-arylethylamines was synthesized as fluorinated analogs of octopamine and noradrenaline with the expectation of bioisosteric OH/F exchanges. The syntheses of these compounds were performed by a Suzuki–Miyaura cross-coupling reaction of 4-(bromodifluoroacetyl)morpholine with aryl boronic acids to produce the intermediate 2,2-difluoro-2-arylacetamides, followed by transformation of difluoroacetamide to difluoroethylamine.

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

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

Structural insights into ligand recognition, activation, and signaling of the α2A adrenergic receptor

The α_2A adrenergic receptor (α_2AAR) is a G protein (heterotrimeric guanine nucleotide-binding protein)-coupled receptor that mediates important physiological functions in response to the endogenous neurotransmitters norepinephrine and epinephrine, as well as numerous chemically distinct drugs. However, the molecular mechanisms of drug actions remain poorly understood. Here, we report the cryo-electron microscopy structures of the human α_2AAR-GoA complex bound to norepinephrine and three imidazoline derivatives (brimonidine, dexmedetomidine, and oxymetazoline). Together with mutagenesis and functional data, these structures provide important insights into the molecular basis of ligand recognition, activation, and signaling at the α_2AAR. Further structural analyses uncover different molecular determinants between α_2AAR and βARs for recognition of norepinephrine and key regions that determine the G protein coupling selectivity. Overall, our studies provide a framework for understanding the signal transduction of the adrenergic system at the atomic level, which will facilitate rational structure-based discovery of safer and more effective medications for α_2AAR.

New Horizons in Structural Biology of Membrane Proteins: Experimental Evaluation of the Role of Conformational Dynamics and Intrinsic Flexibility

A plethora of membrane proteins are found along the cell surface and on the convoluted labyrinth of membranes surrounding organelles. Since the advent of various structural biology techniques, a sub-population of these proteins has become accessible to investigation at near-atomic resolutions. The predominant bona fide methods for structure solution, X-ray crystallography and cryo-EM, provide high resolution in three-dimensional space at the cost of neglecting protein motions through time. Though structures provide various rigid snapshots, only an amorphous mechanistic understanding can be inferred from interpolations between these different static states. In this review, we discuss various techniques that have been utilized in observing dynamic conformational intermediaries that remain elusive from rigid structures. More specifically we discuss the application of structural techniques such as NMR, cryo-EM and X-ray crystallography in studying protein dynamics along with complementation by conformational trapping by specific binders such as antibodies. We finally showcase the strength of various biophysical techniques including FRET, EPR and computational approaches using a multitude of succinct examples from GPCRs, transporters and ion channels.

Glyco-steroidal amphiphiles (GSAs) for membrane protein structural study

Integral membrane proteins pose considerable challenges to high resolution structural analysis. Maintaining membrane proteins in their native state during protein isolation is essential for structural study of these bio-macromolecules. Detergents are the most commonly used amphiphilic compounds for stabilizing membrane proteins in solution outside a lipid bilayer. We previously introduced a glyco-diosgenin (GDN) detergent that was shown to be highly effective at stabilizing a wide range of membrane proteins. This steroidal detergent has additionally gained attention due to its compatibility with membrane protein structure study via cryo-EM. However, synthetic inconvenience limits widespread use of GDN in membrane protein study. To improve its synthetic accessibility and to further enhance detergent efficacy for protein stabilization, we designed a new class of glyco-steroid-based detergents using three steroid units: cholestanol, cholesterol and diosgenin. These new detergents were efficiently prepared and showed marked efficacy for protein stabilization in evaluation with a few model membrane proteins including two G protein-coupled receptors. Some new agents were not only superior to a gold standard detergent, DDM, but were also more effective than the original GDN at preserving protein integrity long term. These agents represent valuable alternatives to GDN, and are likely to facilitate structural determination of challenging membrane proteins.

Magic angle spinning NMR of G protein-coupled receptors

G protein-coupled receptors (GPCRs) have a simple seven transmembrane helix architecture which has evolved to recognize a diverse number of chemical signals. The more than 800 GPCRs encoded in the human genome function as receptors for vision, smell and taste, and mediate key physiological processes. Consequently, these receptors are a major target for pharmaceuticals. Protein crystallography and electron cryo-microscopy have provided high resolution structures of many GPCRs in both active and inactive conformations. However, these structures have not sparked a surge in rational drug design, in part because GPCRs are inherently dynamic and the structural changes induced by ligand or drug binding to stabilize inactive or active conformations are often subtle rearrangements in packing or hydrogen-bonding interactions. NMR spectroscopy provides a sensitive probe of local structure and dynamics at specific sites within these receptors as well as global changes in receptor structure and dynamics. These methods can also capture intermediate states and conformations with low populations that provide insights into the activation pathways. We review the use of solid-state magic angle spinning NMR to address the structure and activation mechanisms of GPCRs. The focus is on the large and diverse class A family of receptors. We highlight three specific class A GPCRs in order to illustrate how solid-state, as well as solution-state, NMR spectroscopy can answer questions in the field involving how different GPCR classes and subfamilies are activated by their associated ligands, and how small molecule drugs can modulate GPCR activation.

Allosteric modulation of dopamine D2L receptor in complex with Gi1 and Gi2 proteins: the effect of subtle structural and stereochemical ligand modifications

Background

Allosteric modulation of G protein-coupled receptors (GPCRs) is nowadays one of the hot topics in drug discovery. In particular, allosteric modulators of D₂ receptor have been proposed as potential modern therapeutics to treat schizophrenia and Parkinson’s disease.

Methods

To address some subtle structural and stereochemical aspects of allosteric modulation of D₂ receptor, we performed extensive in silico studies of both enantiomers of two compounds (compound 1 and compound 2), and one of them (compound 2) was synthesized as a racemate in-house and studied in vitro.

Results

Our molecular dynamics simulations confirmed literature reports that the R enantiomer of compound 1 is a positive allosteric modulator of the D_2L receptor, while its S enantiomer is a negative allosteric modulator. Moreover, based on the principal component analysis (PCA), we hypothesized that both enantiomers of compound 2 behave as silent allosteric modulators, in line with our in vitro studies. PCA calculations suggest that the most pronounced modulator-induced receptor rearrangements occur at the transmembrane helix 7 (TM7). In particular, TM7 bending at the conserved P7.50 and G7.42 was observed. The latter resides next to the Y7.43, which is a significant part of the orthosteric binding site. Moreover, the W7.40 conformation seems to be affected by the presence of the positive allosteric modulator.

Conclusions

Our work reveals that allosteric modulation of the D_2L receptor can be affected by subtle ligand modifications. A change in configuration of a chiral carbon and/or minor structural modulator modifications are solely responsible for the functional outcome of the allosteric modulator.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s43440-021-00352-x.

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.

Molecular mechanism of agonism and inverse agonism in ghrelin receptor

Much effort has been invested in the investigation of the structural basis of G protein-coupled receptors (GPCRs) activation. Inverse agonists, which can inhibit GPCRs with constitutive activity, are considered useful therapeutic agents, but the molecular mechanism of such ligands remains insufficiently understood. Here, we report a crystal structure of the ghrelin receptor bound to the inverse agonist PF-05190457 and a cryo-electron microscopy structure of the active ghrelin receptor-Go complex bound to the endogenous agonist ghrelin. Our structures reveal a distinct binding mode of the inverse agonist PF-05190457 in the ghrelin receptor, different from the binding mode of agonists and neutral antagonists. Combining the structural comparisons and cellular function assays, we find that a polar network and a notable hydrophobic cluster are required for receptor activation and constitutive activity. Together, our study provides insights into the detailed mechanism of ghrelin receptor binding to agonists and inverse agonists, and paves the way to design specific ligands targeting ghrelin receptors.

Ghrelin receptor regulates energy homeostasis through constitutive activity or by the ghrelin. Here the authors report two structures of ghrelin receptor bound to agonist and inverse agonist, providing insights into the mechanism of inverse agonism, which is of interest for specific ligand design.

Allosteric Effect of Nanobody Binding on Ligand-Specific Active States of the β2 Adrenergic Receptor

Nanobody binding stabilizes G-protein-coupled receptors (GPCR) in a fully active state and modulates their affinity for bound ligands. However, the atomic-level basis for this allosteric regulation remains elusive. Here, we investigate the conformational changes induced by the binding of a nanobody (Nb80) on the active-like β2 adrenergic receptor (β2AR) via enhanced sampling molecular dynamics simulations. Dimensionality reduction analysis shows that Nb80 stabilizes structural features of the β2AR with an ∼14 Å outward movement of transmembrane helix 6 and a close proximity of transmembrane (TM) helices 5 and 7, and favors the fully active-like conformation of the receptor, independent of ligand binding, in contrast to the conditions under which no intracellular binding partner is bound, in which case the receptor is only stabilized in an intermediate-active state. This activation is supported by the residues located at hotspots located on TMs 5, 6, and 7, as shown by supervised machine learning methods. Besides, ligand-specific subtle differences in the conformations assumed by intracellular loop 2 and extracellular loop 2 are captured from the trajectories of various ligand-bound receptors in the presence of Nb80. Dynamic network analysis further reveals that Nb80 binding triggers tighter and stronger local communication networks between the Nb80 and the ligand-binding sites, primarily involving residues around ICL2 and the intracellular end of TM3, TM5, TM6, as well as ECL2, ECL3, and the extracellular ends of TM6 and TM7. In particular, we identify unique allosteric signal transmission mechanisms between the Nb80-binding site and the extracellular domains in conformations modulated by a full agonist, BI167107, and a G-protein-biased partial agonist, salmeterol, involving mainly TM1 and TM2, and TM5, respectively. Altogether, our results provide insights into the effect of intracellular binding partners on the GPCR activation mechanism, which should be taken into account in structure-based drug discovery.

Recent Advances in Structure, Function, and Pharmacology of Class A Lipid GPCRs: Opportunities and Challenges for Drug Discovery

Great progress has been made over the past decade in understanding the structural, functional, and pharmacological diversity of lipid GPCRs. From the first determination of the crystal structure of bovine rhodopsin in 2000, much progress has been made in the field of GPCR structural biology. The extraordinary progress in structural biology and pharmacology of GPCRs, coupled with rapid advances in computational approaches to study receptor dynamics and receptor-ligand interactions, has broadened our comprehension of the structural and functional facets of the receptor family members and has helped usher in a modern age of structure-based drug design and development. First, we provide a primer on lipid mediators and lipid GPCRs and their role in physiology and diseases as well as their value as drug targets. Second, we summarize the current advancements in the understanding of structural features of lipid GPCRs, such as the structural variation of their extracellular domains, diversity of their orthosteric and allosteric ligand binding sites, and molecular mechanisms of ligand binding. Third, we close by collating the emerging paradigms and opportunities in targeting lipid GPCRs, including a brief discussion on current strategies, challenges, and the future outlook.

Analysis of L-DOPA and droxidopa binding to human β2-adrenergic receptor

Over the last two decades, an increasing number of studies has been devoted to a deeper understanding of the molecular process involved in the binding of various agonists and antagonists to active and inactive conformations of β2-adrenergic receptor (β2AR). The 3.2 Å x-ray crystal structure of human β2AR active state in combination with the endogenous low affinity agonist adrenaline offers an ideal starting structure for studying the binding of various catecholamines to adrenergic receptors. We show that molecular docking of levodopa (L-DOPA) and droxidopa into rigid and flexible β2AR models leads for both ligands to binding anchor sites comparable to those experimentally reported for adrenaline, namely D113/N312 and S203/S204/S207 side chains. Both ligands have a hydrogen bond network that is extremely similar to those of noradrenaline and dopamine. Interestingly, redocking neutral and protonated versions of adrenaline to rigid and flexible β2AR models results in binding poses that are more energetically stable and distinct from the x-ray crystal structure. Similarly, lowest energy conformations of noradrenaline and dopamine generated by docking into flexible β2AR models had binding free energies lower than those of best poses in rigid receptor models. Furthermore, our findings show that L-DOPA and droxidopa molecules have binding affinities comparable to those predicted for adrenaline, noradrenaline, and dopamine, which are consistent with previous experimental and computational findings and supported by the molecular dynamics simulations of β2AR-ligand complexes performed here.

Identification and characterization of an atypical Gαs-biased β2AR agonist that fails to evoke airway smooth muscle cell tachyphylaxis

We sought β₂AR agonists for treating obstructive lung diseases such as asthma, in which this receptor relaxes airway smooth muscle (ASM) cells and opens airways. Agonists favoring Gs coupling (leads to airway relaxation) compared with activating β-arrestin (limits effectiveness due to receptor desensitization) were pursued in a 40-million-compound screening library. Of several agonists identified, one was apparently biased away from β-arrestin. Agonist–receptor–G protein modeling revealed different receptor interactions compared with other agonists. The favorable effects of the apparent biasing with this agonist were demonstrated in a physiologic system (ASM relaxation). These studies point to a different structural class of β-agonists that might be used to treat obstructive lung diseases without the adverse effects associated with tachyphylaxis.

G protein–coupled receptors display multifunctional signaling, offering the potential for agonist structures to promote conformational selectivity for biased outputs. For β₂-adrenergic receptors (β₂AR), unbiased agonists stabilize conformation(s) that evoke coupling to Gαs (cyclic adenosine monophosphate [cAMP] production/human airway smooth muscle [HASM] cell relaxation) and β-arrestin engagement, the latter acting to quench Gαs signaling, contributing to receptor desensitization/tachyphylaxis. We screened a 40-million-compound scaffold ranking library, revealing unanticipated agonists with dihydroimidazolyl-butyl-cyclic urea scaffolds. The S-stereoisomer of compound C1 shows no detectable β-arrestin engagement/signaling by four methods. However, C1-S retained Gαs signaling—a divergence of the outputs favorable for treating asthma. Functional studies with two models confirmed the biasing: β₂AR-mediated cAMP signaling underwent desensitization to the unbiased agonist albuterol but not to C1-S, and desensitization of HASM cell relaxation was observed with albuterol but not with C1-S. These HASM results indicate biologically pertinent biasing of C1-S, in the context of the relevant physiologic response, in the human cell type of interest. Thus, C1-S was apparently strongly biased away from β-arrestin, in contrast to albuterol and C5-S. C1-S structural modeling and simulations revealed binding differences compared with unbiased epinephrine at transmembrane (TM) segments 3,5,6,7 and ECL2. C1-S (R2 = cyclohexane) was repositioned in the pocket such that it lost a TM6 interaction and gained a TM7 interaction compared with the analogous unbiased C5-S (R2 = benzene group), which appears to contribute to C1-S biasing away from β-arrestin. Thus, an agnostic large chemical-space library identified agonists with receptor interactions that resulted in relevant signal splitting of β₂AR actions favorable for treating obstructive lung disease.

A novel red fluorescence dopamine biosensor selectively detects dopamine in the presence of norepinephrine in vitro

Dopamine (DA) and norepinephrine (NE) are pivotal neuromodulators that regulate a broad range of brain functions, often in concert. Despite their physiological importance, untangling the relationship between DA and NE in the fine control of output function is currently challenging, primarily due to a lack of techniques to allow the observation of spatiotemporal dynamics with sufficiently high selectivity. Although genetically encoded fluorescent biosensors have been developed to detect DA, their poor selectivity prevents distinguishing DA from NE. Here, we report the development of a red fluorescent genetically encoded GPCR (G protein-coupled receptor)-activation reporter for DA termed 'R-GenGAR-DA'. More specifically, a circular permutated red fluorescent protein (cpmApple) was replaced by the third intracellular loop of human DA receptor D1 (DRD1) followed by the screening of mutants within the linkers between DRD1 and cpmApple. We developed two variants: R-GenGAR-DA1.1, which brightened following DA stimulation, and R-GenGAR-DA1.2, which dimmed. R-GenGAR-DA1.2 demonstrated a reasonable dynamic range (ΔF/F0 = - 43%), DA affinity (EC50 = 0.92 µM) and high selectivity for DA over NE (66-fold) in HeLa cells. Taking advantage of the high selectivity of R-GenGAR-DA1.2, we monitored DA in presence of NE using dual-color fluorescence live imaging, combined with the green-NE biosensor GRABNE1m, which has high selectivity for NE over DA (> 350-fold) in HeLa cells and hippocampal neurons grown from primary culture. Thus, this is a first step toward the multiplex imaging of these neurotransmitters in, for example, freely moving animals, which will provide new opportunities to advance our understanding of the high spatiotemporal dynamics of DA and NE in normal and abnormal brain function.

Molecular Mechanisms of Diverse Activation Stimulated by Different Biased Agonists for the β2-Adrenergic Receptor

β2AR is an important drug target protein involving many diseases. Biased drugs induce specific signaling and provide additional clinical utility to optimize β2AR-based therapies. However, the biased signaling mechanism has not been elucidated. Motivated by the issue, we chose four agonists with divergent bias (balanced agonist, G-protein-biased agonist, and β-arrestin-biased agonists) and utilized Gaussian accelerated molecular dynamics simulation coupled with a dynamic network to probe the molecular mechanisms of distinct biased activation induced by the structural differences between the four agonists. Our simulations reveal that the G-protein-biased agonist induces an open conformation with the outward shifts of TM6 and TM7 for the intracellular domain, which will be beneficial to couple G protein. In contrast, the β-arrestin-biased agonists regulate an occluded conformation with a slightly outward movement of TM6 and an inward shift of TM7, which should favor β-arrestin signaling. The balanced agonist does not induce an observable outward shift for TM6 but, along with a slight tilt for TM7, leads to an inactive-like conformation. In addition, our results reveal the first time that ICL3 presents specific conformations with different agonists. The G-protein-biased agonist drives ICL3 to open so that the G protein-binding pocket can be available, while the β-arrestin-biased agonists induce ICL3 to form a closed conformation with a stable local α-helix. MM/PBSA analysis further reveals that the hydroxyl groups in the resorcinol of the G-protein-biased agonist form strong interactions with Y5.38 and S5.42, thus preventing tilting of the TM5 extracellular end. The catechol of the balanced agonist and the β-arrestin-biased ones induces the rearrangement of two hydrophobic residues F6.52 and W6.48. However, different from the balanced agonist, the ethyl substituent of β-arrestin-biased agonists forms additional hydrophobic interactions with W6.48 and F6.51 after the rearrangement, which should contribute to the β-arrestin bias. The shortest pathway analysis further reveals that the three residues Y7.43, N7.45, and N7.49 are crucial for allosterically regulating G-protein-biased signaling, while the two residues W6.48 and F6.44 make an important contribution to regulate β-arrestin-biased signaling. For the balanced agonist NE, the allosteric regulation pathway simultaneously involves the residue associated with G-protein-biased signaling like S5.46 and the residues related to β-arrestin-biased signaling like W6.48 and F6.44, thus producing unbiased signaling. The observations could advance our understanding of the biased activation mechanism on class A GPCRs and provide a useful guideline for the design of biased drugs.

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

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

Optical tools to study the subcellular organization of GPCR neuromodulation

Modulation of neuronal circuit activity is key to information processing in the brain. G protein-coupled receptors (GPCRs), the targets of most neuromodulatory ligands, show extremely diverse expression patterns in neurons and receptors can be localized in various sub-neuronal membrane compartments. Upon activation, GPCRs promote signaling cascades that alter the level of second messengers, drive phosphorylation changes, modulate ion channel function, and influence gene expression, all of which critically impact neuron physiology. Because of its high degree of complexity, this form of interneuronal communication has remained challenging to integrate into our conceptual understanding of brain function. Recent technological advances in fluorescence microscopy and the development of optical biosensors now allow investigating neuromodulation with unprecedented resolution on the level of individual cells. In this review, we will highlight recent imaging techniques that enable determining the precise localization of GPCRs in neurons, with specific focus on the subcellular and nanoscale level. Downstream of receptors, we describe novel conformation-specific biosensors that allow for real-time monitoring of GPCR activation and of distinct signal transduction events in neurons. Applying these new tools has the potential to provide critical insights into the function and organization of GPCRs in neuronal cells and may help decipher the molecular and cellular mechanisms that underlie neuromodulation.

Encoding mu-opioid receptor biased agonism with interaction fingerprints

Opioids are potent painkillers, however, their therapeutic use requires close medical monitoring to diminish the risk of severe adverse effects. The G-protein biased agonists of the μ-opioid receptor (MOR) have shown safer therapeutic profiles than non-biased ligands. In this work, we performed extensive all-atom molecular dynamics simulations of two markedly biased ligands and a balanced reference molecule. From those simulations, we identified a protein-ligand interaction fingerprint that characterizes biased ligands. Then, we built and virtually screened a database containing 68,740 ligands with proven or potential GPCR agonistic activity. Exemplary molecules that fulfill the interacting pattern for biased agonism are showcased, illustrating the usefulness of this work for the search of biased MOR ligands and how this contributes to the understanding of MOR biased signaling.

SYNBIP: synthetic binding proteins for research, diagnosis and therapy

The success of protein engineering and design has extensively expanded the protein space, which presents a promising strategy for creating next-generation proteins of diverse functions. Among these proteins, the synthetic binding proteins (SBPs) are smaller, more stable, less immunogenic, and better of tissue penetration than others, which make the SBP-related data attracting extensive interest from worldwide scientists. However, no database has been developed to systematically provide the valuable information of SBPs yet. In this study, a database named ‘Synthetic Binding Proteins for Research, Diagnosis, and Therapy (SYNBIP)’ was thus introduced. This database is unique in (a) comprehensively describing thousands of SBPs from the perspectives of scaffolds, biophysical & functional properties, etc.; (b) panoramically illustrating the binding targets & the broad application of each SBP and (c) enabling a similarity search against the sequences of all SBPs and their binding targets. Since SBP is a human-made protein that has not been found in nature, the discovery of novel SBPs relied heavily on experimental protein engineering and could be greatly facilitated by in-silico studies (such as AI and computational modeling). Thus, the data provided in SYNBIP could lay a solid foundation for the future development of novel SBPs. The SYNBIP is accessible without login requirement at both official (https://idrblab.org/synbip/) and mirror (http://synbip.idrblab.net/) sites.

Unique Positive Cooperativity Between the β-Arrestin-Biased β-Blocker Carvedilol and a Small Molecule Positive Allosteric Modulator of the β2-Adrenergic Receptor

Among β-blockers that are clinically prescribed for heart failure, carvedilol is a first-choice agent with unique pharmacological properties. Carvedilol is distinct from other β-blockers in its ability to elicit β-arrestin–biased agonism, which has been suggested to underlie its cardioprotective effects. Augmenting the pharmacologic properties of carvedilol thus holds the promise of developing more efficacious and/or biased β-blockers. We recently identified compound-6 (cmpd-6), the first small molecule positive allosteric modulator of the β2-adrenergic receptor (β2AR). Cmpd-6 is positively cooperative with orthosteric agonists at the β2AR and enhances agonist-mediated transducer (G-protein and β-arrestin) signaling in an unbiased manner. Here, we report that cmpd-6, quite unexpectedly, displays strong positive cooperativity only with carvedilol among a panel of structurally diverse β-blockers. Cmpd-6 enhances the binding affinity of carvedilol for the β2AR and augments its ability to competitively antagonize agonist-induced cAMP generation. Cmpd-6 potentiates β-arrestin1– but not Gs-protein–mediated high-affinity binding of carvedilol at the β2AR and β-arrestin–mediated cellular functions in response to carvedilol including extracellular signal-regulated kinase phosphorylation, receptor endocytosis, and trafficking into lysosomes. Importantly, an analog of cmpd-6 that selectively retains positive cooperativity with carvedilol acts as a negative modulator of agonist-stimulated β2AR signaling. These unprecedented cooperative properties of carvedilol and cmpd-6 have implications for fundamental understanding of G-protein–coupled receptor (GPCR) allosteric modulation, as well as for the development of more effective biased beta blockers and other GPCR therapeutics.

Universal Properties and Specificities of the β2-Adrenergic Receptor-Gs Protein Complex Activation Mechanism Revealed by All-Atom Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are transmembrane proteins of high pharmacological relevance. It has been proposed that their activity is linked to structurally distinct, dynamically interconverting functional states and the process of activation relies on an interconnecting network of conformational switches in the transmembrane domain. However, it is yet to be uncovered how ligands with different extents of functional effect exert their actions. According to our recent hypothesis, based on indirect observations and the literature data, the transmission of the external stimulus to the intracellular surface is accompanied by the shift of macroscopic polarization in the transmembrane domain, furnished by concerted movements of highly conserved polar motifs and the rearrangement of polar species. In this follow-up study, we have examined the β2-adrenergic receptor (β2AR) to see if our hypothesis drawn from an extensive study of the μ-opioid receptor (MOP) is fundamental and directly transferable to other class A GPCRs. We have found that there are some general similarities between the two receptors, in agreement with previous studies, and there are some receptor-specific differences that could be associated with different signaling pathways.

A resource of high-quality and versatile nanobodies for drug delivery

Therapeutic and diagnostic efficacies of small biomolecules and chemical compounds are hampered by suboptimal pharmacokinetics. Here, we developed a repertoire of robust and high-affinity antihuman serum albumin nanobodies (Nb_HSA) that can be readily fused to small biologics for half-life extension. We characterized the thermostability, binding kinetics, and cross-species reactivity of Nb_HSAs, mapped their epitopes, and structurally resolved a tetrameric HSA-Nb complex. We parallelly determined the half-lives of a cohort of selected Nb_HSAs in an HSA mouse model by quantitative proteomics. Compared to short-lived control nanobodies, the half-lives of Nb_HSAs were drastically prolonged by 771-fold. Nb_HSAs have distinct and diverse pharmacokinetics, positively correlating with their albumin binding affinities at the endosomal pH. We then generated stable and highly bioactive Nb_HSA-cytokine fusion constructs “Duraleukin” and demonstrated Duraleukin's high preclinical efficacy for cancer treatment in a melanoma model. This high-quality and versatile Nb toolkit will help tailor drug half-life to specific medical needs.

•

We provide a resource of high-affinity and versatile albumin nanobodies for drug delivery

•

We systematically map albumin nanobody epitopes by hybrid structural approaches

•

We parallelly measure the pharmacokinetics of nanobodies in a humanized mouse model

•

We develop nanobody-cytokine conjugates “Duraleukin” for cancer immunotherapy

2D SIFt: a matrix of ligand-receptor interactions

Depicting a ligand-receptor complex via Interaction Fingerprints has been shown to be both a viable data visualization and an analysis tool. The spectrum of its applications ranges from simple visualization of the binding site through analysis of molecular dynamics runs, to the evaluation of the homology models and virtual screening. Here we present a novel tool derived from the Structural Interaction Fingerprints providing a detailed and unique insight into the interactions between receptor and specific regions of the ligand (grouped into pharmacophore features) in the form of a matrix, a 2D-SIFt descriptor. The provided implementation is easy to use and extends the python library, allowing the generation of interaction matrices and their manipulation (reading and writing as well as producing the average 2D-SIFt). The library for handling the interaction matrices is available via repository http://bitbucket.org/zchl/sift2d.

Nanobodies as sensors of GPCR activation and signaling

Nanobodies have emerged as useful tools to study G protein-coupled receptor (GPCR) structure, dynamic, and subcellular localization. Initially, several nanobodies have been developed as chaperones to facilitate GPCR crystallization. To explore their potential as biosensors to monitor receptor activation and dynamics, we here described protocols to characterize nanobody's interaction with GPCRs and their application as probes for protein identification and visualization on the cellular level. We also introduced a chimeric approach to enable a kappa-opioid receptor derived nanobody to bind to other GPCRs, including orphan GPCRs whose endogenous ligand or intracellular transducers are unknown. This approach provides a reporter assay to identify tool molecules to study the function of orphan GPCRs.

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.