Putative Binding Sites,and Pathways to Them,for Amidine and Guanidine Current Inhibitors on Acid‐Sensing Ion Channels (ASIC). A Theoretical Approach with hASIC1a Homology Model

Central inhibition of the acid‐sensing hASIC1a channel, acting upstream of the opiate system, might serve to treat any type of pain, avoiding the unwanted addiction problems of the opioid drugs. To this end, inhibition of hASIC1a channel by PcTx1, a peptide from the Trinidad chevron tarantula, is under development. New inhibitors of the hASIC1a channel are also being sought, in the hope of further modulating the activity, from which antiplasmodial amidine and guanidine phenyl drugs have emerged as promising candidates. However, how such current inhibition takes place remains obscure from the molecular point of view, hindering any further progress in developing drugs. Therefore, the nature of the binding sites, and how they are reached by the amidine‐guanidine drugs, was investigated here via automated docking and molecular dynamics with hASIC1a homology models. This study has revealed that this ion channel is rich in binding sites, and that flexible drugs, such as nafamostat, may penetrate it in a snake‐like elongated conformation. Then, crawling like a snake through temporary holes in the protein, nafamostat either simply flips, or changes to a high‐energy folded conformation to become adapted to the shape of the binding site.     

Computational elucidation,mutational and hot spot-based designing of potential inhibitors against human acid-sensing ion channels (hASIC-1a) to treat various physiological conditions

Acid-sensing ion channels are ligand/proton-gated ion channels belonging to the family of the degenerin/epithelial Na⁺ channel (DEG/ENaC). They function as a sodium-selective pore for Ca²⁺ entry into neuronal cells during pathological conditions. The blocking of this channel has therapeutic importance, because at basal physiological pH (7.2), it is in a closed state and under a more acidic condition, and the ASIC1a ion channel is activated. To investigate the different states of the hASIC1a channel based on mutational analysis, structure-based virtual screening and molecular dynamics simulation studies. The system showed stability after 30 ns (after 1500 frame), and it was stabilized to an average value around 2.2Å. During the simulation, the ion channel residues in persistent contact with toxin PcTx1 were D237, E238, D347, D351, E219 and E355. These residues are important physiologically for the activation of the channel. From in silico alanine scanning, the significant hotspots obtained in hASIC1 are E344, P347, F352, D351, E355 and E219. From the sitemap analysis, it was evident that the sitemap found one of the active sites at the PcTx1 binding site with a site score of 1.086 and a D-score of 1.035 for hASIC1. We obtained a few promising hits and final potential hits from the virtual screening in hASIC1 that made interactions with the residues in the acidic pocket (E344, P347, F352, D351, E355 and E219). Based on these studies, the hits and scaffolds of potential therapeutic interest against various pathological conditions are associated with hASIC1a for future studies.     

On the putative binding site of RFamide-family neuropeptides from the western Atlantic clam Sunray Venus and cephalopods on acid-sensing ion channels. An automated docking and molecular-dynamics study with hASIC1a homology model

