Asymmetric structural motions of the homomeric alpha7 nicotinic receptor ligand binding domain revealed by molecular dynamics simulation

Kinetic Mechanism of Surfactant-Based Molecular Recognition: Selective Permeability across an Oil-Water Interface Regulated by Supramolecular Aggregates

Lipids are known to play a vital role in the molecular organization of all cellular life. Molecular recognition is another fundamental biological process that is generally attributed to biological polymers, such as proteins and nucleic acids. However, there is evidence that aggregates of lipids and lipid-like molecules are also capable of selectively binding to or regulating the partitioning of other molecules. We previously demonstrated that a model two-phase octanol/water system can selectively partition Red 40 and Blue 1 dyes added to an aqueous phase, with the selectivity depending on the surfactant (e.g., cetyltrimethylammonium bromide) dissolved in the organic phase. Here, we elucidate the mechanism of molecular recognition in this system by using quantitative partitioning experiments and molecular dynamics (MD) simulations. Our results indicate that the selectivity for the red dye is thermodynamically favored at all surfactant concentrations, while selectivity for the blue dye is kinetically favored at high surfactant concentrations. The kinetic selectivity for the blue dye correlates with the presence of molecular aggregation at the oil-water interface. Coarse-grained MD simulations elucidate nanoscale supramolecular structures that can preferentially bind one small molecule rather than another at an interface, providing a selectively permeable barrier in the absence of proteins. The results suggest a new supramolecular mechanism for molecular recognition with potential applications in drug delivery, drug discovery, and biosensing.

Therapeutic Targeting of α7 Nicotinic Acetylcholine Receptors

The α7-type nicotinic acetylcholine receptor is one of the most unique and interesting of all the members of the cys-loop superfamily of ligand-gated ion channels. Since it was first identified initially as a binding site for α-bungarotoxin in mammalian brain and later as a functional homomeric receptor with relatively high calcium permeability, it has been pursued as a potential therapeutic target for numerous indications, from Alzheimer disease to asthma. In this review, we discuss the history and state of the art for targeting α7 receptors, beginning with subtype-selective agonists and the basic pharmacophore for the selective activation of α7 receptors. A key feature of α7 receptors is their rapid desensitization by standard "orthosteric" agonist, and we discuss insights into the conformational landscape of α7 receptors that has been gained by the development of ligands binding to allosteric sites. Some of these sites are targeted by positive allosteric modulators that have a wide range of effects on the activation profile of the receptors. Other sites are targeted by direct allosteric agonist or antagonists. We include a perspective on the potential importance of α7 receptors for metabotropic as well as ionotropic signaling. We outline the challenges that exist for future development of drugs to target this important receptor and approaches that may be considered to address those challenges. SIGNIFICANCE STATEMENT: The α7-type nicotinic acetylcholine receptor (nAChR) is acknowledged as a potentially important therapeutic target with functional properties associated with both ionotropic and metabotropic signaling. The functional properties of α7 nAChR can be regulated in diverse ways with the variety of orthosteric and allosteric ligands described in this review.

Stable desensitization of α7 nicotinic acetylcholine receptors by NS6740 requires interaction with S36 in the orthosteric agonist binding site

NS6740 is an α₇ nicotinic acetylcholine receptor-selective partial agonist with low efficacy for channel activation, capable of promoting the stable conversion of the receptors to nonconducting (desensitized) states that can be reactivated with the application of positive allosteric modulators (PAMs). In spite of its low efficacy for channel activation, NS6740 is an effective activator of the cholinergic anti-inflammatory pathway. We observed that the concentration-response relationships for channel activation, both when applied alone and when co-applied with the PAM PNU-120596 are inverted-U shaped with inhibitory/desensitizing activities dominant at high concentrations. We evaluated the potential importance of recently identified binding sites for allosteric activators and tested the hypotheses that the stable desensitization produced by NS6740 may be due to binding to these sites. Our experiments were guided by molecular modeling of NS6740 binding to both the allosteric and orthosteric activation sites on the receptor. Our results indicate that with α₇C190A mutants, which have compromised orthosteric activation sites, NS6740 may work at the allosteric activation sites to promote transient PAM-dependent currents but not the stable desensitization seen with wild-type α₇ receptors. Modeling NS6740 in the orthosteric binding sites identified S36 as an important residue for NS6740 binding and predicted that an S36V mutation would limit NS6740 activity. The efficacy of NS6740 for α₇S36V receptors was reduced to zero, and applications of the compound to α₇S36V receptors failed to induce the desensitization observed with wild-type receptors. The results indicate that the unique properties of NS6740 are due primarily to binding at the sites for orthosteric agonists.

Proline Residues Contribute to Efficient GABAp Receptor Function

GABAp receptors are homomeric pentameric ligand-gated ion channels (pLGICs) and are useful for probing the molecular details of the mechanism of action in this important protein family. Here, we explore the role of proline (Pro) residues by creating mutant receptors, expressing them in HEK293 cells, and using fluorescent membrane potential sensitive dye to monitor receptor activity. The data revealed that 3 of the Pro-to-alanine substitutions resulted in nonfunctional receptors (one in the Cys-loop, one in loop A and one in the β2-β3 loop), 7 resulted in increased EC₅₀ values, and the remaining 13 resulted in receptors with properties similar to WT receptors. Further exploration of the nonfunctional mutants using expression in Xenopus laevis oocytes and whole-cell voltage-clamp electrophysiology, incorporating both canonical and noncanonical amino acids, revealed that the Pro in the Cys-loop had a preference for analogues with a high intrinsic cis bias, the Pro in loop A required a ring, and the β2-β3 loop Pro contributes to expression. The data from the Cys-loop Pro are consistent with those from other pLGICs, while those of the loop A Pro and some of the other Pros surprisingly differ. Thus, overall, the data clarify the roles of many of the Pros in the GABAp receptor and also suggest that caution must be applied in using data from one receptor to understand molecular details of all pLGICs.

A Computational Analysis of the Factors Governing the Dynamics of α7 nAChR and Its Homologs

The α7 nicotinic acetylcholine receptor is a homopentameric ion channel from the Cys-loop receptor superfamily targeted for psychiatric indications and inflammatory pain. Molecular dynamics studies of the receptor have focused on residue mobility and global conformational changes to address receptor function. However, a comparative analysis of α7 with its homologs that cannot trigger channel opening has not been made so far. To identify the residues involved in α7 activation, we ran triplicate 500-ns molecular dynamics simulations with an α7 extracellular domain homology model and two acetylcholine-binding protein homologs. We tested the effect of ligand binding and amino acid sequence on the structure and dynamics of the three proteins. We found that mobile regions identified based on root mean-square deviation and root mean-square fluctuation values are not always consistent among the individual α7 extracellular domain simulations. Comparison of the replica-average properties of the three proteins based on dynamic cross-correlation maps showed that ligand binding affects the coupling between the C-loop and the Cys-loop, vestibular loop, and β1-β2 loops. In addition, the main-immunogenic-region-like domain of α7 went through correlated motions with multiple domains of the receptor. These correlated motions were absent or diminished in α7 homologs, suggesting a unique role in α7 activation.

Many Proline Residues in the Extracellular Domain Contribute to Glycine Receptor Function

