Functional validation of virtual screening for novel agents with general anesthetic action at ligand-gated ion channels

A bupropion modulatory site in the Gloeobacter violaceus ligand-gated ion channel

Bupropion is an atypical antidepressant and smoking cessation drug that causes adverse effects such as insomnia, irritability, and anxiety. Bupropion inhibits dopamine and norepinephrine reuptake transporters and eukaryotic cation-conducting pentameric ligand-gated ion channels, such as nicotinic acetylcholine and serotonin type 3A receptors, at clinically relevant concentrations. Here, we demonstrate that bupropion also inhibits a prokaryotic homolog of pentameric ligand-gated ion channels, the Gloeobacter violaceus ligand-gated ion channel (GLIC). Using the GLIC as a model, we used molecular docking to predict binding sites for bupropion. Bupropion was found to bind to several sites within the transmembrane domain, with the predominant site being localized to the interface between transmembrane segments M1 and M3 of two adjacent subunits. Residues W213, T214, and W217 in the first transmembrane segment, M1, and F267 and I271 in the third transmembrane segment, M3, most frequently reside within a 4 Å distance from bupropion. We then used single amino acid substitutions at these positions and two-electrode voltage-clamp recordings to determine their impact on bupropion inhibitory effects. The substitution T214F alters bupropion potency by shifting the half-maximal inhibitory concentration to a 13-fold higher value compared to wild-type GLIC. Residue T214 is found within a previously identified binding pocket for neurosteroids and lipids in the GLIC. This intersubunit binding pocket is structurally conserved and almost identical to a binding pocket described for neurosteroids in γ-aminobutyric acid type A receptors. Our data thus suggest that the T214 that lines a previously identified lipophilic binding pocket in GLIC and γ-aminobutyric acid type A receptors is also a modulatory site for bupropion interaction with the GLIC.

Identification of a binding site for bupropion in Gloeobacter violaceus ligand-gated ion channel

Bupropion is an atypical antidepressant and smoking cessation drug which causes adverse effects such as insomnia, irritability, and anxiety. Bupropion inhibits dopamine and norepinephrine reuptake transporters and eukaryotic cation-conducting pentameric ligand-gated ion channels (pLGICs), such as nicotinic acetylcholine (nACh) and serotonin type 3A (5-HT3A) receptors, at clinically relevant concentrations. However, the binding sites and binding mechanisms of bupropion are still elusive. To further understand the inhibition of pLGICs by bupropion, in this work, using a prokaryotic homologue of pLGICs as a model, we examined the inhibitory potency of bupropion in Gloeobacter violaceus ligand-gated ion channel (GLIC), a proton-gated ion channel. Bupropion inhibited proton-induced currents in GLIC with an inhibitory potency of 14.9 ± 2.0 μM, comparable to clinically attainable concentrations previously shown to also modulate eukaryotic pLGICs. Using single amino acid substitutions in GLIC and two-electrode voltage-clamp recordings, we further determined a binding site for bupropion in the lower third of the first transmembrane segment M1 at residue T214. The sidechain of M1 T214 together with additional residues of M1 and also of M3 of the adjacent subunit have previously been shown to contribute to binding of other lipophilic molecules like allopregnanolone and pregnanolone.

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.

Physicochemistry shapes bioactivity landscape of pan-ABC transporter modulators: Anchor point for innovative Alzheimer's disease therapeutics

Alzheimer's disease (AD) is a devastating neurological disorder characterized by the pathological accumulation of macromolecular Aβ and tau leading to neuronal death. Drugs approved to treat AD may ameliorate disease symptoms, however, no curative treatment exists. Aβ peptides were discovered to be substrates of adenosine triphosphate-(ATP)-binding cassette (ABC) transporters. Activators of these membrane-bound efflux proteins that promote binding and/or translocation of Aβ could revolutionize AD medicine. The knowledge about ABC transporter activators is very scarce, however, the few molecules that were reported contain substructural features of multitarget (pan-)ABC transporter inhibitors. A cutting-edge strategy to obtain new drug candidates is to explore and potentially exploit the recently proposed multitarget binding site of pan-ABC transporter inhibitors as anchor point for the development of innovative activators to promote Aβ clearance from the brain. Molecular associations between functional bioactivities and physicochemical properties of small-molecules are key to understand these processes. This study provides an analysis of a recently reported unique multitarget dataset for the correlation between multitarget bioactivity and physicochemistry. Six novel pan-ABC transporter inhibitors were validated containing substructural features of ABC transporter activators, which underpins the relevance of the multitarget binding site for the targeted development of novel AD diagnostics and therapeutics.

