Extra-helical binding site of a glucagon receptor antagonist

Identifying and Assessing Putative Allosteric Sites and Modulators for CXCR4 Predicted through Network Modeling and Site Identification by Ligand Competitive Saturation

The chemokine receptor CXCR4 is a critical target for the treatment of several cancer types and HIV-1 infections. While orthosteric and allosteric modulators have been developed targeting its extracellular or transmembrane regions, the intramembrane region of CXCR4 may also include allosteric binding sites suitable for the development of allosteric drugs. To investigate this, we apply the Gaussian Network Model (GNM) to the monomeric and dimeric forms of CXCR4 to identify residues essential for its local and global motions located in the hinge regions of the protein. Residue interaction network (RIN) analysis suggests hub residues that participate in allosteric communication throughout the receptor. Mutual residues from the network models reside in regions with a high capacity to alter receptor dynamics upon ligand binding. We then investigate the druggability of these potential allosteric regions using the site identification by ligand competitive saturation (SILCS) approach, revealing two putative allosteric sites on the monomer and three on the homodimer. Two screening campaigns with Glide and SILCS-Monte Carlo docking using FDA-approved drugs suggest 20 putative hit compounds including antifungal drugs, anticancer agents, HIV protease inhibitors, and antimalarial drugs. In vitro assays considering mAB 12G5 and CXCL12 demonstrate both positive and negative allosteric activities of these compounds, supporting our computational approach. However, in vivo functional assays based on the recruitment of β-arrestin to CXCR4 do not show significant agonism and antagonism at a single compound concentration. The present computational pipeline brings a new perspective to computer-aided drug design by combining conformational dynamics based on network analysis and cosolvent analysis based on the SILCS technology to identify putative allosteric binding sites using CXCR4 as a showcase.

G protein-coupled receptors (GPCRs): advances in structures, mechanisms, and drug discovery

G protein-coupled receptors (GPCRs), the largest family of human membrane proteins and an important class of drug targets, play a role in maintaining numerous physiological processes. Agonist or antagonist, orthosteric effects or allosteric effects, and biased signaling or balanced signaling, characterize the complexity of GPCR dynamic features. In this study, we first review the structural advancements, activation mechanisms, and functional diversity of GPCRs. We then focus on GPCR drug discovery by revealing the detailed drug-target interactions and the underlying mechanisms of orthosteric drugs approved by the US Food and Drug Administration in the past five years. Particularly, an up-to-date analysis is performed on available GPCR structures complexed with synthetic small-molecule allosteric modulators to elucidate key receptor-ligand interactions and allosteric mechanisms. Finally, we highlight how the widespread GPCR-druggable allosteric sites can guide structure- or mechanism-based drug design and propose prospects of designing bitopic ligands for the future therapeutic potential of targeting this receptor family.

Unraveling the Interplay of Extracellular Domain Conformational Changes and Parathyroid Hormone Type 1 Receptor Activation in Class B1 G Protein-Coupled Receptors: Integrating Enhanced Sampling Molecular Dynamics Simulations and Markov State Models

Parathyroid hormone (PTH) type 1 receptor (PTH1R), as a typical class B1 G protein-coupled receptor (GPCR), is responsible for regulating bone turnover and maintaining calcium homeostasis, and its dysregulation has been implicated in the development of several diseases. The extracellular domain (ECD) of PTH1R is crucial for the recognition and binding of ligands, and the receptor may exhibit an autoinhibited state with the closure of the ECD in the absence of ligands. However, the correlation between ECD conformations and PTH1R activation remains unclear. Thus, this study combines enhanced sampling molecular dynamics (MD) simulations and Markov state models (MSMs) to reveal the possible relevance between the ECD conformations and the activation of PTH1R. First, 22 intermediate structures are generated from the autoinhibited state to the active state and conducted for 10 independent 200 ns simulations each. Then, the MSM is constructed based on the cumulative 44 μs simulations with six identified microstates. Finally, the potential interplay between ECD conformational changes and PTH1R activation as well as cryptic allosteric pockets in the intermediate states during receptor activation is revealed. Overall, our findings reveal that the activation of PTH1R has a specific correlation with ECD conformational changes and provide essential insights for GPCR biology and developing novel allosteric modulators targeting cryptic sites.

Molecular determinants of antagonist interactions with chemokine receptors CCR2 and CCR5

By driving monocyte chemotaxis, the chemokine receptor CCR2 shapes inflammatory responses and the formation of tumor microenvironments. This makes it a promising target in inflammation and immuno-oncology; however, despite extensive efforts, there are no FDA-approved CCR2-targeting therapeutics. Cited challenges include the redundancy of the chemokine system, suboptimal properties of compound candidates, and species differences that confound the translation of results from animals to humans. Structure-based drug design can rationalize and accelerate the discovery and optimization of CCR2 antagonists to address these challenges. The prerequisites for such efforts include an atomic-level understanding of the molecular determinants of action of existing antagonists. In this study, using molecular docking and artificial-intelligence-powered compound library screening, we uncover the structural principles of small molecule antagonism and selectivity towards CCR2 and its sister receptor CCR5. CCR2 orthosteric inhibitors are shown to universally occupy an inactive-state-specific tunnel between receptor helices 1 and 7; we also discover an unexpected role for an extra-helical groove accessible through this tunnel, suggesting its potential as a new targetable interface for CCR2 and CCR5 modulation. By contrast, only shape complementarity and limited helix 8 hydrogen bonding govern the binding of various chemotypes of allosteric antagonists. CCR2 residues S1012.63 and V2446.36 are implicated as determinants of CCR2/CCR5 and human/mouse orthosteric and allosteric antagonist selectivity, respectively, and the role of S1012.63 is corroborated through experimental gain-of-function mutagenesis. We establish a critical role of induced fit in antagonist recognition, reveal strong chemotype selectivity of existing structures, and demonstrate the high predictive potential of a new deep-learning-based compound scoring function. Finally, this study expands the available CCR2 structural landscape with computationally generated chemotype-specific models well-suited for structure-based antagonist design.

Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) make up the largest superfamily of human membrane proteins and represent primary targets of ∼1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon the binding of positive and negative allosteric modulators (PAMs and NAMs). The mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), deep learning (DL), and free energy prOfiling Workflow (GLOW). GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence and absence of the modulator. DL and free energy calculations revealed significantly reduced dynamic fluctuations and conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G-protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes. Therefore, GPCR allostery exhibits a dynamic "conformational selection" mechanism. In the absence of available modulator-bound structures as for most current GPCRs, it is critical to use a structural ensemble of representative GPCR conformations rather than a single structure for compound docking ("ensemble docking"), which will potentially improve structure-based design of novel allosteric drugs of GPCRs.

Single-molecule analysis reveals that a glucagon-bound extracellular domain of the glucagon receptor is dynamic