Prolines in signaling proteins are of particular interest because they have a range of unique properties that may be critical for function. Here we show that many proline residues in the extracellular domain (ECD) of the glycine receptor are involved in the correct functioning of this ligand-gated ion channel. We explore their role by creating mutant receptors, expressing them in cells, and using fluorescent membrane potential sensitive dye to monitor receptor activity. We then interpret the changes in receptor parameters using structural information from the open and closed states of the receptor. The data reveal that substitution with alanine of ten of the 13 Pro residues in the ECD alters the function of the receptor: one substitution ablates function, six cause a decrease in the EC₅₀, and three cause an increase. Only three of these mutants result in EC₅₀ values similar to WT. The nonfunctional mutant, Pro30Ala, was further probed in oocytes, and the data suggest a role in both expression and function. Examination of the locations of sensitive Pro residues in the receptor and identification of potential interactions with nearby residues reveal how these residues could contribute to the correct functioning of this typical pentameric ligand-gated ion channel.

Nicotinic receptor pharmacology in silico: Insights and challenges

Nicotinic acetylcholine receptors (nAChR) are homo- or hetero-pentameric ligand-gated ion channels of the Cys-loop superfamily and play important roles in the nervous system and muscles. Studies on nAChR benefit from in silico modeling due to the lack of high-resolution structures for most receptor subtypes and challenges in experiments addressing the complex mechanism of activation involving allosteric sites. Although there is myriad of computational modeling studies on nAChR, the multitude of the methods and parameters used in these studies makes modeling nAChR a daunting task, particularly for the non-experts in the field. To address this problem, the modeling literature on Torpedo nAChR and α7 nAChR were focused on as examples of heteromeric and homomeric nAChR, and the key in silico modeling studies between the years 1995-2019 were concisely reviewed. This was followed by a critical analysis of these studies by comparing the findings with each other and with the emerging experimental and computational data on nAChR. Based on these critical analyses, suggestions were made to guide the future researchers in the field of in silico modeling of nAChR. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

In Silico Modeling of the α7 Nicotinic Acetylcholine Receptor: New Pharmacological Challenges Associated with Multiple Modes of Signaling

The α7 nicotinic acetylcholine receptor is a homopentameric ion-channel of the Cys-loop superfamily characterized by its low probability of opening, high calcium permeability, and rapid desensitization. The α7 receptor has been targeted for the treatment of the cognitive symptoms of schizophrenia, depression, and Alzheimer's disease, but it is also involved in inflammatory modulation as a part of the cholinergic anti-inflammatory pathway. Despite its functional importance, in silico studies of the α7 receptor cannot produce a general model explaining the structural features of receptor activation, nor predict the mode of action for various ligand classes. Two particular problems in modeling the α7 nAChR are the absence of a high-resolution structure and the presence of five potentially nonequivalent orthosteric ligand binding sites. There is wide variability regarding the templates used for homology modeling, types of ligands investigated, simulation methods, and simulation times. However, a systematic survey focusing on the methodological similarities and differences in modeling α7 has not been done. In this work, we make a critical analysis of the modeling literature of α7 nAChR by comparing the findings of computational studies with each other and with experimental studies under the main topics of structural studies, ligand binding studies, and comparisons with other nAChR. In light of our findings, we also summarize current problems in the field and make suggestions for future studies concerning modeling of the α7 receptor.

Structural effects of divalent calcium cations on the α7 nicotinic acetylcholine receptor: A molecular dynamics simulation study

The α7 subtype of neuronal nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel protein that is vital to various neurological functions, including modulation of neurotransmitter release. A relatively high concentration of extracellular Ca²⁺ in the neuronal environment is likely to exert substantial structural and functional influence on nAChRs, which may affect their interactions with agonists and antagonists. In this work, we employed atomistic molecular dynamics (MD) simulations to examine the effects of elevated Ca²⁺ on the structure and dynamics of α7 nAChR embedded in a model phospholipid bilayer. Our results suggest that the presence of Ca²⁺ in the α7 nAChR environment results in closure of loop C-in the extracellular ligand-binding domain, a motion normally associated with agonist binding and receptor activation. Elevated Ca²⁺ also alters the conformation of key regions of the receptor, including the inter-helical loops, pore-lining helices and the "gate" residues, and causes partial channel opening in the absence of an agonist, leading to an attendant reduction in the free energy of Ca²⁺ permeation through the pore as elucidated by umbrella sampling simulations. Overall, the structural and permeability changes in α7 nAChR suggest that elevated Ca²⁺ induces a partially activated receptor state that is distinct from both the resting and the agonist-activated states. These results are consistent with the notion that divalent ions can serve as a potentiator of nAChRs, resulting in a higher rate of receptor activation (and subsequent desensitization) in the presence of agonists, with possible implications for diseases involving calcium dysregulation.

Allosteric Agonism of α7 Nicotinic Acetylcholine Receptors: Receptor Modulation Outside the Orthosteric Site

Nicotinic acetylcholine receptors (nAChRs) are members of the Cys-loop superfamily of ligand-gated ion channels. Typically, channel activation follows the binding of agonists to the orthosteric binding sites of the receptor. α7 nAChRs have a very low probability of channel activation, which can be reversed by the binding of α7 selective positive allosteric modulators (PAMs) to putative sites within the transmembrane domains. Although typical PAMs, like PNU-120596, require coapplication of an orthosteric agonist to produce large channel activations, some, like GAT107 and B-973B [(S)-3-(3,4-difluorophenyl)-N-(1-(6-(4-(pyridin-2-yl)piperazin-1-yl)pyrazin-2-yl)ethyl)propanamide], are characterized as allosteric activating PAMs, which also bind to an allosteric activation (AA) site in the extracellular domain and activate the α7 ion channel by themselves. We had previously characterized N,N-diethyl-N'-phenylpiperazine analogs with various functions. In this work, we docked members of this family to a homology model of the α7 receptor extracellular domain. The compound 1,1-diethyl-4(naphthalene-2-yl)piperazin-1-ium (2NDEP) a weak partial agonist, showed particularly favorable docking and binding energies at the putative AA site of the receptor. We hypothesized that 2NDEP could couple with PAMs through the AA site. This hypothesis was tested with the α7 mutant C190A, which is not activated by orthosteric agonists but is effectively activated by GAT107. The results showed that 2NDEP acts as an allosteric agonist of α7C190A when coapplied with the PAM PNU-120596. Also, the allosteric activity was nearly abolished upon coapplication with the AA site-selective antagonist 2,3,5,6MP-TQS (cis-trans-4-(2,3,5,6-tetramethylphenyl)-3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-8-sulfonamide), consistent with AA site involvement. Overall, our findings show a novel mode of agonism through an allosteric site in the extracellular domain of α7 nAChR.

Macroscopic and Microscopic Activation of α7 Nicotinic Acetylcholine Receptors by the Structurally Unrelated Allosteric Agonist-Positive Allosteric Modulators (ago-PAMs) B-973B and GAT107

B-973 is an efficacious type II positive allosteric modulator (PAM) of α7 nicotinic acetylcholine receptors that, like 4BP-TQS and its active isomer GAT107, can produce direct allosteric activation in addition to potentiation of orthosteric agonist activity, which identifies it as an allosteric activating (ago)-PAM. We compared the properties of B-973B, the active enantiomer of B-973, with those of GAT107 regarding the separation of allosteric potentiation and activation. Both ago-PAMs can strongly activate mutants of α7 that are insensitive to standard orthosteric agonists like acetylcholine. Likewise, the activity of both ago-PAMs is largely eliminated by the M254L mutation in the putative transmembrane PAM-binding site. Allosteric activation by B-973B appeared more protracted than that produced by GAT107, and B-973B responses were relatively insensitive to the noncompetitive antagonist mecamylamine compared with GAT107 responses. Similar differences are also seen in the single-channel currents. The two agents generate unique profiles of full-conductance and subconductance states, with B-973B producing protracted bursts, even in the presence of mecamylamine. Modeling and docking studies suggest that the molecular basis for these effects depends on specific interactions in both the extracellular and transmembrane domains of the receptor.