Elephants in the Dark: Insights and Incongruities in Pentameric Ligand-gated Ion Channel Models

The superfamily of pentameric ligand-gated ion channels (pLGICs) comprises key players in electrochemical signal transduction across evolution, including historic model systems for receptor allostery and targets for drug development. Accordingly, structural studies of these channels have steadily increased, and now approach 250 depositions in the protein data bank. This review contextualizes currently available structures in the pLGIC family, focusing on morphology, ligand binding, and gating in three model subfamilies: the prokaryotic channel GLIC, the cation-selective nicotinic acetylcholine receptor, and the anion-selective glycine receptor. Common themes include the challenging process of capturing and annotating channels in distinct functional states; partially conserved gating mechanisms, including remodeling at the extracellular/transmembrane-domain interface; and diversity beyond the protein level, arising from posttranslational modifications, ligands, lipids, and signaling partners. Interpreting pLGIC structures can be compared to describing an elephant in the dark, relying on touch alone to comprehend the many parts of a monumental beast: each structure represents a snapshot in time under specific experimental conditions, which must be integrated with further structure, function, and simulations data to build a comprehensive model, and understand how one channel may fundamentally differ from another.

Structural Basis for a Bimodal Allosteric Mechanism of General Anesthetic Modulation in Pentameric Ligand-Gated Ion Channels

Ion channel modulation by general anesthetics is a vital pharmacological process with implications for receptor biophysics and drug development. Functional studies have implicated conserved sites of both potentiation and inhibition in pentameric ligand-gated ion channels, but a detailed structural mechanism for these bimodal effects is lacking. The prokaryotic model protein GLIC recapitulates anesthetic modulation of human ion channels, and it is accessible to structure determination in both apparent open and closed states. Here, we report ten X-ray structures and electrophysiological characterization of GLIC variants in the presence and absence of general anesthetics, including the surgical agent propofol. We show that general anesthetics can allosterically favor closed channels by binding in the pore or favor open channels via various subsites in the transmembrane domain. Our results support an integrated, multi-site mechanism for allosteric modulation, and they provide atomic details of both potentiation and inhibition by one of the most common general anesthetics.

High-throughput Screening in Larval Zebrafish Identifies Novel Potent Sedative-hypnotics

WHAT WE ALREADY KNOW ABOUT THIS TOPIC

WHAT THIS ARTICLE TELLS US THAT IS NEW: BACKGROUND:: Many general anesthetics were discovered empirically, but primary screens to find new sedative-hypnotics in drug libraries have not used animals, limiting the types of drugs discovered. The authors hypothesized that a sedative-hypnotic screening approach using zebrafish larvae responses to sensory stimuli would perform comparably to standard assays, and efficiently identify new active compounds.

METHODS

The authors developed a binary outcome photomotor response assay for zebrafish larvae using a computerized system that tracked individual motions of up to 96 animals simultaneously. The assay was validated against tadpole loss of righting reflexes, using sedative-hypnotics of widely varying potencies that affect various molecular targets. A total of 374 representative compounds from a larger library were screened in zebrafish larvae for hypnotic activity at 10 µM. Molecular mechanisms of hits were explored in anesthetic-sensitive ion channels using electrophysiology, or in zebrafish using a specific reversal agent.

RESULTS

Zebrafish larvae assays required far less drug, time, and effort than tadpoles. In validation experiments, zebrafish and tadpole screening for hypnotic activity agreed 100% (n = 11; P = 0.002), and potencies were very similar (Pearson correlation, r > 0.999). Two reversible and potent sedative-hypnotics were discovered in the library subset. CMLD003237 (EC50, ~11 µM) weakly modulated γ-aminobutyric acid type A receptors and inhibited neuronal nicotinic receptors. CMLD006025 (EC50, ~13 µM) inhibited both N-methyl-D-aspartate and neuronal nicotinic receptors.

CONCLUSIONS

Photomotor response assays in zebrafish larvae are a mechanism-independent platform for high-throughput screening to identify novel sedative-hypnotics. The variety of chemotypes producing hypnosis is likely much larger than currently known.

