Structural basis for the cooperative allosteric activation of the free fatty acid receptor GPR40

Influence of the Water Model on the Structure and Interactions of the GPR40 Protein with the Lipid Membrane and the Solvent: Rigid versus Flexible Water Models

G protein-coupled receptors (GPCR) are responsible for modulating various physiological functions and are thus related to the pathophysiology of different diseases. Being potential therapeutic targets, multiple computational methodologies have been developed to analyze their behavior and interactions with other species. The solvent, on the other hand, has received much less attention. In this work, we analyzed the effect of four explicit water models on the structure and interactions of the GPR40 receptor in its apo form. We employed the rigid SPC/E and TIP4P models, and their flexible versions, the FBA/ϵ and TIP4P/ϵflex. We explored the structural changes and their correlation with some bulk dynamic properties of water. Our results showed an adverse effect on the conservation of the secondary structure of the receptor with all the models due to the breaking of the intramolecular hydrogen bond network, being more evident for the TIP4P models. Notably, all four models brought the receptor to states similar to the active one, modifying the intracellular part of the TM5 and TM6 domains in a "hinge" type movement, allowing the opening of the structure. Regarding the dynamic properties, the rigid models showed results comparable to those obtained in other studies on membrane systems. However, flexible models exhibit disparities in the molecular representation of systems. Surprisingly, the FBA/ϵ model improves the molecular picture of several properties, even though their agreement with bulk diffusion is poorer. These findings reinforce our idea that exploring other water models or improving the current ones, to better represent the membrane interface, can lead to a positive impact on the description of the signal transduction mechanisms and the search of new drugs by targeting these receptors.

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

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

Structure-based identification of a G protein-biased allosteric modulator of cannabinoid receptor CB1

Cannabis sativa is known for its therapeutic benefit in various diseases including pain relief by targeting cannabinoid receptors. The primary component of cannabis, Δ9-tetrahydrocannabinol (THC), and other agonists engage the orthosteric site of CB1, activating both Gi and β-arrestin signaling pathways. The activation of diverse pathways could result in on-target side effects and cannabis addiction, which may hinder therapeutic potential. A significant challenge in pharmacology is the design of a ligand that can modulate specific signaling of CB1. By leveraging insights from the structure-function selectivity relationship (SFSR), we have identified Gi signaling-biased agonist-allosteric modulators (ago-BAMs). Further, two cryoelectron microscopy (cryo-EM) structures reveal the binding mode of ago-BAM at the extrahelical allosteric site of CB1. Combining mutagenesis and pharmacological studies, we elucidated the detailed mechanism of ago-BAM-mediated biased signaling. Notably, ago-BAM CB-05 demonstrated analgesic efficacy with fewer side effects, minimal drug toxicity and no cannabis addiction in mouse pain models. In summary, our finding not only suggests that ago-BAMs of CB1 provide a potential nonopioid strategy for pain management but also sheds light on BAM identification for GPCRs.

Cryo-electron microscopy for GPCR research and drug discovery in endocrinology and metabolism

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors, with many GPCRs having crucial roles in endocrinology and metabolism. Cryogenic electron microscopy (cryo-EM) has revolutionized the field of structural biology, particularly regarding GPCRs, over the past decade. Since the first pair of GPCR structures resolved by cryo-EM were published in 2017, the number of GPCR structures resolved by cryo-EM has surpassed the number resolved by X-ray crystallography by 30%, reaching >650, and the number has doubled every ~0.63 years for the past 6 years. At this pace, it is predicted that the structure of 90% of all human GPCRs will be completed within the next 5-7 years. This Review highlights the general structural features and principles that guide GPCR ligand recognition, receptor activation, G protein coupling, arrestin recruitment and regulation by GPCR kinases. The Review also highlights the diversity of GPCR allosteric binding sites and how allosteric ligands could dictate biased signalling that is selective for a G protein pathway or an arrestin pathway. Finally, the authors use the examples of glycoprotein hormone receptors and glucagon-like peptide 1 receptor to illustrate the effect of cryo-EM on understanding GPCR biology in endocrinology and metabolism, as well as on GPCR-related endocrine diseases and drug discovery.

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.

Structural and dynamic insights into the activation of the μ-opioid receptor by an allosteric modulator

G-protein-coupled receptors (GPCRs) play pivotal roles in various physiological processes. These receptors are activated to different extents by diverse orthosteric ligands and allosteric modulators. However, the mechanisms underlying these variations in signaling activity by allosteric modulators remain largely elusive. Here, we determine the three-dimensional structure of the μ-opioid receptor (MOR), a class A GPCR, in complex with the Gi protein and an allosteric modulator, BMS-986122, using cryogenic electron microscopy. Our results reveal that BMS-986122 binding induces changes in the map densities corresponding to R1673.50 and Y2545.58, key residues in the structural motifs conserved among class A GPCRs. Nuclear magnetic resonance analyses of MOR in the absence of the Gi protein reveal that BMS-986122 binding enhances the formation of the interaction between R1673.50 and Y2545.58, thus stabilizing the fully-activated conformation, where the intracellular half of TM6 is outward-shifted to allow for interaction with the Gi protein. These findings illuminate that allosteric modulators like BMS-986122 can potentiate receptor activation through alterations in the conformational dynamics in the core region of GPCRs. Together, our results demonstrate the regulatory mechanisms of GPCRs, providing insights into the rational development of therapeutics targeting GPCRs.

