Stabilization of the GluCl ligand-gated ion channel in the presence and absence of ivermectin

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.

Microsecond-timescale simulations suggest 5-HT-mediated preactivation of the 5-HT3A serotonin receptor

Aided by efforts to improve their speed and efficiency, molecular dynamics (MD) simulations provide an increasingly powerful tool to study the structure-function relationship of pentameric ligand-gated ion channels (pLGICs). However, accurate reporting of the channel state and observation of allosteric regulation by agonist binding with MD remains difficult due to the timescales necessary to equilibrate pLGICs from their artificial and crystalized conformation to a more native, membrane-bound conformation in silico. Here, we perform multiple all-atom MD simulations of the homomeric 5-hydroxytryptamine 3A (5-HT3A) serotonin receptor for 15 to 20 μs to demonstrate that such timescales are critical to observe the equilibration of a pLGIC from its crystalized conformation to a membrane-bound conformation. These timescales, which are an order of magnitude longer than any previous simulation of 5-HT3A, allow us to observe the dynamic binding and unbinding of 5-hydroxytryptamine (5-HT) (i.e., serotonin) to the binding pocket located on the extracellular domain (ECD) and allosteric regulation of the transmembrane domain (TMD) from synergistic 5-HT binding. While these timescales are not long enough to observe complete activation of 5-HT3A, the allosteric regulation of ion gating elements by 5-HT binding is indicative of a preactive state, which provides insight into molecular mechanisms that regulate channel activation from a resting state. This mechanistic insight, enabled by microsecond-timescale MD simulations, will allow a careful examination of the regulation of pLGICs at a molecular level, expanding our understanding of their function and elucidating key structural motifs that can be targeted for therapeutic regulation.

Lipid-Dependent Regulation of Ion Channels and G Protein-Coupled Receptors: Insights from Structures and Simulations

Ion channels and G protein-coupled receptors (GPCRs) are regulated by lipids in their membrane environment. Structural studies combined with biophysical and molecular simulation investigations reveal interaction sites for specific lipids on membrane protein structures. For K channels, PIP₂ plays a key role in regulating Kv and Kir channels. Likewise, several recent cryo-EM structures of TRP channels have revealed bound lipids, including PIP₂ and cholesterol. Among the pentameric ligand-gated ion channel family, structural and biophysical studies suggest the M4 TM helix may act as a lipid sensor, e.g., forming part of the binding sites for neurosteroids on the GABA_A receptor. Structures of GPCRs have revealed multiple cholesterol sites, which may modulate both receptor dynamics and receptor oligomerization. PIP₂ also interacts with GPCRs and may modulate their interactions with G proteins. Overall, it is evident that multiple lipid binding sites exist on channels and receptors that modulate their function allosterically and are potential druggable sites.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid–protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid–protein interactions, and to link lipid–protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid–protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid–protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid–protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid–protein interactions.

Closed-Locked and Apo-Resting State Structures of the Human α7 Nicotinic Receptor: A Computational Study

Nicotinic acetylcholine receptors, belonging to the Cys-loop superfamily of ligand-gated ion channels (LGICs), are membrane proteins present in neurons and at neuromuscular junctions. They are responsible for signal transmission, and their function is regulated by neurotransmitters, agonists, and antagonists drugs. A detailed knowledge of their conformational transition in response to ligand binding is critical to understanding the basis of ligand-receptor interaction, in view of new pharmacological approaches to control receptor activity. However, the scarcity of experimentally derived structures of human channels makes this perspective extremely challenging. To contribute overcoming this issue, we have recently reported structural models for the open and the desensitized states of the human α7 nicotinic receptor. Here, we provide all-atom structural models of the same receptor in two different nonconductive states. The first structure, built via homology modeling and relaxed with extensive Molecular Dynamics simulations, represents the receptor bound to the natural antagonist α-conotoxin ImI. After comparison with available experimental data and computational models of other eukaryotic LGICs, we deem it consistent with the "closed-locked" state. The second model, obtained with simulations from the spontaneous relaxation of the open, agonist-bound α7 structure after ligand removal, recapitulates the characteristics of the apo-resting state of the receptor. These results add to our previous work on the active and desensitized state conformations, contributing to the structural characterization of the conformational landscape of the human α7 receptor and suggesting benchmarks to discriminate among conformations found in experiments or in simulations of LGICs. In particular key interactions at the interface between the extracellular domain and the transmembrane domain are identified, that could be critical to the α7 receptor function.