Dynamic information is vital to understanding the activation mechanism of G protein-coupled receptors (GPCRs). Despite the availability of high-resolution structures of different conformational states, the dynamics of those states at the molecular level are poorly understood. Here, we used total internal reflection fluorescence microscopy to study the extracellular domain (ECD) of the glucagon receptor (GCGR), a class B family GPCR that controls glucose homeostasis. Single-molecule fluorescence resonance energy transfer was used to observe the ECD dynamics of GCGR molecules expressed and purified from mammalian cells. We observed that for apo-GCGR, the ECD is dynamic and spent time predominantly in a closed conformation. In the presence of glucagon, the ECD is wide open and also shows more dynamic behavior than apo-GCGR, a finding that was not previously reported. These results suggest that both apo-GCGR and glucagon-bound GCGRs show reversible opening and closing of the ECD with respect to the seven-transmembrane (7TM) domain. This work demonstrates a molecular approach to visualizing the dynamics of the GCGR ECD and provides a foundation for understanding the conformational changes underlying GPCR activation, which is critical in the development of new therapeutics.

Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are the largest superfamily of human membrane proteins and represent primary targets of ~ 1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon binding of positive and negative allosteric modulators (PAMs and NAMs). Mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), Deep Learning (DL) and free energy prOfiling Workflow (GLOW). A total of 18 available high-resolution experimental structures of allosteric modulator-bound class A and B GPCRs were collected for simulations. A number of 8 computational models were generated to examine selectivity of the modulators by changing their target receptors to different subtypes. All-atom GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence/absence of the modulator. DL and free energy calculations revealed significantly reduced conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes in the computational models. Therefore, comprehensive DL of extensive GaMD simulations has revealed a general dynamic mechanism of GPCR allostery, which will greatly facilitate rational design of selective allosteric drugs of GPCRs.

Allosteric modulation of G protein-coupled receptor signaling

G protein-coupled receptors (GPCRs), the largest family of transmembrane proteins, regulate a wide array of physiological processes in response to extracellular signals. Although these receptors have proven to be the most successful class of drug targets, their complicated signal transduction pathways (including different effector G proteins and β-arrestins) and mediation by orthosteric ligands often cause difficulties for drug development, such as on- or off-target effects. Interestingly, identification of ligands that engage allosteric binding sites, which are different from classic orthosteric sites, can promote pathway-specific effects in cooperation with orthosteric ligands. Such pharmacological properties of allosteric modulators offer new strategies to design safer GPCR-targeted therapeutics for various diseases. Here, we explore recent structural studies of GPCRs bound to allosteric modulators. Our inspection of all GPCR families reveals recognition mechanisms of allosteric regulation. More importantly, this review highlights the diversity of allosteric sites and presents how allosteric modulators control specific GPCR pathways to provide opportunities for the development of new valuable agents.

Deep Learning Dynamic Allostery of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are the largest superfamily of human membrane proteins and represent primary targets of ~1/3 of currently marketed drugs. Allosteric modulators have emerged as more selective drug candidates compared with orthosteric agonists and antagonists. However, many X-ray and cryo-EM structures of GPCRs resolved so far exhibit negligible differences upon binding of positive and negative allosteric modulators (PAMs and NAMs). Mechanism of dynamic allosteric modulation in GPCRs remains unclear. In this work, we have systematically mapped dynamic changes in free energy landscapes of GPCRs upon binding of allosteric modulators using the Gaussian accelerated molecular dynamics (GaMD), Deep Learning (DL) and free energy prOfiling Workflow (GLOW). A total of 18 available high-resolution experimental structures of allosteric modulator-bound class A and B GPCRs were collected for simulations. A number of 8 computational models were generated to examine selectivity of the modulators by changing their target receptors to different subtypes. All-atom GaMD simulations were performed for a total of 66 μs on 44 GPCR systems in the presence/absence of the modulator. DL and free energy calculations revealed significantly reduced conformational space of GPCRs upon modulator binding. While the modulator-free GPCRs often sampled multiple low-energy conformational states, the NAMs and PAMs confined the inactive and active agonist-G protein-bound GPCRs, respectively, to mostly only one specific conformation for signaling. Such cooperative effects were significantly reduced for binding of the selective modulators to "non-cognate" receptor subtypes in the computational models. Therefore, comprehensive DL of extensive GaMD simulations has revealed a general dynamic mechanism of GPCR allostery, which will greatly facilitate rational design of selective allosteric drugs of GPCRs.

A reappraisal of the role of cyclic AMP in the physiological action of glucagon

Effects of the metabolic hormone glucagon can be physiological or supraphysiological, based on agonist concentration and the mediating cellular signal. The threshold concentration (TC) for activating the AC/cAMP signal pathway in liver is ≥ 100 pM. By contrast, mean plasma concentrations are around 20-45 pM, depending on the vascular bed. Accordingly, effects produced at TCs below 100 pM are physiological and mediated by cellular signal pathways other than AC/cAMP. Effects generated at concentrations above 100 pM are supraphysiological, often mediated by simultaneous activation of cAMP-independent and -dependent pathways. Physiological responses, and their established or implicated signal pathways, include stimulation of: glucose mobilization, fatty acid oxidation, and urea synthesis in liver (PLC/IP3/Ca²⁺/CaM); lipolysis in white and brown adipose tissue and oxygen consumption in brown adipose of the rat but not in humans (PLC/IP3/Ca²⁺/CaM); renal potassium and phosphate excretion in rodents and GFR in humans (signal undetermined); and glucose utilization in rat heart (PI3K/akt). Supraphysiological responses involve the AC/cAMP pathway and include: enhanced stimulation of glucose mobilization and stimulation of urea synthesis in liver; further stimulation of white and brown adipose lipolysis and thermogenesis in brown adipose tissue; stimulation of renal Cl^- transport; and increased rat heart contractility. The AC/cAMP pathway is likely recruited when plasma glucagon rises above 100 pM during periods of elevated metabolic stress and systemic glucose demand, such as in the early neonate or strenuously exercising adult. The current cAMP-centered model should therefore be reconsidered and replaced with one that places more emphasis on cAMP-independent pathways.

New insights into the structure and function of class B1 GPCRs

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors. Class B1 GPCRs constitute a subfamily of 15 receptors that characteristically contain large extracellular domains (ECDs) and respond to long polypeptide hormones. Class B1 GPCRs are critical regulators of homeostasis, and as such, many are important drug targets. While most transmembrane proteins, including GPCRs, are recalcitrant to crystallization, recent advances in electron cryo-microscopy (cryo-EM) have facilitated a rapid expansion of the structural understanding of membrane proteins. As a testament to this success, structures for all the class B1 receptors bound to G proteins have been determined by cryo-EM in the past five years. Further advances in cryo-EM have uncovered dynamics of these receptors, ligands, and signalling partners. Here, we examine the recent structural underpinnings of the class B1 GPCRs with an emphasis on structure-function relationships.

Structural basis for recognition of antihistamine drug by human histamine receptor