Computational insights into diverse binding modes of the allosteric modulator and their regulation on dopamine D1 receptor

Allosteric drugs hold the promise of addressing many challenges in the current drug development of GPCRs. However, the molecular mechanism underlying their allosteric modulations remain largely elusive. The dopamine D1 receptor (DRD1), a member of Class A GPCRs, is critical for treating psychiatric disorders, and LY3154207 serves as its promising positive allosteric modulator (PAM). In the work, we utilized extensive Gaussian-accelerated molecular dynamics simulations (a total of 41μs) for the first time probe the diverse binding modes of the allosteric modulator and their regulation effects, based on the DRD1 and LY3154207 as representative. Our simulations identify four binding modes of LY3154207 (one boat mode, two metastable vertical modes and a novel cleft-anchored mode), in which the boat mode is the most stable while there three modes are similar in the stability. However, it is interesting to observed that the most stable boat mode inversely exhibits the weakest positive allosteric effect on influencing the orthosteric ligand binding and maintaining the activity of the transducer binding site. It should result from its induced weaker correlation between the allosteric site and the orthosteric site, and between the orthosteric site and the transducer binding site than the other three binding modes, as well as its weakened interaction between a crucial activation-related residue (S2025.46) and the orthosteric ligand (dopamine). Overall, the work offers atomic-level information to advance our understanding of the complex allosteric regulation on GPCRs, which is beneficial to the allosteric modulator design and development.

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.

Structure of GPR101-Gs enables identification of ligands with rejuvenating potential

GPR101 is an orphan G protein-coupled receptor actively participating in energy homeostasis. Here we report the cryo-electron microscopy structure of GPR101 constitutively coupled to Gs heterotrimer, which reveals unique features of GPR101, including the interaction of extracellular loop 2 within the 7TM bundle, a hydrophobic chain packing-mediated activation mechanism and the structural basis of disease-related mutants. Importantly, a side pocket is identified in GPR101 that facilitates in silico screening to identify four small-molecule agonists, including AA-14. The structure of AA-14-GPR101-Gs provides direct evidence of the AA-14 binding at the side pocket. Functionally, AA-14 partially restores the functions of GH/IGF-1 axis and exhibits several rejuvenating effects in wild-type mice, which are abrogated in Gpr101-deficient mice. In summary, we provide a structural basis for the constitutive activity of GPR101. The structure-facilitated identification of GPR101 agonists and functional analysis suggest that targeting this orphan receptor has rejuvenating potential.

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

Free fatty acid receptors 1 to 4 (FFA1 to FFA4) are class A G protein-coupled receptors (GPCRs). FFA1 to FFA3 share substantial sequence similarity, whereas FFA4 is unrelated. However, FFA1 and FFA4 are activated by long-chain fatty acids, while FFA2 and FFA3 respond to short-chain fatty acids generated by intestinal microbiota. FFA1, FFA2, and FFA4 are potential drug targets for metabolic and inflammatory conditions. Here, we determined the active structures of FFA1 and FFA4 bound to docosahexaenoic acid, FFA4 bound to the synthetic agonist TUG-891, and butyrate-bound FFA2, each complexed with an engineered heterotrimeric Gq protein (miniGq), by cryo-electron microscopy. Together with computational simulations and mutagenesis studies, we elucidated the similarities and differences in the binding modes of fatty acid ligands to their respective GPCRs. Our findings unveiled distinct mechanisms of receptor activation and G protein coupling. We anticipate that these outcomes will facilitate structure-based drug development and underpin future research on this group of GPCRs.

Structural basis of G protein-Coupled receptor CMKLR1 activation and signaling induced by a chemerin-derived agonist

Chemokine-like receptor 1 (CMKLR1), also known as chemerin receptor 23 (ChemR23) or chemerin receptor 1, is a chemoattractant G protein-coupled receptor (GPCR) that responds to the adipokine chemerin and is highly expressed in innate immune cells, including macrophages and neutrophils. The signaling pathways of CMKLR1 can lead to both pro- and anti-inflammatory effects depending on the ligands and physiological contexts. To understand the molecular mechanisms of CMKLR1 signaling, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of the CMKLR1-Gi signaling complex with chemerin9, a nanopeptide agonist derived from chemerin, which induced complex phenotypic changes of macrophages in our assays. The cryo-EM structure, together with molecular dynamics simulations and mutagenesis studies, revealed the molecular basis of CMKLR1 signaling by elucidating the interactions at the ligand-binding pocket and the agonist-induced conformational changes. Our results are expected to facilitate the development of small molecule CMKLR1 agonists that mimic the action of chemerin9 to promote the resolution of inflammation.

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.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Predicting GPR40 Agonists with A Deep Learning-Based Ensemble Model

