Allosteric Ligand Binding Site Research Articles

Molecular insights into disease-associated glutamate transporter (EAAT1 / SLC1A3) variants using in silico and in vitro approaches

Glutamate is an essential excitatory neurotransmitter and an intermediate for energy metabolism. Depending on the tumor site, cancer cells have increased or decreased expression of excitatory amino acid transporter 1 or 2 (EAAT1/2, SLC1A3/2) to regulate glutamate uptake for the benefit of tumor growth. Thus, EAAT1/2 may be an attractive target for therapeutic intervention in oncology. Genetic variation of EAAT1 has been associated with rare cases of episodic ataxia, but the occurrence and functional contribution of EAAT1 mutants in other diseases, such as cancer, is poorly understood. Here, 105 unique somatic EAAT1 mutations were identified in cancer patients from the Genomic Data Commons dataset. Using EAAT1 crystal structures and in silico studies, eight mutations were selected based on their close proximity to the orthosteric or allosteric ligand binding sites and the predicted change in ligand binding affinity. In vitro functional assessment in a live-cell, impedance-based phenotypic assay demonstrated that these mutants differentially affect L-glutamate and L-aspartate transport, as well as the inhibitory potency of an orthosteric (TFB-TBOA) and allosteric (UCPH-101) inhibitor. Moreover, two episodic ataxia-related mutants displayed functional responses that were in line with literature, which confirmed the validity of our assay. Of note, ataxia-related mutant M128R displayed inhibitor-induced functional responses never described before. Finally, molecular dynamics (MD) simulations were performed to gain mechanistic insights into the observed functional effects. Taken together, the results in this work demonstrate 1) the suitability of the label-free phenotypic method to assess functional variation of EAAT1 mutants and 2) the opportunity and challenges of using in silico techniques to rationalize the in vitro phenotype of disease-relevant mutants.

Conservation of Allosteric Ligand Binding Sites in G-Protein Coupled Receptors.

Despite the growing number of G protein-coupled receptor (GPCR) structures, only 39 structures have been cocrystallized with allosteric inhibitors. These structures have been studied by protein mapping using the FTMap server, which determines the clustering of small organic probe molecules distributed on the protein surface. The method has found druggable sites overlapping with the cocrystallized allosteric ligands in 21 GPCR structures. Mapping of Alphafold2 generated models of these proteins confirms that the same sites can be identified without the presence of bound ligands. We then mapped the 394 GPCR X-ray structures available at the time of the analysis (September 2020). Results show that for each of the 21 structures with bound ligands there exist many other GPCRs that have a strong binding hot spot at the same location, suggesting potential allosteric sites in a large variety of GPCRs. These sites cluster at nine distinct locations, and each can be found in many different proteins. However, ligands binding at the same location generally show little or no similarity, and the amino acid residues interacting with these ligands also differ. Results confirm the possibility of specifically targeting these sites across GPCRs for allosteric modulation and help to identify the most likely binding sites among the limited number of potential locations. The FTMap server is available free of charge for academic and governmental use at https://ftmap.bu.edu/.

Validation of an Allosteric Binding Site of Src Kinase Identified by Unbiased Ligand Binding Simulations

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often “cryptic” without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a V-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein–ligand systems.

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

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

Covalent Occlusion of the RORγt Ligand Binding Pocket Allows Unambiguous Targeting of an Allosteric Site.

The nuclear receptor RORγt is a key positive regulator in the differentiation and proliferation of T helper 17 (Th17) cells and the production of proinflammatory cytokines like IL-17a. Dysregulation of this pathway can result in the development of various autoimmune diseases, and inhibition of RORγt with small molecules thus holds great potential as a therapeutic strategy. RORγt has a unique allosteric ligand binding site in the ligand binding domain, which is distinct from the canonical, orthosteric binding site. Allosteric modulation of RORγt shows high potential, but the targeted discovery of novel allosteric ligands is highly challenging via currently available methods. Here, we introduce covalent, orthosteric chemical probes for RORγt that occlude the binding of canonical, orthosteric ligands but still allow allosteric ligand binding. Ultimately, these probes could be used to underpin screening approaches for the unambiguous and rapid identification of novel allosteric RORγt ligands.

