Size matters in activation/inhibition of ligand-gated ion channels

Structural Biology-Guided Design, Synthesis, and Biological Evaluation of Novel Insect Nicotinic Acetylcholine Receptor Orthosteric Modulators

The development of novel and safe insecticides remains an important need for a growing world population to protect crops and animal and human health. New chemotypes modulating the insect nicotinic acetylcholine receptors have been recently brought to the agricultural market, yet with limited understanding of their molecular interactions at their target receptor. Herein, we disclose the first crystal structures of these insecticides, namely, sulfoxaflor, flupyradifurone, triflumezopyrim, flupyrimin, and the experimental compound, dicloromezotiaz, in a double-mutated acetylcholine-binding protein which mimics the insect-ion-channel orthosteric site. Enabled by these findings, we discovered novel pharmacophores with a related mode of action, and we describe herein their design, synthesis, and biological evaluation.

Signalling assemblies: the odds of symmetry

The assembly of proteins into complexes is fundamental to nearly all biological signalling processes. Symmetry is a dominant feature of the structures of experimentally determined protein complexes, observed in the vast majority of homomers and many heteromers. However, some asymmetric structures exist, and asymmetry also often forms transiently, intractable to traditional structure determination methods. Here, we explore the role of protein complex symmetry and asymmetry in cellular signalling, focusing on receptors, transcription factors and transmembrane channels, among other signalling assemblies. We highlight a recurrent tendency for asymmetry to be crucial for signalling function, often being associated with activated states. We conclude with a discussion of how consideration of protein complex symmetry and asymmetry has significant potential implications and applications for pharmacology and human disease.

Gating Motions and Stationary Gating Properties of Ionotropic Glutamate Receptors: Computation Meets Electrophysiology

Ionotropic glutamate receptors (iGluRs) are tetrameric ligand-gated ion channels essential to all aspects of brain function, including higher order processes such as learning and memory. For decades, electrophysiology was the primary means for characterizing the function of iGluRs and gaining mechanistic insight. Since the turn of the century, structures of isolated water-soluble domains and transmembrane-domain-containing constructs have provided the basis for formulating mechanistic hypotheses. Because these structures only represent sparse, often incomplete snapshots during iGluR activation, significant gaps in knowledge remain regarding structures, energetics, and dynamics of key substates along the functional processes. Some of these gaps have recently been filled by molecular dynamics simulations and theoretical modeling. In this Account, I describe our work in the latter arena toward characterizing iGluR gating motions and developing a formalism for calculating thermodynamic and kinetic properties of stationary gating. The structures of iGluR subunits have a highly modular architecture, in which the ligand-binding domain and the transmembrane domain are well separated and connected by flexible linkers. The ligand-binding domain in turn is composed of two subdomains. During activation, agonist binding induces the closure of the intersubdomain cleft. The cleft closure leads to the outward pulling of a linker tethered to the extracellular terminus of the major pore-lining helix of the transmembrane domain, thereby opening the channel. This activation model based on molecular dynamics simulations was validated by residue-specific information from electrophysiological data on cysteine mutants. A further critical test was made through introducing glycine insertions in the linker. Molecular dynamics simulations showed that, with lengthening by glycine insertions, the linker became less effective in pulling the pore-lining helix, leading to weaker stabilization of the channel-open state. In full agreement, single-channel recordings showed that the channel open probability decreased progressively as the linker was lengthened by glycine insertions. Crystal structures of ligand-binding domains showing different degrees of cleft closure between full and partial agonists suggested a simple mechanism for one subtype of iGluRs, but mysteries surrounded a second subtype, where the ligand-binding domains open to similar degrees when bound with either full or partial agonists. Our free energy simulations now suggest that broadening of the free energy basin for cleft closure is a plausible solution. A theoretical basis for these mechanistic hypotheses on partial agonisms was provided by a model for the free energy surface of a full receptor, where the stabilization by cleft closure is transmitted via the linker to the channel-open state. This model can be implemented by molecular dynamics simulations to predict thermodynamic and kinetics properties of stationary gating that are amenable to direct test by single-channel recordings. Close integration between computation and electrophysiology holds great promises in revealing the conformations of key substates in functional processes and the mechanisms of disease-associated mutations.

Molecular blueprint of allosteric binding sites in a homologue of the agonist-binding domain of the α7 nicotinic acetylcholine receptor