Allosteric potentiation of a ligand-gated ion channel is mediated by access to a deep membrane-facing cavity

Theories of general anesthesia have shifted in focus from bulk lipid effects to specific interactions with membrane proteins. Target receptors include several subtypes of pentameric ligand-gated ion channels; however, structures of physiologically relevant proteins in this family have yet to define anesthetic binding at high resolution. Recent cocrystal structures of the bacterial protein GLIC provide snapshots of state-dependent binding sites for the common surgical agent propofol (PFL), offering a detailed model system for anesthetic modulation. Here, we combine molecular dynamics and oocyte electrophysiology to reveal differential motion and modulation upon modification of a transmembrane binding site within each GLIC subunit. WT channels exhibited net inhibition by PFL, and a contraction of the cavity away from the pore-lining M2 helix in the absence of drug. Conversely, in GLIC variants exhibiting net PFL potentiation, the cavity was persistently expanded and proximal to M2. Mutations designed to favor this deepened site enabled sensitivity even to subclinical concentrations of PFL, and a uniquely prolonged mode of potentiation evident up to ∼30 min after washout. Dependence of these prolonged effects on exposure time implicated the membrane as a reservoir for a lipid-accessible binding site. However, at the highest measured concentrations, potentiation appeared to be masked by an acute inhibitory effect, consistent with the presence of a discrete, water-accessible site of inhibition. These results support a multisite model of transmembrane allosteric modulation, including a possible link between lipid- and receptor-based theories that could inform the development of new anesthetics.

From gene networks to drugs: systems pharmacology approaches for AUD

The alcohol research field has amassed an impressive number of gene expression datasets spanning key brain areas for addiction, species (humans as well as multiple animal models), and stages in the addiction cycle (binge/intoxication, withdrawal/negative effect, and preoccupation/anticipation). These data have improved our understanding of the molecular adaptations that eventually lead to dysregulation of brain function and the chronic, relapsing disorder of addiction. Identification of new medications to treat alcohol use disorder (AUD) will likely benefit from the integration of genetic, genomic, and behavioral information included in these important datasets. Systems pharmacology considers drug effects as the outcome of the complex network of interactions a drug has rather than a single drug-molecule interaction. Computational strategies based on this principle that integrate gene expression signatures of pharmaceuticals and disease states have shown promise for identifying treatments that ameliorate disease symptoms (called in silico gene mapping or connectivity mapping). In this review, we suggest that gene expression profiling for in silico mapping is critical to improve drug repurposing and discovery for AUD and other psychiatric illnesses. We highlight studies that successfully apply gene mapping computational approaches to identify or repurpose pharmaceutical treatments for psychiatric illnesses. Furthermore, we address important challenges that must be overcome to maximize the potential of these strategies to translate to the clinic and improve healthcare outcomes.

Permeating disciplines: Overcoming barriers between molecular simulations and classical structure-function approaches in biological ion transport

Ion translocation across biological barriers is a fundamental requirement for life. In many cases, controlling this process-for example with neuroactive drugs-demands an understanding of rapid and reversible structural changes in membrane-embedded proteins, including ion channels and transporters. Classical approaches to electrophysiology and structural biology have provided valuable insights into several such proteins over macroscopic, often discontinuous scales of space and time. Integrating these observations into meaningful mechanistic models now relies increasingly on computational methods, particularly molecular dynamics simulations, while surfacing important challenges in data management and conceptual alignment. Here, we seek to provide contemporary context, concrete examples, and a look to the future for bridging disciplinary gaps in biological ion transport. This article is part of a Special Issue entitled: Beyond the Structure-Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.

Rational screening of peroxisome proliferator-activated receptor-γ agonists from natural products: potential therapeutics for heart failure

CONTEXT

Peroxisome proliferator-activated receptor-γ (PPARγ) is a member of the nuclear hormone receptor superfamily of ligand-activated transcription factors. Activation of PPARγ pathway has been shown to enhance fatty acid oxidation, improve endothelial cell function, and decrease myocardial fibrosis in heart failure. Thus, the protein has been raised as an attractive target for heart failure therapy.

OBJECTIVE

This work attempted to discover new and potent PPARγ agonists from natural products using a synthetic strategy of computer virtual screening and transactivation reporter assay.

MATERIALS AND METHODS