Multiple regions in the extracellular domain of the glycine receptor determine receptor activity

Glycine receptors (GlyRs) are Cys-loop receptors that mediate fast synaptic inhibition in the brain stem and spinal cord. They are involved in the generation of motor rhythm, reflex circuit coordination, and sensory signal processing and therefore represent targets for therapeutic interventions. The extracellular domains (ECDs) of Cys-loop receptors typically contain many aromatic amino acids, but only those in the receptor binding pocket have been extensively studied. Here, we show that many Phe residues in the ECD that are not located in the binding pocket are also involved in GlyR function. We examined these Phe residues by creating several GlyR variants, characterizing these variants with the two-electrode voltage clamp technique in Xenopus oocytes, and interpreting changes in receptor parameters by using currently available structural information on the open and closed states of the GlyR. Substitution of six of the eight Phe residues in the ECD with Ala resulted in loss of function or significantly increased the EC₅₀ and also altered the maximal response to the partial GlyR agonist taurine compared with glycine in those receptor variants that were functional. Substitutions with other amino acids, combined with examination of nearby residues that could potentially interact with these Phe residues, suggested interactions that could be important for GlyR function, and possibly similar interactions could contribute to the function of other members of the Cys-loop receptor family. Overall, our results suggest that many ECD regions are important for GlyR function and that these regions could inform the design of therapeutic agents targeting GlyR activity.

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.

Interactions of acetylcholine binding site residues contributing to nicotinic acetylcholine receptor gating: role of residues Y93, Y190, K145 and D200

The nicotinic acetylcholine receptor exhibits multiple conformational states, resting (channel closed), active (channel open) and desensitized (channel closed). The resting state may be distinguished from the active and desensitized states by the orientation of loop C in the extracellular ligand binding domain (LBD). Homology modeling was used to generate structures of the Torpedo californica α2βδγ nAChR that initially represent the resting state (loop C open) and the desensitized state (loop C closed). Molecular dynamics (MD) simulations were performed on the extracellular LBD on each nAChR conformational state, with and without the agonist anabaseine present in each binding site (the αγ and the αδ sites). Three MD simulations of 10ns each were performed for each of the four conditions. Comparison of dynamics revealed that in the presence of agonist, loop C was drawn inward and attains a more stable conformation. Examination of side-chain interactions revealed that residue αY190 exhibited hydrogen-bonding interactions either with residue αY93 in the ligand binding site or with residue αK145 proximal to the binding site. αK145 also exhibited side chain (salt bridge) interactions with αD200 and main chain interactions with αY93. Residues αW149, αY198, γY116/δT119, γL118/δL121 and γL108/δL111 appear to play the role of stabilizing ligand in the binding site. In MD simulations for the desensitized state, the effect of ligand upon the interactions among αK145, αY190, and αY93 as well as ligand-hydrogen-bonding to αW149 were more pronounced at the αγ interface than at the αδ interface. Differences in affinity for the desensitized state were determined experimentally to be 10-fold. The changes in side chain interactions observed for the two conformations and induced by ligand support a model wherein hydrogen bond interactions between αD200 and αY93 are broken and rearrange to form a salt-bridge between αK145 and αD200 and hydrogen bond interactions between αY93 and αY190 and between αK145 and αY190.

End-plate acetylcholine receptor: structure, mechanism, pharmacology, and disease

The synapse is a localized neurohumoral contact between a neuron and an effector cell and may be considered the quantum of fast intercellular communication. Analogously, the postsynaptic neurotransmitter receptor may be considered the quantum of fast chemical to electrical transduction. Our understanding of postsynaptic receptors began to develop about a hundred years ago with the demonstration that electrical stimulation of the vagus nerve released acetylcholine and slowed the heart beat. During the past 50 years, advances in understanding postsynaptic receptors increased at a rapid pace, owing largely to studies of the acetylcholine receptor (AChR) at the motor endplate. The endplate AChR belongs to a large superfamily of neurotransmitter receptors, called Cys-loop receptors, and has served as an exemplar receptor for probing fundamental structures and mechanisms that underlie fast synaptic transmission in the central and peripheral nervous systems. Recent studies provide an increasingly detailed picture of the structure of the AChR and the symphony of molecular motions that underpin its remarkably fast and efficient chemoelectrical transduction.

Structural link between γ-aminobutyric acid type A (GABAA) receptor agonist binding site and inner β-sheet governs channel activation and allosteric drug modulation

Rapid opening and closing of pentameric ligand-gated ion channels (pLGICs) regulate information flow throughout the brain. For pLGICs, it is postulated that neurotransmitter-induced movements in the extracellular inner β-sheet trigger channel activation. Homology modeling reveals that the β4-β5 linker physically connects the neurotransmitter binding site to the inner β-sheet. Inserting 1, 2, 4, and 8 glycines in this region of the GABA(A) receptor β-subunit progressively decreases GABA activation and converts the competitive antagonist SR-95531 into a partial agonist, demonstrating that this linker is a key element whose length and flexibility are optimized for efficient signal propagation. Insertions in the α- and γ-subunits have little effect on GABA or SR-95531 actions, suggesting that asymmetric motions in the extracellular domain power pLGIC gating. The effects of insertions on allosteric modulator actions, pentobarbital, and benzodiazepines, have different subunit dependences, indicating that modulator-induced signaling is distinct from agonist gating.

Blockade of neuronal α7-nAChR by α-conotoxin ImI explained by computational scanning and energy calculations

α-Conotoxins potently inhibit isoforms of nicotinic acetylcholine receptors (nAChRs), which are essential for neuronal and neuromuscular transmission. They are also used as neurochemical tools to study nAChR physiology and are being evaluated as drug leads to treat various neuronal disorders. A number of experimental studies have been performed to investigate the structure-activity relationships of conotoxin/nAChR complexes. However, the structural determinants of their binding interactions are still ambiguous in the absence of experimental structures of conotoxin-receptor complexes. In this study, the binding modes of α-conotoxin ImI to the α7-nAChR, currently the best-studied system experimentally, were investigated using comparative modeling and molecular dynamics simulations. The structures of more than 30 single point mutants of either the conotoxin or the receptor were modeled and analyzed. The models were used to explain qualitatively the change of affinities measured experimentally, including some nAChR positions located outside the binding site. Mutational energies were calculated using different methods that combine a conformational refinement procedure (minimization with a distance dependent dielectric constant or explicit water, or molecular dynamics using five restraint strategies) and a binding energy function (MM-GB/SA or MM-PB/SA). The protocol using explicit water energy minimization and MM-GB/SA gave the best correlations with experimental binding affinities, with an R² value of 0.74. The van der Waals and non-polar desolvation components were found to be the main driving force for binding of the conotoxin to the nAChR. The electrostatic component was responsible for the selectivity of the various ImI mutants. Overall, this study provides novel insights into the binding mechanism of α-conotoxins to nAChRs and the methodological developments reported here open avenues for computational scanning studies of a rapidly expanding range of wild-type and chemically modified α-conotoxins.

