The Glycine Receptor Allosteric Ligands Library (GRALL)

A millisecond coarse-grained simulation approach to decipher allosteric cannabinoid binding at the glycine receptor α1

Glycine receptors (GlyR) are regulated by small-molecule binding at several allosteric sites. Cannabinoids like tetrahydrocannabinol (THC) and N-arachidonyl-ethanol-amide (AEA) potentiate the GlyR response but their mechanism of action is not fully established. By combining millisecond coarse-grained (CG) MD simulations powered by Martini 3 with backmapping to all-atom representations, we have characterized the cannabinoid-binding site(s) at the zebrafish GlyR-α1 active state with atomic resolution. Based on hundreds of thousand ligand-binding events, we find that cannabinoids bind to the transmembrane domain of the receptor at both intrasubunit and intersubunit sites. For THC, the intrasubunit binding mode predicted in simulation is in excellent agreement with recent cryo-EM structures, while intersubunit binding recapitulates in full previous mutagenesis experiments. Intriguingly, AEA is predicted to bind at the same intersubunit site despite the strikingly different chemistry. Statistical analyses of the ligand-receptor interactions highlight potentially relevant residues for GlyR potentiation, offering experimentally testable predictions. The predictions for AEA have been validated by electrophysiology recordings of rationally designed mutants. The results highlight the existence of multiple cannabinoid-binding sites for the allosteric regulation of GlyR and put forward an effective strategy for the identification and structural characterization of allosteric binding sites.

Allosteric Cholesterol Site in Glycine Receptors Characterized through Molecular Simulations

Glycine receptors are pentameric ligand-gated ion channels that conduct chloride ions across postsynaptic membranes to facilitate fast inhibitory neurotransmission. In addition to gating by the glycine agonist, interactions with lipids and other compounds in the surrounding membrane environment modulate their function, but molecular details of these interactions remain unclear, in particular, for cholesterol. Here, we report coarse-grained simulations in a model neuronal membrane for three zebrafish glycine receptor structures representing apparent resting, open, and desensitized states. We then converted the systems to all-atom models to examine detailed lipid interactions. Cholesterol bound to the receptor at an outer-leaflet intersubunit site, with a preference for the open and desensitized versus resting states, indicating that it can bias receptor function. Finally, we used short atomistic simulations and iterative amino acid perturbations to identify residues that may mediate allosteric gating transitions. Frequent cholesterol contacts in atomistic simulations clustered with residues identified by perturbation analysis and overlapped with mutations influencing channel function and pathology. Cholesterol binding at this site was also observed in a recently reported pig heteromeric glycine receptor. These results indicate state-dependent lipid interactions relevant to allosteric transitions of glycine receptors, including specific amino acid contacts applicable to biophysical modeling and pharmaceutical design.

Machine learning approaches in predicting allosteric sites

Allosteric regulation is a fundamental biological mechanism that can control critical cellular processes via allosteric modulator binding to protein distal functional sites. The advantages of allosteric modulators over orthosteric ones have sparked the development of numerous computational approaches, such as the identification of allosteric binding sites, to facilitate allosteric drug discovery. Building on the success of machine learning (ML) models for solving complex problems in biology and chemistry, several ML models for predicting allosteric sites have been developed. In this review, we provide an overview of these models and discuss future perspectives powered by the field of artificial intelligence such as protein language models.

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.

Allosteric modulation of SHP2: Quest from known to unknown

Src homology-2 domain-containing protein tyrosine phosphatase-2 (SHP2) is a key regulatory factor in the cell cycle and its activating mutations play an important role in the development of various cancers, making it an important target for antitumor drugs. Due to the highly conserved amino acid sequence and positively charged nature of the active site of SHP2, it is difficult to discover inhibitors with high affinity for the catalytic site of SHP2 and sufficient cell permeability, making it considered an "undruggable" target. However, the discovery of allosteric regulation mechanisms provides new opportunities for transforming undruggable targets into druggable ones. Given the limitations of orthosteric inhibitors, SHP2 allosteric inhibitors have become a more selective and safer research direction. In this review, we elucidate the oncogenic mechanism of SHP2 and summarize the discovery methods of SHP2 allosteric inhibitors, providing new strategies for the design and improvement of SHP2 allosteric inhibitors.

Novel Phenylene Lipids That Are Positive Allosteric Modulators of Glycine Receptors and Inhibitors of Glycine Transporter 2