The neural γ2α1β2α1β2 gamma amino butyric acid ion channel receptor: structural analysis of the effects of the ivermectin molecule and disulfide bridges

While ~30% of the human genome encodes membrane proteins, only a handful of structures of membrane proteins have been resolved to high resolution. Here, we studied the structure of a member of the Cys-loop ligand gated ion channel protein superfamily of receptors, human type A γ₂α₁β₂α₁β₂ gamma amino butyric acid receptor complex in a lipid bilayer environment. Studying the correlation between the structure and function of the gamma amino butyric acid receptor may enhance our understanding of the molecular basis of ion channel dysfunctions linked with epilepsy, ataxia, migraine, schizophrenia and other neurodegenerative diseases. The structure of human γ₂α₁β₂α₁β₂ has been modeled based on the X-ray structure of the Caenorhabditis elegans glutamate-gated chloride channel via homology modeling. The template provided the first inhibitory channel structure for the Cys-loop superfamily of ligand-gated ion channels. The only available template structure before this glutamate-gated chloride channel was a cation selective channel which had very low sequence identity with gamma aminobutyric acid receptor. Here, our aim was to study the effect of structural corrections originating from modeling on a more reliable template structure. The homology model was analyzed for structural properties via a 100 ns molecular dynamics (MD) study. Due to the structural shifts and the removal of an open channel potentiator molecule, ivermectin, from the template structure, helical packing changes were observed in the transmembrane segment. Namely removal of ivermectin molecule caused a closure around the Leu 9 position along the ion channel. In terms of the structural shifts, there are three potential disulfide bridges between the M1 and M3 helices of the γ₂ and 2 α₁ subunits in the model. The effect of these disulfide bridges was investigated via monitoring the differences in root mean square fluctuations (RMSF) of individual amino acids and principal component analysis of the MD trajectory of the two homology models-one with the disulfide bridge and one with protonated Cys residues. In all subunit types, RMSF of the transmembrane domain helices are reduced in the presence of disulfide bridges. Additionally, loop A, loop F and loop C fluctuations were affected in the extracellular domain. In cross-correlation analysis of the trajectory, the two model structures displayed different coupling in between the M2-M3 linker region, protruding from the membrane, and the β1-β2/D loop and cys-loop regions in the extracellular domain. Correlations of the C loop, which collapses directly over the bound ligand molecule, were also affected by differences in the packing of transmembrane helices. Finally, more localized correlations were observed in the transmembrane helices when disulfide bridges were present in the model. The differences observed in this study suggest that dynamic coupling at the interface of extracellular and ion channel domains differs from the coupling introduced by disulfide bridges in the transmembrane region. We hope that this hypothesis will be tested experimentally in the near future.

An intersubunit electrostatic interaction in the GABAA receptor facilitates its responses to benzodiazepines