Investigated here are interactions of C-terminal amidated peptides with the hASIC1a acid-sensing ion channel. The peptides comprise endogenous FMRFa, present in the western Atlantic clam Sunray Venus, and FIRFa, present in cephalopods, as well as non-endogenous ones for comparison. The interaction is investigated by automated docking. The resulting key hASIC1a-FMRFa complex, set in a lipidic POPC (=1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membrane surrounded by H(2) O and Na(+) -neutralized, was also investigated by molecular dynamics. It was observed that all investigated peptides become encapsulated into the ion channel, on one side by the thumb and finger of a subunit, and, on the opposite side, by the knuckle and β-ball of a second subunit. The third subunit is not involved. This is much the same binding site that was disclosed previously by both a similar computational approach, and electrophysiological and binding experiments for the hASIC1a ion channel-blocker tarantula toxin PCTX1. This paves the way to a better understanding of the role of these peptides in invertebrates.     

A flexible GAS belt responds to pore mutations changing the ion selectivity of proton-gated channels

Proton-gated ion channels conduct mainly Na⁺ to induce postsynaptic membrane depolarization. Finding the determinants of ion selectivity requires knowledge of the pore structure in the open conformation, but such information is not yet available. Here, the open conformation of the hASIC1a channel was computationally modeled, and functional effects of pore mutations were analyzed in light of the predicted structures. The open pore structure shows two constrictions of similar diameter formed by the backbone of the GAS belt and, right beneath it, by the side chains of H28 from the reentrant loop. Models of nonselective mutant channels, but not those that maintain ion selectivity, predict enlargement of the GAS belt, suggesting that this motif is quite flexible and that the loss of stabilizing interactions in the central pore leads to changes in size/shape of the belt. Our results are consistent with the “close-fit” mechanism governing selectivity of hASIC1a, wherein the backbone of the GAS substitutes at least part of the hydration shell of a permeant ion to enable crossing the pore constriction.     

Potentiation of acid-sensing ion channels by sulfhydryl compounds

The acid-sensing ion channels (ASICs) are voltage-independent ion channels activated by acidic extracellular pH. ASICs play a role in sensory transduction, behavior, and acidotoxic neuronal death, which occurs during stroke and ischemia. During these conditions, the extracellular concentration of sulfhydryl reducing agents increases. We used perforated patch-clamp technique to analyze the impact of sulfhydryls on H⁺-gated currents from Chinese hamster ovary (CHO) cells expressing human ASIC1a (hASIC1a). We found that hASIC1a currents activated by pH 6.5 were increased almost twofold by the sulfhydryl-containing reducing agents dithiothreitol (DTT) and glutathione. DTT shifted the pH-dose response of hASIC1a toward a more neutral pH (pH_0.5 from 6.54 to 6.69) and slowed channel desensitization. The effect of reducing agents on native mouse hippocampal neurons and transfected mouse ASIC1a was similar. We found that the effect of DTT on hASIC1a was mimicked by the metal chelator TPEN, and mutant hASIC1a channels with reduced TPEN potentiation showed reduced DTT potentiation. Furthermore, the addition of DTT in the presence of TPEN did not result in further increases in current amplitude. These results suggest that the effect of DTT on hASIC1a is due to relief of tonic inhibition by transition metal ions. We found that all ASICs examined remained potentiated following the removal of DTT. This effect was reversed by the oxidizing agent DTNB in hASIC1a, supporting the hypothesis that DTT also impacts ASICs via a redox-sensitive site. Thus sulfhydryl compounds potentiate H⁺-gated currents via two mechanisms, metal chelation and redox modulation of target amino acids. glutathione; DTT; redox; zinc     

An arginine residue in the outer segment of hASIC1a TM1 affects both proton affinity and channel desensitization

Acid-sensing ion channels (ASICs) respond to changes in pH in the central and peripheral nervous systems and participate in synaptic plasticity and pain perception. Understanding the proton-mediated gating mechanism remains elusive despite the of their structures in various conformational states. We report here that R64, an arginine located in the outer segment of the first transmembrane domain of all three isoforms of mammalian ASICs, markedly impacts the apparent proton affinity of activation and the degree of desensitization from the open and preopen states. Rosetta calculations of free energy changes predict that substitutions of R64 in hASIC1a by aromatic residues destabilize the closed conformation while stabilizing the open conformation. Accordingly, F64 enhances the efficacy of proton-mediated gating of hASIC1a, which increases the apparent pH₅₀ and facilitates channel opening when only one or two subunits are activated. F64 also lengthens the duration of opening events, thus keeping channels open for extended periods of time and diminishing low pH-induced desensitization. Our results indicate that activation of a proton sensor(s) with pH₅₀ equal to or greater than pH 7.2–7.1 opens F64hASIC1a, whereas it induces steady-state desensitization in wildtype channels due to the high energy of activation imposed by R64, which prevents opening of the pore. Together, these findings suggest that activation of a high-affinity proton-sensor(s) and a common gating mechanism may mediate the processes of activation and steady-state desensitization of hASIC1a.     

High-throughput characterization of photocrosslinker-bearing ion channel variants to map residues critical for function and pharmacology

Incorporation of noncanonical amino acids (ncAAs) can endow proteins with novel functionalities, such as crosslinking or fluorescence. In ion channels, the function of these variants can be studied with great precision using standard electrophysiology, but this approach is typically labor intensive and low throughput. Here, we establish a high-throughput protocol to conduct functional and pharmacological investigations of ncAA-containing human acid-sensing ion channel 1a (hASIC1a) variants in transiently transfected mammalian cells. We introduce 3 different photocrosslinking ncAAs into 103 positions and assess the function of the resulting 309 variants with automated patch clamp (APC). We demonstrate that the approach is efficient and versatile, as it is amenable to assessing even complex pharmacological modulation by peptides. The data show that the acidic pocket is a major determinant for current decay, and live-cell crosslinking provides insight into the hASIC1a–psalmotoxin 1 (PcTx1) interaction. Further, we provide evidence that the protocol can be applied to other ion channels, such as P2X2 and GluA2 receptors. We therefore anticipate the approach to enable future APC-based studies of ncAA-containing ion channels in mammalian cells.

This study describes a method to rapidly screen hundreds of ion channel variants containing non-canonical amino acids. A proof-of-principle introducing photocrosslinking non-canonical amino acids into the human ion channel hASIC1a shows how this approach can provide insights into function and pharmacology.     

Substrate recognition and translocation by polyspecific organic cation transporters

Organic cation transporters (OCTs) of the SLC22 family play a pivotal role in distribution and excretion of cationic drugs. They mediate electrogenic translocation of cations in both directions. OCTs are polyspecific transporters. During substrate translocation they perform a series of conformational changes involving an outward-facing conformation, an occluded state and an inward-facing conformation. Mutagenesis of OCT1 in combination with homology modeling showed that identical amino acids form the innermost parts of the outward-open and inward-open binding clefts. In addition to low affinity substrate binding sites, OCT1 contains high affinity substrate binding sites that can mediate inhibition via non-transported compounds.     

Calcitonin gene-related peptide receptor independently stimulates 3′,5′-cyclic adenosine monophosphate and Ca2+ signaling pathways

The inhibition of ion transporting ATPases (Na⁺,K⁺-ATPase, Ca²⁺,Mg²⁺- and Ca²⁺-ATPase) by two amphiphilic drugs e.g. chlorpromazine (antipsychotic) and chloroquine (antimalarial) are found to be competitive in nature in vitro with respect to the substrate. Two binding sites - high and low affinity are found to exist on all the three ATPases toward these drugs as evident from the plot of F/F₀ vs. different drug concentrations of tryptophan fluorescence of the enzymes. Circular dichroism analysis suggest that binding of drugs to the high affinity site does not involve any change in conformation of ATPase molecules which occur only when drug binds to the low affinity sites. The drug binding sites and possible effect on conformational change of ATPase molecules of these two drugs have been described in this report.     

General anesthetics predicted to block the GLIC pore with micromolar affinity

Although general anesthetics are known to modulate the activity of ligand-gated ion channels in the Cys-loop superfamily, there is at present neither consensus on the underlying mechanisms, nor predictive models of this modulation. Viable models need to offer quantitative assessment of the relative importance of several identified anesthetic binding sites. However, to date, precise affinity data for individual sites has been challenging to obtain by biophysical means. Here, the likely role of pore block inhibition by the general anesthetics isoflurane and propofol of the prokaryotic pentameric channel GLIC is investigated by molecular simulations. Microscopic affinities are calculated for both single and double occupancy binding of isoflurane and propofol to the GLIC pore. Computations are carried out for an open-pore conformation in which the pore is restrained to crystallographic radius, and a closed-pore conformation that results from unrestrained molecular dynamics equilibration of the structure. The GLIC pore is predicted to be blocked at the micromolar concentrations for which inhibition by isofluorane and propofol is observed experimentally. Calculated affinities suggest that pore block by propofol occurs at signifcantly lower concentrations than those for which inhibition is observed: we argue that this discrepancy may result from binding of propofol to an allosteric site recently identified by X-ray crystallography, which may cause a competing gain-of-function effect. Affinities of isoflurane and propofol to the allosteric site are also calculated, and shown to be 3 mM for isoflurane and 10 μM for propofol; both anesthetics have a lower affinity for the allosteric site than for the unoccupied pore.     

Anesthetic Binding in a Pentameric Ligand-Gated Ion Channel: GLIC

Cys-loop receptors are molecular targets of general anesthetics, but the knowledge of anesthetic binding to these proteins remains limited. Here we investigate anesthetic binding to the bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC), a structural homolog of cys-loop receptors, using an experimental and computational hybrid approach. Tryptophan fluorescence quenching experiments showed halothane and thiopental binding at three tryptophan-associated sites in the extracellular (EC) domain, transmembrane (TM) domain, and EC-TM interface of GLIC. An additional binding site at the EC-TM interface was predicted by docking analysis and validated by quenching experiments on the N200W GLIC mutant. The binding affinities (K_D) of 2.3 ± 0.1 mM and 0.10 ± 0.01 mM were derived from the fluorescence quenching data of halothane and thiopental, respectively. Docking these anesthetics to the original GLIC crystal structure and the structures relaxed by molecular dynamics simulations revealed intrasubunit sites for most halothane binding and intersubunit sites for thiopental binding. Tryptophans were within reach of both intra- and intersubunit binding sites. Multiple molecular dynamics simulations on GLIC in the presence of halothane at different sites suggested that anesthetic binding at the EC-TM interface disrupted the critical interactions for channel gating, altered motion of the TM23 linker, and destabilized the open-channel conformation that can lead to inhibition of GLIC channel current. The study has not only provided insights into anesthetic binding in GLIC, but also demonstrated a successful fusion of experiments and computations for understanding anesthetic actions in complex proteins.     

