Lysine scanning mutagenesis delineates structural model of the nicotinic receptor ligand binding domain

The Surface of Protein λ6-85 Can Act as a Template for Recurring Poly(ethylene glycol) Structure

PEGylated proteins play an increasingly important role in pharmaceutical drug delivery. We recently showed that short poly(ethylene glycol) (PEG) chains can affect protein structure, even when they are not making extensive contact with the protein surface. In contrast, PEG is generally thought to form a relatively unstructured coil, and its compactness depends on solvent conditions. Here we test whether a host protein could allow PEG to form recurrent structural motifs while the PEG chain is in contact with the protein surface. We link a PEG oligomer (n = 45) to one of two nearly opposite locations on the small α-helical protein λ_6-85 to investigate this question. We first demonstrate experimentally that in these particular positions, PEG does not significantly affect the thermodynamic stability or folding kinetics of λ_6-85. We then use several all-atom molecular dynamics (MD) simulations 1 μs in duration to show how PEG equilibrates between states extending into the solvent and states packed onto the protein surface. The packing reveals recurring structures, including persistent hydrogen bond and hydrophobic contact patterns that appear multiple times. Some interactions of PEG with surface lysines are best described as an "intermittent slithering" motion of the PEG around the side chain, as seen in short MD movies. Thus, PEG achieves a variety of metastable organized structures on the protein surface, somewhere between a random globule and true folding. We also investigated the PEG-protein interaction in the unfolded state of the protein. We find that PEG has a propensity to stabilize certain helices of λ_6-85, no matter which of the two positions it was attached to. Thus, sufficiently long PEG chains are organized by the protein surface and in turn interact with certain elements of protein structure more than others, even when PEG is attached to very different sites.

Investigation of Congenital Myasthenia Reveals Functional Asymmetry of Invariant Acetylcholine Receptor (AChR) Cys-loop Aspartates

We identify two heteroallelic mutations in the acetylcholine receptor δ-subunit from a patient with severe myasthenic symptoms since birth: a novel δD140N mutation in the signature Cys-loop and a mutation in intron 7 of the δ-subunit gene that disrupts splicing of exon 8. The mutated Asp residue, which determines the disease phenotype, is conserved in all eukaryotic members of the Cys-loop receptor superfamily. Studies of the mutant acetylcholine receptor expressed in HEK 293 cells reveal that δD140N attenuates cell surface expression and apparent channel gating, predicting a reduced magnitude and an accelerated decay of the synaptic response, thus reducing the safety margin for neuromuscular transmission. Substituting Asn for Asp at equivalent positions in the α-, β-, and ϵ-subunits also suppresses apparent channel gating, but the suppression is much greater in the α-subunit. Mutant cycle analysis applied to single and pairwise mutations reveals that αAsp-138 is energetically coupled to αArg-209 in the neighboring pre-M1 domain. Our findings suggest that the conserved αAsp-138 and αArg-209 contribute to a principal pathway that functionally links the ligand binding and pore domains.

Crosslinking constraints and computational models as complementary tools in modeling the extracellular domain of the glycine receptor

The glycine receptor (GlyR), a member of the pentameric ligand-gated ion channel superfamily, is the major inhibitory neurotransmitter-gated receptor in the spinal cord and brainstem. In these receptors, the extracellular domain binds agonists, antagonists and various other modulatory ligands that act allosterically to modulate receptor function. The structures of homologous receptors and binding proteins provide templates for modeling of the ligand-binding domain of GlyR, but limitations in sequence homology and structure resolution impact on modeling studies. The determination of distance constraints via chemical crosslinking studies coupled with mass spectrometry can provide additional structural information to aid in model refinement, however it is critical to be able to distinguish between intra- and inter-subunit constraints. In this report we model the structure of GlyBP, a structural and functional homolog of the extracellular domain of human homomeric α1 GlyR. We then show that intra- and intersubunit Lys-Lys crosslinks in trypsinized samples of purified monomeric and oligomeric protein bands from SDS-polyacrylamide gels may be identified and differentiated by MALDI-TOF MS studies of limited resolution. Thus, broadly available MS platforms are capable of providing distance constraints that may be utilized in characterizing large complexes that may be less amenable to NMR and crystallographic studies. Systematic studies of state-dependent chemical crosslinking and mass spectrometric identification of crosslinked sites has the potential to complement computational modeling efforts by providing constraints that can validate and refine allosteric models.

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.

Accommodation of a central arginine in a transmembrane peptide by changing the placement of anchor residues

Both Trp and Arg in transmembrane protein domains make important interactions with lipids at the membrane/water interface, but at different depths. Derivatives of the designed peptide GWALP23, acetyl-GGALW(5)LALALALALALALW(19)LAGA-amide, with single Trp anchors, have proven useful for characterizing such interactions. Indeed, previous work revealed quite different effects emanating from Arg substitutions at positions 12 and 14 within GWALP23, with the R12 peptide exhibiting multiple positions and orientations with respect to DOPC bilayer membranes (Vostrikov et al. J. Am. Chem. Soc. 2010, 132, 5803-5811). To gain further understanding of the multistate behavior, we moved the Trp "anchor" residues to more outer positions 3 and 21 in GWALP23 itself, and in the R12 and R14 derivatives. The locations and orientations of the peptides with respect to lipid bilayer membranes of differing thickness were investigated by means of solid-state (2)H NMR spectroscopy, using labeled alanines, and coarse-grained molecular dynamics simulations. Interestingly, relatively intense and narrow (2)H resonances from selected backbone C(α) deuterons were observed over quite narrow ranges of frequency and sample orientation. The backbone resonances reflect dynamic complexities and at the same time provide important contributions for the analysis of peptide transmembrane orientation. With the Trp(3,21) anchors relatively far from the peptide and bilayer center, the results indicate significantly large apparent tilt angles, for example, close to 30° for the new R12 and R14 peptides with respect to the bilayer normal of DLPC membranes. The R12 side chain indeed is "rescued" to a stable position, where it is accommodated within the transmembrane helix, when the Trp anchors are moved outward and to another face of the helix. At the same time, the R14 side chain of transmembrane GW(3,21)ALP23 also retains a stable favored position.