Conotoxins are peptide toxins extracted from the venom of carnivorous marine cone snails. Members of the α-conotoxin subfamily potently block nicotinic acetylcholine receptors (nAChRs), which are involved in signal transmission between two neurons or between neurons and muscle fibers. nAChRs are important pharmacological targets due to their involvement in the transmission of pain stimuli and also in numerous neurone diseases and disorders. Their potency and specificity have led to the development of α-conotoxins as drug leads, and also to their use in the investigation of the role of nAChRs in various physiological processes. The most studied conotoxin/nAChR system, ImI/α7, was modeled in this study, and several computational methods were tested for their ability to explain the perturbations observed experimentally after introducing single point mutations into either ImI or the α7 receptor. The aim of this study was to establish a theoretical basis to rapidly identify new α-conotoxin mutants that might have improved specificity and affinity for a given receptor subtype. Furthermore, hundreds of thousands of conotoxins are predicted to exist, and computational methods are needed to help streamline the discovery of their molecular targets.

Decrypting the sequence of structural events during the gating transition of pentameric ligand-gated ion channels based on an interpolated elastic network model

Despite many experimental and computational studies of the gating transition of pentameric ligand-gated ion channels (pLGICs), the structural basis of how ligand binding couples to channel gating remains unknown. By using a newly developed interpolated elastic network model (iENM), we have attempted to compute a likely transition pathway from the closed- to the open-channel conformation of pLGICs as captured by the crystal structures of two prokaryotic pLGICs. The iENM pathway predicts a sequence of structural events that begins at the ligand-binding loops and is followed by the displacements of two key loops (loop 2 and loop 7) at the interface between the extracellular and transmembrane domain, the tilting/bending of the pore-lining M2 helix, and subsequent movements of M4, M3 and M1 helices in the transmembrane domain. The predicted order of structural events is in broad agreement with the Φ-value analysis of α subunit of nicotinic acetylcholine receptor mutants, which supports a conserved core mechanism for ligand-gated channel opening in pLGICs. Further perturbation analysis has supported the critical role of certain intra-subunit and inter-subunit interactions in dictating the above sequence of events.

Unresponsive correlated motion in alpha7 nAChR to halothane binding explains its functional insensitivity to volatile anesthetics

Neuronal nicotinic acetylcholine receptors (nAChRs) have been implicated as targets for general anesthetics, but the functional responses to anesthetic modulation vary considerably among different subtypes of nAChRs. Inhaled general anesthetics, such as halothane, could effectively inhibit the channel activity of the alpha4beta2 nAChR but not the homologous alpha7 nAChR. To understand why alpha7 is insensitive to inhaled general anesthetics, we performed multiple sets of 20 ns molecular dynamics (MD) simulations on the closed- and open-channel alpha7 in the absence and presence of halothane and critically compared the results with those from our studies on the alpha4beta2 nAChR (Liu et al. J. Phys. Chem. B 2009, 113, 12581 and Liu et al. J. Phys. Chem. B 2010, 114, 626). Several halothane binding sites with fairly high binding affinities were identified in both closed- and open-channel alpha7, suggesting that a lack of sensitive functional responses of the alpha7 nAChR to halothane in the previous experiments was unlikely due to a lack of halothane interaction with alpha7. The binding affinities of halothane in alpha7 seemed to be protein conformation-dependent. Overall, halothane affinity was higher in the closed-channel alpha7. Halothane binding to alpha7 did not induce profound changes in alpha7 structure and dynamics that could be related to the channel function. In contrast, correlated motion of the open-channel alpha4beta2 was reduced substantially in the presence of halothane, primarily due to the more susceptible nature of beta2 to anesthetic modulation. The amphiphilic extracellular and transmembrane domain interface of the beta2 subunit is attractive to halothane and is susceptible to halothane perturbation, which may be why alpha4beta2 is functionally more sensitive to halothane than alpha7.

Mapping spatial relationships between residues in the ligand-binding domain of the 5-HT3 receptor using a molecular ruler

The serotonin 5-HT(3) receptor (5-HT(3)R) is a member of the Cys-loop ligand-gated ion channel family. We used a combination of site-directed mutagenesis, homology modeling, and ligand-docking simulations to analyze antagonist-receptor interactions. Mutation of E236, which is near loop C of the binding site, to aspartate prevents expression of the receptor on the cell surface, and no specific ligand binding can be detected. On the other hand, mutation to glutamine, asparagine, or alanine produces receptors that are expressed on the cell surface, but decreases receptor affinity for the competitive antagonist d-tubocurarine (dTC) 5-35-fold. The results of a double-mutant cycle analysis employing a panel of dTC analogs to identify specific points of interactions between the dTC analogs and E236 are consistent with E236 making a direct physical interaction with the 12 -OH of dTC. dTC is a rigid molecule of known three-dimensional structure. Together with previous studies linking other regions of dTC to specific residues in the binding site, these data allow us to define the relative spatial arrangement of three different residues in the ligand-binding site: R92 (loop D), N128 (loop A), and E236 (near loop C). Molecular modeling employing these distance constraints followed by molecular-dynamics simulations produced a dTC/receptor complex consistent with the experimental data. The use of the rigid ligands as molecular rulers in conjunction with double-mutant cycle analysis provides a means of mapping the relative positions of various residues in the ligand-binding site of any ligand-receptor complex, and thus is a useful tool for delineating the architecture of the binding site.

In silico characterization of cytisinoids docked into an acetylcholine binding protein

Homology models of nicotinic acetylcholine receptors (nAChRs) suggest that subtype specificity is due to non-conserved residues in the complementary subunit of the ligand-binding pocket. Cytisine and its derivatives generally show a strong preference for heteromeric alpha4beta2* nAChRs over the homomeric alpha7 subtype, and the structural modifications studied do not cause large changes in their nAChR subtype selectivity. In the present work we docked cytisine, N-methylcytisine, and several pyridone ring-substituted cytisinoids into the crystallographic structure of the Lymnaea stagnalis acetylcholine binding protein (AChBP) co-crystallized with nicotine (1UW6). The graphical analysis of the best poses showed that cytisinoids have weak interactions with the side chains of the non-conserved amino acids in the complementary subunit justifying the use of PDB 1UWB as a surrogate for nAChR. Furthermore, we found a high correlation (R(2)=0.96) between the experimental pIC(50) values at alpha4beta2* nAChR and docking energy (S) of the best cytisinoid poses within the AChBP. Due to the quality of the correlation we suggest that this equation might be used as a predictive model to propose new cytisine-derived nAChRs ligands. Our docking results also suggest that further structural modifications of these cytisinoids will not greatly alter their alpha4beta2*/alpha7 selectivity.

Higher susceptibility to halothane modulation in open- than in closed-channel alpha4beta2 nAChR revealed by molecular dynamics simulations

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is a potential molecular target for general anesthetics. It is unclear, however, whether anesthetic action produces the same effect on the open and closed channels. Computations parallel to our previous open channel study (J. Phys. Chem. B 2009, 113, 12581) were performed on the closed-channel alpha4beta2 nAChR to investigate the conformation-dependent anesthetic effects on channel structures and dynamics. Flexible ligand docking and over 20 ns molecular dynamics simulations revealed similar halothane-binding sites in the closed and open channels. The sites with relatively high binding affinities (approximately -6.0 kcal/mol) were identified at the interface of extracellular (EC) and transmembrane (TM) domains or at the interface between alpha4 and beta2 subunits. Despite similar sites for halothane binding, the closed-channel conformation showed much less sensitivity than the open channel to the structural and dynamical perturbations from halothane. Compared to the systems without anesthetics, the amount of water inside the pore decreased by 22% in the presence of halothane in the open channel but only by 6% in the closed channel. Comparison of the nonbonded interactions at the EC/TM interfaces suggested that the beta2 subunits were more prone than the alpha4 subunits to halothane binding. In addition, our data support the notion that halothane exerts its effect by disturbing the quaternary structure and dynamics of the channel. The study concludes that sensitivity and global dynamics responsiveness of alpha4beta2 nAChR to halothane are conformation dependent. The effect of halothane on the global dynamics of the open-channel conformation might also account for the action of other inhaled general anesthetics.