Orthosteric and Allosteric Dual Targeting of the Nuclear Receptor RORγt with a Bitopic Ligand

The RORγt nuclear receptor (NR) is of critical importance for the differentiation and proliferation of T helper 17 (Th17) cells and their production of the pro-inflammatory cytokine IL-17a. Dysregulation of RORγt has been linked to various autoimmune diseases, and small molecule inhibition of RORγt is therefore an attractive strategy to treat these diseases. RORγt is a unique NR in that it contains both a canonical, orthosteric and a second, allosteric ligand binding site in its ligand binding domain (LBD). Hence, dual targeting of both binding pockets constitutes an attractive alternative molecular entry for pharmacological modulation. Here, we report a chemical biology approach to develop a bitopic ligand for the RORγt NR, enabling concomitant engagement of both binding pockets. Three candidate bitopic ligands, Bit-L15, Bit-L9, and Bit-L4, comprising an orthosteric and allosteric RORγt pharmacophore linked via a polyethylene glycol (PEG) linker, were designed, synthesized, and evaluated to examine the influence of linker length on the RORγt binding mode. Bit-L15 and Bit-L9 show convincing evidence of concomitant engagement of both RORγt binding pockets, while the shorter Bit-L4 does not show this evidence, as was anticipated during the ligand design. As the most potent bitopic RORγt ligand, Bit-L15, antagonizes RORγt function in a potent manner in both a biochemical and cellular context. Furthermore, Bit-L15 displays an increased selectivity for RORγt over RORα and PPARγ compared to the purely orthosteric and allosteric parent compounds. Combined, these results highlight potential advantages of bitopic NR modulation over monovalent targeting strategies.

Discovery of an Allosteric Ligand Binding Site in SMYD3 Lysine Methyltransferase.

SMYD3 is a multifunctional epigenetic enzyme with lysine methyltransferase activity and various interaction partners. It is implicated in the pathophysiology of cancers but with an unclear mechanism. To discover tool compounds for clarifying its biochemistry and potential as a therapeutic target, a set of drug‐like compounds was screened in a biosensor‐based competition assay. Diperodon was identified as an allosteric ligand; its R and S enantiomers were isolated, and their affinities to SMYD3 were determined (K D=42 and 84 μM, respectively). Co‐crystallization revealed that both enantiomers bind to a previously unidentified allosteric site in the C‐terminal protein binding domain, consistent with its weak inhibitory effect. No competition between diperodon and HSP90 (a known SMYD3 interaction partner) was observed although SMYD3–HSP90 binding was confirmed (K D=13 μM). Diperodon clearly represents a novel starting point for the design of tool compounds interacting with a druggable allosteric site, suitable for the exploration of noncatalytic SMYD3 functions and therapeutics with new mechanisms of action.

Cooperativity between the orthosteric and allosteric ligand binding sites of RORγt

Cooperative ligand binding is an important phenomenon in biological systems where ligand binding influences the binding of another ligand at an alternative site of the protein via an intramolecular network of interactions. The underlying mechanisms behind cooperative binding remain poorly understood, primarily due to the lack of structural data of these ternary complexes. Using time-resolved fluorescence resonance energy transfer (TR-FRET) studies, we show that cooperative ligand binding occurs for RORγt, a nuclear receptor associated with the pathogenesis of autoimmune diseases. To provide the crucial structural insights, we solved 12 crystal structures of RORγt simultaneously bound to various orthosteric and allosteric ligands. The presence of the orthosteric ligand induces a clamping motion of the allosteric pocket via helices 4 to 5. Additional molecular dynamics simulations revealed the unusual mechanism behind this clamping motion, with Ala355 shifting between helix 4 and 5. The orthosteric RORγt agonists regulate the conformation of Ala355, thereby stabilizing the conformation of the allosteric pocket and cooperatively enhancing the affinity of the allosteric inverse agonists.