Benzodiazepines are positive allosteric modulators of the GABA_A receptor (GABA_AR), acting at the α-γ subunit interface to enhance GABA_AR function. GABA or benzodiazepine binding induces distinct conformational changes in the GABA_AR. The molecular rearrangements in the GABA_AR following benzodiazepine binding remain to be fully elucidated. Using two molecular models of the GABA_AR, we identified electrostatic interactions between specific amino acids at the α-γ subunit interface that were broken by, or formed after, benzodiazepine binding. Using two-electrode voltage clamp electrophysiology in Xenopus laevis oocytes, we investigated these interactions by substituting one or both amino acids of each potential pair. We found that Lys¹⁰⁴ in the α₁ subunit forms an electrostatic bond with Asp⁷⁵ of the γ₂ subunit after benzodiazepine binding and that this bond stabilizes the positively modified state of the receptor. Substitution of these two residues to cysteine and subsequent covalent linkage between them increased the receptor's sensitivity to low GABA concentrations and decreased its response to benzodiazepines, producing a GABA_AR that resembles a benzodiazepine-bound WT GABA_AR. Breaking this bond restored sensitivity to GABA to WT levels and increased the receptor's response to benzodiazepines. The α₁ Lys¹⁰⁴ and γ₂ Asp⁷⁵ interaction did not play a role in ethanol or neurosteroid modulation of GABA_AR, suggesting that different modulators induce different conformational changes in the receptor. These findings may help explain the additive or synergistic effects of modulators acting at the GABA_AR.

Un-gating and allosteric modulation of a pentameric ligand-gated ion channel captured by molecular dynamics

Pentameric ligand-gated ion channels (pLGICs) mediate intercellular communication at synapses through the opening of an ion pore in response to the binding of a neurotransmitter. Despite the increasing availability of high-resolution structures of pLGICs, a detailed understanding of the functional isomerization from closed to open (gating) and back is currently missing. Here, we provide the first atomistic description of the transition from open to closed (un-gating) in the glutamate-gated chloride channel (GluCl) from Caenorhabditis Elegans. Starting with the active-state structure solved in complex with the neurotransmitter L-glutamate and the positive allosteric modulator (PAM) ivermectin, we analyze the spontaneous relaxation of the channel upon removal of ivermectin by explicit solvent/membrane Molecular Dynamics (MD) simulations. The μs-long trajectories support the conclusion that ion-channel deactivation is mediated by two distinct quaternary transitions, i.e. a global receptor twisting followed by the radial expansion (or blooming) of the extracellular domain. At variance with previous models, we show that pore closing is exclusively regulated by the global twisting, which controls the position of the β1-β2 loop relative to the M2-M3 loop at the EC/TM domain interface. Additional simulations with L-glutamate restrained to the crystallographic binding mode and ivermectin removed indicate that the same twisting isomerization is regulated by agonist binding at the orthosteric site. These results provide a structural model for gating in pLGICs and suggest a plausible mechanism for the pharmacological action of PAMs in this neurotransmitter receptor family. The simulated un-gating converges to the X-ray structure of GluCl resting state both globally and locally, demonstrating the predictive character of state-of-art MD simulations.

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.

Identification of an Inhibitory Alcohol Binding Site in GABAA ρ1 Receptors