Structural model of rho1 GABAC receptor based on evolutionary analysis: Testing of predicted protein-protein interactions involved in receptor assembly and function

The homopentameric rho1 GABA(C) receptor is a ligand-gated ion channel with a binding pocket for gamma-aminobutyric acid (GABA) at the interfaces of N-terminal extracellular domains. We combined evolutionary analysis, structural modeling, and experimental testing to study determinants of GABA(C) receptor assembly and channel gating. We estimated the posterior probability of selection pressure at amino acid residue sites measured as omega-values and built a comparative structural model, which identified several polar residues under strong selection pressure at the subunit interfaces that may form intersubunit hydrogen bonds or salt bridges. At three selected sites (R111, T151, and E55), mutations disrupting intersubunit interactions had strong effects on receptor folding, assembly, and function. We next examined the role of a predicted intersubunit salt bridge for residue pair R158-D204. The mutant R158D, where the positively charged residue is replaced by a negatively charged aspartate, yielded a partially degraded receptor and lacked membrane surface expression. The membrane surface expression was rescued by the double mutant R158D-D204R, where positive and negative charges are switched, although the mutant receptor was inactive. The single mutants R158A, D204R, and D204A exhibited diminished activities and altered kinetic profiles with fast recovery kinetics, suggesting that R158-D204 salt bridge perhaps stabilizes the open state of the GABA(C) receptor. Our results emphasize the functional importance of highly conserved polar residues at the protein-protein interfaces in GABA(C) rho1 receptors and demonstrate how the integration of computational and experimental approaches can aid discovery of functionally important interactions.

Main immunogenic region structure promotes binding of conformation-dependent myasthenia gravis autoantibodies, nicotinic acetylcholine receptor conformation maturation, and agonist sensitivity

The main immunogenic region (MIR) is a conformation-dependent region at the extracellular apex of alpha1 subunits of muscle nicotinic acetylcholine receptor (AChR) that is the target of half or more of the autoantibodies to muscle AChRs in human myasthenia gravis and rat experimental autoimmune myasthenia gravis. By making chimeras of human alpha1 subunits with alpha7 subunits, both MIR epitopes recognized by rat mAbs and by the patient-derived human mAb 637 to the MIR were determined to consist of two discontiguous sequences, which are adjacent only in the native conformation. The MIR, including loop alpha1 67-76 in combination with the N-terminal alpha helix alpha1 1-14, conferred high-affinity binding for most rat mAbs to the MIR. However, an additional sequence corresponding to alpha1 15-32 was required for high-affinity binding of human mAb 637. A water soluble chimera of Aplysia acetylcholine binding protein with the same alpha1 MIR sequences substituted was recognized by a majority of human, feline, and canine myasthenia gravis sera. The presence of the alpha1 MIR sequences in alpha1/alpha7 chimeras greatly promoted AChR expression and significantly altered the sensitivity to activation. This reveals a structural and functional, as well as antigenic, significance of the MIR.

Molecular-dynamics simulations of ELIC-a prokaryotic homologue of the nicotinic acetylcholine receptor

The ligand-gated ion channel from Erwinia chrysanthemi (ELIC) is a prokaryotic homolog of the eukaryotic nicotinic acetylcholine receptor (nAChR) that responds to the binding of neurotransmitter acetylcholine and mediates fast signal transmission. ELIC is similar to the nAChR in its primary sequence and overall subunit organization, but despite their structural similarity, it is not clear whether these two ligand-gated ion channels operate in a similar manner. Further, it is not known to what extent mechanistic insights gleaned from the ELIC structure translate to eukaryotic counterparts such as the nAChR. Here we use molecular-dynamics simulations to probe the conformational dynamics and hydration of the transmembrane pore of ELIC. The results are compared with those from our previous simulation of the human alpha7 nAChR. Overall, ELIC displays increased stability compared to the nAChR, whereas the two proteins exhibit remarkable similarity in their global motion and flexibility patterns. The majority of the increased stability of ELIC does not stem from the deficiency of the models used in the simulations, and but rather seems to have a structural basis. Slightly altered dynamical correlation features are also observed among several loops within the membrane region. In sharp contrast to the nAChR, ELIC is completely dehydrated from the pore center to the extracellular end throughout the simulation. Finally, the simulation of an ELIC mutant substantiates the important role of F246 on the stability, hydration and possibly function of the ELIC channel.

Role of acetylcholine receptor domains in ion selectivity

The nicotinic acetylcholine receptor (nAChR) is a ligand gated ion channel protein, composed of three domains: a transmembrane domain (TM-domain), extracellular domain (EC-domain), and intracellular domain (IC-domain). Due to its biological importance, much experimental and theoretical research has been carried out to explore its mechanisms of gating and selectivity, but there are still many unresolved issues, especially on the ion selectivity. Moreover, most of the previous theoretical work has concentrated on the TM-domain or EC-domain of nAChR, which may be insufficient to understand the entire structure-function relation. In this work, we perform molecular dynamics, Brownian dynamics simulations and continuum electrostatic calculations to investigate the role of different nAChR domains in ion conduction and selectivity. The results show that although both the EC and IC domains contain strong negative charges that create large cation concentrations at either end of the pore, this alone is not sufficient to create the observed cation selectivity and may play a greater role in determining the channel conductance. The presence of cations in the wide regions of the pore can screen out the protein charge allowing anions to enter, meaning that local regions of the TM-domain are most likely responsible for discriminating between ions. These new results complement our understanding about the ion conduction and selectivity mechanism of nAChR.

Binding of long-chain alpha-neurotoxin would stabilize the resting state of nAChR: a comparative study with alpha-conotoxin

Background

The details of interaction in a complex between potent antagonists such as long chain α-neurotoxins and α-conotoxins with nicotinic acetylcholine receptor (nAChR), and conformational changes induced by these antagonists, are not yet clear.

Modeling

In order to uncover some of these critical structural features, we conducted a docking simulation and a molecular dynamics simulation (MD) of a model of the ligand binding domain of nAChR in complex with a long-chain α-neurotoxin and an α-conotoxin.

Results

Our docking results confirm the claim that T.nAChR is in the basal or resting state, which favors binding to the alpha-neurotoxins. Moreover, more correct "hits" for the α/γ interface upon docking for conotoxin-nAChR confirm the preference of conotoxin GI for the α/γ interface. More importantly, upon binding of α-neurotoxin, ligand-bonded nAChR is less dynamic in certain domains than the apo form of the conotoxin-AChR complex. Some critical interactions in the binding site such as the salt bridge formed between K145/D200 in the neurotoxin-nAChR complex is further stabilized during the MD simulation, while it is obviously more labile in the apo form.

Conclusion

These observations could support the claim that alpha neurotoxins stabilize the nAChR resting state.

Nicotinic receptors containing the alpha7 subunit: a model for rational drug design

The neuronal nicotinic receptor has gained considerable recognition as a target, not just for combating drug addiction but also for treating a number of illnesses ranging from neurodegenerative diseases to psychotic disorders like schizophrenia. This recognition has led to a burgeoning field examining the receptor at all levels. A class of nicotinic receptors that contains the alpha7 gene product, apparently as a homomer, illustrates this multidisciplinary approach. Here, we review recent progress in our understanding of this class of receptors based on data from molecular, structural, physiological and patho-physiological studies. These studies have set the stage for rational drug design to combat disorders of the central nervous system. The studies also exemplify the cautious approach needed in developing CNS therapies and the importance of physiology in tempering drug design.