Chronic pain is a complex condition that remains resistant to current therapeutics. We previously synthesized a series of N-acyl amino acids (NAAAs) that inhibit the glycine transporter, GlyT2, some of which are also positive allosteric modulators of glycine receptors (GlyRs). In this study, we have synthesized a library of NAAAs that contain a phenylene ring within the acyl tail with the objective of improving efficacy at both GlyT2 and GlyRs and also identifying compounds that are efficacious as dual-acting modulators to enhance glycine neurotransmission. The most efficacious positive allosteric modulator of GlyRs was 2-[8-(2-octylphenyl)octanoylamino]acetic acid (8-8 OPGly) which potentiates the EC5 for glycine activation of GlyRα1 by 1500% with an EC50 of 664 nM. Phenylene-containing NAAAs with a lysine headgroup were the most potent inhibitors of GlyT2 with (2S)-6-amino-2-[8-(3-octylphenyl)octanoylamino]hexanoic acid (8-8 MPLys) inhibiting GlyT2 with an IC50 of 32 nM. The optimal modulator across both proteins was (2S)-6-amino-2-[8-(2-octylphenyl)octanoylamino]hexanoic acid (8-8 OPLys), which inhibits GlyT2 with an IC50 of 192 nM and potentiates GlyRs by up to 335% at 1 μM. When tested in a dual GlyT2/GlyRα1 expression system, 8-8 OPLys caused the greatest reductions in the EC50 for glycine. This suggests that the synergistic effects of a dual-acting modulator cause greater enhancements in glycinergic activity compared to single-target modulators and may provide an alternate approach to the development of new non-opioid analgesics for the treatment of chronic pain.

Illumination of a progressive allosteric mechanism mediating the glycine receptor activation

Pentameric ligand-gated ion channel mediate signal transduction at chemical synapses by transiting between resting and open states upon neurotransmitter binding. Here, we investigate the gating mechanism of the glycine receptor fluorescently labeled at the extracellular-transmembrane interface by voltage-clamp fluorometry (VCF). Fluorescence reports a glycine-elicited conformational change that precedes pore opening. Low concentrations of glycine, partial agonists or specific mixtures of glycine and strychnine trigger the full fluorescence signal while weakly activating the channel. Molecular dynamic simulations of a partial agonist bound-closed Cryo-EM structure show a highly dynamic nature: a marked structural flexibility at both the extracellular-transmembrane interface and the orthosteric site, generating docking properties that recapitulate VCF data. This work illuminates a progressive propagating transition towards channel opening, highlighting structural plasticity within the mechanism of action of allosteric effectors.

Positive Allosteric Modulators of Glycine Receptors and Their Potential Use in Pain Therapies

Glycine receptors are ligand-gated ion channels that mediate synaptic inhibition throughout the mammalian spinal cord, brainstem, and higher brain regions. They have recently emerged as promising targets for novel pain therapies due to their ability to produce antinociception by inhibiting nociceptive signals within the dorsal horn of the spinal cord. This has greatly enhanced the interest in developing positive allosteric modulators of glycine receptors. Several pharmaceutical companies and research facilities have attempted to identify new therapeutic leads by conducting large-scale screens of compound libraries, screening new derivatives from natural sources, or synthesizing novel compounds that mimic endogenous compounds with antinociceptive activity. Advances in structural techniques have also led to the publication of multiple high-resolution structures of the receptor, highlighting novel allosteric binding sites and providing additional information for previously identified binding sites. This has greatly enhanced our understanding of the functional properties of glycine receptors and expanded the structure activity relationships of novel pharmacophores. Despite this, glycine receptors are yet to be used as drug targets due to the difficulties in obtaining potent, selective modulators with favorable pharmacokinetic profiles that are devoid of side effects. This review presents a summary of the structural basis for how current compounds cause positive allosteric modulation of glycine receptors and discusses their therapeutic potential as analgesics. SIGNIFICANCE STATEMENT: Chronic pain is a major cause of disability, and in Western societies, this will only increase as the population ages. Despite the high level of prevalence and enormous socioeconomic burden incurred, treatment of chronic pain remains limited as it is often refractory to current analgesics, such as opioids. The National Institute for Drug Abuse has set finding effective, safe, nonaddictive strategies to manage chronic pain as their top priority. Positive allosteric modulators of glycine receptors may provide a therapeutic option.

Structure and Mechanism of Glycine Receptor Elucidated by Cryo-Electron Microscopy

Glycine receptors (GlyRs) are pentameric ion channels that mediate fast inhibitory neurotransmission. GlyRs are found in the central nervous system including the spinal cord, brain stem, and cerebellum, as well as in the retina, sperm, macrophages, hippocampus, cochlea, and liver. Due to their crucial roles in counter-balancing excitatory signals and pain signal transmission, GlyR dysfunction can lead to severe diseases, and as a result, compounds that modify GlyR activity may have tremendous therapeutic potential. Despite this potential, the development of GlyR-specific small-molecule ligands is lacking. Over the past few years, high-resolution structures of both homomeric and heteromeric GlyRs structures in various conformations have provided unprecedented details defining the pharmacology of ligand binding, subunit composition, and mechanisms of channel gating. These high-quality structures will undoubtedly help with the development of GlyR-targeted therapies.