Modelling topoisomerase I inhibition by minor groove binders

Topoisomerase inhibition is an extremely useful target for anticancer and antimicrobial drugs, and an undesirable side effect of some drugs targeting other proteins. Published modelling studies are sparse, and have used small data sets with relatively low molecular diversity. Given the important role of minor groove binding in the mechanism of topoisomerase I inhibition, we have conducted the first 3D QSAR study of topoisomerase I inhibition of a large, diverse set of minor groove binders using the minor groove binding conformation as the alignment template. The highly significant QSAR models resulting from this alignment identify the roles played by molecular features, most importantly the hydrogen bond donor properties.     

Affinity of calcium channel inhibitors, benzodiazepines, and other vasoactive compounds for the nucleoside transport system

There is evidence to suggest that several different groups of drugs including the so-called coronary vasodilators, benzodiazepines, and calcium channel inhibitors may owe their vasoactivity, in part, to the potentiation of the vasorelaxant effects of endogenous adenosine. To measure the affinity of some of these agents for the membrane-located nucleoside transport system, competition binding assays have been performed using the high-affinity radioligand [3H]nitrobenzylthioinosine (NBMPR). Experiments were performed on human erythrocytes and cardiac membranes from guinea pigs and rats. Recognized nucleoside transport inhibitors had high affinity (less than 50 nM) for NBMPR recognition sites associated with the nucleoside transporter complex in human erythrocytes, whereas calcium channel inhibitors and benzodiazepines had predominantly low affinity (greater than 1 microM). Although some recognized transport inhibitors, such as dipyridamole, show marked differences in affinity for NBMPR sites in guinea pig and rat tissues, benzodiazepines and calcium channel blockers displayed no such species selectivity and had low affinity (greater than 1 microM) for NBMPR sites in both guinea pig and rat cardiac membranes. Consequently, it is unlikely that agents such as benzodiazepines and calcium channel inhibitors cause significant inhibition of adenosine transport, and hence potentiate adenosine actions, at the concentrations required to induce effects through occupation of their respective, specific high-affinity sites.     