Highly fatal fast-channel syndrome caused by AChR ε subunit mutation at the agonist binding site

OBJECTIVE

To characterize the molecular basis of a novel fast-channel congenital myasthenic syndrome.

METHODS

We used the candidate gene approach to identify the pathogenic mutation in the acetylcholine receptor (AChR) ε subunit, genetically engineered the mutant AChR into HEK cells, and evaluated the level of expression and kinetic properties of the mutant receptor.

RESULTS

An 8-year-old boy born to consanguineous parents had severe myasthenic symptoms since birth. He is wheelchair bound and pyridostigmine therapy enables him to take only a few steps. Three similarly affected siblings died in infancy. He carries a homozygous p.W55R mutation at the α/ε subunit interface of the AChR agonist binding site. The mutant protein expresses well in HEK cells. Patch-clamp analysis of the mutant receptor expressed in HEK cells reveals 30-fold reduced apparent agonist affinity, 75-fold reduced apparent gating efficiency, and strikingly attenuated channel opening probability (P(open)) over a range agonist concentrations.

CONCLUSION

Introduction of a cationic Arg into the anionic environment of α/ε AChR binding site hinders stabilization of cationic ACh by aromatic residues and accounts for the markedly perturbed kinetic properties of the receptor. The very low P(open) explains the poor response to pyridostigmine and the high fatality of the disease.

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.

Receptor signaling and the cell biology of synaptic transmission

This volume describes a series of psychiatric and neuropsychiatric disorders, connects some aspects of somatic and psychiatric medicine, and describes various current and emerging therapies. The purpose of this chapter is to set the stage for the volume by developing the theoretical basis of synaptic transmission and introducing the various neurotransmitters and their receptors involved in the process. The intent is to provide not only a historical context through which to understand neurotransmitters, but a current contextual basis for understanding neuronal signal transduction and applying this knowledge to facilitate treatment of maladies of the brain and mind.

Intramembrane proton binding site linked to activation of bacterial pentameric ion channel

Prokaryotic orthologs of eukaryotic Cys-loop receptor channels recently emerged as structural and mechanistic surrogates to investigate this superfamily of intercellular signaling proteins. Here, we examine proton activation of the prokaryotic ortholog GLIC using patch clamp electrophysiology, mutagenesis, and molecular dynamics (MD) simulations. Whole-cell current recordings from human embryonic kidney (HEK) 293 cells expressing GLIC show half-maximal activation at pH 6, close to the pK(a) of histidine, implicating the three native His residues in proton sensing linked to activation. The mutation H235F abolishes proton activation, H277Y is without effect, and all nine mutations of His-127 prevent expression on the cell surface. In the GLIC crystal structure, His-235 on transmembrane (TM) α-helix 2, hydrogen bonds to the main chain carbonyl oxygen of Ile-259 on TM α-helix 3. MD simulations show that when His-235 is protonated, the hydrogen bond persists, and the channel remains in the open conformation, whereas when His-235 is deprotonated, the hydrogen bond dissociates, and the channel closes. Mutations of the proximal Tyr-263, which also links TM α-helices 2 and 3 via a hydrogen bond, alter proton sensitivity over a 1.5 pH unit range. MD simulations show that mutations of Tyr-263 alter the hydrogen bonding capacity of His-235. The overall findings show that His-235 in the TM region of GLIC is a novel proton binding site linked to channel activation.

The use of site-directed mutagenesis to study GPCRs

G protein coupled receptors (GPCRs) are highly flexible and dynamic proteins, which are able to interact with diverse ligands, effectors, and regulatory proteins. Site-directed mutagenesis (SDM) is a powerful tool for providing insight into how these proteins actually work, both in its own right and when used in conjunction with information provided by other techniques such as crystallography or molecular modelling. Mutagenesis has been used to identify and characterise a myriad of functionally important residues, motifs and domains within the GPCR architecture, and to identify aspects of similarity and differences between the major families of GPCRs. This chapter presents the necessary information for undertaking informative SDM of these proteins. Whilst this is relevant to protein structure/function studies in -general, specific pitfalls and protocols suited to investigating GPCRs in particular will be highlighted.

The positive allosteric modulator morantel binds at noncanonical subunit interfaces of neuronal nicotinic acetylcholine receptors

We are interested in the positive allosteric modulation of neuronal nicotinic acetylcholine (ACh) receptors and have recently shown that the anthelmintic compound morantel potentiates by enhancing channel gating of the alpha3beta2 subtype. Based on the demonstration that morantel-elicited currents were inhibited by the classic ACh competitor dihydro-beta-erythroidine in a noncompetitive manner and that morantel still potentiates at saturating concentrations of agonist (Wu et al., 2008), we hypothesized that morantel binds at the noncanonical beta2(+)/alpha3(-) subunit interface. In the present study, we created seven cysteine-substituted subunits by site-directed mutagenesis, choosing residues in the putative morantel binding site with the aid of structural homology models. We coexpressed the mutant subunits and their respective wild-type partners in Xenopus oocytes and characterized the morantel potentiation of ACh-evoked currents, as well as morantel-evoked currents, before and after treatment with a variety of methanethiosulfonate (MTS)-based compounds, using voltage-clamp recordings. The properties of four of the seven mutants, two residues on each side of the interface, were changed by MTS treatments. Coapplication with ACh enhanced the extent of MTS modification for alpha3A106Cbeta2 and alpha3beta2S192C receptors. The activities of two mutants, alpha3T115Cbeta2 and alpha3beta2T150C, were dramatically altered by MTS modification. For alpha3beta2T150C, while peak current amplitudes were reduced, potentiation was enhanced. For alpha3T115Cbeta2, both current amplitudes and potentiation were reduced. MTS modification and morantel were mutually inhibitory: MTS treatment decreased morantel-evoked currents and morantel decreased the rate of MTS modification. We conclude that the four residues showing MTS effects contribute to the morantel binding site.