The histamine receptors belong to the G protein-coupled receptor (GPCR) superfamily, and play important roles in the regulation of histamine and other neurotransmitters in the central nervous system, as potential targets for the treatment of neurologic and psychiatric disorders. Here we report the crystal structure of human histamine receptor H₃R bound to an antagonist PF-03654746 at 2.6 Å resolution. Combined with the computational and functional assays, our structure reveals binding modes of the antagonist and allosteric cholesterol. Molecular dynamic simulations and molecular docking of different antihistamines further elucidate the conserved ligand-binding modes. These findings are therefore expected to facilitate the structure-based design of novel antihistamines.

Intermolecular Interactions in G Protein-Coupled Receptor Allosteric Sites at the Membrane Interface from Molecular Dynamics Simulations and Quantum Chemical Calculations

Allosteric modulators are called promising candidates in G protein-coupled receptor (GPCR) drug development by displaying subtype selectivity and more specific receptor modulation. Among the allosteric sites known to date, cavities at the receptor-lipid interface represent an uncharacteristic binding location that raises many questions about the ligand interactions and stability, the binding site structure, and how all of these are affected by lipid molecules. In this work, we analyze interactions in the allosteric sites of the PAR2, C5aR1, and GCGR receptors in three lipid compositions using molecular dynamics simulations. In addition, we performed quantum chemical calculations involving the symmetry-adapted perturbation theory (SAPT) and the natural population analysis to quantify the strength of intermolecular interactions. We show that besides classical hydrogen bonds, weak polar interactions such as O-HC, O-Br, and long-range electrostatics with the backbone amides contribute to the stability of allosteric modulators at the receptor-lipid interface. The allosteric cavities are detectable in various membrane compositions. The availability of polar atoms for interactions in such cavities can be assessed by water molecules from simulations. Although ligand-lipid interactions are weak, lipid tails play a role in ligand binding pose stability and the size of allosteric cavities. We discuss physicochemical aspects of ligand binding at the receptor-lipid interface and suggest a compound library enriched by weak donor groups for ligand search in such sites.

Hepatocyte cholesterol content modulates glucagon receptor signalling

OBJECTIVE

To determine whether glucagon receptor (GCGR) actions are modulated by cellular cholesterol levels.

METHODS

We determined the effects of experimental cholesterol depletion and loading on glucagon-mediated cAMP production, ligand internalisation and glucose production in human hepatoma cells, mouse and human hepatocytes. GCGR interactions with lipid bilayers were explored using coarse-grained molecular dynamic simulations. Glucagon responsiveness was measured in mice fed a high cholesterol diet with or without simvastatin to modulate hepatocyte cholesterol content.

RESULTS

GCGR cAMP signalling was reduced by higher cholesterol levels across different cellular models. Ex vivo glucagon-induced glucose output from mouse hepatocytes was enhanced by simvastatin treatment. Mice fed a high cholesterol diet had increased hepatic cholesterol and a blunted hyperglycaemic response to glucagon, both of which were partially reversed by simvastatin. Simulations identified likely membrane-exposed cholesterol binding sites on the GCGR, including a site where cholesterol is a putative negative allosteric modulator.

CONCLUSIONS

Our results indicate that cellular cholesterol content influences glucagon sensitivity and indicate a potential molecular basis for this phenomenon. This could be relevant to the pathogenesis of non-alcoholic fatty liver disease, which is associated with both hepatic cholesterol accumulation and glucagon resistance.

Antagonistic Effect and In Vitro Activity of Dauricine on Glucagon Receptor

Abnormal increases in glucagon (GCG) are the primary cause of type II diabetes mellitus. When GCG interacts with a glucagon receptor (GCGR), GCG can increase the blood glucose level. In this paper, a compound that could interfere with the binding of GCG and GCGR to inhibit the increase of blood glucose was investigated. First, molecular docking was used to conduct preliminary screening of compounds whose active components could combine with GCGR by AutoDock Vina. The binding of the receptor-ligand complex was analyzed by PyMOL. Results showed that dauricine could tightly bind to the receptor pocket. Second, the plasmid pcDNA3.1(+)-GCGR containing the target gene was transfected into HEK293 cells for expression, which was the cell model established to screen GCGR antagonist. Dauricine, the lead compound of glucagon receptor antagonist (GRA), was screened using the GRA screening model in vitro. Finally, using [Des-His1, Glu9]-Glucagon amide as the positive control, flow cytometry was used to express the antagonistic effect of the compound. Consequently, dauricine can antagonize the GCGR.

Targeting VIP and PACAP Receptor Signaling: New Insights into Designing Drugs for the PACAP Subfamily of Receptors

Pituitary Adenylate Cyclase-Activating Peptide (PACAP) and Vasoactive Intestinal Peptide (VIP) are neuropeptides involved in a diverse array of physiological and pathological processes through activating the PACAP subfamily of class B1 G protein-coupled receptors (GPCRs): VIP receptor 1 (VPAC1R), VIP receptor 2 (VPAC2R), and PACAP type I receptor (PAC1R). VIP and PACAP share nearly 70% amino acid sequence identity, while their receptors PAC1R, VPAC1R, and VPAC2R share 60% homology in the transmembrane regions of the receptor. PACAP binds with high affinity to all three receptors, while VIP binds with high affinity to VPAC1R and VPAC2R, and has a thousand-fold lower affinity for PAC1R compared to PACAP. Due to the wide distribution of VIP and PACAP receptors in the body, potential therapeutic applications of drugs targeting these receptors, as well as expected undesired side effects, are numerous. Designing selective therapeutics targeting these receptors remains challenging due to their structural similarities. This review discusses recent discoveries on the molecular mechanisms involved in the selectivity and signaling of the PACAP subfamily of receptors, and future considerations for therapeutic targeting.

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

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

Probe Confined Dynamic Mapping for G Protein-Coupled Receptor Allosteric Site Prediction

Targeting G protein-coupled receptors (GPCRs) through allosteric sites offers advantages over orthosteric sites in identifying drugs with increased selectivity and potentially reduced side effects. In this study, we developed a probe confined dynamic mapping protocol that allows the prediction of allosteric sites at both the GPCR extracellular and intracellular sides, as well as at the receptor-lipid interface. The applied harmonic wall potential enhanced sampling of probe molecules in a selected area of a GPCR while preventing membrane distortion in molecular dynamics simulations. The specific probes derived from GPCR allosteric ligand structures performed better in allosteric site mapping compared to commonly used cosolvents. The M₂ muscarinic, β₂ adrenergic, and P₂Y₁ purinergic receptors were selected for the protocol's retrospective validation. The protocol was next validated prospectively to locate the binding site of [5-fluoro-4-(hydroxymethyl)-2-methoxyphenyl]-(4-fluoro-1H-indol-1-yl)methanone at the D₂ dopamine receptor, and subsequent mutagenesis confirmed the prediction. The protocol provides fast and efficient prediction of key amino acid residues surrounding allosteric sites in membrane proteins and facilitates the structure-based design of allosteric modulators.