Myasthenia gravis and the tops and bottoms of AChRs: antigenic structure of the MIR and specific immunosuppression of EAMG using AChR cytoplasmic domains

The main immunogenic region (MIR), against which half or more of the autoantibodies to acetylcholine receptors (AChRs) in myasthenia gravis (MG) or experimental autoimmune MG (EAMG) are directed, is located at the extracellular end of alpha1 subunits. Rat monoclonal antibodies (mAbs) to the MIR efficiently compete with MG patient autoantibodies for binding to human muscle AChRs. Antibodies bound to the MIR do not interfere with cholinergic ligand binding or AChR function, but target complement and trigger antigenic modulation. Rat mAbs to the MIR also bind to human ganglionic AChR alpha3 subunits, but MG patient antibodies do not. By making chimeras of alpha1 subunits with alpha7 subunits or ACh binding protein, the structure of the MIR and its functional effects are being investigated. Many mAbs to the MIR bind only to the native conformation of alpha1 subunits because they bind to sequences that are adjacent only in the native structure. The MIR epitopes recognized by these mAbs are not recognized by most patient antibodies whose epitopes must be nearby. The presence of the MIR epitopes in alpha1/alpha7 chimeras greatly promotes AChR expression and sensitivity to activation. EAMG can be suppressed by treatment with denatured, bacterially expressed mixtures of extracellular and cytoplasmic domains of human alpha1, beta1, gamma, delta, and epsilon subunits. A mixture of only the cytoplasmic domains not only avoids the potential liability of provoking formation antibodies to pathologically significant epitopes on the extracellular surface, but also potently suppresses the development of EAMG.

Non-agonist-binding subunit interfaces confer distinct functional signatures to the alternate stoichiometries of the alpha4beta2 nicotinic receptor: an alpha4-alpha4 interface is required for Zn2+ potentiation

The alpha4beta2 subtype is the most abundant nicotinic acetylcholine receptor (nAChR) in the brain and possesses the high-affinity binding site for nicotine. The alpha4 and beta2 nAChR subunits assemble into two alternate stoichiometries, (alpha4)(2)(beta2)(3) and (alpha4)(3)(beta2)(2), which differ in their functional properties and sensitivity to chronic exposure to nicotine. Here, we investigated the sensitivity of both receptor stoichiometries to modulation by Zn2+. We show that Zn2+ exerts an inhibitory modulatory effect on (alpha4)(2)(beta2)(3) receptors, whereas it potentiates or inhibits, depending on its concentration, the function of (alpha4)(3)(beta2)(2) receptors. Furthermore, Zn2+ inhibition on (alpha4)(2)(beta2)(3) nAChRs is voltage-dependent, whereas it is not on (alpha4)(3)(beta2)(2) receptors. We used molecular modeling in conjunction with alanine substitution and functional studies to identify two distinct sets of residues that determine these effects and may coordinate Zn(2+). Zn(2+) inhibition is mediated by a site located on the beta2(+)/alpha4(-) subunit interfaces on both receptor stoichiometries. alpha4(H195) and beta2(D218) are key determinants of this site. Zn2+ potentiation on (alpha4)(3)(beta2)(2) nAChRs is exerted by a site that resides on the alpha4(+)/alpha4(-) of this receptor stoichiometry. alpha4(H195) on the (-) side of the ACh-binding alpha4 subunit and alpha4(E224) on the (+) side of the non-ACh-binding alpha4 subunit critically contribute to this site. We also identified residues within the beta2 subunit that confer voltage dependency to Zn2+ inhibition on (alpha4)(2)(beta2)(3), but not on (alpha4)(3)(beta2)(2) nAChRs.

Spontaneous conformational change and toxin binding in alpha7 acetylcholine receptor: insight into channel activation and inhibition

Nicotinic AChRs (nAChRs) represent a paradigm for ligand-gated ion channels. Despite intensive studies over many years, our understanding of the mechanisms of activation and inhibition for nAChRs is still incomplete. Here, we present molecular dynamics (MD) simulations of the alpha7 nAChR ligand-binding domain, both in apo form and in alpha-Cobratoxin-bound form, starting from the respective homology models built on crystal structures of the acetylcholine-binding protein. The toxin-bound form was relatively stable, and its structure was validated by calculating mutational effects on the toxin-binding affinity. However, in the apo form, one subunit spontaneously moved away from the conformation of the other four subunits. This motion resembles what has been proposed for leading to channel opening. At the top, the C loop and the adjacent beta7-beta8 loop swing downward and inward, whereas at the bottom, the F loop and the C terminus of beta10 swing in the opposite direction. These swings appear to tilt the whole subunit clockwise. The resulting changes in solvent accessibility show strong correlation with experimental results by the substituted cysteine accessibility method upon addition of acetylcholine. Our MD simulation results suggest a mechanistic model in which the apo form, although predominantly sampling the "closed" state, can make excursions into the "open" state. The open state has high affinity for agonists, leading to channel activation, whereas the closed state upon distortion has high affinity for antagonists, leading to inhibition.

Structural answers and persistent questions about how nicotinic receptors work

The electron diffraction structure of nicotinic acetylcholine receptor (nAChR) from Torpedo marmorata and the X-ray crystallographic structure of acetylcholine binding protein (AChBP) are providing new answers to persistent questions about how nAChRs function as biophysical machines and as participants in cellular and systems physiology. New high-resolution information about nAChR structures might come from advances in crystallography and NMR, from extracellular domain nAChRs as high fidelity models, and from prokaryotic nicotinoid proteins. At the level of biophysics, structures of different nAChRs with different pharmacological profiles and kinetics will help describe how agonists and antagonists bind to orthosteric binding sites, how allosteric modulators affect function by binding outside these sites, how nAChRs control ion flow, and how large cytoplasmic domains affect function. At the level of cellular and systems physiology, structures of nAChRs will help characterize interactions with other cellular components, including lipids and trafficking and signaling proteins, and contribute to understanding the roles of nAChRs in addiction, neurodegeneration, and mental illness. Understanding nAChRs at an atomic level will be important for designing interventions for these pathologies.

A molecular basis for agonist and antagonist actions at GABA(C) receptors

We modelled the N-terminal ligand-binding domain of the rho1 GABA(C) receptor based on the Lymnaea stagnalis acetylcholine-binding protein (L-AChBP) crystal structure using comparative modelling and validated using flexible docking guided by known mutagenesis studies. A range of known rho1 GABA(C) receptor ligands comprising seven full agonists, 10 partial agonists, 43 antagonists and 12 inactive molecules were used to evaluate and validate the models. Of the 50 models identified, six models that allowed flexible ligand docking in accordance with the experimental data were selected and used to study detailed receptor-ligand interactions. The most refined model to accommodate all known active ligands featured a cavity comprising of a volume of 488 A(3). A detailed analysis of the interaction between the rho1 GABA(C) receptor model and the docked ligands revealed possible H-bonds and cation-pi interactions between the different ligands and binding site residues. Based on quantum mechanical/molecular mechanical (QM/MM) calculations, the model showed distinctive conformations of loop C that provided a molecular basis for agonist and antagonist actions. Agonists elicit loop C closure, while a more open loop C was observed upon antagonist binding. The model differentiates the role for key residues known to be involved in either binding and/or gating.

Dynamics of heteropentameric nicotinic acetylcholine receptor: implications of the gating mechanism