Role of intracellular domains in the function of the herg potassium channel

The functional role of the large intracellular regions (which include the cyclic nucleotide binding domain, cNBD, and the Per-Arnt-Sim domain, PAS) in the herg channel is not well understood. We have studied possible interactions of the cNBD with other parts of the channel protein using lysine mutations to disrupt such interactions. Some lysine mutations caused significant right shifts in the voltage dependence of inactivation; almost all the mutants caused speeding up of deactivation time course. In a homology model of the cNBD, lysine mutations that affected both inactivation and deactivation lie in a hydrophobic band on the surface of the structure of this domain. Some known mutations in the Long QT Syndrome type 2, with effects on deactivation, are located at residues close to hydrophobic bands on the cNBD and the PAS domains. Such bands of residues in these intracellular domains may play an important part in channel function.

Tryptophan scanning of the acetylcholine receptor's betaM4 transmembrane domain: decoding allosteric linkage at the lipid-protein interface with ion-channel gating

The nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel protein that mediates fast excitatory synaptic transmission in the peripheral and central nervous systems. Changes in the structure and function of the AChR can lead to serious impairment of physiological processes. In this study, we combined site-directed mutagenesis, radioligand binding assays, electrophysiological recordings and Fourier analyses to characterize the functional role and structural aspects of the betaM4 transmembrane domain of the Torpedo AChR. We performed tryptophan replacements, from residues L438 through F455, along the betaM4 transmembrane domain. Expression levels of mutants F439W-G450W and F452W-I454W produced peak currents similar to or lower than those in wild-type (WT). Tryptophan substitutions at positions L438 and T451 led to a deficiency in either subunit expression or receptor assembly. Mutations L440W, V442W, C447W and S453W produced a gain-of-function response. Mutation F455W produced a loss of ion channel function. The periodicity profile of the normalized expression level (closed state) and EC(50) (open state) revealed a minor conformational change of 0.4 residues/turn of the betaM4 domain. These findings suggest that a minor movement of the betaM4 domain occurs during channel activation.

Multiple modes for conferring surface expression of homomeric beta1 GABAA receptors

The gamma-aminobutyric acid type A (GABA(A)) receptor assembles from individual subunits to form ligand-gated ion channels. Human (h) beta3 subunits assemble to form homomeric surface receptors in somatic cells, but hbeta1 subunits do not. We have identified three distinct sets of amino acid residues in the N-terminal extracellular domain of the hbeta1 subunit, which when mutated to the homologous residue in hbeta3 allow expression as a functional homomeric receptor. The three sets likely result in three modes of assembly. Mode 1 expression results from a single amino acid change at residue hbeta1 Asp-37. Mode 2 expression results from mutations of residues between positions 44 and 73 together with residues between positions 169 and 173. Finally, mode 3 results from the mutations A45V and K196R. Examination of homology-based structural models indicates that many of the residues are unlikely to be involved in physical inter-subunit interactions, suggesting that a major alteration is stabilization of an assembly competent form of the subunit. These mutations do not, however, have a major effect on the surface expression of heteromeric receptors which include the alpha1 subunit.

Control of cation permeation through the nicotinic receptor channel

We used molecular dynamics (MD) simulations to explore the transport of single cations through the channel of the muscle nicotinic acetylcholine receptor (nAChR). Four MD simulations of 16 ns were performed at physiological and hyperpolarized membrane potentials, with and without restraints of the structure, but all without bound agonist. With the structure unrestrained and a potential of −100 mV, one cation traversed the channel during a transient period of channel hydration; at −200 mV, the channel was continuously hydrated and two cations traversed the channel. With the structure restrained, however, cations did not traverse the channel at either membrane potential, even though the channel was continuously hydrated. The overall results show that cation selective transport through the nAChR channel is governed by electrostatic interactions to achieve charge selectivity, but ion translocation relies on channel hydration, facilitated by a trans-membrane field, coupled with dynamic fluctuations of the channel structure.

Communication between a cell and its environment relies on channel-forming proteins to provide a low energy pathway for ions to move in and out. Although channel-forming proteins are essential to all life forms, the atomic-scale mechanisms that enable ions to pass through the channel remain elusive due to the lack of experimental approaches to monitor the protein and ion in real time and at atomic resolution. A powerful alternative approach is molecular dynamics (MD) simulation based on the laws of physics applied to the increasing body of protein structures resolved at atomic resolution. Here we present all-atom MD simulations applied to the nicotinic acetylcholine receptor (nAChR) that initiates voluntary movement in skeletal muscle. By focusing on individual permeant cations, we find that selective cation translocation occurs in stages: cations are first selected through a series of oppositely charged residues within the protein vestibule leading to a narrow hydrophobic constriction, but then hydration of the narrow region and dynamic fluctuations of the protein enable the cation to pass through. The findings provide a general framework for understanding how ions are selected for transport based on charge, and how the dynamic interplay between water, the ion, and the channel protein enable rapid ion translocation through the broad class of channel-forming proteins with hydrophobic barriers.

The kinetics of competitive antagonism of nicotinic acetylcholine receptors at physiological temperature

Detailed information about the ligand-binding site of nicotinic acetylcholine receptors has emerged from structural and mutagenesis experiments. However, these approaches provide only static images of ligand-receptor interactions. Kinetic measurements of changes in protein function are needed to develop a more dynamic picture. Previously, we measured association and dissociation rate constants for competitive inhibition of current through embryonic muscle acetylcholine receptor channels at 25 degrees C. Little is known about competitive antagonism at physiological temperatures. Here, we performed measurements at 37 degrees C and used thermodynamics to estimate the energetics of antagonism. We used rapid solution exchange protocols to determine equilibrium and kinetics of inhibition of acetylcholine-activated currents in outside-out patches by (+)-tubocurarine, pancuronium and cisatracurium. Kinetic rates as high as 600 s(-1) were resolved by this technique. Binding was primarily enthalpy driven. The 12 degrees C increase in temperature decreased equilibrium antagonist binding by 1.7- to 1.9-fold. In contrast, association and dissociation rate constants increased 1.9- to 6.0-fold. Activation energies for dissociation were 90 +/- 6, 106 +/- 8 and 116 +/- 10 kJ mol(-1) for cisatracurium, (+)-tubocurarine and pancuronium, respectively. The corresponding apparent activation energies for association were 38 +/- 6, 85 +/- 6 and 107 +/- 13 kJ mol(-1). The higher activation energy for association of (+)-tubocurarine and pancuronium compared with cisatracurium is notable. This may arise from either a more superficial binding site for the large antagonist cisatracurium compared to the other ligands, or from a change in receptor conformation upon binding of (+)-tubocurarine and pancuronium but not cisatracurium. Differences in ligand desolvation and ligand conformation are not likely to be important.