Molecular and in vivo phenotyping of missense variants of the human glucagon receptor

Naturally occurring missense variants of G protein–coupled receptors with loss of function have been linked to metabolic disease in case studies and in animal experiments. The glucagon receptor, one such G protein–coupled receptor, is involved in maintaining blood glucose and amino acid homeostasis; however, loss-of-function mutations of this receptor have not been systematically characterized. Here, we observed fewer glucagon receptor missense variants than expected, as well as lower allele diversity and fewer variants with trait associations as compared with other class B1 receptors. We performed molecular pharmacological phenotyping of 38 missense variants located in the receptor extracellular domain, at the glucagon interface, or with previously suggested clinical implications. These variants were characterized in terms of cAMP accumulation to assess glucagon-induced Gα_s coupling, and of recruitment of β-arrestin-1/2. Fifteen variants were impaired in at least one of these downstream functions, with six variants affected in both cAMP accumulation and β-arrestin-1/2 recruitment. For the eight variants with decreased Gα_s signaling (D63^ECDN, P86^ECDS, V96^ECDE, G125^ECDC, R225^3.30H, R308^5.40W, V368^6.59M, and R378^7.35C) binding experiments revealed preserved glucagon affinity, although with significantly reduced binding capacity. Finally, using the UK Biobank, we found that variants with wildtype-like Gα_s signaling did not associate with metabolic phenotypes, whereas carriers of cAMP accumulation-impairing variants displayed a tendency toward increased risk of obesity and increased body mass and blood pressure. These observations are in line with the essential role of the glucagon system in metabolism and support that Gα_s is the main signaling pathway effecting the physiological roles of the glucagon receptor.

Synthesis and anti-diabetic activity of novel biphenylsulfonamides as glucagon receptor antagonists

Type 2 diabetes is characterized by chronic hyperglycemia. Insulin, a hormone secreted from pancreatic β-cells, decreases blood glucose levels, and glucagon, a hormone secreted from pancreatic α-cells, increases blood glucose levels by counterregulation of insulin through stimulation of hepatic glucose production. In diabetic patients, dysregulation of glucagon secretion contributes to hyperglycemia. Thus, inhibition of the glucagon receptor is one strategy for the treatment of hyperglycemia in type 2 diabetes. In this paper, we report a series of biphenylsulfonamide derivatives that were designed, synthesized, and then evaluated by cAMP and hepatic glucose production assays as glucagon receptor antagonists. Of these, compound 7aB-3 decreased glucagon-induced cAMP production and glucagon-induced glucose production in the in vitro assays. Glucagon challenge tests and glucose tolerance tests showed that compound 7aB-3 significantly inhibited glucagon-induced glucose increases and improved glucose tolerance. These results suggest that compound 7aB-3 has therapeutic potential for the treatment of type 2 diabetes.

The Evolutionary History of Vertebrate Adhesion GPCRs and Its Implication on Their Classification

Adhesion G protein-coupled receptors (aGPCRs) form a structurally separate class of GPCRs with an unresolved evolutionary history and classification. Based on phylogenetic relations of human aGPCRs, nine families (A-G, L, V) were distinguished. Taking advantage of available genome data, we determined the aGPCR repertoires in all vertebrate classes. Although most aGPCR families show a high numerical stability in vertebrate genomes, the full repertoire of family E, F, and G members appeared only after the fish-tetrapod split. We did not find any evidence for new aGPCR families in vertebrates which are not present in the human genome. Based on ortholog sequence alignments, selection analysis clearly indicated two types of tetrapod aGPCRs: (i) aGPCR under strong purifying selection in tetrapod evolution (families A, B, D, L, V); and (ii) aGPCR with signatures of positive selection in some tetrapod linages (families C, E, G, F). The alignments of aGPCRs also allowed for a revised definition of reference positions within the seven-transmembrane-helix domain (relative position numbering scheme). Based on our phylogenetic cluster analysis, we suggest a revised nomenclature of aGPCRs including their transcript variants. Herein, the former families E and L are combined to one family (L) and GPR128/ADGRG7 forms a separate family (E). Furthermore, our analyses provide valuable information about the (patho)physiological relevance of individual aGPCR members.

Identification of a Novel Allosteric Site at the M5 Muscarinic Acetylcholine Receptor

The M₅ muscarinic acetylcholine receptor (mAChR) has emerged as an exciting therapeutic target for the treatment of addiction and behavioral disorders. This has been in part due to promising preclinical studies with the M₅ mAChR selective negative allosteric modulator (NAM), ML375. The binding site of ML375 remains unknown, however, making it difficult to develop improved M₅ mAChR selective modulators. To determine the possible location of the ML375 binding site, we used radioligand binding and functional assays to show that ML375 does not interact with the well-characterized "common" mAChR allosteric site located in the receptor's extracellular vestibule, nor a previously proposed second allosteric site recognized by the modulator, amiodarone. Molecular docking was used to predict potential allosteric sites within the transmembrane (TM) domain of the M₅ mAChR. These predicted sites were assessed using M₅-M₂ mAChR receptor chimeras and further targeted with site-directed mutagenesis, which enabled the identification of a putative binding site for ML375 at the interface of TMs 2-4. Collectively, these results identify a third allosteric site at the M₅ mAChR and highlight the ability of allosteric modulators to selectively target highly conserved proteins.

Ligand binding at the protein-lipid interface: strategic considerations for drug design

Many drug targets are embedded within the phospholipid bilayer of cellular membranes, including G protein-coupled receptors, ion channels, transporters and membrane-bound enzymes. Increasing evidence from biophysical and structural studies suggests that many small-molecule drugs commonly associate with these targets at binding sites at the protein-phospholipid interface. Without a direct path from bulk solvent to a binding site, a drug must first partition in the phospholipid membrane before interacting with the protein target. This membrane access mechanism necessarily affects the interpretation of potency data, structure-activity relationships, pharmacokinetics and physicochemical properties for drugs that target these sites. With an increasing number of small-molecule intramembrane binding sites revealed through X-ray crystallography and cryogenic electron microscopy, we suggest that ligand-lipid interactions likely play a larger role in small-molecule drug action than commonly appreciated. This Perspective introduces key concepts and drug design considerations to aid discovery teams operating within this target space, and discusses challenges and future opportunities in the field.

Structural insights into hormone recognition by the human glucose-dependent insulinotropic polypeptide receptor

Glucose-dependent insulinotropic polypeptide (GIP) is a peptide hormone that exerts crucial metabolic functions by binding and activating its cognate receptor, GIPR. As an important therapeutic target, GIPR has been subjected to intensive structural studies without success. Here, we report the cryo-EM structure of the human GIPR in complex with GIP and a Gs heterotrimer at a global resolution of 2.9 Å. GIP adopts a single straight helix with its N terminus dipped into the receptor transmembrane domain (TMD), while the C terminus is closely associated with the extracellular domain and extracellular loop 1. GIPR employs conserved residues in the lower half of the TMD pocket to recognize the common segments shared by GIP homologous peptides, while uses non-conserved residues in the upper half of the TMD pocket to interact with residues specific for GIP. These results provide a structural framework of hormone recognition and GIPR activation.