Recent studies have identified G protein-coupled receptor 40 (GPR40) as a promising target for treating type 2 diabetes mellitus, and GPR40 agonists have several superior effects over other hypoglycemic drugs, including cardiovascular protection and suppression of glucagon levels. In this study, we constructed an up-to-date GPR40 ligand dataset for training models and performed a systematic optimization of the ensemble model, resulting in a powerful ensemble model (ROC AUC: 0.9496) for distinguishing GPR40 agonists and non-agonists. The ensemble model is divided into three layers, and the optimization process is carried out in each layer. We believe that these results will prove helpful for both the development of GPR40 agonists and ensemble models. All the data and models are available on GitHub. (https://github.com/Jiamin-Yang/ensemble_model).

An inverse agonist of orphan receptor GPR61 acts by a G protein-competitive allosteric mechanism

GPR61 is an orphan GPCR related to biogenic amine receptors. Its association with phenotypes relating to appetite makes it of interest as a druggable target to treat disorders of metabolism and body weight, such as obesity and cachexia. To date, the lack of structural information or a known biological ligand or tool compound has hindered comprehensive efforts to study GPR61 structure and function. Here, we report a structural characterization of GPR61, in both its active-like complex with heterotrimeric G protein and in its inactive state. Moreover, we report the discovery of a potent and selective small-molecule inverse agonist against GPR61 and structural elucidation of its allosteric binding site and mode of action. These findings offer mechanistic insights into an orphan GPCR while providing both a structural framework and tool compound to support further studies of GPR61 function and modulation.

Biased allosteric activation of ketone body receptor HCAR2 suppresses inflammation

Hydroxycarboxylic acid receptor 2 (HCAR2), modulated by endogenous ketone body β-hydroxybutyrate and exogenous niacin, is a promising therapeutic target for inflammation-related diseases. HCAR2 mediates distinct pathophysiological events by activating Gi/o protein or β-arrestin effectors. Here, we characterize compound 9n as a Gi-biased allosteric modulator (BAM) of HCAR2 and exhibit anti-inflammatory efficacy in RAW264.7 macrophages via a specific HCAR2-Gi pathway. Furthermore, four structures of HCAR2-Gi complex bound to orthosteric agonists (niacin or monomethyl fumarate), compound 9n, and niacin together with compound 9n simultaneously reveal a common orthosteric site and a unique allosteric site. Combined with functional studies, we decipher the action framework of biased allosteric modulation of compound 9n on the orthosteric site. Moreover, co-administration of compound 9n with orthosteric agonists could enhance anti-inflammatory effects in the mouse model of colitis. Together, our study provides insight to understand the molecular pharmacology of the BAM and facilitates exploring the therapeutic potential of the BAM with orthosteric drugs.

Dynamic allosteric networks drive adenosine A1 receptor activation and G-protein coupling

G-protein coupled receptors (GPCRs) present specific activation pathways and signaling among receptor subtypes. Hence, an extensive knowledge of the structural dynamics of the receptor is critical for the development of therapeutics. Here, we target the adenosine A1 receptor (A1R), for which a negligible number of drugs have been approved. We combine molecular dynamics simulations, enhanced sampling techniques, network theory and pocket detection to decipher the activation pathway of A1R, decode the allosteric networks and identify transient pockets. The A1R activation pathway reveal hidden intermediate and pre-active states together with the inactive and fully-active states observed experimentally. The protein energy networks computed throughout these conformational states successfully unravel the extra and intracellular allosteric centers and the communication pathways that couples them. We observe that the allosteric networks are dynamic, being increased along activation and fine-tuned in presence of the trimeric G-proteins. Overlap of transient pockets and energy networks uncover how the allosteric coupling between pockets and distinct functional regions of the receptor is altered along activation. By an in-depth analysis of the bridge between activation pathway, energy networks and transient pockets, we provide a further understanding of A1R. This information can be useful to ease the design of allosteric modulators for A1R.

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

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

A non-canonical mechanism of GPCR activation

The goal of designing safer, more effective drugs has led to tremendous interest in molecular mechanisms through which ligands can precisely manipulate signaling of G-protein-coupled receptors (GPCRs), the largest class of drug targets. Decades of research have led to the widely accepted view that all agonists-ligands that trigger GPCR activation-function by causing rearrangement of the GPCR's transmembrane helices, opening an intracellular pocket for binding of transducer proteins. Here we demonstrate that certain agonists instead trigger activation of free fatty acid receptor 1 by directly rearranging an intracellular loop that interacts with transducers. We validate the predictions of our atomic-level simulations by targeted mutagenesis; specific mutations which disrupt interactions with the intracellular loop convert these agonists into inverse agonists. Further analysis suggests that allosteric ligands could regulate signaling of many other GPCRs via a similar mechanism, offering rich possibilities for precise control of pharmaceutically important targets.

Targeting the endocannabinoid system: Structural determinants and molecular mechanism of allosteric modulation

Although drugs targeting the orthosteric binding site of cannabinoid receptors (CBRs) have several therapeutic effects on human physiological and pathological conditions, they can also cause serious adverse effects. Only a few orthosteric ligands have successfully passed clinical trials. Recently, allosteric modulation has become a novel option for drug discovery, with fewer adverse effects and the potential to avoid drug overdose. In this review, we highlight novel findings related to the drug discovery of allosteric modulators (AMs) targeting CBRs. We summarize newly synthesized AMs and the reported/predicted allosteric binding sites. We also discuss the structural determinants of the AMs binding as well as the molecular mechanism of CBR allostery.

A Free Fatty Acid Synthetic Agonist Accelerates Wound Healing and Improves Scar Quality in Mice

BACKGROUND

Impaired wound healing is a health problem around the world, and the search for a novel product to repair wounded skin is a major topic in the field. GW9508 is a synthetic molecule described as a selective agonist of free fatty acid receptors (FFARs) 1 and 4, and there is evidence of its anti-inflammatory effects on several organs of the body.

PURPOSE

Here, we aimed to evaluate the effects of topical GW9508 on wound healing in mice.

RESEARCH DESIGN

First, we used bioinformatic methods to determine the expression of FFAR1 and FFAR4 mRNA in the skin from a human cell atlas assembled with single-cell transcriptomes. Next, we employed 6-week-old C57BL6J mice with 2 wounds inflicted in the back. The mice were randomly divided into 2 groups, a control group, which received topical vehicle, and a treatment group, which received GW9508, for 12 days. The wound was monitored by photographic documentation every 2 days, and samples were collected at day 6 and 12 post injury for RT-PCR, western blot and histology analyses.

RESULTS

FFAR1 and FFAR4 mRNA are expressed in skin cells in similar amounts to those in other tissues. Topical GW9508 accelerated wound healing and decreased gene expression of IL-10 and metalloproteinase 9 on days 6 and 12 post injury. It increased the quantity of Collagen I and improved the organization of collagen fibres. Conclusions: Our results show that GW9508 could be an attractive drug treatment for wounded skin. Future studies need to be performed to assess the impact of GW9508 in chronic wound models.

Trans-palmitoleic acid, a dairy fat biomarker, stimulates insulin secretion and activates G protein-coupled receptors with a different mechanism from the cis isomer

Dietary trans-palmitoleic acid (trans 16:1n-7, tPOA), a biomarker for high-fat dairy product intake, has been associated with a lower risk of type 2 diabetes mellitus (T2DM) in some cross-sectional and prospective epidemiological studies. Here, we investigated the insulin secretion-promoting activity of tPOA and compared them with the effects evoked by the cis-POA isomer (cPOA), an endogenous lipokine biosynthesized in the liver and adipose tissue, and found in some natural food sources. The debate about the positive and negative relationships of those two POA isomers with metabolic risk factors and the underlying mechanisms is still going on. Therefore, we examined the potency of both POA isomers to potentiate insulin secretion in murine and human pancreatic β cell lines. We also investigated whether POA isomers activate G protein-coupled receptors proposed as potential targets for T2DM treatment. We show that tPOA and cPOA augment glucose-stimulated insulin secretion (GSIS) to a similar extent; however, their insulin secretagogue activity is associated with different signaling pathways. We also performed ligand docking and molecular dynamics simulations to predict the preferred orientation of POA isomers and the strength of association between those two fatty acids and GPR40, GPR55, GPR119, and GPR120 receptors. Overall, this study provides insight into the bioactivity of tPOA and cPOA toward selected GPCR functions, indicating them as targets responsible for the insulin secretagogue action of POA isomers. It reveals that both tPOA and cPOA may promote insulin secretion and subsequently regulate glucose homeostasis.

Structural basis of CMKLR1 signaling induced by chemerin9

Chemokine-like receptor 1 (CMKLR1), also known as chemerin receptor 23 (ChemR23) or chemerin receptor 1, is a chemoattractant G protein-coupled receptor (GPCR) that responds to the adipokine chemerin and is highly expressed in innate immune cells, including macrophages and neutrophils. The signaling pathways of CMKLR1 can lead to both pro- and anti-inflammatory effects depending on the ligands and physiological contexts. To understand the molecular mechanisms of CMKLR1 signaling, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of the CMKLR1-Gi signaling complex with chemerin9, a nanopeptide agonist derived from chemerin, which induced complex phenotypic changes of macrophages in our assays. The cryo-EM structure, together with molecular dynamics simulations and mutagenesis studies, revealed the molecular basis of CMKLR1 signaling by elucidating the interactions at the ligand-binding pocket and the agonist-induced conformational changes. Our results are expected to facilitate the development of small molecule CMKLR1 agonists that mimic the action of chemerin9 to promote the resolution of inflammation.

Molecular mechanism of fatty acid activation of FFAR1

FFAR1 is a G-protein-coupled receptor (GPCR) that responds to circulating free fatty acids to enhance glucose-stimulated insulin secretion and release of incretin hormones. Due to the glucose-lowering effect of FFAR1 activation, potent agonists for this receptor have been developed for the treatment of diabetes. Previous structural and biochemical studies of FFAR1 showed multiple sites of ligand binding to the inactive state but left the mechanism of fatty acid interaction and receptor activation unknown. We used cryo-electron microscopy to elucidate structures of activated FFAR1 bound to a Gq mimetic, which were induced either by the endogenous FFA ligand docosahexaenoic acid or γ-linolenic acid and the agonist drug TAK-875. Our data identify the orthosteric pocket for fatty acids and show how both endogenous hormones and synthetic agonists induce changes in helical packing along the outside of the receptor that propagate to exposure of the G-protein-coupling site. These structures show how FFAR1 functions without the highly conserved "DRY" and "NPXXY" motifs of class A GPCRs and also illustrate how the orthosteric site of a receptor can be bypassed by membrane-embedded drugs to confer full activation of G protein signaling.

Unsaturated bond recognition leads to biased signal in a fatty acid receptor

Individual free fatty acids (FAs) play important roles in metabolic homeostasis, many through engagement with more than 40G protein-coupled receptors. Searching for receptors to sense beneficial omega-3 FAs of fish oil enabled the identification of GPR120, which is involved in a spectrum of metabolic diseases. Here, we report six cryo-electron microscopy structures of GPR120 in complex with FA hormones or TUG891 and Gi or Giq trimers. Aromatic residues inside the GPR120 ligand pocket were responsible for recognizing different double-bond positions of these FAs and connect ligand recognition to distinct effector coupling. We also investigated synthetic ligand selectivity and the structural basis of missense single-nucleotide polymorphisms. We reveal how GPR120 differentiates rigid double bonds and flexible single bonds. The knowledge gleaned here may facilitate rational drug design targeting to GPR120.

Neurotensin Receptor Allosterism Revealed in Complex with a Biased Allosteric Modulator

The NTSR1 neurotensin receptor (NTSR1) is a G protein-coupled receptor (GPCR) found in the brain and peripheral tissues with neurotensin (NTS) being its endogenous peptide ligand. In the brain, NTS modulates dopamine neuronal activity, induces opioid-independent analgesia, and regulates food intake. Recent studies indicate that biasing NTSR1 toward β-arrestin signaling can attenuate the actions of psychostimulants and other drugs of abuse. Here, we provide the cryoEM structures of NTSR1 ternary complexes with heterotrimeric Gq and GoA with and without the brain-penetrant small-molecule SBI-553. In functional studies, we discovered that SBI-553 displays complex allosteric actions exemplified by negative allosteric modulation for G proteins that are Gα subunit selective and positive allosteric modulation and agonism for β-arrestin translocation at NTSR1. Detailed structural analysis of the allosteric binding site illuminated the structural determinants for biased allosteric modulation of SBI-553 on NTSR1.

Ligand recognition and allosteric modulation of the human MRGPRX1 receptor

The human MAS-related G protein-coupled receptor X1 (MRGPRX1) is preferentially expressed in the small-diameter primary sensory neurons and involved in the mediation of nociception and pruritus. Central activation of MRGPRX1 by the endogenous opioid peptide fragment BAM8-22 and its positive allosteric modulator ML382 has been shown to effectively inhibit persistent pain, making MRGPRX1 a promising target for non-opioid pain treatment. However, the activation mechanism of MRGPRX1 is still largely unknown. Here we report three high-resolution cryogenic electron microscopy structures of MRGPRX1-Gαq in complex with BAM8-22 alone, with BAM8-22 and ML382 simultaneously as well as with a synthetic agonist compound-16. These structures reveal the agonist binding mode for MRGPRX1 and illuminate the structural requirements for positive allosteric modulation. Collectively, our findings provide a molecular understanding of the activation and allosteric modulation of the MRGPRX1 receptor, which could facilitate the structure-based design of non-opioid pain-relieving drugs.

Optimization of Physicochemical Properties of Pyrrolidine GPR40 AgoPAMs Results in a Differentiated Profile with Improved Pharmacokinetics and Reduced Off-Target Activities

GPR40 AgoPAMs are highly effective antidiabetic agents that have a dual mechanism of action, stimulating both glucose-dependent insulin and GLP-1 secretion. The early lipophilic, aromatic pyrrolidine and dihydropyrazole GPR40 AgoPAMs from our laboratory were highly efficacious in lowering plasma glucose levels in rodents but possessed off-target activities and triggered rebound hyperglycemia in rats at high doses. A focus on increasing molecular complexity through saturation and chirality in combination with reducing polarity for the pyrrolidine AgoPAM chemotype resulted in the discovery of compound 46, which shows significantly reduced off-target activities as well as improved aqueous solubility, rapid absorption, and linear PK. In vivo, compound 46 significantly lowers plasma glucose levels in rats during an oral glucose challenge yet does not demonstrate the reactive hyperglycemia effect at high doses that was observed with earlier GPR40 AgoPAMs.

The impact of cryo-EM on determining allosteric modulator-bound structures of G protein-coupled receptors

G-protein coupled receptors (GPCRs) are important therapeutic targets for the treatment of human disease. Although GPCRs are highly successful drug targets, there are many challenges associated with the discovery and translation of small molecule ligands that target the endogenous ligand-binding site for GPCRs. Allosteric modulators are a class of ligands that target alternative binding sites known as allosteric sites and offer fresh opportunities for the development of new therapeutics. However, only a few allosteric modulators have been approved as drugs. Advances in GPCR structural biology enabled by the cryogenic electron microscopy (cryo-EM) revolution have provided new insights into the molecular mechanism and binding location of small molecule allosteric modulators. This review highlights the latest findings from allosteric modulator-bound structures of Class A, B, and C GPCRs with a focus on small molecule ligands. Emerging methods that will facilitate cryo-EM structures of more difficult ligand-bound GPCR complexes are also discussed. The results of these studies are anticipated to aid future structure-based drug discovery efforts across many different GPCRs.

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.

Comparative Study of Receptor-, Receptor State-, and Membrane-Dependent Cholesterol Binding Sites in A2A and A1 Adenosine Receptors Using Coarse-Grained Molecular Dynamics Simulations

We used coarse-grained molecular dynamics (CG MD) simulations to study protein-cholesterol interactions for different activation states of the A_2A adenosine receptor (A_2AR) and the A₁ adenosine receptor (A₁R) and predict new cholesterol binding sites indicating amino acid residues with a high residence time in three biologically relevant membranes. Compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-cholesterol and POPC-phosphatidylinositol-bisphosphate (PIP₂)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A_2AR and revealed the binding sites with long-lasting amino acid residues. We observed that using the plasma mimetic membrane and plotting residues with cholesterol residence time ≥2 μs, our CG MD simulations captured most obviously the cholesterol-protein interactions. For the inactive A_2AR, we identified one more binding site in which cholesterol is bound to residues with a long residence time compared to the previously detected, for the active A₁R, three binding sites, and for the inactive A₁R, two binding sites. We calculated that for the active states, cholesterol binds to residues with a much longer residence time compared to the inactive state for both A_2AR and A₁R. The stability of the identified binding sites to A₁R or A_2AR with CG MD simulations was additionally investigated with potential of mean force calculations using umbrella sampling. We observed that the binding sites with residues to which cholesterol has a long residence time in A_2AR have shallow binding free energy minima compared to the related binding sites in A₁R, suggesting a stronger binding for cholesterol to A₁R. The differences in binding sites in which cholesterol is stabilized and interacts with residues with a long residence time between active and inactive states of A₁R and A_2AR can be important for differences in functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators, which can bind through lipid pathways. We observed a stronger binding for cholesterol to A₁R (i.e., generally higher association rates) compared to A_2AR, which remains to be demonstrated. For the active states, cholesterol binds to residues with much longer residence times compared to the inactive state for both A_2AR and A₁R. Taken together, binding sites of active A₁R may be considered as promising allosteric targets.

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.

Learn from failures and stay hopeful to GPR40, a GPCR target with robust efficacy, for therapy of metabolic disorders

GPR40 is a class A G-protein coupled receptor (GPCR) mainly expressed in pancreas, intestine, and brain. Its endogenous ligand is long-chain fatty acids, which activate GPR40 after meal ingestion to induce secretion of incretins in the gut, including GLP-1, GIP, and PYY, the latter control appetite and glucose metabolism. For its involvement in satiety regulation and metabolic homeostasis, partial and AgoPAM (Positive Allosteric Modulation agonist) GPR40 agonists had been developed for type 2 diabetes (T2D) by many pharmaceutical companies. The proof-of-concept of GPR40 for control of hyperglycemia was achieved by clinical trials of partial GPR40 agonist, TAK-875, demonstrating a robust decrease in HbA1c (-1.12%) after chronic treatment in T2D. The development of TAK-875, however, was terminated due to liver toxicity in 2.7% patients with more than 3-fold increase of ALT in phase II and III clinical trials. Different mechanisms had since been proposed to explain the drug-induced liver injury, including acyl glucuronidation, inhibition of mitochondrial respiration and hepatobiliary transporters, ROS generation, etc. In addition, activation of GPR40 by AgoPAM agonists in pancreas was also linked to β-cell damage in rats. Notwithstanding the multiple safety concerns on the development of small-molecule GPR40 agonists for T2D, some partial and AgoPAM GPR40 agonists are still under clinical development. Here we review the most recent progress of GPR40 agonists development and the possible mechanisms of the side effects in different organs, and discuss the possibility of developing novel strategies that retain the robust efficacy of GPR40 agonists for metabolic disorders while avoid toxicities caused by off-target and on-target mechanisms.

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.

Structural identification of lysophosphatidylcholines as activating ligands for orphan receptor GPR119

Lysophosphatidylcholine (LPC) is an essential mediator in human lipid metabolism and is associated with a variety of diseases, but the exact identity of LPC receptors remains controversial. Through extensive biochemical and structural analyses, we have identified the orphan receptor GPR119 as the receptor for LPC. The structure of the GPR119-G-protein complex without any added ligands reveals a density map that fits well with LPC, which is further confirmed by mass spectrometry and functional studies. As LPCs are abundant on the cell membrane, their preoccupancy in the receptor may lead to 'constitutive activity' of GPR119. The structure of GPR119 bound to APD668, a clinical drug candidate for type 2 diabetes, reveals an exceedingly similar binding mode to LPC. Together, these data highlight structural evidence for LPC function in regulating glucose-dependent insulin secretion through direct binding and activation of GPR119, and provide structural templates for drug design targeting GPR119.

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

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

Mechanistic Studies on the Stereoselectivity of FFAR1 Modulators

Free fatty acid receptor 1 (FFAR1) is a potential therapeutic target for the treatment of type 2 diabetes (T2D). It has been validated that agonists targeting FFAR1 can achieve the initial therapeutic endpoints of T2D, and the epimer agonists (R,S) AM-8596 can activate FFAR1 differently, with one acting as a partial agonist and the other as a full agonist. Up to now, the origin of the stereoselectivity of FFAR1 agonists remains elusive. In this work, we used molecular simulation methods to elucidate the mechanism of the stereoselectivity of the FFAR1 agonists (R)-AM-8596 and (S)-AM-8596. We found that the full agonist (R)-AM-8596 disrupts the residue interaction network around the receptor binding pocket and promotes the opening of the binding site for the G-protein, thereby resulting in the full activation of FFAR1. In contrast, the partial agonist (S)-AM-8596 forms stable electrostatic interactions with FFAR1, which stabilizes the residue network and hinders the conformational transition of the receptor. Our work thus clarifies the selectivity and underlying molecular activation mechanism of FFAR1 agonists.

The pocketome of G-protein-coupled receptors reveals previously untargeted allosteric sites

G-protein-coupled receptors do not only feature the orthosteric pockets, where most endogenous agonists bind, but also a multitude of other allosteric pockets that have come into the focus as potential binding sites for synthetic modulators. Here, to better characterise such pockets, we investigate 557 GPCR structures by exhaustively docking small molecular probes in silico and converting the ensemble of binding locations to pocket-defining volumes. Our analysis confirms all previously identified pockets and reveals nine previously untargeted sites. In order to test for the feasibility of functional modulation of receptors through binding of a ligand to such sites, we mutate residues in two sites, in two model receptors, the muscarinic acetylcholine receptor M₃ and β₂-adrenergic receptor. Moreover, we analyse the correlation of inter-residue contacts with the activation states of receptors and show that contact patterns closely correlating with activation indeed coincide with these sites.

Activation mechanism of the μ-opioid receptor by an allosteric modulator

SignificanceThe allosteric modulators, which bind to nonorthosteric sites to enhance the signaling activities of G-protein-coupled receptors (GPCRs), are new candidates for GPCR-targeting drugs. Our solution NMR analyses of the μ-opioid receptor (MOR) revealed that the MOR activity was determined by a conformational equilibrium between three conformations. Interestingly, an allosteric modulator shifted the equilibrium toward a conformation with the highest activity to a level that cannot be reached by orthosteric ligands alone, leading to the increased activity of MOR. Our NMR analyses also identified the binding site of the allosteric modulator, including the residues contributing to the regulation of the equilibrium. These findings provide insights into the rational developments of novel allosteric modulators.

GLOW: A Workflow Integrating Gaussian-Accelerated Molecular Dynamics and Deep Learning for Free Energy Profiling

We introduce a Gaussian-accelerated molecular dynamics (GaMD), deep learning (DL), and free energy profiling workflow (GLOW) to predict molecular determinants and map free energy landscapes of biomolecules. All-atom GaMD-enhanced sampling simulations are first performed on biomolecules of interest. Structural contact maps are then calculated from GaMD simulation frames and transformed into images for building DL models using a convolutional neural network. Important structural contacts are further determined from DL models of attention maps of the structural contact gradients, which allow us to identify the system reaction coordinates. Finally, free energy profiles are calculated for the selected reaction coordinates through energetic reweighting of the GaMD simulations. We have also successfully demonstrated GLOW for the characterization of activation and allosteric modulation of a G protein-coupled receptor, using the adenosine A₁ receptor (A₁AR) as a model system. GLOW findings are highly consistent with previous experimental and computational studies of the A₁AR, while also providing further mechanistic insights into the receptor function. In summary, GLOW provides a systematic approach to mapping free energy landscapes of biomolecules. The GLOW workflow and its user manual can be downloaded at http://miaolab.org/GLOW.

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.

Mapping the binding sites of UDP and prostaglandin E2 glyceryl ester in the nucleotide receptor P2Y6

Cyclooxygenase‐2 catalyzes the biosynthesis of prostaglandins from arachidonic acid and the biosynthesis of prostaglandin glycerol esters (PG‐Gs) from 2‐arachidonoylglycerol. PG‐Gs are mediators of several biological actions such as macrophage activation, hyperalgesia, synaptic plasticity, and intraocular pressure. Recently, the human UDP receptor P2Y₆ was identified as a target for the prostaglandin E2 glycerol ester (PGE₂‐G). Here, we show that UDP and PGE₂‐G are evolutionary conserved endogenous agonists at vertebrate P2Y₆ orthologs. Using sequence comparison of P2Y₆ orthologs, homology modeling, and ligand docking studies, we proposed several receptor positions participating in agonist binding. Site‐directed mutagenesis and functional analysis of these P2Y₆ mutants revealed that both UDP and PGE₂‐G share in parts one ligand‐binding site. Thus, the convergent signaling of these two chemically very different agonists has already been manifested in the evolutionary design of the ligand‐binding pocket.

In Silico Prediction and Validation of CB2 Allosteric Binding Sites to Aid the Design of Allosteric Modulators

Although the 3D structures of active and inactive cannabinoid receptors type 2 (CB2) are available, neither the X-ray crystal nor the cryo-EM structure of CB2-orthosteric ligand-modulator has been resolved, prohibiting the drug discovery and development of CB2 allosteric modulators (AMs). In the present work, we mainly focused on investigating the potential allosteric binding site(s) of CB2. We applied different algorithms or tools to predict the potential allosteric binding sites of CB2 with the existing agonists. Seven potential allosteric sites can be observed for either CB2-CP55940 or CB2-WIN 55,212-2 complex, among which sites B, C, G and K are supported by the reported 3D structures of Class A GPCRs coupled with AMs. Applying our novel algorithm toolset-MCCS, we docked three known AMs of CB2 including Ec2la (C-2), trans-β-caryophyllene (TBC) and cannabidiol (CBD) to each site for further comparisons and quantified the potential binding residues in each allosteric binding site. Sequentially, we selected the most promising binding pose of C-2 in five allosteric sites to conduct the molecular dynamics (MD) simulations. Based on the results of docking studies and MD simulations, we suggest that site H is the most promising allosteric binding site. We plan to conduct bio-assay validations in the future.

A FFAR1 full agonist restores islet function in models of impaired glucose-stimulated insulin secretion and diabetic non-human primates

The free fatty acid receptor 1 (FFAR1/GPR40) mediates fatty acid-induced insulin secretion from pancreatic β-cells. At least 3 distinct binding sites exist on the FFAR1 receptor and numerous synthetic ligands have been investigated for their anti-diabetic actions. Fasiglifam, binds to site-1 and stimulates intra-cellular calcium release and improves glycemic control in diabetic patients. Recently, small molecule FFAR1 agonists were discovered which bind to site-3, stimulating both intra-cellular calcium and cAMP, resulting in insulin and glucagon-like peptide-1 (GLP-1) secretion. The ability of our site-3 FFAR1 agonist (compound A) to control blood glucose was evaluated in spontaneously diabetic cynomolgus monkeys during an oral glucose tolerance test. In type-2 diabetic (T2D) animals, significant reductions in blood glucose and insulin were noted. To better understand the mechanism of these in vivo findings, we evaluated the effect of compound A in islets under several conditions of dysfunction. First, healthy human and non-human primate islets were treated with compound A and showed potentiation of insulin and glucagon secretion from both species. Next, we determined glucose-responsive insulin secretion under gluco-lipotoxic conditions and from islets isolated from type-2 diabetic humans. Despite a dysfunctional phenotype that failed to secrete insulin in response to glucose, site-3 FFAR1 agonism not only enhanced insulin secretion, but restored glucose responsiveness across a range of glucose concentrations. Lastly, we treated ex vivo human islets chronically with a sulfonylurea to induce secondary beta-cell failure. Again, this model showed reduced glucose-responsive insulin secretion that was restored and potentiated by site-3 FFAR1 agonism. Together these data suggest a mechanism for FFAR1 where agonists have direct effects on islet hormone secretion that can overcome a dysfunctional T2D phenotype. These unique characteristics of FFAR1 site-3 agonists make them an appealing potential therapy to treat type-2 diabetes.