Roles of amino acids and subunits in determining the inhibition of nicotinic acetylcholine receptors by competitive antagonists

BACKGROUND

Binding sites for agonists and competitive antagonists (nondepolarizing neuromuscular blocking agents) are located at the alpha-delta and alpha-epsilon subunit interfaces of adult nicotinic acetylcholine receptors. Most information about the amino acids that participate in antagonist binding comes from binding studies with (+)-tubocurarine and metocurine. These bind selectively to the alpha-epsilon interface but are differentially sensitive to mutations. To test the generality of this observation, the authors measured current inhibition by five competitive antagonists on wild-type and mutant acetylcholine receptors.

METHODS

HEK293 cells were transfected with wild-type or mutant (alphaY198F, epsilonD59A, epsilonD59N, epsilonD173A, epsilonD173N, deltaD180K) mouse muscle acetylcholine receptor complementary DNA. Outside-out patches were excised and perfused with acetylcholine in the absence and presence of antagonist. Concentration-response curves were constructed to determine antagonist IC50. An antagonist-removal protocol was used to determine dissociation and association rates.

RESULTS

Effects of mutations were antagonist specific. alphaY198F decreased the IC50 of (+)-tubocurarine 10-fold, increased the IC50 of vecuronium 5-fold, and had smaller effects on other antagonists. (+)-Tubocurarine was the most sensitive antagonist to epsilonD173 mutations. epsilonD59 mutations had large effects on metocurine and cisatracurium. deltaD180K decreased inhibition by pancuronium, vecuronium, and cisatracurium. Inhibition by these antagonists was increased for receptors containing two delta subunits but no epsilon subunit. Differences in IC50 arose from differences in both dissociation and association rates.

CONCLUSION

Competitive antagonists exhibited different patterns of sensitivity to mutations. Except for pancuronium, the antagonists were sensitive to mutations at the alpha-epsilon interface. Pancuronium, vecuronium, and cisatracurium were selective for the alpha-delta interface. This suggests the possibility of synergistic inhibition by pairs of antagonists.

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.

Agonist-driven conformational changes in the inner beta-sheet of alpha7 nicotinic receptors

Cys-loop ligand-gated ion channels assemble as pentameric proteins, and each monomer contributes two structural elements: an extracellular ligand-binding domain (LBD) and a transmembrane ion channel domain. Models of receptor activation include rotational movements of subunits leading to opening of the ion channel. We tested this idea using substituted cysteine accessibility to track conformational changes in the inner beta sheet of the LBD. Using a nondesensitizing chick alpha7 background (L(247)T), we constructed 18 consecutive cysteine replacement mutants (Leu(36) to Ile(53)) and tested each for expression of acetylcholine (ACh)-evoked currents and functional sensitivity to thiol modification. We measured rates of modification in the presence and absence of ACh to identify conformational changes associated with receptor activation. Resting modification rates of eight substituted cysteines in the beta1 and beta2 strands and the sequence between them (loop 2) varied over several orders of magnitude, suggesting substantial differences in the accessibility or electrostatic environment of individual side chains. These differences were in general agreement with structural models of the LBD. Eight of 18 cysteine replacements displayed ACh-dependent changes in modification rates, indicating a change in the accessibility or electrostatic environment of the introduced cysteine during activation. We were surprised that the effects of agonist exposure were difficult to reconcile with rotational models of activation. Acetylcholine reduced the modification rate of M(40)C but increased it at N(52)C despite the close physical proximity of these residues. Our results suggest that models that depend strictly on rigid-body rotation of the LBD may provide an incomplete description of receptor activation.

Computational evidence for the ligand selectivity to the alpha4beta2 and alpha3beta4 nicotinic acetylcholine receptors

The homology models of the alpha4beta2 and alpha3beta4 nicotinic acetylcholine receptors (nAChRs) suggest that the two nAChR subtypes are different in their ligand-binding pockets due to the non-conserved residues in the beta-subunits. The docking of nicotine, epibatidine, A-84543, and the two analogs of A-84543 ligands 1 and 2 to the homology models of alpha4beta2 and alpha3beta4 is presented. It is found that the protonated amino groups of these ligands bind to the alpha-subunits, whereas the remaining parts of the ligands bind to the beta-subunits. The two non-conserved amino acids Lys77 and Phe117 in the beta2-subunit corresponding to Ile77 and Gln117 in the beta4-subunit are identified to be the key players determining the binding modes of the ligands. We demonstrate how the increase in the number of the atoms connecting the pyrrolidine and pyridine rings in A-84543, 1, and 2, and an introduction of the alkynyl substituent in the pyridine ring affect the binding and shift the selectivity of these ligands toward the beta2-containing receptors. Further improvement in affinity and selectivity in this and other series of the ligands may be achieved by designing molecules that would specifically target the non-conserved regions in nAChRs.

Computed pore potentials of the nicotinic acetylcholine receptor