Discovery of Novel Allosteric Modulators Targeting an Extra-Helical Binding Site of GLP-1R Using Structure- and Ligand-Based Virtual Screening

Allosteric modulators have emerged with many potential pharmacological advantages as they do not compete the binding of agonist or antagonist to the orthosteric sites but ultimately affect downstream signaling. To identify allosteric modulators targeting an extra-helical binding site of the glucagon-like peptide-1 receptor (GLP-1R) within the membrane environment, the following two computational approaches were applied: structure-based virtual screening with consideration of lipid contacts and ligand-based virtual screening with the maintenance of specific allosteric pocket residue interactions. Verified by radiolabeled ligand binding and cAMP accumulation experiments, two negative allosteric modulators and seven positive allosteric modulators were discovered using structure-based and ligand-based virtual screening methods, respectively. The computational approach presented here could possibly be used to discover allosteric modulators of other G protein-coupled receptors.

Crystal structure of dopamine D1 receptor in complex with G protein and a non-catechol agonist

Dopamine D1 receptor (D1R) is an important drug target implicated in many psychiatric and neurological disorders. Selective agonism of D1R are sought to be the therapeutic strategy for these disorders. Most selective D1R agonists share a dopamine-like catechol moiety in their molecular structure, and their therapeutic potential is therefore limited by poor pharmacological properties in vivo. Recently, a class of non-catechol D1R selective agonists with a distinct scaffold and pharmacological properties were reported. Here, we report the crystal structure of D1R in complex with stimulatory G protein (Gs) and a non-catechol agonist Compound 1 at 3.8 Å resolution. The structure reveals the ligand bound to D1R in an extended conformation, spanning from the orthosteric site to extracellular loop 2 (ECL2). Structural analysis reveals that the unique features of D1R ligand binding pocket explains the remarkable selectivity of this scaffold for D1R over other aminergic receptors, and sheds light on the mechanism for D1R activation by the non-catechol agonist.

MolAICal: a soft tool for 3D drug design of protein targets by artificial intelligence and classical algorithm

Deep learning is an important branch of artificial intelligence that has been successfully applied into medicine and two-dimensional ligand design. The three-dimensional (3D) ligand generation in the 3D pocket of protein target is an interesting and challenging issue for drug design by deep learning. Here, the MolAICal software is introduced to supply a way for generating 3D drugs in the 3D pocket of protein targets by combining with merits of deep learning model and classical algorithm. The MolAICal software mainly contains two modules for 3D drug design. In the first module of MolAICal, it employs the genetic algorithm, deep learning model trained by FDA-approved drug fragments and Vinardo score fitting on the basis of PDBbind database for drug design. In the second module, it uses deep learning generative model trained by drug-like molecules of ZINC database and molecular docking invoked by Autodock Vina automatically. Besides, the Lipinski's rule of five, Pan-assay interference compounds (PAINS), synthetic accessibility (SA) and other user-defined rules are introduced for filtering out unwanted ligands in MolAICal. To show the drug design modules of MolAICal, the membrane protein glucagon receptor and non-membrane protein SARS-CoV-2 main protease are chosen as the investigative drug targets. The results show MolAICal can generate the various and novel ligands with good binding scores and appropriate XLOGP values. We believe that MolAICal can use the advantages of deep learning model and classical programming for designing 3D drugs in protein pocket. MolAICal is freely for any nonprofit purpose and accessible at https://molaical.github.io.

Ligand-Receptor Interactions and Machine Learning in GCGR and GLP-1R Drug Discovery

The large amount of data that has been collected so far for G protein-coupled receptors requires machine learning (ML) approaches to fully exploit its potential. Our previous ML model based on gradient boosting used for prediction of drug affinity and selectivity for a receptor subtype was compared with explicit information on ligand-receptor interactions from induced-fit docking. Both methods have proved their usefulness in drug response predictions. Yet, their successful combination still requires allosteric/orthosteric assignment of ligands from datasets. Our ligand datasets included activities of two members of the secretin receptor family: GCGR and GLP-1R. Simultaneous activation of two or three receptors of this family by dual or triple agonists is not a typical kind of information included in compound databases. A precise allosteric/orthosteric ligand assignment requires a continuous update based on new structural and biological data. This data incompleteness remains the main obstacle for current ML methods applied to class B GPCR drug discovery. Even so, for these two class B receptors, our ligand-based ML model demonstrated high accuracy (5-fold cross-validation Q² > 0.63 and Q² > 0.67 for GLP-1R and GCGR, respectively). In addition, we performed a ligand annotation using recent cryogenic-electron microscopy (cryo-EM) and X-ray crystallographic data on small-molecule complexes of GCGR and GLP-1R. As a result, we assigned GLP-1R and GCGR actives deposited in ChEMBL to four small-molecule binding sites occupied by positive and negative allosteric modulators and a full agonist. Annotated compounds were added to our recently released repository of GPCR data.

Optimization of Truncated Glucagon Peptides to Achieve Selective, High Potency, Full Antagonists

Antagonism of glucagon's biological action is a proven strategy for decreasing glucose in diabetic animals and patients. To achieve full, potent, and selective suppression, we chemically optimized N-terminally truncated glucagon fragments for the identification and establishment of the minimum sequence peptide, [Glu9]glucagon(6-29) amide (11) as a full antagonist in cellular signaling and receptor binding (IC₅₀ = 36 nM). Substitution of Phe6 with l-3-phenyllactic acid (Pla) produced [Pla6, Glu9]glucagon(6-29) amide (21), resulting in a 3-fold improvement in receptor binding (IC₅₀ = 12 nM) and enhanced antagonist potency. Further substitution of Glu9 and Asn28 with aspartic acid yielded [Pla6, Asp28]glucagon amide (26), which demonstrated a further increase in inhibitory potency (IC₅₀ = 9 nM), and improved aqueous solubility. Peptide 26 and a palmitoylated analogue, [Pla6, Lys10(γGluγGlu-C16), Asp28]glucagon(6-29) amide (31), displayed sustained duration in vivo action that successfully reversed glucagon-induced glucose elevation in mice.

Evidence for Glucagon Secretion and Function Within the Human Gut