Novel dynamic residue network analysis approaches to study allosteric modulation: SARS-CoV-2 Mpro and its evolutionary mutations as a case study.

The rational search for allosteric modulators and the allosteric mechanisms of these modulators in the presence of mutations is a relatively unexplored field. Here, we established novel in silico approaches and applied them to SARS-CoV-2 main protease (Mpro) as a case study. First, we identified six potential allosteric modulators. Then, we focused on understanding the allosteric effects of these modulators on each of its protomers. We introduced a new combinatorial approach and dynamic residue network (DRN) analysis algorithms to examine patterns of change and conservation of critical nodes, according to five independent criteria of network centrality. We observed highly conserved network hubs for each averaged DRN metric on the basis of their existence in both protomers in the absence and presence of all ligands (persistent hubs). We also detected ligand specific signal changes. Using eigencentrality (EC) persistent hubs and ligand introduced hubs we identified a residue communication path connecting the allosteric binding site to the catalytic site. Finally, we examined the effects of the mutations on the behavior of the protein in the presence of selected potential allosteric modulators and investigated the ligand stability. One crucial outcome was to show that EC centrality hubs form an allosteric communication path between the allosteric ligand binding site to the active site going through the interface residues of domains I and II; and this path was either weakened or lost in the presence of some of the mutations. Overall, the results revealed crucial aspects that need to be considered in rational computational drug discovery.

Soluble CD40L activates soluble and cell-surface integrin αvβ3, α5β1, and α4β1 by binding to the allosteric ligand-binding site (site 2)

CD40L is a member of the TNF superfamily that participates in immune cell activation. It binds to and signals through several integrins, including αvβ3 and α5β1, which bind to the trimeric interface of CD40L. We previously showed that several integrin ligands can bind to the allosteric site (site 2), which is distinct from the classical ligand-binding site (site 1), raising the question of if CD40L activates integrins. In our explorations of this question, we determined that integrin α4β1, which is prevalently expressed on the same CD4+ T cells as CD40L, is another receptor for CD40L. Soluble (s)CD40L activated soluble integrins αvβ3, α5β1, and α4β1 in cell-free conditions, indicating that this activation does not require inside-out signaling. Moreover, sCD40L activated cell-surface integrins in CHO cells that do not express CD40. To learn more about the mechanism of binding, we determined that sCD40L bound to a cyclic peptide from site 2. Docking simulations predicted that the residues of CD40L that bind to site 2 are located outside of the CD40L trimer interface, at a site where four HIGM1 (hyper-IgM syndrome type 1) mutations are clustered. We tested the effect of these mutations, finding that the K143T and G144E mutants were the most defective in integrin activation, providing support that this region interacts with site 2. We propose that allosteric integrin activation by CD40L also plays a role in CD40L signaling, and defective site 2 binding may be related to the impaired CD40L signaling functions of these HIGM1 mutants.

Discovery of an Inhibitor for the TREK-1 Channel Targeting an Intermediate Transition State of Channel Gating.

Modulators can be designed to stabilize the inactive and active states of ion channels, but whether intermediate (IM) states of channel gating are druggable remains underexplored. In this study, using molecular dynamics simulations of the TWIK-related potassium channel 1 (TREK-1) channel, a two-pore domain potassium channel, we captured an IM state during the transition from the down (inactive) state to the up (active-like) state. The IM state contained a druggable allosteric pocket that was not present in the down or up state. Drug design targeting the pocket led to the identification of the TKIM compound as an inhibitor of TREK-1. Using integrated methods, we verified that TKIM binds to the pocket of the IM state of TREK-1, which differs from the binding of common inhibitors, which bind to channels in the inactive state. Overall, this study identified an allosteric ligand-binding site and a new mechanistic inhibitor for TREK-1, suggesting that IM states of ion channels may be promising druggable targets for use in discovering allosteric modulators.