Nicotinic receptors: From protein allostery to computational neuropharmacology

We propose an extension and further development of the Monod-Wyman-Changeux model for allosteric transitions of regulatory proteins to brain communications and specifically to neurotransmitters receptors, with the nicotinic acetylcholine receptor (nAChR) as a model of ligand-gated ion channels. The present development offers an expression of the change of the gating isomerization constant caused by pharmacological ligand binding in terms of its value in the absence of ligands and several "modulation factors", which vary with orthosteric ligand binding (agonists/antagonists), allosteric ligand binding (positive allosteric modulators/negative allosteric modulators) and receptor desensitization. The new - explicit - formulation of such "modulation factors", provides expressions for the pharmacological attributes of potency, efficacy, and selectivity for the modulatory ligands (including endogenous neurotransmitters) in terms of their binding affinity for the active, resting, and desensitized states of the receptor. The current formulation provides ways to design neuroactive compounds with a controlled pharmacological profile, opening the field of computational neuro-pharmacology.

Attenuation of clinical and immunological outcomes during SARS-CoV-2 infection by ivermectin

The devastating pandemic due to SARS-CoV-2 and the emergence of antigenic variants that jeopardize the efficacy of current vaccines create an urgent need for a comprehensive understanding of the pathophysiology of COVID-19, including the contribution of inflammation to disease. It also warrants for the search of immunomodulatory drugs that could improve disease outcome. Here, we show that standard doses of ivermectin (IVM), an anti-parasitic drug with potential immunomodulatory activities through the cholinergic anti-inflammatory pathway, prevent clinical deterioration, reduce olfactory deficit, and limit the inflammation of the upper and lower respiratory tracts in SARS-CoV-2-infected hamsters. Whereas it has no effect on viral load in the airways of infected animals, transcriptomic analyses of infected lungs reveal that IVM dampens type I interferon responses and modulates several other inflammatory pathways. In particular, IVM dramatically reduces the Il-6/Il-10 ratio in lung tissue and promotes macrophage M2 polarization, which might account for the more favorable clinical presentation of IVM-treated animals. Altogether, this study supports the use of immunomodulatory drugs such as IVM, to improve the clinical condition of SARS-CoV-2-infected patients.

Identification of N-acyl amino acids that are positive allosteric modulators of glycine receptors

Glycine receptors (GlyRs) mediate inhibitory neurotransmission within the spinal cord and play a crucial role in nociceptive signalling. This makes them primary targets for the development of novel chronic pain therapies. Endogenous lipids have previously been shown to modulate glycine receptors and produce analgesia in pain models, however little is known about what chemical features mediate these effects. In this study, we characterised lipid modulation of GlyRs by screening a library of N-acyl amino acids across all receptor subtypes and determined chemical features crucial for their activity. Acyl-glycine's with a C18 carbon tail were found to produce the greatest potentiation, and require a cis double bond within the central region of the carbon tail (ω6 - ω9) to be active. At 1 µM, C18 ω6,9 glycine potentiated glycine induced currents in α₃ and α₃β receptors by over 50%, and α₁, α₂, α₁β and α₂β receptors by over 100%. C18 ω9 glycine (N-oleoyl glycine) significantly enhance glycine induced peak currents and cause a dose-dependent shift in the glycine concentration response. In the presence of 3 µM C18 ω9 glycine, the EC_5o of glycine at the α₁ receptor was reduced from 17 µM to 10 µM. This study has identified several acyl-amino acids which are positive allosteric modulators of GlyRs and make promising lead compounds for the development of novel chronic pain therapies.

Modulators of the Inhibitory Glycine Receptor

The inhibitory glycine receptor is a member of the Cys-loop superfamily of ligand-gated ion channels. It is the principal mediator of rapid synaptic inhibition in the spinal cord and brainstem and plays an important role in the modulation of higher brain functions including vision, hearing, and pain signaling. Glycine receptor function is controlled by only a few agonists, while the number of antagonists and positive or biphasic modulators is steadily increasing. These modulators are important for the study of receptor activation and regulation and have found clinical interest as potential analgesics and anticonvulsants. High-resolution structures of the receptor have become available recently, adding to our understanding of structure-function relationships and revealing agonistic, inhibitory, and modulatory sites on the receptor protein. This Review presents an overview of compounds that activate, inhibit, or modulate glycine receptor function in vitro and in vivo.