Glucagon is secreted by pancreatic α cells in response to hypoglycemia and increases hepatic glucose output through hepatic glucagon receptors (GCGRs). There is evidence supporting the notion of extrapancreatic glucagon but its source and physiological functions remain elusive. Intestinal tissue samples were obtained from patients undergoing surgical resection of cancer. Mass spectrometry analysis was used to detect glucagon from mucosal lysate. Static incubations of mucosal tissue were performed to assess glucagon secretory response. Glucagon concentration was quantitated using a highly specific sandwich enzyme-linked immunosorbent assay. A cholesterol uptake assay and an isolated murine colonic motility assay were used to assess the physiological functions of intestinal GCGRs. Fully processed glucagon was detected by mass spectrometry in human intestinal mucosal lysate. High glucose evoked significant glucagon secretion from human ileal tissue independent of sodium glucose cotransporter and KATP channels, contrasting glucose-induced glucagon-like peptide 1 (GLP-1) secretion. The GLP-1 receptor agonist Exendin-4 attenuated glucose-induced glucagon secretion from the human ileum. GCGR blockade significantly increased cholesterol uptake in human ileal crypt culture and markedly slowed ex vivo colonic motility. Our findings describe the human gut as a potential source of extrapancreatic glucagon and demonstrate a novel enteric glucagon/GCGR circuit with important physiological functions beyond glycemic regulation.

Molecular Basis of Class B GPCR Selectivity for the Neuropeptides PACAP and VIP

The related neuropeptides PACAP and VIP, and their shared PAC1, VPAC1 and VPAC2 receptors, regulate a large array of physiological activities in the central and peripheral nervous systems. However, the lack of comparative and molecular mechanistic investigations hinder further understanding of their preferred binding selectivity and function. PACAP and VIP have comparable affinity at the VPAC1 and VPAC2 receptor, but PACAP is 400-1,000 fold more potent than VIP at the PAC1 receptor. A molecular understanding of the differing neuropeptide-receptor interactions and the details underlying the receptor transitions leading to receptor activation are much needed for the rational design of selective ligands. To these ends, we have combined structural information and advanced simulation techniques to study PACAP/VIP binding selectivity, full-length receptor conformation ensembles and transitions of the PACAP/VIP receptor variants and subtypes, and a few key interactions in the orthosteric-binding pocket. Our results reveal differential peptide-receptor interactions (at the atomistic detail) important for PAC1, VPAC1 and VPAC2 receptor ligand selectivity. Using microsecond-long molecular dynamics simulations and the Markov State Models, we have also identified diverse receptor conformational ensembles and microstate transition paths for each receptor, the potential mechanisms underlying receptor open and closed states, and the interactions and dynamics at the transmembrane orthosteric pocket for receptor activation. These analyses reveal important features in class B GPCR structure-dynamics-function relationships, which provide novel insights for structure-based drug discovery.

Engineering of Challenging G Protein-Coupled Receptors for Structure Determination and Biophysical Studies

Membrane proteins such as G protein-coupled receptors (GPCRs) exert fundamental biological functions and are involved in a multitude of physiological responses, making these receptors ideal drug targets. Drug discovery programs targeting GPCRs have been greatly facilitated by the emergence of high-resolution structures and the resulting opportunities to identify new chemical entities through structure-based drug design. To enable the determination of high-resolution structures of GPCRs, most receptors have to be engineered to overcome intrinsic hurdles such as their poor stability and low expression levels. In recent years, multiple engineering approaches have been developed to specifically address the technical difficulties of working with GPCRs, which are now beginning to make more challenging receptors accessible to detailed studies. Importantly, successfully engineered GPCRs are not only valuable in X-ray crystallography, but further enable biophysical studies with nuclear magnetic resonance spectroscopy, surface plasmon resonance, native mass spectrometry, and fluorescence anisotropy measurements, all of which are important for the detailed mechanistic understanding, which is the prerequisite for successful drug design. Here, we summarize engineering strategies based on directed evolution to reduce workload and enable biophysical experiments of particularly challenging GPCRs.

Molecular mechanism of an antagonistic antibody against glucose-dependent insulinotropic polypeptide receptor

Glucose-dependent insulinotropic polypeptide (GIP) is an incretin hormone involved in regulating glucose and lipid metabolism. GIP receptor (GIPR) antagonism is believed to offer therapeutic potential for various metabolic diseases. Pharmacological intervention of GIPR, however, has limited success due to lack of effective antagonistic reagents. Previously we reported the discovery of two mouse anti-murine GIPR monoclonal antibodies (mAbs) with distinctive properties in rodent models. Here, we report the detailed structural and biochemical characterization of these two antibodies, mAb1 and mAb2. In vitro and in vivo characterizations demonstrated mAb2 is a full GIPR antagonistic antibody and mAb1 is a non-neutralizing GIPR binder. To understand the molecular basis of these two antibodies, we determined the co-crystal structures of GIPR extracellular domain in complex with mAb1 and with mAb2 at resolutions of 2.1 and 2.6 Å, respectively. While the non-neutralizing mAb1 binds to GIPR without competing with the ligand peptide, mAb2 not only partially occludes the ligand peptide binding, but also recognizes the GIPR C-terminal stalk region in a helical conformation that acts as a molecular mimic of the ligand peptide and locks GIPR in a novel auto-inhibited state. Furthermore, administration of mAb2 in diet-induced obesity mice for 7 weeks leads to both reduction in body weight gain and improvement of metabolic profiles. In contrast, mAb1 has no effect on body weight or other metabolic improvement. Together, our studies reveal the unique molecular mechanism of action underlying the superior antagonistic activity of mAb2 and signify the promising therapeutic potential of effective GIPR antagonism for the treatment of metabolic disorders.

G protein-coupled receptors: structure- and function-based drug discovery

As one of the most successful therapeutic target families, G protein-coupled receptors (GPCRs) have experienced a transformation from random ligand screening to knowledge-driven drug design. We are eye-witnessing tremendous progresses made recently in the understanding of their structure-function relationships that facilitated drug development at an unprecedented pace. This article intends to provide a comprehensive overview of this important field to a broader readership that shares some common interests in drug discovery.

The Role of ICL1 and H8 in Class B1 GPCRs; Implications for Receptor Activation

The first intracellular loop (ICL1) of G protein-coupled receptors (GPCRs) has received little attention, although there is evidence that, with the 8^th helix (H8), it is involved in early conformational changes following receptor activation as well as contacting the G protein β subunit. In class B1 GPCRs, the distal part of ICL1 contains a conserved R^12.48KLRCxR^2.46b motif that extends into the base of the second transmembrane helix; this is weakly conserved as a [R/H]^12.48KL[R/H] motif in class A GPCRs. In the current study, the role of ICL1 and H8 in signaling through cAMP, _iCa²⁺ and ERK1/2 has been examined in two class B1 GPCRs, using mutagenesis and molecular dynamics. Mutations throughout ICL1 can either enhance or disrupt cAMP production by CGRP at the CGRP receptor. Alanine mutagenesis identified subtle differences with regard elevation of _iCa²⁺, with the distal end of the loop being particularly sensitive. ERK1/2 activation displayed little sensitivity to ICL1 mutation. A broadly similar pattern was observed with the glucagon receptor, although there were differences in significance of individual residues. Extending the study revealed that at the CRF1 receptor, an insertion in ICL1 switched signaling bias between _iCa²⁺ and cAMP. Molecular dynamics suggested that changes in ICL1 altered the conformation of ICL2 and the H8/TM7 junction (ICL4). For H8, alanine mutagenesis showed the importance of E390^8.49b for all three signal transduction pathways, for the CGRP receptor, but mutations of other residues largely just altered ERK1/2 activation. Thus, ICL1 may modulate GPCR bias via interactions with ICL2, ICL4 and the Gβ subunit.