The dynamics characteristics of the currently available structure of Torpedo nicotinic acetylcholine receptor (nAChR), including the extracellular, transmembrane, and intracellular domains (ICDs), were analyzed using the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM). We found that a symmetric quaternary twist motion, reported previously in the literature in a homopentameric receptor (Cheng et al. J Mol Biol 2006;355:310-324; Taly et al. Biophys J 2005;88:3954-3965), occurred also in the heteropentameric Torpedo nAChR. We believe, however, that the symmetric twist alone is not sufficient to explain a large body of experimental data indicating asymmetry and subunit nonequivalence during gating. Here we report our results supporting the hypothesis that a combination of symmetric and asymmetric motions opens the gate. We show that the asymmetric motion involves tilting of the TM2 helices. Furthermore, our study reveals three additional aspects of channel dynamics: (1) loop A serves as an allosteric mediator between the ligand binding loops and those at the domain interface, particularly the linker between TM2 and TM3; (2) the ICD can modulate the pore dynamics and thus should not be neglected in gating studies; and (3) the F loops, which are peculiarly longer and poorly-conserved in non-alpha-subunits, have important dynamical implications.

Homology modeling and molecular dynamics simulations of the alpha1 glycine receptor reveals different states of the channel

Homology modeling is used to build initial models of the transmembrane domain of the human alpha1 glycine receptor (GlyR) based on the most recently published refined structure of nAChR (PDB ID: 2BG9). Six preliminary GlyR models are constructed using two different approaches. In one approach, five different homopentamers are built by symmetric assembly of alpha1 GlyR subunits using only one of the five unique chains of nAChR as a template. In a second approach, each nAChR subunit serves as a template for an alpha1 GlyR subunit. All six initial GlyR constructs are then embedded into a hydrated POPC lipid bilayer and subjected to molecular dynamics simulation for at least six nanoseconds. Each model is stable throughout the simulation, and the final models fall into three distinct categories. Homopentameric GlyR bundles using a single alpha nAChR subunit as a template appear to be in an open conformation. Under an applied external potential, permeation of Cl(-) ions is observed within several ns in a channel built on an alpha chain. Model channels built on non-alpha chains have a constriction either near the intracellular mouth or more centrally located in the pore domain, both of which may be narrow enough to close the channel and whose locations correspond to putative gates observed in nicotinicoid receptors. The differences between these three general models suggest that channel closure may be effected by either rotation or tangential tilting of TM2.

Nanosecond-timescale conformational dynamics of the human alpha7 nicotinic acetylcholine receptor

We explore the conformational dynamics of a homology model of the human alpha7 nicotinic acetylcholine receptor using molecular dynamics simulation and analyses of root mean-square fluctuations, block partitioning of segmental motion, and principal component analysis. The results reveal flexible regions and concerted global motions of the subunits encompassing extracellular and transmembrane domains of the subunits. The most relevant motions comprise a bending, hinged at the beta10-M1 region, accompanied by concerted tilting of the M2 helices that widens the intracellular end of the channel. Despite the nanosecond timescale, the observations suggest that tilting of the M2 helices may initiate opening of the pore. The results also reveal direct coupling between a twisting motion of the extracellular domain and dynamic changes of M2. Covariance analysis of interresidue motions shows that this coupling arises through a network of residues within the Cys and M2-M3 loops where Phe135 is stabilized within a hydrophobic pocket formed by Leu270 and Ile271. The resulting concerted motion causes a downward shift of the M2 helices that disrupts a hydrophobic girdle formed by 9' and 13' residues.

Central Nicotinic Receptors: Structure, Function, Ligands, and Therapeutic Potential

The growing interest in nicotinic receptors, because of their wide expression in neuronal and non-neuronal tissues and their involvement in several important CNS pathologies, has stimulated the synthesis of a high number of ligands able to modulate their function. These membrane proteins appear to be highly heterogeneous, and still only incomplete information is available on their structure, subunit composition, and stoichiometry. This is due to the lack of selective ligands to study the role of nAChR under physiological or pathological conditions; so far, only compounds showing selectivity between alpha4beta2 and alpha7 receptors have been obtained. The nicotinic receptor ligands have been designed starting from lead compounds from natural sources such as nicotine, cytisine, or epibatidine, and, more recently, through the high-throughput screening of chemical libraries. This review focuses on the structure of the new agonists, antagonists, and allosteric ligands of nicotinic receptors, it highlights the current knowledge on the binding site models as a molecular modeling approach to design new compounds, and it discusses the nAChR modulators which have entered clinical trials.

Homology modeling and molecular dynamics simulations of the glycine receptor ligand binding domain

We present a homology based model of the ligand binding domain (LBD) of the homopentameric alpha1 glycine receptor (GlyR). The model is based on multiple sequence alignment with other members of the nicotinicoid ligand gated ion channel superfamily and two homologous acetylcholine binding proteins (AChBP) from the freshwater (Lymnaea stagnalis) and saltwater (Aplysia californica) snails with known high resolution structure. Using two template proteins with known structure to model three dimensional structure of a target protein is especially advantageous for sequences with low homology as in the case presented in this paper. The final model was cross-validated by critical evaluation of experimental and published mutagenesis, functional and other biochemical studies. In addition, a complex structure with strychnine antagonist in the putative binding site is proposed based on docking simulation using Autodock program. Molecular dynamics (MD) simulations with simulated annealing protocol are reported on the proposed LBD of GlyR, which is stable in 5 ns simulation in water, as well as for a deformed LBD structure modeled on the corresponding domain determined in low-resolution cryomicroscopy structure of the alpha subunit of the full-length acetylcholine receptor (AChR). Our simulations demonstrate that the beta-sandwich central core of the protein monomer is fairly rigid in the simulations and resistant to deformations in water.

An intersubunit trigger of channel gating in the muscle nicotinic receptor

Binding of neurotransmitter triggers gating of synaptic receptor channels, but our understanding of the structures that link the binding site to the channel is just beginning to develop. Here, we identify an intersubunit triggering element required for rapid and efficient gating of muscle nicotinic receptors using a structural model of the Torpedo receptor at 4 A resolution, recordings of currents through single receptor channels, measurements of inter-residue energetic coupling, and functional consequences of disulfide trapping. Mutation of the conserved residues, alphaTyr 127, epsilonAsn 39, and deltaAsn 41, located at the two subunit interfaces that form the agonist binding sites, markedly attenuates acetylcholine-elicited channel gating; mutant cycle analyses based on changes in the channel gating equilibrium constant reveal strong energetic coupling among these residues. After each residue is substituted with Cys, oxidizing conditions that promote disulfide bond formation attenuate gating of mutant, but not wild-type receptors. Gating is similarly attenuated when the Cys substitutions are confined to either of the binding-site interfaces, but can be restored by reducing conditions that promote disulfide bond breakage. Thus, the Tyr-Asn pair is an intersubunit trigger of rapid and efficient gating of muscle nicotinic receptors.

Chemical-scale studies on the role of a conserved aspartate in preorganizing the agonist binding site of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor and related Cys-loop receptors are ligand-gated ion channels that mediate fast synaptic transmission throughout the central and peripheral nervous system. A highly conserved aspartate residue (D89) that is near the agonist binding site but does not directly contact the ligand plays a critical part in receptor function. Here we probe the role of D89 using unnatural amino acid mutagenesis coupled with electrophysiology. Homology modeling implicates several hydrogen bonds involving D89. We find that no single hydrogen bond is essential to proper receptor function. Apparently, the side chain of D89 establishes a redundant network of hydrogen bonds; these bonds preorganize the agonist binding site by positioning a critical tryptophan residue that directly contacts the ligand. Earlier studies of the D89N mutant led to the proposal that a negative charge at this position is essential for receptor function. However, we find that receptors with neutral side chains at position 89 can function well, if the side chain is less perturbing than the amide of asparagine (nitro or keto groups allow function) or if a compensating backbone mutation is introduced to relieve unfavorable electrostatics.