Electrostatic surface potentials in the vestibule of the nicotinic acetylcholine receptor (nAChR) were computed from structural models using the University of Houston Brownian Dynamics program to determine their effect on ion conduction and ionic selectivity. To further determine whether computed potentials accurately reflect the electrostatic environment of the channel, the potentials were used to predict the rate constants for diffusion-enhanced fluorescence energy transfer; the calculated energy transfer rates are directly comparable with those determined experimentally (see companion article by Meltzer et al. in this issue). To include any effects on the local potentials by the bound acceptor fluorophore crystal violet, its binding site was first localized within the pore by fluorescence energy transfer measurements from dansyl-C6-choline bound to the agonist sites and also by simulations of binding using Autodock. To compare the computed potentials with those determined experimentally, we used the predicted energy transfer rates from Tb3+ chelates of varying charge to calculate an expected potential using the Boltzmann relationship. This expected potential (from -20 to -40 mV) overestimates the values determined experimentally (from -10 to -25 mV) by two- to fourfold at similar conditions of ionic strength. Although the results indicate a basic discrepancy between experimental and computed surface potentials, both methods demonstrate that the vestibular potential has a relatively small effect on conduction and selectivity.

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.

Recent advances in Cys-loop receptor structure and function

Throughout the nervous system, moment-to-moment communication relies on postsynaptic receptors to detect neurotransmitters and change the membrane potential. For the Cys-loop superfamily of receptors, recent structural data have catalysed a leap in our understanding of the three steps of chemical-to-electrical transduction: neurotransmitter binding, communication between the binding site and the barrier to ions, and opening and closing of the barrier. The emerging insights might be expected to explain how mutations of receptors cause neurological disease, but the opposite is generally true. Namely, analyses of disease-causing mutations have clarified receptor structure-function relationships as well as mechanisms governing the postsynaptic response.

Activation of the Torpedo nicotinic acetylcholine receptor. The contribution of residues alphaArg55 and gammaGlu93

The Torpedo nicotinic acetylcholine receptor is a heteropentamer (alpha2betagammadelta) in which structurally homologous subunits assemble to form a central ion pore. Viewed from the synaptic cleft, the likely arrangement of these subunits is alpha-gamma-alpha-delta-beta lying in an anticlockwise orientation. High affinity binding sites for agonists and competitive antagonists have been localized to the alpha-gamma and alpha-delta subunit interfaces. We investigated the involvement of amino acids lying at an adjacent interface (gamma-alpha) in receptor properties. Recombinant Torpedo receptors, expressed in Xenopus oocytes, were used to investigate the consequences of mutating alphaArg55 and gammaGlu93, residues that are conserved in most species of the peripheral nicotinic receptors. Based on homology modeling, these residues are predicted to lie in close proximity to one another and it has been suggested that they may form a salt bridge in the receptor's three-dimensional structure (Sine et al. 2002 J Biol Chem277, 29 210-29 223). Although substitution of alphaR55 by phenylalanine or tryptophan resulted in approximately a six-fold increase in the EC50 value for acetylcholine activation, the charge reversal mutation (alphaR55E) had no significant effect. In contrast, the replacement of gammaE93 by an arginine conferred an eight-fold increase in the potency for acetylcholine-induced receptor activation. In the receptor carrying the double mutations, alphaR55E-gammaE93R or alphaR55F-gammaE93R, the potency for acetylcholine activation was partially restored to that of the wild-type. The results suggest that, although individually these residues influence receptor activation, direct interactions between them are unlikely to play a major role in the stabilization of different conformational states of the receptor.

Acetylcholine nicotinic receptors: finding the putative binding site of allosteric modulators using the "blind docking" approach

Allosteric potentiation of acetylcholine nicotinic receptors is considered to be one of the most promising approaches for the treatment of Alzheimer's disease. However, the exact localization of the allosteric binding site and the potentiation mechanism at the molecular level are presently unknown. We have performed the "blind docking" of three known allosteric modulators (galanthamine, codeine and eserine) with the Acetylcholine Binding Protein and models of human alpha7, alpha3beta4 and alpha4beta2 nicotinic receptors, created by homology modeling. Three putative binding sites were identified in the channel pore, each one showing different affinities for the ligands. One of these sites is localized opposite to the agonist binding site and is probably implicated in the potentiation process. On the basis of these results, a possible mechanism for nicotinic acetylcholine receptor (nAChRs) activation is proposed. The present findings may represent an important advance for understanding the allosteric modulation mechanism of nAChRs. [Figure: see text].

Modeling multiple species of nicotine and deschloroepibatidine interacting with alpha4beta2 nicotinic acetylcholine receptor: from microscopic binding to phenomenological binding affinity

A variety of molecular modeling, molecular docking, and first-principles electronic structure calculations were performed to study how the alpha4beta2 nicotinic acetylcholine receptor (nAChR) binds with different species of two typical agonists, (S)-(-)-nicotine and (R)-(-)-deschloroepibatidine, each of which is distinguished by different free bases and protonation states. On the basis of these results, predictions were made regarding the corresponding microscopic binding free energies. Hydrogen-bonding and cation-pi interactions between the receptor and the respective ligands were found to be the dominant factors differentiating the binding strengths of different microscopic binding species. The calculated results and analyses demonstrate that, for each agonist, all the species are interchangeable and can quickly achieve a thermodynamic equilibrium in solution and at the nAChR binding site. This allows quantitation of the equilibrium concentration distributions of the free ligand species and the corresponding microscopic ligand-receptor binding species, their pH dependence, and their contributions to the phenomenological binding affinity. The predicted equilibrium concentration distributions, pK(a) values, absolute phenomenological binding affinities, and their pH dependence are all in good agreement with available experimental data, suggesting that the computational strategy from the microscopic binding species and affinities to the phenomenological binding affinity is reliable for studying alpha4beta2 nAChR-ligand binding. This should provide valuable information for future rational design of drugs targeting nAChRs. The general strategy of the "from-microscopic-to-phenomenological" approach for studying interactions of alpha4beta2 nAChRs with (S)-(-)-nicotine and (R)-(-)-deschloroepibatidine may also be useful in studying other types of ligand-protein interactions involving multiple molecular species of a ligand and in associated rational drug design.

Docking of 6-chloropyridazin-3-yl derivatives active on nicotinic acetylcholine receptors into molluscan acetylcholine binding protein (AChBP)