Drug Repositioning For Allosteric Modulation of VIP and PACAP Receptors

Vasoactive intestinal peptide (VIP) and pituitary adenylate cyclase-activating polypeptide (PACAP) are two neuropeptides that contribute to the regulation of intestinal motility and secretion, exocrine and endocrine secretions, and homeostasis of the immune system. Their biological effects are mediated by three receptors named VPAC1, VPAC2 and PAC1 that belong to class B GPCRs. VIP and PACAP receptors have been identified as potential therapeutic targets for the treatment of chronic inflammation, neurodegenerative diseases and cancer. However, pharmacological use of endogenous ligands for these receptors is limited by their lack of specificity (PACAP binds with high affinity to VPAC1, VPAC2 and PAC1 receptors while VIP recognizes both VPAC1 and VPAC2 receptors), their poor oral bioavailability (VIP and PACAP are 27- to 38-amino acid peptides) and their short half-life. Therefore, the development of non-peptidic small molecules or specific stabilized peptidic ligands is of high interest. Structural similarities between VIP and PACAP receptors are major causes of difficulties in the design of efficient and selective compounds that could be used as therapeutics. In this study we performed structure-based virtual screening against the subset of the ZINC15 drug library. This drug repositioning screen provided new applications for a known drug: ticagrelor, a P2Y12 purinergic receptor antagonist. Ticagrelor inhibits both VPAC1 and VPAC2 receptors which was confirmed in VIP-binding and calcium mobilization assays. A following analysis of detailed ticagrelor binding modes to all three VIP and PACAP receptors with molecular dynamics revealed its allosteric mechanism of action. Using a validated homology model of inactive VPAC1 and a recently released cryo-EM structure of active VPAC1 we described how ticagrelor could block conformational changes in the region of 'tyrosine toggle switch' required for the receptor activation. We also discuss possible modifications of ticagrelor comparing other P2Y12 antagonist - cangrelor, closely related to ticagrelor but not active for VPAC1/VPAC2. This comparison with inactive cangrelor could lead to further improvement of the ticagrelor activity and selectivity for VIP and PACAP receptor sub-types.

Design and preparation of the class B G protein-coupled receptors GLP-1R and GCGR for 19 F-NMR studies in solution

The human glucagon-like peptide-1 receptor (GLP-1R) and the glucagon receptor (GCGR) are class B G protein-coupled receptors (GPCRs) that are activated by interactions with, respectively, the glucagon-like peptide-1 (GLP-1) and glucagon (GCG). These polypeptide hormones are involved in the regulation of lipid and cholic acid metabolism, and thus play an important role in the pathogenesis of glucose metabolism and diabetes mellitus, which attracts keen interest of these GPCRs as drug targets. GLP-1R and GCGR have therefore been extensively investigated by X-ray crystallography and cryo-electron microscopy (cryo-EM), so that their structures are well known. Here, we present the groundwork for using nuclear magnetic resonance (NMR) spectroscopy in solution to complement the molecular architectures with information on intramolecular dynamics and on the thermodynamics and kinetics of interactions with physiological ligands and extrinsic drug candidates. This includes the generation of novel, near-wild-type constructs of GLP-1R and GCGR, optimization of the solution conditions for NMR studies in detergent micelles and in nanodiscs, post-translational chemical introduction of fluorine-19 NMR probes, and sequence-specific assignments of the ¹⁹ F-labels attached to indigenous cysteines. Addition of the negative allosteric modulator (NAM) NNC0640 was critically important for obtaining the long-time stability needed for our NMR experiments, and we report on novel insights into the allosteric effects arising from binding of NNC0640 to the transmembrane domain of GLP-1R (GLP-1R[TMD]).

Dynamical Correlations Reveal Allosteric Sites in G Protein-Coupled Receptors

G protein-coupled Receptors (GPCRs) play a central role in many physiological processes and, consequently, constitute important drug targets. In particular, the search for allosteric drugs has recently drawn attention, since they could be more selective and lead to fewer side effects. Accordingly, computational tools have been used to estimate the druggability of allosteric sites in these receptors. In spite of many successful results, the problem is still challenging, particularly the prediction of hydrophobic sites in the interface between the protein and the membrane. In this work, we propose a complementary approach, based on dynamical correlations. Our basic hypothesis was that allosteric sites are strongly coupled to regions of the receptor that undergo important conformational changes upon activation. Therefore, using ensembles of experimental structures, normal mode analysis and molecular dynamics simulations we calculated correlations between internal fluctuations of different sites and a collective variable describing the activation state of the receptor. Then, we ranked the sites based on the strength of their coupling to the collective dynamics. In the β2 adrenergic (β2AR), glucagon (GCGR) and M2 muscarinic receptors, this procedure allowed us to correctly identify known allosteric sites, suggesting it has predictive value. Our results indicate that this dynamics-based approach can be a complementary tool to the existing toolbox to characterize allosteric sites in GPCRs.

Protease-activated receptor-2 ligands reveal orthosteric and allosteric mechanisms of receptor inhibition

Protease-activated receptor-2 (PAR2) has been implicated in multiple pathophysiologies but drug discovery is challenging due to low small molecule tractability and a complex activation mechanism. Here we report the pharmacological profiling of a potent new agonist, suggested by molecular modelling to bind in the putative orthosteric site, and two novel PAR2 antagonists with distinctly different mechanisms of inhibition. We identify coupling between different PAR2 binding sites. One antagonist is a competitive inhibitor that binds to the orthosteric site, while a second antagonist is a negative allosteric modulator that binds at a remote site. The allosteric modulator shows probe dependence, more effectively inhibiting peptide than protease activation of PAR2 signalling. Importantly, both antagonists are active in vivo, inhibiting PAR2 agonist-induced acute paw inflammation in rats and preventing activation of mast cells and neutrophils. These results highlight two distinct mechanisms of inhibition that potentially could be targeted for future development of drugs that modulate PAR2.

Kennedy et al. report the pharmacological and in vivo profiling of two small molecule PAR2 inhibitors and an agonist. They conclude that while the small molecule agonist and one of the inhibitors bind to the orthosteric PAR2 binding site, the other inhibitor is a negative allosteric modulator, highlighting two distinct mechanisms of inhibition that could be targeted for future development of drugs that modulate PAR2.

Allosteric binding sites at the receptor-lipid bilayer interface: novel targets for GPCR drug discovery