Exploring Volatile General Anesthetic Binding to a Closed Membrane-Bound Bacterial Voltage-Gated Sodium Channel via Computation

Despite the clinical ubiquity of anesthesia, the molecular basis of anesthetic action is poorly understood. Amongst the many molecular targets proposed to contribute to anesthetic effects, the voltage gated sodium channels (VGSCs) should also be considered relevant, as they have been shown to be sensitive to all general anesthetics tested thus far. However, binding sites for VGSCs have not been identified. Moreover, the mechanism of inhibition is still largely unknown. The recently reported atomic structures of several members of the bacterial VGSC family offer the opportunity to shed light on the mechanism of action of anesthetics on these important ion channels. To this end, we have performed a molecular dynamics “flooding” simulation on a membrane-bound structural model of the archetypal bacterial VGSC, NaChBac in a closed pore conformation. This computation allowed us to identify binding sites and access pathways for the commonly used volatile general anesthetic, isoflurane. Three sites have been characterized with binding affinities in a physiologically relevant range. Interestingly, one of the most favorable sites is in the pore of the channel, suggesting that the binding sites of local and general anesthetics may overlap. Surprisingly, even though the activation gate of the channel is closed, and therefore the pore and the aqueous compartment at the intracellular side are disconnected, we observe binding of isoflurane in the central cavity. Several sampled association and dissociation events in the central cavity provide consistent support to the hypothesis that the “fenestrations” present in the membrane-embedded region of the channel act as the long-hypothesized hydrophobic drug access pathway.     

Two neighboring residues of loop A of the alpha1 subunit point towards the benzodiazepine binding site of GABAA receptors

Benzodiazepines are widely used drugs exerting sedative, anxiolytic, muscle relaxant, and anticonvulsant effects by acting through specific high affinity binding sites on some GABA(A) receptors. It is important to understand how these ligands are positioned in this binding site. We are especially interested here in the conformation of loop A of the alpha(1)beta(2)gamma(2) GABA(A) receptor containing a key residue for the interaction of benzodiazepines: alpha(1)H101. We describe a direct interaction of alpha(1)N102 with a diazepam- and an imidazobenzodiazepine-derivative. Our observations help to better understand the conformation of this region of the benzodiazepine pocket in GABA(A) receptor.     