The α7 nicotinic acetylcholine receptor (nAChR) belongs to the family of pentameric ligand-gated ion channels and is involved in fast synaptic signaling. In this study, we take advantage of a recently identified chimera of the extracellular domain of the native α7 nicotinic acetylcholine receptor and acetylcholine binding protein, termed α7-AChBP. This chimeric receptor was used to conduct an innovative fragment-library screening in combination with X-ray crystallography to identify allosteric binding sites. One allosteric site is surface-exposed and is located near the N-terminal α-helix of the extracellular domain. Ligand binding at this site causes a conformational change of the α-helix as the fragment wedges between the α-helix and a loop homologous to the main immunogenic region of the muscle α1 subunit. A second site is located in the vestibule of the receptor, in a preexisting intrasubunit pocket opposite the agonist binding site and corresponds to a previously identified site involved in positive allosteric modulation of the bacterial homolog ELIC. A third site is located at a pocket right below the agonist binding site. Using electrophysiological recordings on the human α7 nAChR we demonstrate that the identified fragments, which bind at these sites, can modulate receptor activation. This work presents a structural framework for different allosteric binding sites in the α7 nAChR and paves the way for future development of novel allosteric modulators with therapeutic potential.

Mechanism-Based Mathematical Model for Gating of Ionotropic Glutamate Receptors