As a superfamily of membrane receptors, G-protein-coupled receptors (GPCRs) have significant roles in human physiological processes, including cell proliferation, metabolism, and neuromodulation. GPCRs are vital targets of therapeutic drugs, and their allosteric regulation represents a novel direction for drug discovery. Given the numerous breakthroughs in structural biology, diverse allosteric sites on GPCRs have been identified within the extracellular and intracellular loops, and the seven core transmembrane helices. However, a unique type of allosteric site has also been discovered at the interface of the receptor-lipid bilayer, similar to the β₂-adrenergic receptor. Here, we review recent identifications of these allosteric sites and the detailed modulator-target interactions within the interface for each modulator to highlight the role of lipids in GPCR allosteric drug discovery.

Covalent allosteric modulation: An emerging strategy for GPCRs drug discovery

Designing covalent allosteric modulators brings new opportunities to the field of drug discovery towards G-protein-coupled receptors (GPCRs). Targeting an allosteric binding pocket can allow a modulator to have protein subtype selectivity and low drug resistance. Utilizing covalent warheads further enables the modulator to increase the binding potency and extend the duration of action. This review starts with GPCR allosteric modulation to discuss the structural biology of allosteric binding pockets, the different types of allosteric modulators, as well as the advantages of employing allosteric modulation. This is followed by a discussion on covalent modulators to clarify how covalent ligands can benefit the receptor modulation and to illustrate moieties that can commonly be used as covalent warheads. Finally, case studies are presented on designing class A, B, and C GPCR covalent allosteric modulators to demonstrate successful stories on combining allosteric modulation and covalent binding. Limitations and future perspectives are also covered.

Glucagon and Glucagon-like Peptide-1 Receptors: Promising Therapeutic Targets for an Effective Management of Diabetes Mellitus

G-protein-coupled receptors (GPCRs) are membrane-bound proteins, which are responsible for the detection of extracellular stimuli and the origination of intracellular responses. Both glucagon and glucagon-like peptide-1 (GLP-1) receptors belong to G protein-coupled receptor (GPCR) superfamily. Along with insulin, glucagon and GLP-1 are critical hormones for maintaining normal serum glucose within the human body. Glucagon generally plays its role in the liver through cyclic adenosine monophosphate (cAMP), where it compensates for the action of insulin. GLP-1 is secreted by the L-cells of the small intestine to stimulate insulin secretion and inhibit glucagon action. Despite extensive research efforts and the multiple approaches adopted, the glycemic control in the case of type-2 diabetes mellitus remains a major challenge. Therefore, a deep understanding of the structure-function relationship of these receptors will have great implications for future therapies in order to maintain a normal glucose level for an extended period of time. The antagonists of glucagon receptors that can effectively block the hepatic glucose production, as a result of glucagon action, are highly desirable for the tuning of the hyperglycemic state in type 2 diabetes mellitus. In the same manner, GLP-1R agonists act as important treatment modalities, thanks to their multiple anti-diabetic actions to attain normal glucose levels. In this review article, the structural diversity of glucagon and GLP-1 receptors along with their signaling pathways, site-directed mutations and significance in drug discovery against type-2 diabetes are illustrated. Moreover, the promising non-peptide antagonists of glucagon receptor and agonists of GLP-1 receptor, for the management of diabetes are presented with elaboration on the structure-activity relationship (SAR).

The putative binding site and SAR rationalization of small molecules against glucagon-like peptide-1 receptor using homology model and crystal structures: a comparative study

Glucagon-like peptide-1 (GLP-1) is involved in glucose-stimulated insulin secretion and weight regulating actions through the activation of the GLP-1 receptor (GLP-1R). Clinical effectiveness of GLP-1 mimetics is effective in improving glucose control in patients. Thus, identifying and developing orally active small-molecule agonists are highly desirable. This study summarizes the structure-function relationship of hGLP-1R through computational approaches and search of small molecule GLP-1R agonists. We carried out mutation guided data-driven study, for developing the GLP-1R model to explore and validate the putative site for quinoxaline analogues. The developed GLP-1R homology model was subjected to 500 ns MD simulation for validation. Various snapshots were considered to identify the best structure of GLP-1R based on correlation between experimental pEC₅₀ and various theoretical parameters (docking score, MM-GBSA ΔG bind, WM/MM ΔG bind). The putative binding site (Sitemap and WaterMap has been predicted and it matched well with the available data. Excellent correlation (R² =0.94), between pEC50 and WM/MM ΔG bind for the snapshot at 350 ns was observed after including induced-fit docking results of the most potent molecule. Enrichment calculation indicates better AUC (=0.75) for predicted complex structure. A comparison of the developed GLP-1R model with the available crystal structure shows excellent similarities and it was used for virtual screening to find small molecule agonists. The good correlation of our model with crystal structures of GLP-1R may help to understand the structure-function relationship of other secretin families.

The rapid "teabag" method for high-end purification of membrane proteins

Overproduction and purification of membrane proteins are generally challenging and time-consuming procedures due to low expression levels, misfolding, and low stability once extracted from the membrane. Reducing processing steps and shortening the timespan for purification represent attractive approaches to overcome some of these challenges. We have therefore compared a fast “teabag” purification method with conventional purification for five different membrane proteins (MraY, AQP10, ClC-1, PAR2 and KCC2). Notably, this new approach reduces the purification time significantly, and the quality of the purified membrane proteins is equal to or exceeds conventional methods as assessed by size exclusion chromatography, SDS-PAGE and downstream applications such as ITC, crystallization and cryo-EM. Furthermore, the method is scalable, applicable to a range of affinity resins and allows for parallelization. Consequently, the technique has the potential to substantially simplify purification efforts of membrane proteins in basic and applied sciences.

A combined activation mechanism for the glucagon receptor

We report on a combined activation mechanism for a class B G-protein-coupled receptor (GPCR), the glucagon receptor. By computing the conformational free-energy landscape associated with the activation of the receptor-agonist complex and comparing it with that obtained with the ternary complex (receptor-agonist-G protein) we show that the agonist stabilizes the receptor in a preactivated complex, which is then fully activated upon binding of the G protein. The proposed mechanism contrasts with the generally assumed GPCR activation mechanism, which proceeds through an opening of the intracellular region allosterically elicited by the binding of the agonist. The mechanism found here is consistent with electron cryo-microscopy structural data and might be general for class B GPCRs. It also helps us to understand the mode of action of the numerous allosteric antagonists of this important drug target.

Tuning Down the Pain - An Overview of Allosteric Modulation of Opioid Receptors: Mechanisms of Modulation, Allosteric Sites, Modulator Syntheses

Opioid signaling plays a central role in pain perception. As such, it remains the main target in the development of antinociceptive agents, despite serious side effects involved. In recent years, hopes for improved opioid painkillers are rising, together with our understanding of allosterism and biased signaling mechanisms. In this review, we focus on recently discovered allosteric modulators of opioid receptors, insights into phenomena underlying their action, as well as on how they extend our understanding of mechanisms of previously known compounds. A brief overlook of their synthesis is also presented.