The crystal structure of Acetylcholine Binding Protein (AChBP), homolog of the ligand binding domain of nAChR, has been used as model for computational investigations on the ligand-receptor interactions of derivatives of 6-chloropyridazine substituted at C3 with 3,8-diazabicyclo[3.2.1]octane, 2,5-diazabicyclo[2.2.1]heptane and with piperazine and homopiperazine, substituted or not at N4. The ligand-receptor complexes have been analyzed by docking techniques using the binding site of HEPES complexed with AChBP as template. The good relationship between the observed binding affinity and the calculated docking energy confirms that this model provides a good starting point for understanding the binding domain of neuronal nicotinic receptors. An analysis of the possible factors significant for the ligand recognition has evidenced, besides the cation-pi interaction, the distance between the chlorine atom of the pyridazinyl group and the carbonylic oxygen of Leu B112 as an important parameter in the modulation of the binding energy.

Initial coupling of binding to gating mediated by conserved residues in the muscle nicotinic receptor

We examined functional consequences of intrasubunit contacts in the nicotinic receptor alpha subunit using single channel kinetic analysis, site-directed mutagenesis, and structural modeling. At the periphery of the ACh binding site, our structural model shows that side chains of the conserved residues alphaK145, alphaD200, and alphaY190 converge to form putative electrostatic interactions. Structurally conservative mutations of each residue profoundly impair gating of the receptor channel, primarily by slowing the rate of channel opening. The combined mutations alphaD200N and alphaK145Q impair channel gating to the same extent as either single mutation, while alphaK145E counteracts the impaired gating due to alphaD200K, further suggesting electrostatic interaction between these residues. Interpreted in light of the crystal structure of acetylcholine binding protein (AChBP) with bound carbamylcholine (CCh), the results suggest in the absence of ACh, alphaK145 and alphaD200 form a salt bridge associated with the closed state of the channel. When ACh binds, alphaY190 moves toward the center of the binding cleft to stabilize the agonist, and its aromatic hydroxyl group approaches alphaK145, which in turn loosens its contact with alphaD200. The positional changes of alphaK145 and alphaD200 are proposed to initiate the cascade of perturbations that opens the receptor channel: the first perturbation is of beta-strand 7, which harbors alphaK145 and is part of the signature Cys-loop, and the second is of beta-strand 10, which harbors alphaD200 and connects to the M1 domain. Thus, interplay between these three conserved residues relays the initial conformational change from the ACh binding site toward the ion channel.

Crystal structure of nicotinic acetylcholine receptor homolog AChBP in complex with an α-conotoxin PnIA variant

Conotoxins (Ctx) form a large family of peptide toxins from cone snail venoms that act on a broad spectrum of ion channels and receptors. The subgroup alpha-Ctx specifically and selectively binds to subtypes of nicotinic acetylcholine receptors (nAChRs), which are targets for treatment of several neurological disorders. Here we present the structure at a resolution of 2.4 A of alpha-Ctx PnIA (A10L D14K), a potent blocker of the alpha(7)-nAChR, bound with high affinity to acetylcholine binding protein (AChBP), the prototype for the ligand-binding domains of the nAChR superfamily. Alpha-Ctx is buried deep within the ligand-binding site and interacts with residues on both faces of adjacent subunits. The toxin itself does not change conformation, but displaces the C loop of AChBP and induces a rigid-body subunit movement. Knowledge of these contacts could facilitate the rational design of drug leads using the Ctx framework and may lead to compounds with increased receptor subtype selectivity.

Theoretical studies of the M2 transmembrane segment of the glycine receptor: models of the open pore structure and current-voltage characteristics

The pentameric glycine receptor (GlyR), a member of the nicotinicoid superfamily of ligand-gated ion channels, is an inhibitory Cl(-) channel that is gated by glycine. Using recently published NMR data of the second transmembrane segment (M2) of the human alpha1 GlyR, structural models of pentameric assemblies embedded in a lipid bilayer were constructed using a combination of experimentally determined constraints coupled with all-atom energy minimization. Based on this structure of the pentameric M2 "pore", Brownian dynamics simulations of ion permeation through this putative conducting open state of the channel were carried out. Simulated I-V curves were in good agreement with published experimental current-voltage curves and the anion/cation permeability ratio, suggesting that our open-state model may be representative of the conducting channel of the full-length receptor. These studies also predicted regions of chloride occupancy and suggested residues critical to anion permeation. Calculations of the conductance of the cation-selective mutant A251E channel are also consistent with experimental data. In addition, both rotation and untilting of the pore helices of our model were found to be broadly consistent with closing of the channel, albeit at distinct regions that may reflect alternate gates of the receptor.

Mutation of single murine acetylcholine receptor subunits reveals differential contribution of P121 to acetylcholine binding and channel opening

The nicotinic acetylcholine receptor (AChR) is a heteropentameric, ligand-gated ion channel at the neuromuscular junction, where it is responsible for signal transduction between the motorneuron and the muscle. Point mutations in the subunits of the receptor change the channel's electrophysiological properties and underlie inherited forms of muscle weakness, the congenital myasthenic syndromes. One point mutation (P121L) has been identified in the epsilon-subunit of patients suffering from the fast-channel congenital myasthenic syndrome, which is evoked by reduced AChR openings. We introduced the P121L mutation into all murine AChR subunits and performed electrophysiological studies in Xenopus laevis oocytes. The P121L mutation in the epsilon-subunit of the adult mouse AChR affected ligand binding and channel gating in a manner similar to that described for human AChR. At equivalent positions in the alpha- and beta-subunits, the mutation caused only minor electrophysiological changes. Mutation of the delta-subunit had similar, but less pronounced functional consequences compared to epsilonP121L, reflecting the asymmetry of the acetylcholine binding sites and the dominant effect of the alpha-epsilon site on channel opening.

Invariant aspartic Acid in muscle nicotinic receptor contributes selectively to the kinetics of agonist binding