A large library of structurally diverse, drug-like natural products was compiled, from which those with unsatisfactory pharmacokinetic profile and/or structurally redundant compounds were excluded. The binding mode of remaining candidates to PPARγ ligand-binding domain (LBD) was computationally modelled using molecular docking and their relative binding potency was ranked by an empirical scoring scheme. Consequently, eight commercially available hits with top scores were selected and their biological activity was determined using a cell-based reporter-gene assay.

RESULTS

Four natural product compounds, namely ZINC13408172, ZINC4292805, ZINC44179 and ZINC901461, were identified to have high or moderate agonistic potency against human PPARγ with EC₅₀ values of 0.084, 2.1, 0.35 and 5.6 μM, respectively, which are comparable to or even better than that of the approved PPARγ full agonists pioglitazone (EC₅₀=0.16 μM) and rosiglitazone (EC₅₀=0.034 μM). Hydrophobic interactions and van der Waals contacts are the primary chemical forces to stabilize the complex architecture of PPARγ LBD domain with these agonist ligands, while few hydrogen bonds, salt bridges and/or π-π stacking at the complex interfaces confer selectivity and specificity for the domain-agonist recognition.

DISCUSSION AND CONCLUSION

The integrated in vitro-in silico screening strategy can be successfully applied to rational discovery of biologically active compounds. The newly identified natural products with PPARγ agonistic potency are considered as promising lead scaffolds to develop novel chemical therapeutics for heart failure.

Docking Methodologies and Recent Advances

Docking, a molecular modelling method, has wide applications in identification and optimization in modern drug discovery. This chapter addresses the recent advances in the docking methodologies like fragment docking, covalent docking, inverse docking, post processing, hybrid techniques, homology modeling etc. and its protocol like searching and scoring functions. Advances in scoring functions for e.g. consensus scoring, quantum mechanics methods, clustering and entropy based methods, fingerprinting, etc. are used to overcome the limitations of the commonly used force-field, empirical and knowledge based scoring functions. It will cover crucial necessities and different algorithms of docking and scoring. Further different aspects like protein flexibility, ligand sampling and flexibility, and the performance of scoring function will be discussed.

Inter- and Intra-Subunit Butanol/Isoflurane Sites of Action in the Human Glycine Receptor

Glycine receptors (GlyRs) mediate inhibitory neurotransmission and are targets for alcohols and anesthetics in brain. GlyR transmembrane (TM) domains contain critical residues for alcohol/anesthetic action: amino acid A288 in TM3 forms crosslinks with TM1 (I229) in the adjacent subunit as well as TM2 (S267) and TM4 (Y406, W407, I409, Y410) in the same subunit. We hypothesized that these residues may participate in intra-subunit and inter-subunit sites of alcohol/anesthetic action. The following double and triple mutants of GLRA1 cDNA (encoding human glycine receptor alpha 1 subunit) were injected into Xenopus laevis oocytes: I229C/A288C, I229C/A288C/C290S, A288C/Y406C, A288C/W407C, A288C/I409C, and A288C/Y410C along with the corresponding single mutants and wild-type GLRA1. Butanol (22 mM) or isoflurane (0.6 mM) potentiation of GlyR-mediated currents before and after application of the cysteine crosslinking agent HgCl₂ (10 μM) was measured using two-electrode voltage clamp electrophysiology. Crosslinking nearly abolished butanol and isoflurane potentiation in the I229C/A288C and I229C/A288C/C290S mutants but had no effect in single mutants or wild-type. Crosslinking also inhibited butanol and isoflurane potentiation in the TM3-4 mutants (A288C/Y406C, A288C/W407C, A288C/I409C, A288C/Y410C) with no effect in single mutants or wild-type. We extracted proteins from oocytes expressing I229C/288C, A288C/Y410C, or wild-type GlyRs, used mass spectrometry to verify their expression and possible inter-subunit dimerization, plus immunoblotting to investigate the biochemical features of proposed crosslinks. Wild-type GlyR subunits measured about 50 kDa; after crosslinking, the dimeric/monomeric 100:50 kDa band ratio was significantly increased in I229C/288C but not A288C/Y410C mutants or wild-type, providing support for TM1-3 inter-subunit and TM3-4 intra-subunit crosslinking. A GlyR homology model based on the GluCl template provides further evidence for a multi-site model for alcohol/anesthetic interaction with human GLRA1.

Functional characterization of neurotransmitter activation and modulation in a nematode model ligand-gated ion channel