Modeling Subtype-Selective Agonists Binding with α4β2 and α7 Nicotinic Acetylcholine Receptors: Effects of Local Binding and Long-Range Electrostatic Interactions

The subtype-selective binding of 14 representative agonists with alpha4beta2 and alpha7 nicotinic acetylcholine receptors (nAChRs) has been studied by performing homology modeling, molecular docking, geometry optimizations, and microscopic and phenomenological binding free energy calculations. All of the computational results demonstrate that the subtype selectivity of the agonists binding with alpha4beta2 and alpha7 7 nAChRs is affected by both local binding and long-range electrostatic interactions between the receptors and the protonated structures of the agonists. The effects of the long-range electrostatic interactions are mainly due to the distinct difference in the net charge of the ligand-binding domain between the two nAChR subtypes. For the alpha4beta2-selective agonists examined, the microscopic binding modes with the alpha4beta2 nAChR are very similar to the corresponding modes with the alpha7 nAChR, and therefore, the subtype selectivity of these agonists binding with alpha4beta2 and alpha7 nAChRs is dominated by the long-range electrostatic interactions. For the alpha7-selective agonists, their microscopic binding modes with the alpha7 nAChR are remarkably different from those with the alpha4beta2 nAChR so that the local binding (including the hydrogen bonding and cation-pi interactions) with the alpha7 nAChR is much stronger than that with the alpha4beta2 nAChR. The calculated phenomenological binding free energies are in good agreement with available experimental data for the relative binding free energies concerning the subtype selectivity of agonists binding with the two different nAChR subtypes. The fundamental insights obtained in the present study should be valuable for future rational design of potential therapeutic agents targeted to specific nAChR subtypes.

Blocking of the Nicotinic Acetylcholine Receptor Ion Channel by Chlorpromazine, a Noncompetitive Inhibitor: A Molecular Dynamics Simulation Study

A large series of pharmacological agents, distinct from the typical competitive antagonists, block in a noncompetitive manner the permeability response of the nicotinic acetylcholine receptor (nAChR) to the neurotransmitter acetylcholine. Taking the neuroleptic chlorpromazine (CPZ) as an example of such agents, the blocking mechanism of noncompetitive inhibitors to the ion channel pore of the nAChR has been explored at the atomic level using both conventional and steered molecular dynamics (MD) simulations. Repeated steered MD simulations have permitted calculation of the free energy (approximately 36 kJ/mol) of CPZ binding and identification of the optimal site in the region of the serine and leucine rings, at approximately 4 A from the pore entrance. Coulomb and the Lennard-Jones interactions between CPZ and the ion channel as well as the conformational fluctuations of CPZ were examined to assess the contribution of each to the binding of CPZ to the nAChR. The MD simulations disclose a dynamic interaction of CPZ binding to the nAChR ionic channel. The cationic ammonium head of CPZ forms strong hydrogen bonds with Glu262 (alpha), Asp268 (beta), Glu272 (beta), Ser276 (beta), Glu280 (delta), Gln271 (gamma), Glu275 (gamma), and Asn279 (gamma) nAChR residues. Finally, the conventional MD simulation of CPZ at its identified binding site demonstrates that the binding of CPZ not only blocks ion transport through the channel but also markedly inhibits the conformational transitions of the channel, necessary for nAChR to carry out its biological function.

Solution NMR of Acetylcholine Binding Protein Reveals Agonist-Mediated Conformational Change of the C-Loop

Previous X-ray crystallography, molecular dynamics simulation, fluorescence spectroscopy, and deuterium-hydrogen exchange of acetylcholine binding protein (AChBP) suggest that after binding of the agonist, the C-loop at the periphery of the binding site draws inward to cap the site and envelop the agonist. In this study, we use high-resolution solution NMR to monitor changes in the chemical environment of the C-loop without and with acetylcholine (ACh) bound. Substitution of [15N]cysteine for the native cysteines 123, 136, 187, and 188 provided intrinsic monitors of the chemical environments of the Cys- and C-loops, respectively. Two-dimensional transverse relaxation-optimized spectroscopy 15N-1H HSQC spectroscopy of apo-AChBP revealed seven well resolved cross-peaks for the group of cysteines. The spectrum of AChBP with Ser substituted for Cys 187 and 188 shows only two main cross-peaks, corresponding to Cys 123 and 136 from the Cys-loop, enabling resonance assignments. After binding of ACh, the five cross-peaks associated with cysteines from the C-loop condense into two predominant cross-peaks not observed in the spectrum from the apo protein, indicating a restricted range of conformations and change in chemical environment of the C-loop. The results show that isotopic cysteine can be incorporated into specified positions of AChBP expressed from a eukaryotic source, that the C-loop assumes multiple conformations without ACh, but that its conformation becomes restricted with ACh bound. The collective findings suggest a structural mechanism for agonist recognition in AChBP and related Cys-loop receptors.

Fully flexible docking models of the complex between α7 nicotinic receptor and a potent heptapeptide inhibitor of the β-amyloid peptide binding

The heptapeptide IQTTWSR (IQ), recently reported as inhibitor of the beta-amyloid (Abeta) binding to nicotinic acetylcholine receptors (nAChrs), was docked to the homology model of the alpha7 nicotinic acetylcholine receptor. The most representative models were further subjected to molecular dynamics simulations. The data obtained here suggest that Abeta needs highly specific structural motifs to bind to the alpha7nAChR. These structural motifs are located principally in the upper and lower surroundings of loop C, including loop F and sheets beta1, beta2, beta6, beta9, and beta10 of the receptor. Overall, these results suggest that IQ can be mimicked by more bioavailable, stable compounds that would be helpful for the understanding of the Abeta binding site and its dynamics, and for the design of novel agents to be used as an effective alternative against Alzheimer's disease.

Modulating inhibitory ligand-gated ion channels

The glycine and gamma-aminobutyric acid receptors (GlyR and GABA(A)R, respectively) are the major inhibitory neurotransmitter-gated receptors in the central nervous system of animals. Given the important role of these receptors in neuronal inhibition, they are prime targets of many therapeutic agents and are the object of intense studies aimed at correlating their structure and function. In this review, the structure and dynamics of these and other homologous members of the nicotinicoid superfamily are described. The modulatory actions of the major biological macromolecules that bind and allosterically affect these receptors are also discussed.

Potentials of mean force for acetylcholine unbinding from the alpha7 nicotinic acetylcholine receptor ligand-binding domain

The nicotinic acetylcholine receptor is a prototype ligand-gated ion channel that mediates signal transduction in the neuromuscular junctions and other cholinergic synapses. The molecular basis for the energetics of ligand binding and unbinding is critical to our understanding of the pharmacology of this class of receptors. Here, we used steered molecular dynamics to investigate the unbinding of acetylcholine from the ligand-binding domain of human alpha7 nicotinic acetylcholine receptor along four different predetermined pathways. Pulling forces were found to correlate well with interactions between acetylcholine and residues in the binding site during the unbinding process. From multiple trajectories along these unbinding pathways, we calculated the potentials of mean force for acetylcholine unbinding. Four available methods based on Jarzynski's equality were used and compared for their efficiencies. The most probable pathway was identified to be along a direction approximately parallel to the membrane. The derived binding energy for acetylcholine was in good agreement with that derived from the experimental binding constant for acetylcholine binding protein, but significantly higher than that for the complete human alpha7 nicotinic acetylcholine receptor. In addition, it is likely that several intermediate states exist along the unbinding pathways.