We examined functional contributions of interdomain contacts within the nicotinic receptor ligand binding site using single channel kinetic analyses, site-directed mutagenesis, and a homology model of the major extracellular region. At the principal face of the binding site, the invariant αD89 forms a highly conserved interdomain contact near αT148, αW149, and αT150. Patch-clamp recordings show that the mutation αD89N markedly slows acetylcholine (ACh) binding to receptors in the resting closed state, but does not affect rates of channel opening and closing. Neither αT148L, αT150A, nor mutations at both positions substantially affects the kinetics of receptor activation, showing that hydroxyl side chains at these positions are not hydrogen bond donors for the strong acceptor αD89. However substituting a negative charge at αT148, but not at αT150, counteracts the effect of αD89N, demonstrating that a negative charge in the region of interdomain contact confers rapid association of ACh. Interpreted within the structural framework of ACh binding protein and a homology model of the receptor ligand binding site, these results implicate main chain amide groups in the domain harboring αW149 as principal hydrogen bond donors for αD89. The specific effect of αD89N on ACh association suggests that interdomain hydrogen bonding positions αW149 for optimal interaction with ACh.

Agonist-mediated Conformational Changes in Acetylcholine-binding Protein Revealed by Simulation and Intrinsic Tryptophan Fluorescence

We delineated acetylcholine (ACh)-dependent conformational changes in a prototype of the nicotinic receptor ligand binding domain by molecular dynamics simulation and changes in intrinsic tryptophan (Trp) fluorescence. Prolonged molecular dynamics simulation of ACh-binding protein showed that binding of ACh establishes close register of Trps from adjacent subunits, Trp(143) and Trp(53), and draws the peripheral C-loop inward to occlude the entrance to the binding cavity. Close register of Trp(143) and Trp(53) was demonstrated by ACh-mediated quenching of intrinsic Trp fluorescence, elimination of quenching by mutation of one or both Trps to Phe, and decreased lifetime of Trp fluorescence by bound ACh. Occlusion of the binding cavity by the C-loop was demonstrated by restricted access of an extrinsic quencher of binding site Trp fluorescence by ACh. The collective findings showed that ACh initially establishes close register of conserved Trps from adjacent subunits and then draws the C-loop inward to occlude the entrance to the binding cavity.

Refined structure of the nicotinic acetylcholine receptor at 4A resolution

We present a refined model of the membrane-associated Torpedo acetylcholine (ACh) receptor at 4A resolution. An improved experimental density map was obtained from 342 electron images of helical tubes, and the refined structure was derived to an R-factor of 36.7% (R(free) 37.9%) by standard crystallographic methods, after placing the densities corresponding to a single molecule into an artificial unit cell. The agreement between experimental and calculated phases along the helical layer-lines was used to monitor progress in the refinement and to give an independent measure of the accuracy. The atomic model allowed a detailed description of the whole receptor in the closed-channel form, including the ligand-binding and intracellular domains, which have not previously been interpreted at a chemical level. We confirm that the two ligand-binding alpha subunits have a different extended conformation from the three other subunits in the closed channel, and identify several interactions on both pairs of subunit interfaces, and within the alpha subunits, which may be responsible for their "distorted" structures. The ACh-coordinating amino acid side-chains of the alpha subunits are far apart in the closed channel, indicating that a localised rearrangement, involving closure of loops B and C around the bound ACh molecule, occurs upon activation. A comparison of the structure of the alpha subunit with that of AChBP having ligand present, suggests how the localised rearrangement overcomes the distortions and initiates the rotational movements associated with opening of the channel. Both vestibules of the channel are strongly electronegative, providing a cation-stabilising environment at either entrance of the membrane pore. Access to the pore on the intracellular side is further influenced by narrow lateral windows, which would be expected to screen out electrostatically ions of the wrong charge and size.

Ligand-induced conformational change in the alpha7 nicotinic receptor ligand binding domain

Molecular dynamics simulations of a homology model of the ligand binding domain of the alpha7 nicotinic receptor are conducted with a range of bound ligands to induce different conformational states. Four simulations of 15 ns each are run with no ligand, antagonist d-tubocurarine (dTC), agonist acetylcholine (ACh), and agonist ACh with potentiator Ca(2+), to give insight into the conformations of the active and inactive states of the receptor and suggest the mechanism for conformational change. The main structural factor distinguishing the active and inactive states is that a more open, symmetric arrangement of the five subunits arises for the two agonist simulations, whereas a more closed and asymmetric arrangement results for the apo and dTC cases. Most of the difference arises in the lower portion of the ligand binding domain near its connection to the adjacent transmembrane domain. The transfer of the more open state to the transmembrane domain could then promote ion flow through the channel. Variation in how subunits pack together with no ligand bound appears to give rise to asymmetry in the apo case. The presence of dTC expands the receptor but induces rotations in alternate directions in adjacent subunits that lead to an asymmetric arrangement as in the apo case. Ca(2+) appears to promote a slightly greater expansion in the subunits than ACh alone by stabilizing the C-loop and ACh positions. Although the simulations are unlikely to be long enough to view the full conformational changes between open and closed states, a collection of different motions at a range of length scales are observed that are likely to participate in the conformational change.

Structural basis for epibatidine selectivity at desensitized nicotinic receptors

The agonist binding sites of the fetal muscle nicotinic acetylcholine receptor are formed at the interfaces of alpha-subunits and neighboring gamma- and delta-subunits. When the receptor is in the nonconducting desensitized state, the alpha-gamma site binds the agonist epibatidine 200-fold more tightly than does the alpha-delta site. To determine the structural basis for this selectivity, we constructed gamma/delta-subunit chimeras, coexpressed them with complementary wild-type subunits in HEK 293 cells, and determined epibatidine affinity of the resulting complexes. The results reveal three determinants of epibatidine selectivity: gamma104-117/delta106-delta119, gamma164-171/delta166-177, and gammaPro190/deltaAla196. Point mutations reveal that three sequence differences within the gamma104-117/delta106-delta119 region are determinants of epibatidine selectivity: gammaLys104/deltaTyr106, gammaSer111/deltaTyr113, and gammaTyr117/deltaTyr119. In the delta-subunit, simultaneous mutation of these residues to their gamma equivalent produces high affinity, gamma-like epibatidine binding. However, converting gamma to delta affinity requires replacement of the gamma104-117 segment with delta sequence, suggesting interplay of residues in this region. The structural basis for epibatidine selectivity is explained by computational docking of epibatidine to a homology model of the alpha-gamma binding site.