Enhancing the Signaling of GPCRs via Orthosteric Ions.

G protein-coupled receptors play essential roles in cellular processes such as neuronal signaling, vision, olfaction, tasting, and metabolism. As GPCRs are the most important drug targets, understanding their interactions with ligands is of utmost importance for discovering related new medicines. In many GPCRs, an allosteric sodium ion next to the highly conserved residue D2.50 has been proposed to stabilize the inactive receptor state by mediating interactions between transmembrane helices. Here, we probed the existence of internal and functionally important sodium ions in the dopamine D2 receptor, using molecular dynamics simulations. Besides a new sodium ion at the allosteric ligand binding site, we discovered an additional sodium ion, located close to the orthosteric ligand binding site. Through cell-based activation assays, the signaling of D2 receptor with site-specific mutations was tested against a series of chemically modified agonists. We concluded an important structural role of this newly discovered orthosteric sodium ion in modulating the receptor signaling: It enables the coordination of a polar residue in the ligand binding site with an appropriately designed agonist molecule. An identical interaction was also observed in a recently released high-resolution crystal structure of mu-opioid receptor, which was reresolved in this work. Probably because of similar interactions, various metal ions have been found to increase the signaling of many other GPCRs. This unique principle and strategy could be used to optimize the drug activity of GPCR. Our findings open a new mechanistic opportunity of GPCR signaling and help design the next generation of drugs targeting GPCRs.

How to make an undruggable enzyme druggable: lessons from ras proteins.

Significant advances have been made toward discovering allosteric inhibitors for challenging drug targets such as the Ras family of membrane-associated signaling proteins. Malfunction of Ras proteins due to somatic mutations is associated with up to a quarter of all human cancers. Computational techniques have played critical roles in identifying and characterizing allosteric ligand-binding sites on these proteins, and to screen ligand libraries against those sites. These efforts, combined with a wide range of biophysical, structural, biochemical and cell biological experiments, are beginning to yield promising inhibitors to treat malignancies associated with mutated Ras proteins. In this chapter, we discuss some of these developments and how the lessons learned from Ras might be applied to similar other challenging drug targets.

Pharmacology and function of the orphan GPR139 G protein-coupled receptor.

G protein‐coupled receptors (GPCRs) constitute the largest family of receptors and membrane proteins in the human genome with ~800 members of which half are olfactory. GPCRs are activated by a very broad range of endogenous signalling molecules and are involved in a plethora of physiological functions. All GPCRs contain a transmembrane domain, consisting of a bundle of seven α‐helices spanning the cell membrane, and forming the majority of the known ortho‐ or allosteric ligand binding sites. Due to their many physiological functions and the accessible and druggable transmembrane pocket, GPCRs constitute the largest family of drug targets mediating the actions of 34% of currently marketed drugs. GPCRs activate one or more of the four G protein families (Gq/11, Gi/o, Gs and G12/13) and/or ß‐arrestin. About a third of the non‐olfactory GPCRs are referred to as orphan receptors which means that their endogenous agonist(s) have not yet been found or firmly established. In this MiniReview, we focus on the orphan GPR139 receptor, for which the aromatic amino acids L‐Trp and L‐Phe as well as ACTH/α‐MSH‐related peptides have been proposed as endogenous agonists. GPR139 has been reported to activate several G protein pathways of which Gq/11 is the primary one. The receptor shows the highest expression in the striatum, thalamus, hypothalamus, pituitary and habenula of the human, rat and mouse CNS. We review the surrogate agonists and antagonists that have been published as well as the agonist pharmacophore and binding site. Finally, the putative physiological functions and therapeutic potential are outlined.

CavityPlus: a web server for protein cavity detection with pharmacophore modelling, allosteric site identification and covalent ligand binding ability prediction.