Alcohols inhibit γ-aminobutyric acid type A ρ1 receptor function. After introducing mutations in several positions of the second transmembrane helix in ρ1, we studied the effects of ethanol and hexanol on GABA responses using two-electrode voltage clamp electrophysiology in Xenopus laevis oocytes. The 6' mutations produced the following effects on ethanol and hexanol responses: small increase or no change (T6'M), increased inhibition (T6'V), and small potentiation (T6'Y and T6'F). The 5' mutations produced mainly increases in hexanol inhibition. Other mutations produced small (3' and 9') or no changes (2' and L277 in the first transmembrane domain) in alcohol effects. These results suggest an inhibitory alcohol binding site near the 6' position. Homology models of ρ1 receptors based on the X-ray structure of GluCl showed that the 2', 5', 6', and 9' residues were easily accessible from the ion pore, with 5' and 6' residues from neighboring subunits facing each other; L3' and L277 also faced the neighboring subunit. We tested ethanol through octanol on single and double mutated ρ1 receptors [ρ1(I15'S), ρ1(T6'Y), and ρ1(T6'Y,I15'S)] to further characterize the inhibitory alcohol pocket in the wild-type ρ1 receptor. The pocket can only bind relatively short-chain alcohols and is eliminated by introducing Y in the 6' position. Replacing the bulky 15' residue with a smaller side chain introduced a potentiating binding site, more sensitive to long-chain than to short-chain alcohols. In conclusion, the net alcohol effect on the ρ1 receptor is determined by the sum of its actions on inhibitory and potentiating sites.

Subunit stoichiometry and arrangement in a heteromeric glutamate-gated chloride channel

The invertebrate glutamate-gated chloride-selective receptors (GluClRs) are ion channels serving as targets for ivermectin (IVM), a broad-spectrum anthelmintic drug used to treat human parasitic diseases like river blindness and lymphatic filariasis. The native GluClR is a heteropentamer consisting of α and β subunit types, with yet unknown subunit stoichiometry and arrangement. Based on the recent crystal structure of a homomeric GluClαR, we introduced mutations at the intersubunit interfaces where Glu (the neurotransmitter) binds. By electrophysiological characterization of these mutants, we found heteromeric assemblies with two equivalent Glu-binding sites at β/α intersubunit interfaces, where the GluClβ and GluClα subunits, respectively, contribute the "principal" and "complementary" components of the putative Glu-binding pockets. We identified a mutation in the IVM-binding site (far away from the Glu-binding sites), which significantly increased the sensitivity of the heteromeric mutant receptor to both Glu and IVM, and improved the receptor subunits' cooperativity. We further characterized this heteromeric GluClR mutant as a receptor having a third Glu-binding site at an α/α intersubunit interface. Altogether, our data unveil heteromeric GluClR assemblies having three α and two β subunits arranged in a counterclockwise β-α-β-α-α fashion, as viewed from the extracellular side, with either two or three Glu-binding site interfaces.

Direct-Space Corrections Enable Fast and Accurate Lorentz-Berthelot Combination Rule Lennard-Jones Lattice Summation

Long-range lattice summation techniques such as the particle-mesh Ewald (PME) algorithm for electrostatics have been revolutionary to the precision and accuracy of molecular simulations in general. Despite the performance penalty associated with lattice summation electrostatics, few biomolecular simulations today are performed without it. There are increasingly strong arguments for moving in the same direction for Lennard-Jones (LJ) interactions, and by using geometric approximations of the combination rules in reciprocal space, we have been able to make a very high-performance implementation available in GROMACS. Here, we present a new way to correct for these approximations to achieve exact treatment of Lorentz-Berthelot combination rules within the cutoff, and only a very small approximation error remains outside the cutoff (a part that would be completely ignored without LJ-PME). This not only improves accuracy by almost an order of magnitude but also achieves absolute biomolecular simulation performance that is an order of magnitude faster than any other available lattice summation technique for LJ interactions. The implementation includes both CPU and GPU acceleration, and its combination with improved scaling LJ-PME simulations now provides performance close to the truncated potential methods in GROMACS but with much higher accuracy.

Interactions of L-3,5,3'-Triiodothyronine [corrected], Allopregnanolone, and Ivermectin with the GABAA Receptor: Evidence for Overlapping Intersubunit Binding Modes

Structural mechanisms of modulation of γ-aminobutyric acid (GABA) type A receptors by neurosteroids and hormones remain unclear. The thyroid hormone L-3,5,3’-triiodothyronine (T3) inhibits GABA_A receptors at micromolar concentrations and has common features with neurosteroids such as allopregnanolone (ALLOP). Here we use functional experiments on α₂β₁γ₂ GABA_A receptors expressed in Xenopus oocytes to detect competitive interactions between T3 and an agonist (ivermectin, IVM) with a crystallographically determined binding site at subunit interfaces in the transmembrane domain of a homologous receptor (glutamate-gated chloride channel, GluCl). T3 and ALLOP also show competitive effects, supporting the presence of both a T3 and ALLOP binding site at one or more subunit interfaces. Molecular dynamics (MD) simulations over 200 ns are used to investigate the dynamics and energetics of T3 in the identified intersubunit sites. In these simulations, T3 molecules occupying all intersubunit sites (with the exception of the α-β interface) display numerous energetically favorable conformations with multiple hydrogen bonding partners, including previously implicated polar/acidic sidechains and a structurally conserved deformation in the M1 backbone.

Functional characterization of ivermectin binding sites in α1β2γ2L GABA(A) receptors

GABA_A receptors (GABA_ARs) are the major inhibitory neurotransmitter receptors in the brain and are therapeutic targets for many indications including sedation, anesthesia and anxiolysis. There is, however, considerable scope for the development of new therapeutics with improved beneficial effects and reduced side-effect profiles. The anthelminthic drug, ivermectin, activates the GABA_AR although its binding site is not known. The molecular site of action of ivermectin has, however, been defined by crystallography in the homologous glutamate-gated chloride channel. Resolving the molecular mechanisms of ivermectin binding to α1β2γ2L GABA_ARs may provide insights into the design of improved therapeutics. Given that ivermectin binds to subunit interfaces, we sought to define (1) which subunit interface sites it binds to, (2) whether these sites are equivalent in terms of ivermectin sensitivity or efficacy, and (3) how many must be occupied for maximal efficacy. Our approach involved precluding ivermectin from binding to particular interfaces by introducing bulky M3 domain 36′F sidechains to the “+” side of those interfaces. We thereby demonstrated that ivermectin produces irreversible channel activation only when it binds to the single γ2L-β2 interface site. When it binds to α1-β2 sites it elicits potentiation of GABA-gated currents but has no irreversible activating effect. Ivermectin cannot bind to the β2-α1 interface site due to its endogenous bulky 36′ methionine. Replacing this with an alanine creates a functional site at this interface, but surprisingly it is inhibitory. Molecular docking simulations reveal that the γ2L-β2 interface forms more contacts with ivermectin than the other interfaces, possibly explaining why ivermectin appears to bind irreversibly at this interface. This study demonstrates unexpectedly stark pharmacological differences among GABA_AR ivermectin binding sites.

Conformational gating dynamics in the GluCl anion-selective chloride channel

Cys-loop receptors are central to propagation of signals in the nervous system. The gating of the membrane-spanning pore is triggered by structural rearrangements in the agonist-binding site, located some 50 Å away from the pore. A sequential conformational change, propagating from the ligand-binding site to the pore, has been proposed to govern gating in all Cys-loop receptors. Here, we identify structural and dynamic components of the conformational gating in the eukaryotic glutamate-gated chloride channel (GluCl) by means of molecular dynamics (MD) simulations with and without the l-glutamate agonist bound. A significant increase in pore opening and accompanying hydration is observed in the presence of glutamate. Potential of mean force calculations reveal that the barrier for ion passage drops from 15 kcal/mol to 5-10 kcal/mol with the agonist bound. This appears to be explained by agonist binding that leads to significant changes in the intersubunit hydrogen-bonding pattern, which induce a slight tilt of the extracellular domain relative to the transmembrane domain in the simulations. This rearrangement is subtle, but correspond to the direction of the quaternary twist observed as a key difference between open and closed X-ray structures. While the full reversible gating is still a much slower process, the observed structural dynamics sheds new light on the early stages of how the agonist influences the extracellular domain, how the extracellular domain interacts with the transmembrane domain, and how changes in the transmembrane domain alter the free energy of ion passage.

A Structural Model of the Human α7 Nicotinic Receptor in an Open Conformation

Nicotinic acetylcholine receptors (nAchRs) are ligand-gated ion channels that regulate chemical transmission at the neuromuscular junction. Structural information is available at low resolution from open and closed forms of an eukaryotic receptor, and at high resolution from other members of the same structural family, two prokaryotic orthologs and an eukaryotic GluCl channel. Structures of human channels however are still lacking. Homology modeling and Molecular Dynamics simulations are valuable tools to predict structures of unknown proteins, however, for the case of human nAchRs, they have been unsuccessful in providing a stable open structure so far. This is due to different problems with the template structures: on one side the homology with prokaryotic species is too low, while on the other the open eukaryotic GluCl proved itself unstable in several MD studies and collapsed to a dehydrated, non-conductive conformation, even when bound to an agonist. Aim of this work is to obtain, by a mixing of state-of-the-art homology and simulation techniques, a plausible prediction of the structure (still unknown) of the open state of human α7 nAChR complexed with epibatidine, from which it is possible to start structural and functional test studies. To prevent channel closure we employ a restraint that keeps the transmembrane pore open, and obtain in this way a stable, hydrated conformation. To further validate this conformation, we run four long, unbiased simulations starting from configurations chosen at random along the restrained trajectory. The channel remains stable and hydrated over the whole runs. This allows to assess the stability of the putative open conformation over a cumulative time of 1 μs, 800 ns of which are of unbiased simulation. Mostly based on the analysis of pore hydration and size, we suggest that the obtained structure has reasonable chances to be (at least one of the possible) structures of the channel in the open conformation.

The atypical cation-conduction and gating properties of ELIC underscore the marked functional versatility of the pentameric ligand-gated ion-channel fold

The superfamily of pentameric ligand-gated ion channels (pLGICs) is unique among ionotropic receptors in that the same overall structure has evolved to generate multiple members with different combinations of agonist specificities and permeant-ion charge selectivities. However, aside from these differences, pLGICs have been typically regarded as having several invariant functional properties. These include pore blockade by extracellular quaternary-ammonium cations in the micromolar-to-millimolar concentration range (in the case of the cation-selective members), and a gain-of-function phenotype, which manifests as a slower deactivation time course, as a result of mutations that reduce the hydrophobicity of the transmembrane pore lining. Here, we tested this notion on three distantly related cation-selective members of the pLGIC superfamily: the mouse muscle nicotinic acetylcholine receptor (nAChR), and the bacterial GLIC and ELIC channels. Remarkably, we found that, whereas low millimolar concentrations of TMA(+) and TEA(+) block the nAChR and GLIC, neither of these two quaternary-ammonium cations blocks ELIC at such concentrations; instead, both carry measurable inward currents when present as the only cations on the extracellular side. Also, we found that, whereas lidocaine binding speeds up the current-decay time courses of the nAChR and GLIC in the presence of saturating concentrations of agonists, the binding of lidocaine to ELIC slows this time course down. Furthermore, whereas mutations that reduce the hydrophobicity of the side chains at position 9' of the M2 α-helices greatly slowed the deactivation time course of the nAChR and GLIC, these mutations had little effect--or even sped up deactivation--when engineered in ELIC. Our data indicate that caution should be exercised when generalizing results obtained with ELIC to the rest of the pLGICs, but more intriguingly, they hint at the possibility that ELIC is a representative of a novel branch of the superfamily with markedly divergent pore properties despite a well-conserved three-dimensional architecture.

A predicted binding site for cholesterol on the GABAA receptor

Modulation of the GABA type A receptor (GABAAR) function by cholesterol and other steroids is documented at the functional level, yet its structural basis is largely unknown. Current data on structurally related modulators suggest that cholesterol binds to subunit interfaces between transmembrane domains of the GABAAR. We construct homology models of a human GABAAR based on the structure of the glutamate-gated chloride channel GluCl of Caenorhabditis elegans. The models show the possibility of previously unreported disulfide bridges linking the M1 and M3 transmembrane helices in the α and γ subunits. We discuss the biological relevance of such disulfide bridges. Using our models, we investigate cholesterol binding to intersubunit cavities of the GABAAR transmembrane domain. We find that very similar binding modes are predicted independently by three approaches: analogy with ivermectin in the GluCl crystal structure, automated docking by AutoDock, and spontaneous rebinding events in unbiased molecular dynamics simulations. Taken together, the models and atomistic simulations suggest a somewhat flexible binding mode, with several possible orientations. Finally, we explore the possibility that cholesterol promotes pore opening through a wedge mechanism.

X-ray structures of GluCl in apo states reveal a gating mechanism of Cys-loop receptors

Cys-loop receptors are neurotransmitter-gated ion channels that are essential mediators of fast chemical neurotransmission and are associated with a large number of neurological diseases and disorders, as well as parasitic infections. Members of this ion channel superfamily mediate excitatory or inhibitory neurotransmission depending on their ligand and ion selectivity. Structural information for Cys-loop receptors comes from several sources including electron microscopic studies of the nicotinic acetylcholine receptor, high-resolution X-ray structures of extracellular domains and X-ray structures of bacterial orthologues. In 2011 our group published structures of the Caenorhabditis elegans glutamate-gated chloride channel (GluCl) in complex with the allosteric partial agonist ivermectin, which provided insights into the structure of a possibly open state of a eukaryotic Cys-loop receptor, the basis for anion selectivity and channel block, and the mechanism by which ivermectin and related molecules stabilize the open state and potentiate neurotransmitter binding. However, there remain unanswered questions about the mechanism of channel opening and closing, the location and nature of the shut ion channel gate, the transitions between the closed/resting, open/activated and closed/desensitized states, and the mechanism by which conformational changes are coupled between the extracellular, orthosteric agonist binding domain and the transmembrane, ion channel domain. Here we present two conformationally distinct structures of C. elegans GluCl in the absence of ivermectin. Structural comparisons reveal a quaternary activation mechanism arising from rigid-body movements between the extracellular and transmembrane domains and a mechanism for modulation of the receptor by phospholipids.

Agonist and antagonist binding in human glycine receptors

The human glycine receptor (hGlyR) is an anion-permeable ligand-gated channel that is part of a larger superfamily of receptors called the Cys-loop family. hGlyRs are particularly amenable to single-channel recordings, thus making them a model experimental system for understanding the Cys-loop receptor family in general. Understanding the relationship between agonist binding and efficacy in Cys-loop receptors should improve our future prospects for making specific agonists or antagonists. However, at present, there is no high-resolution structure for the complete hGlyR, and thus, modeling is needed to provide a physical framework on which to interpret single-channel data. The structure of the glutamate-gated chloride channel from Caenorhabditis elegans shows a much higher level of sequence identity to human hGlyR than previous templates such as AChBP or the bacterial channels, GLIC and ELIC. Thus, we constructed a model of the hGlyR and validated it against previously reported mutagenesis data. We used molecular dynamics to refine the model and to explore binding of both an agonist (glycine) and an antagonist (strychnine). The model shows excellent agreement with previous data but also suggests some unique features: (i) a water molecule that forms part of the binding site and allows us to account for some previous results that were difficult to reconcile, (ii) an interaction of the glycine agonist with S129, and (iii) an effect from E211, both of which we confirmed with new site-directed mutagenesis and patch clamp recordings. Finally, examination of the simulations suggests that strychnine binding induces movement to a conformational state distinct from the glycine-bound or apo state, not only within the ligand-binding domain but also in the transmembrane domain.

Pentameric Ligand-gated Ion Channels : Insights from Computation

Pentameric ligand-gated ion channels (pLGICs) conduct upon the binding of an agonist and are fundamental to neurotransmission. New insights into the complex mechanisms underlying pLGIC gating, ion selectivity, and modulation have recently been gained via a series of crystal structures in prokaryotes and C .elegans, as well as computational studies relying on these structures. Here we review contributions from a variety of computational approaches, including normal mode analysis, automated docking, and fully atomistic molecular dynamics simulation. Examples from our own research, particularly concerning interactions with general anesthetics and lipids, are used to illustrate predictive results complementary to crystallographic studies.