Coupling of agonist binding to channel gating in an ACh-binding protein linked to an ion channel

Neurotransmitter receptors from the Cys-loop superfamily couple the binding of agonist to the opening of an intrinsic ion pore in the final step in rapid synaptic transmission. Although atomic resolution structural data have recently emerged for individual binding and pore domains, how they are linked into a functional unit remains unknown. Here we identify structural requirements for functionally coupling the two domains by combining acetylcholine (ACh)-binding protein, whose structure was determined at atomic resolution, with the pore domain from the serotonin type-3A (5-HT3A) receptor. Only when amino-acid sequences of three loops in ACh-binding protein are changed to their 5-HT3A counterparts does ACh bind with low affinity characteristic of activatable receptors, and trigger opening of the ion pore. Thus functional coupling requires structural compatibility at the interface of the binding and pore domains. Structural modelling reveals a network of interacting loops between binding and pore domains that mediates this allosteric coupling process.

Computational approaches to understand alpha-conotoxin interactions at neuronal nicotinic receptors

Recent and increasing use of computational tools in the field of nicotinic receptors has led to the publication of several models of ligand-receptor interactions. These models are all based on the crystal structure at 2.7 A resolution of a protein related to the extracellular N-terminus of nicotinic acetylcholine receptors (nAChRs), the acetylcholine binding protein. In the absence of any X-ray or NMR information on nAChRs, this new structure has provided a reliable alternative to study the nAChR structure. We are now able to build homology models of the binding domain of any nAChR subtype and fit in different ligands using docking programs. This strategy has already been performed successfully for the docking of several nAChR agonists and antagonists. This minireview focuses on the interaction of alpha-conotoxins with neuronal nicotinic receptors in light of our new understanding of the receptor structure. Computational tools are expected to reveal the molecular recognition mechanisms that govern the interaction between alpha-conotoxins and neuronal nAChRs at the molecular level. An accurate determination of their binding modes on the neuronal nAChR may allow the rational design of alpha-conotoxin-based ligands with novel nAChR selectivity.

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

A homology model of the ligand binding domain of the alpha7 nicotinic receptor is constructed based on the acetylcholine-binding protein crystal structure. This structure is refined in a 10 ns molecular dynamics simulation. The modeled structure proves fairly resilient, with no significant changes at the secondary or tertiary structural levels. The hypothesis that the acetylcholine-binding protein template is in the activated or desensitized state, and the absence of a bound agonist in the simulation suggests that the structure may also be relaxing from this state to the activatable state. Candidate motions that take place involve not only the side chains of residues lining the binding sites, but also the subunit positions that determine the overall shape of the receptor. In particular, two nonadjacent subunits move outward, whereas their partners counterclockwise to them move inward, leading to a marginally wider interface between themselves and an overall asymmetric structure. This in turn affects the binding sites, producing two that are more open and characterized by distinct side-chain conformations of W54 and L118, although motions of the side chains of all residues in every binding site still contribute to a reduction in binding site size, especially the outward motion of W148, which hinders acetylcholine binding. The Cys loop at the membrane interface also displays some flexibility. Although the short simulation timescale is unlikely to sample adequately all the conformational states, the pattern of observed motions suggests how ligand binding may correlate with larger-scale subunit motions that would connect with the transmembrane region that controls the passage of ions. Furthermore, the shape of the asymmetry with binding sites of differing affinity for acetylcholine, characteristic of other nicotinic receptors, may be a natural property of the relaxed, activatable state of alpha7.

Base docking model of the homomeric α7 nicotinic receptor–β-amyloid1–42 complex

A homology model of the human alpha7 nicotinic receptor was constructed based on the acetylcholine-binding protein crystal structure. Subsequently, the three-dimensional structure of the complex between the alpha7 nicotinic receptor and the 42-amino acid beta-amyloid peptide was obtained for the first time with the aid of the ESCHER program, a well-known method for protein-protein docking. The final complex showed that the most important interactions occur between the residues V12-K28 from the peptide and the loop C of the receptor. The model agrees with many experimental data, and may be used as a base model for further detailed studies in order to gain insight into the binding and dynamics of the complex at molecular level and their correlation with the memory impairments characteristic of the Alzheimer's disease.

Acrylodan-conjugated Cysteine Side Chains Reveal Conformational State and Ligand Site Locations of the Acetylcholine-binding Protein

We undertook cysteine substitution mutagenesis and fluorophore conjugation at selected residue positions to map sites of ligand binding and changes in solvent exposure of the acetylcholine-binding protein from Lymnaea stagnalis, a nicotinic receptor surrogate. Acrylodan fluorescence emission is highly sensitive to its local environment, and when bound to protein, exhibits changes in both intensity and emission wavelength that are reflected in the degree of solvent exclusion and the effective dielectric constant of the environment of the fluorophore. Hence, cysteine mutants were generated based on the acetylcholine-binding protein crystal structure and predicted ligand binding sites, and fluorescence parameters were assayed on the acrylodan-conjugated proteins. This approach allows one to analyze the environment around the conjugated fluorophore side chain and the changes induced by bound ligand. Introduction of an acrylodan-cysteine conjugate at position 178 yields a large blue shift with alpha-bungarotoxin association, whereas the agonists and alkaloid antagonists induce red shifts reflecting solvent exposure at this position. Such residue-selective changes in fluorescence parameters suggest that certain ligands can induce distinct conformational states of the binding protein, and that mutually exclusive binding results from disparate portals of entry to and orientations of the bound alpha-toxin and smaller acetylcholine congeners at the binding pocket. Labeling at other residue positions around the predicted binding pocket also reveals distinctive spectral changes for alpha-bungarotoxin, agonists, and alkaloid antagonists.