CavityPlus is a web server that offers protein cavity detection and various functional analyses. Using protein three-dimensional structural information as the input, CavityPlus applies CAVITY to detect potential binding sites on the surface of a given protein structure and rank them based on ligandability and druggability scores. These potential binding sites can be further analysed using three submodules, CavPharmer, CorrSite, and CovCys. CavPharmer uses a receptor-based pharmacophore modelling program, Pocket, to automatically extract pharmacophore features within cavities. CorrSite identifies potential allosteric ligand-binding sites based on motion correlation analyses between cavities. CovCys automatically detects druggable cysteine residues, which is especially useful to identify novel binding sites for designing covalent allosteric ligands. Overall, CavityPlus provides an integrated platform for analysing comprehensive properties of protein binding cavities. Such analyses are useful for many aspects of drug design and discovery, including target selection and identification, virtual screening, de novo drug design, and allosteric and covalent-binding drug design. The CavityPlus web server is freely available at http://repharma.pku.edu.cn/cavityplus or http://www.pkumdl.cn/cavityplus.

Stromal cell-derived factor-1 (CXCL12) activates integrins by direct binding to an allosteric ligand-binding site (site 2) of integrins without CXCR4.

Abstract Leukocyte arrest on the endothelial cell surface during leukocyte extravasation is induced by rapid integrin activation by chemokines. We recently reported that fractalkine induces integrin activation without its receptor CX3CR1 through binding to the allosteric site (site 2) of integrins. Peptides from site 2 bound to fractalkine and suppressed integrin activation by fractalkine. We hypothesized that this is not limited to membrane-bound fractalkine. We studied if stromal cell-derived factor 1 (SDF1), another chemokine that plays a critical role in leukocyte arrest, activates integrins through binding to site 2. We describe here that (1) SDF1 activated soluble integrin αvβ3 in cell-free conditions, suggesting that SDF1 can activate αvβ3 without CXCR4; (2) site 2 peptide bound to SDF1, suggesting that SDF1 binds to site 2; (3) SDF1 activated integrins αvβ3, α4β1, and α5β1 on CHO cells (CXCR4-negative) and site 2 peptide suppressed the activation; (4) A CXCR4 antagonist AMD3100 or pertussis toxin did not affect the site 2-mediated integrin activation by SDF1; (5) Cell-surface integrins were fully activated in 1 min (much faster than activation of soluble αvβ3) and the activation lasted at least for 1 h. We propose that the binding of SDF1 to cell-surface proteoglycan facilitates the allosteric activation process; (6) Mutations in the predicted site 2-binding site in SDF1 suppressed integrin activation. These results suggest that SDF1 (e.g., presented on proteoglycans) can rapidly activate integrins in an allosteric manner by binding to site 2 in the absence of CXCR4 or other GPCR. The allosteric integrin activation by SDF1 is a novel target for drug discovery.

Stromal cell-derived factor-1 (CXCL12) activates integrins by direct binding to an allosteric ligand-binding site (site 2) of integrins without CXCR4.

Leukocyte arrest on the endothelial cell surface during leukocyte extravasation is induced by rapid integrin activation by chemokines. We recently reported that fractalkine induces integrin activation without its receptor CX3CR1 through binding to the allosteric site (site 2) of integrins. Peptides from site 2 bound to fractalkine and suppressed integrin activation by fractalkine. We hypothesized that this is not limited to membrane-bound fractalkine. We studied whether stromal cell-derived factor-1 (SDF1), another chemokine that plays a critical role in leukocyte arrest, activates integrins through binding to site 2. We describe here that (1) SDF1 activated soluble integrin αvβ3 in cell-free conditions, suggesting that SDF1 can activate αvβ3 without CXCR4; (2) site 2 peptide bound to SDF1, suggesting that SDF1 binds to site 2; (3) SDF1 activated integrins αvβ3, α4β1, and α5β1 on CHO cells (CXCR4-negative) and site 2 peptide suppressed the activation; (4) A CXCR4 antagonist AMD3100 did not affect the site 2-mediated integrin activation by SDF1; (5) Cell-surface integrins were fully activated in 1 min (much faster than activation of soluble αvβ3) and the activation lasted at least for 1 h. We propose that the binding of SDF1 to cell-surface proteoglycan facilitates the allosteric activation process; (6) Mutations in the predicted site 2-binding site in SDF1 suppressed integrin activation. These results suggest that SDF1 (e.g. presented on proteoglycans) can rapidly activate integrins in an allosteric manner by binding to site 2 in the absence of CXCR4. The allosteric integrin activation by SDF1 is a novel target for drug discovery.