The superfamily of pentameric ligand-gated ion channels includes neurotransmitter receptors that mediate fast synaptic transmission in vertebrates, and are targets for drugs including alcohols, anesthetics, benzodiazepines, and anticonvulsants. However, the mechanisms of ion channel opening, gating, and modulation in these receptors leave many open questions, despite their pharmacological importance. Subtle conformational changes in both the extracellular and transmembrane domains are likely to influence channel opening, but have been difficult to characterize given the limited structural data available for human membrane proteins. Recent crystal structures of a modified Caenorhabditis elegans glutamate-gated chloride channel (GluCl) in multiple states offer an appealing model system for structure-function studies. However, the pharmacology of the crystallographic GluCl construct is not well established. To establish the functional relevance of this system, we used two-electrode voltage-clamp electrophysiology in Xenopus oocytes to characterize activation of crystallographic and native-like GluCl constructs by L-glutamate and ivermectin. We also tested modulation by ethanol and other anesthetic agents, and used site-directed mutagenesis to explore the role of a region of Loop F which was implicated in ligand gating by molecular dynamics simulations. Our findings indicate that the crystallographic construct functionally models concentration-dependent agonism and allosteric modulation of pharmacologically relevant receptors. Specific substitutions at residue Leu174 in loop F altered direct L-glutamate activation, consistent with computational evidence for this region's role in ligand binding. These insights demonstrate conservation of activation and modulation properties in this receptor family, and establish a framework for GluCl as a model system, including new possibilities for drug discovery. In this study, we elucidate the validity of a modified glutamate-gated chloride channel (GluClcryst ) as a structurally accessible model for GABAA receptors. In contrast to native-like controls, GluClcryst exhibits classical activation by its neurotransmitter ligand L-glutamate. The modified channel is also sensitive to allosteric modulators associated with human GABAA receptors, and to site-directed mutations predicted to alter channel opening.

Docking Screens for Novel Ligands Conferring New Biology

It is now plausible to dock libraries of 10 million molecules against targets over several days or weeks. When the molecules screened are commercially available, they may be rapidly tested to find new leads. Although docking retains important liabilities (it cannot calculate affinities accurately nor even reliably rank order high-scoring molecules), it can often can distinguish likely from unlikely ligands, often with hit rates above 10%. Here we summarize the improvements in libraries, target quality, and methods that have supported these advances, and the open access resources that make docking accessible. Recent docking screens for new ligands are sketched, as are the binding, crystallographic, and in vivo assays that support them. Like any technique, controls are crucial, and key experimental ones are reviewed. With such controls, docking campaigns can find ligands with new chemotypes, often revealing the new biology that may be docking's greatest impact over the next few years.

Discovery of a novel general anesthetic chemotype using high-throughput screening

BACKGROUND

The development of novel anesthetics has historically been a process of combined serendipity and empiricism, with most recent new anesthetics developed via modification of existing anesthetic structures.

METHODS

Using a novel high-throughput screen employing the fluorescent anesthetic 1-aminoanthracene and apoferritin as a surrogate for on-pathway anesthetic protein target(s), we screened a 350,000 compound library for competition with 1-aminoanthracene-apoferritin binding. Hit compounds meeting structural criteria had their binding affinities for apoferritin quantified with isothermal titration calorimetry and were tested for γ-aminobutyric acid type A receptor binding using a flunitrazepam binding assay. Chemotypes with a strong presence in the top 700 and exhibiting activity via isothermal titration calorimetry were selected for medicinal chemistry optimization including testing for anesthetic potency and toxicity in an in vivo Xenopus laevis tadpole assay. Compounds with low toxicity and high potency were tested for anesthetic potency in mice.

RESULTS

From an initial chemical library of more than 350,000 compounds, we identified 2,600 compounds that potently inhibited 1-aminoanthracene binding to apoferritin. A subset of compounds chosen by structural criteria (700) was successfully reconfirmed using the initial assay. Based on a strong presence in both the initial and secondary screens the 6-phenylpyridazin-3(2H)-one chemotype was assessed for anesthetic activity in tadpoles. Medicinal chemistry efforts identified four compounds with high potency and low toxicity in tadpoles, two were found to be effective novel anesthetics in mice.

CONCLUSION

The authors demonstrate the first use of a high-throughput screen to successfully identify a novel anesthetic chemotype and show mammalian anesthetic activity for members of that chemotype.