Molecular Cloning and Regional Distribution of a Human Proton Receptor Subunit with Biphasic Functional Properties

Abstract : Small changes of extracellular pH activate depolarizing inward currents in most nociceptive neurons. It has been recently proposed that acid sensitivity of sensory as well as central neurons is mediated by a family of proton-gated cation channels structurally related to Caenorhabditis elegans degenerins and mammalian epithelial sodium channels. We describe here the molecular cloning of a novel human proton receptor, hASIC3, a 531-amino acid-long subunit homologous to rat DRASIC. Expression of homomeric hASIC3 channels in Xenopus oocytes generated biphasic inward currents elicited at pH <5, providing the first functional evidence of a human proton-gated ion channel. Contrary to the DRASIC current phenotype, the fast desensitizing early component and the slow sustained late component differed both by their cationic selectivity and by their response to the antagonist amiloride, but not by their pH sensitivity (pH₅₀ = 3.66 vs. 3.82). Using RT-PCR and mRNA blot hybridization, we detected hASIC3 mRNA in sensory ganglia, brain, and many internal tissues including lung and testis, so hASIC3 gene expression was not restricted to peripheral sensory neurons. These functional and anatomical data strongly suggest that hASIC3 plays a major role in persistent proton-induced currents occurring in physiological and pathological conditions of pH changes, likely through a tissue-specific heteropolymerization with other members of the proton-gated channel family.     

Structural basis for topoisomerase I inhibition by nucleoside analogs

Nucleoside analogs such as 1-beta-D-arabinofuranosyl cytidine (AraC) and 2',2'-difluoro deoxycytidine (dFdC) are important components of the anticancer chemotherapeutic arsenal and are among the most effective anticancer drugs currently available. Although both AraCTP and dFdCTP impede DNA replication through pausing of DNA polymerases, both nucleoside analogs are ultimately incorporated into replicated DNA and interfere in DNA-mediated processes. Our laboratories are investigating the structural basis for the poisoning of topoisomerase I (top1) due to antipyrimidine incorporation into duplex DNA. We recently reported that both AraC and dFdC induce formation of top1 cleavage complexes, and poisoning of top1 contributes to the anticancer activities of both these drugs. Recent NMR and thermodynamic studies from our laboratories provide insight into the mechanism by which AraC and dFdC poison top1. NMR studies from our laboratories have revealed that the arabinosyl sugar of AraC adopted a C2'-endo conformation. Although this is a B-type sugar pucker characteristic of duplex DNA, the conformation is rigid, and this lack of flexibility probably contributes to inhibition of the religation step of the top1 reaction. In contrast to AraC, NMR studies revealed dFdC adopted a C3' endo sugar pucker characteristic of RNA, rather than DNA duplexes. dFdC substitution enhanced formation of top1 cleavage complexes, but did not inhibit religation. The enhancement of top1 cleavage complexes most likely results from a combination of conformational and electrostatic effects. The structural effects of dFdC and AraC are being further investigated in duplex DNA with well-defined top1 cleavage sites to analyze more specifically how these structural perturbations lead to enzyme poisoning.     

Molecular mapping of general anesthetic sites in a voltage-gated ion channel

Several voltage-gated ion channels are modulated by clinically relevant doses of general anesthetics. However, the structural basis of this modulation is not well understood. Previous work suggested that n-alcohols and inhaled anesthetics stabilize the closed state of the Shaw2 voltage-gated (Kv) channel (K-Shaw2) by directly interacting with a discrete channel site. We hypothesize that the inhibition of K-Shaw2 channels by general anesthetics is governed by interactions between binding and effector sites involving components of the channel's activation gate. To investigate this hypothesis, we applied Ala/Val scanning mutagenesis to the S4-S5 linker and the post-PVP S6 segment, and conducted electrophysiological analysis to evaluate the energetic impact of the mutations on the inhibition of the K-Shaw2 channel by 1-butanol and halothane. These analyses identified residues that determine an apparent binding cooperativity and residue pairs that act in concert to modulate gating upon anesthetic binding. In some instances, due to their critical location, key residues also influence channel gating. Complementing these results, molecular dynamics simulations and in silico docking experiments helped us visualize possible anesthetic sites and interactions. We conclude that the inhibition of K-Shaw2 by general anesthetics results from allosteric interactions between distinct but contiguous binding and effector sites involving inter- and intrasubunit interfaces.