In silico mapping of allosteric ligand binding sites in type-1 cannabinoid receptor.

The recent resolution of the crystal structure of type-1 cannabinoid receptor (CB1 ) and the discovery of novel modulators for this target open the way to the possibility of elucidating the structural requirements for CB1 binding, and thereby facilitate a rational drug design. Compounds that target the orthosteric site of CB1 in some cases have shown side effects. Allosteric modulators could potentially avoid these side effects by influencing binding and/or efficacy of orthosteric ligands. Here, we summarize and compare previous data on different putative allosteric binding sites observed in CB1 homology models with an in silico docking study of the recently published crystal structure of the same receptor on endogenous and natural hydrophobic ligands that act as positive allosteric modulators and negative allosteric modulators of CB1 . In particular, a lipid-exposed pocket targeted by most of the tested molecules is reported and discussed.

An allosteric ligand-binding site in the extracellular cap of K2P channels

Two-pore domain potassium (K2P) channels generate leak currents that are responsible for the maintenance of the resting membrane potential, and they are thus potential drug targets for treating diseases. Here, we identify N-(4-cholorphenyl)-N-(2-(3,4-dihydrosioquinolin-2(1H)-yl)-2-oxoethyl)methanesulfonamide (TKDC) as an inhibitor of the TREK subfamily, including TREK-1, TREK-2 and TRAAK channels. Using TKDC as a chemical probe, a study combining computations, mutagenesis and electrophysiology reveals a K2P allosteric ligand-binding site located in the extracellular cap of the channels. Molecular dynamics simulations suggest that ligand-induced allosteric conformational transitions lead to blockage of the ion conductive pathway. Using virtual screening approach, we identify other inhibitors targeting the extracellular allosteric ligand-binding site of these channels. Overall, our results suggest that the allosteric site at the extracellular cap of the K2P channels might be a promising drug target for these membrane proteins.

Structural and biochemical characterization of ligand recognition by CysB, the master regulator of sulfate metabolism

CysB, a member of LysR-type transcriptional regulators, up-regulates the expression of genes associated with sulfate metabolism and cysteine biosynthesis. CysB is activated under sulfur limiting conditions by O-acetylserine (OAS) and N-acetylserine (NAS), but the activation mechanism of CysB remain unknown. Here, we report four crystal structures of ligand binding domains of CysB (CysB-LBD) in apo form and in complex with sulfate, OAS, and NAS. Our results show that CysB has two distinct allosteric ligand binding sites; a sulfate and NAS specific site-1 and a second, NAS and OAS specific site-2. All three ligands bind through the induced-fit mechanism. Surprisingly, OAS remodels the site-1 by binding to site-2, suggesting that site-1 and site-2 are coupled allosterically. Using DNA binding and site-directed mutagenesis approach, we show that OAS enhances NAS mediated activation and mutation at site-1 has no effect on site-2 mediated OAS activation. Results indicate that inducer binding triggered signals from OAS-Specific site-2 are relayed to DBD through site-1. Together, results presented here suggest that induced-fit binding and allosteric coupling between two ligand binding sites and DBD underlie the key feature of CysB activation. Further, this study provides first structural glimpse into recognition of inducer ligands by CysB and provides a general framework to understand how LTTR family regulators respond to dual activators.