We present a mathematical model for ionotropic glutamate receptors (iGluR's) that is built on mechanistic understanding and yields a number of thermodynamic and kinetic properties of channel gating. iGluR's are ligand-gated ion channels responsible for the vast majority of fast excitatory neurotransmission in the central nervous system. The effects of agonist-induced closure of the ligand-binding domain (LBD) are transmitted to the transmembrane channel (TMC) via interdomain linkers. Our model demonstrates that, relative to full agonists, partial agonists may reduce either the degree of LBD closure or the curvature of the LBD free energy basin, leading to less stabilization of the channel open state and hence lower channel open probability. A rigorous relation is derived between the channel closed-to-open free energy difference and the tension within the linker. Finally, by treating LBD closure and TMC opening as diffusive motions, we obtain gating trajectories that resemble stochastic current traces from single-channel recordings and calculate the rate constants for transitions between the channel open and closed states. Our model can be implemented by molecular dynamics simulations to realistically depict iGluR gating and may guide functional experiments in gaining deeper insight into this essential family of channel proteins.

Reduced curvature of ligand-binding domain free-energy surface underlies partial agonism at NMDA receptors

NMDA receptors are ligand-gated ion channels that mediate excitatory synaptic transmission in the central nervous system. Partial agonists elicit submaximal channel activation, but crystal structures of the ligand-binding domains (LBDs) bound with partial and full agonists show little difference. To uncover the molecular mechanism for partial agonism, here we computed the free-energy surfaces of the GluN1 (an obligatory subunit of NMDA receptors) LBD bound with a variety of ligands. The free-energy minima are similarly positioned for full and partial agonists, but the curvatures are significantly reduced in the latter case, indicating higher probabilities for sampling conformations with a not fully closed domain cleft. The free-energy surfaces for antagonists have both shifted minima and further reduced curvatures. Reduced curvature of free-energy surface appears to explain well the partial agonism at NMDA receptors and may present a unique paradigm in producing graded responses for receptors in general.

Study, by use of coarse-grained models, of the functionally crucial residues and allosteric pathway of anesthetic regulation of the Gloeobacter violaceus ligand-gated ion channel

Although pentameric ligand-gated ion channels (pLGICs) have been found to be the targets of general anesthetics, the mechanism of the effects of anesthetics on pLGICs remains elusive. pLGICs from Gloeobacter violaceus (GLIC) can be inhibited by the anesthetic ketamine. X-ray crystallography has shown that the ketamine binding site is distant from the channel gate of the GLIC. It is still not clear how ketamine controls the function of the GLIC by long-range allosteric regulation. In this work, the functionally crucial residues and allosteric pathway of anesthetic regulation of the GLIC were identified by use of a coarse-grained thermodynamic method developed by our group. In our method, the functionally crucial sites were identified as the residues thermodynamically coupled with binding of ketamine. The results from calculation were highly consistent with experimental data. Our study aids understanding of the mechanism of the anesthetic action of ketamine on the GLIC by long-range allosteric modulation.

Comparative pharmacology of flatworm and roundworm glutamate-gated chloride channels: Implications for potential anthelmintics

•

Flatworm and roundworm glutamate-gated chloride channels (GluCls) were compared.

•

Several glutamate analogues activated both GluCls in the millimolar range.

•

Quisqualate selectively activated the flatworm GluCl.

•

Propofol and thymol inhibited both GluCls in the micromolar range.

Pharmacological targeting of glutamate-gated chloride channels (GluCls) is a potent anthelmintic strategy, evidenced by macrocyclic lactones that eliminate numerous roundworm infections by activating roundworm GluCls. Given the recent identification of flatworm GluCls and the urgent need for drugs against schistosomiasis, flatworm GluCls should be evaluated as potential anthelmintic targets. This study sought to identify agonists or modulators of one such GluCl, SmGluCl-2 from the parasitic flatworm Schistosoma mansoni. The effects of nine glutamate-like compounds and three monoterpenoid ion channel modulators were measured by electrophysiology at SmGluCl-2 recombinantly expressed in Xenopus laevis oocytes. For comparison with an established anthelmintic target, experiments were also performed on the AVR-14B GluCl from the parasitic roundworm Haemonchus contortus. l-Glutamate was the most potent agonist at both GluCls, but l-2-aminoadipate, d-glutamate and d-2-aminoadipate activated SmGluCl-2 (EC₅₀ 1.0 ± 0.1 mM, 2.4 ± 0.4 mM, 3.6 ± 0.7 mM, respectively) more potently than AVR-14B. Quisqualate activated only SmGluCl-2 whereas l-aspartate activated only AVR-14B GluCls. Regarding the monoterpenoids, both GluCls were inhibited by propofol, thymol and menthol, SmGluCl-2 most potently by thymol (IC₅₀ 484 ± 85 μM) and least potently by menthol (IC₅₀ > 3 mM). Computational docking suggested that agonist and inhibitor potency is attributable to particular interactions with extracellular or membrane-spanning amino acid residues. These results reveal that flatworm GluCls are pharmacologically susceptible to numerous agonists and modulators and indicate that changes to the glutamate γ-carboxyl or to the propofol 6-isopropyl group can alter the differential pharmacology at flatworm and roundworm GluCls. This should inform the development of more potent compounds and in turn lead to novel anthelmintics.

Involvement of ectodomain Leu 214 in ATP binding and channel desensitization of the P2X4 receptor

P2X receptors are trimeric ATP-gated cation permeable ion channels. When ATP binds, the extracellular head and dorsal fin domains are predicted to move closer to each other. However, there are scant functional data corroborating the role of the dorsal fin in ligand binding. Here using site-directed mutagenesis and electrophysiology, we show that a dorsal fin leucine, L214, contributes to ATP binding. Mutant receptors containing a single substitution of alanine, serine, glutamic acid, or phenylalanine at L214 of the rat P2X4 receptor exhibited markedly reduced sensitivities to ATP. Mutation of other dorsal fin side chains, S216, T223, and D224, did not significantly alter ATP sensitivity. Exposure of L214C to sodium (2-sulfonatoethyl) methanethiosulfonate (MTSES(-)) or (2-aminoethyl) methanethiosulfonate hydrobromide in the absence of ATP blocked responses evoked by subsequent ATP application. In contrast, when MTSES(-) was applied in the presence of ATP, no current inhibition was observed. Furthermore, L214A also slightly reduced the inhibitory effect of the antagonist 2',3'-O-(2,4,6-trinitrophenyl)-ATP, and the blockade was more rapidly reversible after washout. Certain L214 mutants also showed effects on current desensitization in the continued presence of ATP. L214I exhibited an accelerated current decline, whereas L214M exhibited a slower rate. Taken together, these data reveal that position L214 participates in both ATP binding and conformational changes accompanying channel opening and desensitization, providing compelling evidence that the dorsal fin domain indeed has functional properties that are similar to those previously reported for the body domains.

Differential modulation of GABA(A) receptor function by aryl pyrazoles

Several aryl pyrazoles characterized by a different molecular structure (flexible vs constrained), but chemically related to rimonabant and AM251, were tested for their ability to modulate the function of recombinant α1β2γ2L GABAA receptors expressed in Xenopus laevis oocytes. The effects of 6Bio-R, 14Bio-R, NESS 0327, GP1a and GP2a (0.3-30 μM) were evaluated using a two-electrode voltage-clamp technique. 6Bio-R and 14Bio-R potentiated GABA-evoked Cl(-) currents. NESS 0327, GP1a and GP2a did not affect the GABAA receptor function, but they acted as antagonists of 6Bio-R. Moreover, NESS 0327 inhibited the potentiation of the GABAA receptor function induced by rimonabant. The benzodiazepine site seems to participate in the action of these compounds. In fact, flumazenil antagonized the potentiation of the GABAA receptor induced by 6Bio-R, and NESS 0327 reduced the action of lorazepam and zolpidem. On the contrary, NESS 0327 did not antagonize the action of "classic" GABAergic modulators (propanol, anesthetics, barbiturates or steroids). In α1β2 receptors 6Bio-R potentiated the GABAergic function, but flumazenil was still able to antagonize the potentiation induced by 6Bio-R. Aryl pyrazole derivatives activity at the GABAA receptor depends on their molecular structure. These compounds bind to both an αβγ binding site, and to an α/β site which do not require the γ subunit and that may provide structural leads for drugs with potential anticonvulsant effects.

Molecular recognition of the neurotransmitter acetylcholine by an acetylcholine binding protein reveals determinants of binding to nicotinic acetylcholine receptors

Despite extensive studies on nicotinic acetylcholine receptors (nAChRs) and homologues, details of acetylcholine binding are not completely resolved. Here, we report the crystal structure of acetylcholine bound to the receptor homologue acetylcholine binding protein from Lymnaea stagnalis. This is the first structure of acetylcholine in a binding pocket containing all five aromatic residues conserved in all mammalian nAChRs. The ligand-protein interactions are characterized by contacts to the aromatic box formed primarily by residues on the principal side of the intersubunit binding interface (residues Tyr89, Trp143 and Tyr185). Besides these interactions on the principal side, we observe a cation-π interaction between acetylcholine and Trp53 on the complementary side and a water-mediated hydrogen bond from acetylcholine to backbone atoms of Leu102 and Met114, both of importance for anchoring acetylcholine to the complementary side. To further study the role of Trp53, we mutated the corresponding tryptophan in the two different acetylcholine-binding interfaces of the widespread α4β2 nAChR, i.e. the interfaces α4(+)β2(−) and α4(+)α4(−). Mutation to alanine (W82A on the β2 subunit or W88A on the α4 subunit) significantly altered the response to acetylcholine measured by oocyte voltage-clamp electrophysiology in both interfaces. This shows that the conserved tryptophan residue is important for the effects of ACh at α4β2 nAChRs, as also indicated by the crystal structure. The results add important details to the understanding of how this neurotransmitter exerts its action and improves the foundation for rational drug design targeting these receptors.

Heteromeric assembly of P2X subunits

Transcripts and/or proteins of P2X receptor (P2XR) subunits have been found in virtually all mammalian tissues. Generally more than one of the seven known P2X subunits have been identified in a given cell type. Six of the seven cloned P2X subunits can efficiently form functional homotrimeric ion channels in recombinant expression systems. This is in contrast to other ligand-gated ion channel families, such as the Cys-loop or glutamate receptors, where homomeric assemblies seem to represent the exception rather than the rule. P2XR mediated responses recorded from native tissues rarely match exactly the biophysical and pharmacological properties of heterologously expressed homomeric P2XRs. Heterotrimerization of P2X subunits is likely to account for this observed diversity. While the existence of heterotrimeric P2X2/3Rs and their role in physiological processes is well established, the composition of most other P2XR heteromers and/or the interplay between distinct trimeric receptor complexes in native tissues is not clear. After a description of P2XR assembly and the structure of the intersubunit ATP-binding site, this review summarizes the distribution of P2XR subunits in selected mammalian cell types and the biochemically and/or functionally characterized heteromeric P2XRs that have been observed upon heterologous co-expression of P2XR subunits. We further provide examples where the postulated heteromeric P2XRs have been suggested to occur in native tissues and an overview of the currently available pharmacological tools that have been used to discriminate between homo- and heteromeric P2XRs.

Moving through the gate in ATP-activated P2X receptors

P2X receptors are nonselective cation channels gated by extracellular ATP. They represent new therapeutic targets, and they form channels with a unique trimeric architecture. In 2009, the first crystal structure of a P2X receptor was reported, in which the receptor was in an ATP-free, closed channel state. However, our view recently changed when a second crystal structure was reported, in which a P2X receptor was bound to ATP and resolved in an open channel conformation. This remarkable structure not only confirms many key experimental data, including the recent mechanisms of ATP binding and ion permeation, but also reveals unanticipated mechanisms. Certainly, this new information will accelerate our understanding of P2X receptor function and pharmacology at the atomic level.