Dynamics of the acetylcholine receptor pore at the gating transition state

Mechanism of hydrophobic gating in the acetylcholine receptor channel pore

Neuromuscular acetylcholine receptors (AChRs) are hetero-pentameric, ligand-gated ion channels. The binding of the neurotransmitter acetylcholine (ACh) to two target sites promotes a global conformational change of the receptor that opens the channel and allows ion conduction through the channel pore. Here, by measuring free-energy changes from single-channel current recordings and using molecular dynamics simulations, we elucidate how a constricted hydrophobic region acts as a "gate" to regulate the channel opening in the pore of AChRs. Mutations of gate residues, including those implicated in congenital myasthenia syndrome, lower the permeation barrier of the channel substantially and increase the unliganded gating equilibrium constant (constitutive channel openings). Correlations between hydrophobicity and the observed free-energy changes, supported by calculations of water densities in the wild-type versus mutant channel pores, provide evidence for hydrophobic wetting-dewetting transition at the gate. The analysis of a coupled interaction network provides insight into the molecular mechanism of closed- versus open-state conformational changes at the gate. Studies of the transition state by "phi"(φ)-value analysis indicate that agonist binding serves to stabilize both the transition and the open state. Intersubunit interaction energy measurements and molecular dynamics simulations suggest that channel opening involves tilting of the pore-lining M2 helices, asymmetric outward rotation of amino acid side chains, and wetting transition of the gate region that lowers the barrier to ion permeation and stabilizes the channel open conformation. Our work provides new insight into the hydrophobic gate opening and shows why the gate mutations result in constitutive AChR channel activity.

Structural mechanism of muscle nicotinic receptor desensitization and block by curare

Binding of the neurotransmitter acetylcholine to its receptors on muscle fibers depolarizes the membrane and thereby triggers muscle contraction. We sought to understand at the level of three-dimensional structure how agonists and antagonists alter nicotinic acetylcholine receptor conformation. We used the muscle-type receptor from the Torpedo ray to first define the structure of the receptor in a resting, activatable state. We then determined the receptor structure bound to the agonist carbachol, which stabilizes an asymmetric, closed channel desensitized state. We find conformational changes in a peripheral membrane helix are tied to recovery from desensitization. To probe mechanisms of antagonism, we obtained receptor structures with the active component of curare, a poison arrow toxin and precursor to modern muscle relaxants. d-Tubocurarine stabilizes the receptor in a desensitized-like state in the presence and absence of agonist. These findings define the transitions between resting and desensitized states and reveal divergent means by which antagonists block channel activity of the muscle-type nicotinic receptor.

The desensitization pathway of GABAA receptors, one subunit at a time

GABAA receptors mediate most inhibitory synaptic transmission in the brain of vertebrates. Following GABA binding and fast activation, these receptors undergo a slower desensitization, the conformational pathway of which remains largely elusive. To explore the mechanism of desensitization, we used concatemeric α1β2γ2 GABAA receptors to selectively introduce gain-of-desensitization mutations one subunit at a time. A library of twenty-six mutant combinations was generated and their bi-exponential macroscopic desensitization rates measured. Introducing mutations at the different subunits shows a strongly asymmetric pattern with a key contribution of the γ2 subunit, and combining mutations results in marked synergistic effects indicating a non-concerted mechanism. Kinetic modelling indeed suggests a pathway where subunits move independently, the desensitization of two subunits being required to occlude the pore. Our work thus hints towards a very diverse and labile conformational landscape during desensitization, with potential implications in physiology and pharmacology.

Electric discharge of electrocytes: Modelling, analysis and simulation

In this paper, we investigate the electric discharge of electrocytes by extending our previous work on the generation of electric potential. We first give a complete formulation of a single cell unit consisting of an electrocyte and a resistor, based on a Poisson-Nernst-Planck (PNP) system with various membrane currents as interfacial conditions for the electrocyte and a Maxwell's model for the resistor. Our previous work can be treated as a special case with an infinite resistor (or open circuit). Using asymptotic analysis, we simplify our PNP system and reduce it to an ordinary differential equation (ODE) based model. Unlike the case of an infinite resistor, our numerical simulations of the new model reveal several distinct features. A finite current is generated, which leads to non-constant electric potentials in the bulk of intracellular and extracellular regions. Furthermore, the current induces an additional action potential (AP) at the non-innervated membrane, contrary to the case of an open circuit where an AP is generated only at the innervated membrane. The voltage drop inside the electrocyte is caused by an internal resistance due to mobile ions. We show that our single cell model can be used as the basis for a system with stacked electrocytes and the total current during the discharge of an electric eel can be estimated by using our model.

Electric potential generation of electrocytes: Modelling, analysis, and computation

In this paper, we developed a one-dimensional model for electric potential generation of electrocytes in electric eels. The model is based on the Poisson-Nernst-Planck system for ion transport coupled with membrane fluxes including the Hodgkin-Huxley type. Using asymptotic analysis, we derived a simplified zero-dimensional model, which we denote as the membrane model in this paper, as a leading order approximation. Our analysis provides justification for the assumption in membrane models that electric potential is constant in the intracellular space. This is essential to explain the superposition of two membrane potentials that leads to a significant transcellular potential. Numerical simulations are also carried out to support our analytical findings.

Timing of CFTR pore opening and structure of its transition state

In CFTR, the chloride ion channel mutated in cystic fibrosis (CF) patients, pore opening is coupled to ATP-binding-induced dimerization of two cytosolic nucleotide binding domains (NBDs) and closure to dimer disruption following ATP hydrolysis. CFTR opening rate, unusually slow because of its high-energy transition state, is further slowed by CF mutation ΔF508. Here, we exploit equilibrium gating of hydrolysis-deficient CFTR mutant D1370N and apply rate-equilibrium free-energy relationship analysis to estimate relative timing of opening movements in distinct protein regions. We find clear directionality of motion along the longitudinal protein axis and identify an opening transition-state structure with the NBD dimer formed but the pore still closed. Thus, strain at the NBD/pore-domain interface, the ΔF508 mutation locus, underlies the energetic barrier for opening. Our findings suggest a therapeutic opportunity to stabilize this transition-state structure pharmacologically in ΔF508-CFTR to correct its opening defect, an essential step toward restoring CFTR function.

Energetic contributions to channel gating of residues in the muscle nicotinic receptor β1 subunit

In the pentameric ligand-gated ion channel family, transmitter binds in the extracellular domain and conformational changes result in channel opening in the transmembrane domain. In the muscle nicotinic receptor and other heteromeric members of the family one subunit does not contribute to the canonical agonist binding site for transmitter. A fundamental question is whether conformational changes occur in this subunit. We used records of single channel activity and rate-equilibrium free energy relationships to examine the β1 (non-ACh-binding) subunit of the muscle nicotinic receptor. Mutations to residues in the extracellular domain have minimal effects on the gating equilibrium constant. Positions in the channel lining (M2 transmembrane) domain contribute strongly and relatively late during gating. Positions thought to be important in other subunits in coupling the transmitter-binding to the channel domains have minimal effects on gating. We conclude that the conformational changes involved in channel gating propagate from the binding-site to the channel in the ACh-binding subunits and subsequently spread to the non-binding subunit.

Nicotinic acetylcholine receptor and the structural basis of neuromuscular transmission: insights from Torpedo postsynaptic membranes

The nicotinic acetylcholine (ACh) receptor, at the neuromuscular junction, is a neurotransmitter-gated ion channel that has been fine-tuned through evolution to transduce a chemical signal into an electrical signal with maximum efficiency and speed. It is composed from three similar and two identical polypeptide chains, arranged in a ring around a narrow membrane pore. Central to the design of this assembly is a hydrophobic gate in the pore, more than 50 Å away from sites in the extracellular domain where ACh binds. Although the molecular properties of the receptor have been explored intensively over the last few decades, only recently have structures emerged revealing its complex architecture and illuminating how ACh entering the binding sites opens the distant gate. Postsynaptic membranes isolated from the (muscle-derived) electric organ of the Torpedo ray have underpinned most of the structural studies: the membranes form tubular vesicles having receptors arranged on a regular surface lattice, which can be imaged directly in frozen physiological solutions. Advances in electron crystallographic techniques have also been important, enabling analysis of the closed- and open-channel forms of the receptor in unreacted tubes or tubes reacted briefly with ACh. The structural differences between these two forms show that all five subunits participate in a concerted conformational change communicating the effect of ACh binding to the gate, but that three of them (α_γ, β and δ) play a dominant role. Flexing of oppositely facing pore-lining α-helices is the principal motion determining the closed/open state of the gate. These results together with the findings of biochemical, biophysical and other structural studies allow an integrated description of the receptor and of its mode of action at the synapse.

Energy for wild-type acetylcholine receptor channel gating from different choline derivatives

Agonists, including the neurotransmitter acetylcholine (ACh), bind at two sites in the neuromuscular ACh receptor channel (AChR) to promote a reversible, global change in protein conformation that regulates the flow of ions across the muscle cell membrane. In the synaptic cleft, ACh is hydrolyzed to acetate and choline. Replacement of the transmitter's ester acetyl group with a hydroxyl (ACh→choline) results in a + 1.8 kcal/mol reduction in the energy for gating generated by each agonist molecule from a low- to high-affinity change of the transmitter binding site (ΔG(B)). To understand the distinct actions of structurally related agonist molecules, we measured ΔG(B) for 10 related choline derivatives. Replacing the hydroxyl group of choline with different substituents, such as hydrogen, chloride, methyl, or amine, increased the energy for gating (i.e., it made ΔG(B) more negative relative to choline). Extending the ethyl hydroxide tail of choline to propyl and butyl hydroxide also increased this energy. Our findings reveal the amount of energy that is available for the AChR conformational change provided by different, structurally related agonists. We speculate that a hydrogen bond between the choline hydroxyl and the backbone carbonyl of αW149 positions this agonist's quaternary ammonium group so as to reduce the cation-π interaction between this moiety and the aromatic groups at the binding site.

Function of interfacial prolines at the transmitter-binding sites of the neuromuscular acetylcholine receptor

The neuromuscular acetylcholine (ACh) receptor has two conserved prolines in loop D of the complementary subunit at each of its two transmitter-binding sites (α-ε and α-δ). We used single-channel electrophysiology to estimate the energy changes caused by mutations of these prolines with regard to unliganded gating (ΔG0) and the affinity change for ACh that increases the open channel probability (ΔGB). The effects of mutations of ProD2 (εPro-121/δPro-123) were greater than those of its neighbor (εPro-120/δPro-122) and were greater at α-ε versus α-δ. The main consequence of the congenital myasthenic syndrome mutation εProD2-L was to impair the establishment of a high affinity for ACh and thus make ΔGB less favorable. At both binding sites, most ProD2 mutations decreased constitutive activity (increased ΔG0). LRYHQG and RL substitutions reduced substantially the net binding energy (made ΔGB(ACh) less favorable) by ≥2 kcal/mol at α-ε and α-δ, respectively. Mutant cycle analyses were used to estimate energy coupling between the two ProD2 residues and between each ProD2 and glycine residues (αGly-147 and αGly-153) on the primary (α subunit) side of each binding pocket. The distant binding site prolines interact weakly. ProD2 interacts strongly with αGly-147 but only at α-ε and only when ACh is present. The results suggest that in the low to-high affinity change there is a concerted inter-subunit strain in the backbones at εProD2 and αGly-147. It is possible to engineer receptors having a single functional binding site by using a α-ε or α-δ ProD2-R knock-out mutation. In adult-type ACh receptors, the energy from the affinity change for ACh is approximately the same at the two binding sites (approximately -5 kcal/mol).

Loop C and the mechanism of acetylcholine receptor-channel gating

Agonist molecules at the two neuromuscular acetylcholine (ACh) receptor (AChR) transmitter-binding sites increase the probability of channel opening. In one hypothesis for AChR activation (“priming”), the capping of loop C at each binding site transfers energy independently to the distant gate over a discrete structural pathway. We used single-channel analyses to examine the experimental support for this proposal with regard to brief unliganded openings, the effects of loop-C modifications, the effects of mutations to residues either on or off the putative pathway, and state models for describing currents at low [ACh]. The results show that (a) diliganded and brief unliganded openings are generated by the same essential, global transition; (b) the radical manipulation of loop C does not prevent channel opening but impairs agonist binding; (c) both on- and off-pathway mutations alter gating by changing the relative stability of the open-channel conformation by local interactions rather than by perturbing a specific site–gate communication link; and (d) it is possible to estimate directly the rate constants for agonist dissociation from and association to both the low and high affinity forms of the AChR-binding site by using a cyclic kinetic model. We conclude that the mechanism of energy transfer between the binding sites and the gate remains an open question.

The energy and work of a ligand-gated ion channel

Ligand-gated ion channels are allosteric membrane proteins that isomerize between C(losed) and O(pen) conformations. A difference in affinity for ligands in the two states influences the C↔O "gating" equilibrium constant. The energies associated with adult-type mouse neuromuscular nicotinic acetylcholine receptor (AChR) channel gating have been measured by using single-channel electrophysiology. Without ligands, the free energy, enthalpy and entropy of gating are ΔG0=+8.4, ΔH0=+10.9 and TΔS0=+2.5kcal/mol (-100mV, 23°C). Many mutations throughout the protein change ΔG0, including natural ones that cause disease. Agonists and most mutations change approximately independently the ground-state energy difference; thus, it is possible to forecast and engineer AChR responses simply by combining perturbations. The free energy of the low↔high affinity change for the neurotransmitter at each of two functionally equivalent binding sites is ΔGB(ACh)=-5.1kcal/mol. ΔGB(ACh) is set mainly by interactions of ACh with just three binding site aromatic groups. For a series of structurally related agonists, there is a correlation between the energies of low- and high-affinity binding, which implies that gating commences with the formation of the low-affinity complex. Brief, intermediate states in binding and gating have been detected. Several proposals for the nature of the gating transition-state energy landscape and the isomerization mechanism are discussed.

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.

The energetic consequences of loop 9 gating motions in acetylcholine receptor-channels

Acetylcholine receptor-channels (AChRs) mediate fast synaptic transmission between nerve and muscle. In order to better-understand the mechanism by which this protein assembles and isomerizes between closed- and open-channel conformations we measured changes in the diliganded gating equilibrium constant (E(2)) consequent to mutations of residues at the C-terminus of loop 9 (L9) in the α and ε subunits of mouse neuromuscular AChRs. These amino acids are close to two interesting interfaces, between the extracellular and transmembrane domain within a subunit (E–T interface) and between primary and complementary subunits (P–C interface). Most α subunit mutations modestly decreased E(2) (mainly by slowing the channel-opening rate constant) and sometimes produced AChRs that had heterogeneous gating kinetic properties. Mutations in the ε subunit had a larger effect and could either increase or decrease E(2), but did not induce kinetic heterogeneity. There are broad-but-weak energetic interactions between αL9 residues and others at the αE–T interface, as well as between the εL9 residue and others at the P–C interface (in particular, the M2–M3 linker). These interactions serve, in part, to maintain the structural integrity of the AChR assembly at the E–T interface. Overall, the energy changes of L9 residues are significant but smaller than in other regions of the protein.

Novel mouse model reveals distinct activity-dependent and -independent contributions to synapse development

The balanced action of both pre- and postsynaptic organizers regulates the formation of neuromuscular junctions (NMJ). The precise mechanisms that control the regional specialization of acetylcholine receptor (AChR) aggregation, guide ingrowing axons and contribute to correct synaptic patterning are unknown. Synaptic activity is of central importance and to understand synaptogenesis, it is necessary to distinguish between activity-dependent and activity-independent processes. By engineering a mutated fetal AChR subunit, we used homologous recombination to develop a mouse line that expresses AChR with massively reduced open probability during embryonic development. Through histological and immunochemical methods as well as electrophysiological techniques, we observed that endplate anatomy and distribution are severely aberrant and innervation patterns are completely disrupted. Nonetheless, in the absence of activity AChRs form postsynaptic specializations attracting motor axons and permitting generation of multiple nerve/muscle contacts on individual fibers. This process is not restricted to a specialized central zone of the diaphragm and proceeds throughout embryonic development. Phenotypes can be attributed to separate activity-dependent and -independent pathways. The correct patterning of synaptic connections, prevention of multiple contacts and control of nerve growth require AChR-mediated activity. In contrast, myotube survival and acetylcholine-mediated dispersal of AChRs are maintained even in the absence of AChR-mediated activity. Because mouse models in which acetylcholine is entirely absent do not display similar effects, we conclude that acetylcholine binding to the AChR initiates activity-dependent and activity-independent pathways whereby the AChR modulates formation of the NMJ.

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.

Glycine hinges with opposing actions at the acetylcholine receptor-channel transmitter binding site

The extent to which agonists activate synaptic receptor-channels depends on both the intrinsic tendency of the unliganded receptor to open and the amount of agonist binding energy realized in the channel-opening process. We examined mutations of the nicotinic acetylcholine receptor transmitter binding site (α subunit loop B) with regard to both of these parameters. αGly147 is an "activation" hinge where backbone flexibility maintains high values for intrinsic gating, the affinity of the resting conformation for agonists and net ligand binding energy. αGly153 is a "deactivation" hinge that maintains low values for these parameters. αTrp149 (between these two glycines) serves mainly to provide ligand binding energy for gating. We propose that a concerted motion of the two glycine hinges (plus other structural elements at the binding site) positions αTrp149 so that it provides physiologically optimal binding and gating function at the nerve-muscle synapse.

Linking the acetylcholine receptor-channel agonist-binding sites with the gate

The gating isomerization of neuromuscular acetylcholine receptors links the rearrangements of atoms at two transmitter-binding sites with those at a distant gate region in the pore. To explore the mechanism of this reversible process, we estimated the gating rate and equilibrium constants for receptors with point mutations of alpha-subunit residues located between the binding sites and the membrane domain (N95, A96, Y127, and I49). The maximum energy change caused by a side-chain substitution at alphaA96 was huge (approximately 8.6 kcal/mol, the largest value measured so far for any alpha-subunit amino acid). A Phi-value analysis suggests that alphaA96 experiences its change in energy (structure) approximately synchronously with residues alphaY127 and alphaI49, but after the agonist molecule and other residues in loop A. Double mutant-cycle experiments show that the energy changes at alphaA96 are strongly coupled with those of alphaY127 and alphaI49. We identify a column of mutation-sensitive residues in the alpha-subunit that may be a pathway for energy transfer through the extracellular domain in the gating isomerization.

Subunit symmetry at the extracellular domain-transmembrane domain interface in acetylcholine receptor channel gating

Transmitter molecules bind to synaptic acetylcholine receptor channels (AChRs) to promote a global channel-opening conformational change. Although the detailed mechanism that links ligand binding and channel gating is uncertain, the energy changes caused by mutations appear to be more symmetrical between subunits in the transmembrane domain compared with the extracellular domain. The only covalent connection between these domains is the pre-M1 linker, a stretch of five amino acids that joins strand β10 with the M1 helix. In each subunit, this linker has a central Arg (Arg(3')), which only in the non-α-subunits is flanked by positively charged residues. Previous studies showed that mutations of Arg(3') in the α-subunit alter the gating equilibrium constant and reduce channel expression. We recorded single-channel currents and estimated the gating rate and equilibrium constants of adult mouse AChRs with mutations at the pre-M1 linker and the nearby residue Glu(45) in non-α-subunits. In all subunits, mutations of Arg(3') had similar effects as in the α-subunit. In the ε-subunit, mutations of the flanking residues and Glu(45) had only small effects, and there was no energy coupling between εGlu(45) and εArg(3'). The non-α-subunit Arg(3') residues had Φ-values that were similar to those for the α-subunit. The results suggest that there is a general symmetry between the AChR subunits during gating isomerization in this linker and that the central Arg is involved in expression more so than gating. The energy transfer through the AChR during gating appears to mainly involve Glu(45), but only in the α-subunits.

Acetylcholine receptor channels activated by a single agonist molecule

The neuromuscular acetylcholine receptor (AChR) is an allosteric protein that alternatively adopts inactive versus active conformations (R<-->R). The R shape has a higher agonist affinity and ionic conductance than R. To understand how agonists trigger this gating isomerization, we examined single-channel currents from adult mouse muscle AChRs that isomerize normally without agonists but have only a single site able to use agonist binding energy to motivate gating. We estimated the monoliganded gating equilibrium constant E(1) and the energy change associated with the R versus R change in affinity for agonists. AChRs with only one operational binding site gave rise to a single population of currents, indicating that the two transmitter binding sites have approximately the same affinity for the transmitter ACh. The results indicated that E(1) approximately 4.3 x 10(-3) with ACh, and approximately 1.7 x 10(-4) with the partial-agonist choline. From these values and the diliganded gating equilibrium constants, we estimate that the unliganded AChR gating constant is E(0) approximately 6.5 x 10(-7). Gating changes the stability of the ligand-protein complex by approximately 5.2 kcal/mol for ACh and approximately 3.3 kcal/mol for choline.

Energetics of gating at the apo-acetylcholine receptor transmitter binding site

Acetylcholine receptor channels switch between conformations that have a low versus high affinity for the transmitter and conductance for ions (R↔R*; gating). The forward isomerization, which begins at the transmitter binding sites and propagates ∼50 Å to the narrow region of the pore, occurs by approximately the same sequence of molecular events with or without agonists present at the binding sites. To pinpoint the forces that govern the R versus R* agonist affinity ratio, we measured single-channel activation parameters for apo-receptors having combinations of mutations of 10 transmitter binding site residues in the α (Y93, G147, W149, G153, Y190, C192, and Y198), ε (W55 and P121), or δ (W57) subunit. Gating energy changes were largest for the tryptophan residues. The αW149 energy changes were coupled with those of the other aromatic amino acids. Mutating the aromatic residues to Phe reduces the R/R* equilibrium dissociation constant ratio, with αY190 and αW149 being the most sensitive positions. Most of the mutations eliminated long-lived spontaneous openings. The results provide a foundation for understanding how ligands trigger protein conformational change.

Disruption of an intersubunit electrostatic bond is a critical step in glycine receptor activation

Proper regulation of neurotransmission requires that ligand-activated ion channels remain closed until agonist binds. How channels then open remains poorly understood. Glycine receptor (GlyR) gating is initiated by agonist binding at interfaces between adjacent subunits in the extracellular domain. Aspartate-97, located at the alpha1 GlyR interface, is a conserved residue in the cys-loop receptor superfamily. The mutation of D97 to arginine (D97R) causes spontaneous channel opening, with open and closed dwell times similar to those of maximally activated WT GlyR. Using a model of the N-terminal domain of the alpha1 GlyR, we hypothesized that an arginine-119 residue was forming intersubunit electrostatic bonds with D97. The D97R/R119E charge reversal restored this interaction, stabilizing channels in their closed states. Cysteine substitution shows that this link occurs between adjacent subunits. This intersubunit electrostatic interaction among GlyR subunits thus contributes to the stabilization of the closed channel state, and its disruption represents a critical step in GlyR activation.

Single-molecule imaging of a fluorescent unnatural amino acid incorporated into nicotinic receptors

We report on the first, to our knowledge, successful detection of a fluorescent unnatural amino acid (fUAA), Lys(BODIPYFL), incorporated into a membrane protein (the muscle nicotinic acetylcholine receptor, nAChR) in a living cell. Xenopus oocytes were injected with a frameshift-suppressor tRNA, amino-acylated with Lys(BODIPYFL) and nAChR (alpha/beta19'GGGU/gamma/delta) mRNAs. We measured fluorescence from oocytes expressing nAChR beta19'Lys(BODIPYFL), using time-resolved total internal reflection fluorescence microscopy. Under conditions of relatively low receptor density (<0.1 receptors/microm(2)), we observed puncta with diffraction-limited profiles that were consistent with the point-spread function of our microscope. Furthermore, diffraction-limited puncta displayed step decreases in fluorescence intensity, consistent with single-molecule photobleaching. The puncta densities agreed with macroscopic ACh-induced current densities, showing that the fUAA was incorporated, and that receptors were functional. Dose-response relations for the nAChR beta19'Lys(BODIPYFL) receptors were similar to those for wild-type receptors. We also studied nAChR beta19'Lys(BODIPYFL) receptors labeled with alpha-bungarotoxin monoconjugated with Alexa488 (alphaBtxAlexa488). The nAChR has two alphaBtx binding sites, and puncta containing the Lys(BODIPYFL) labeled with alphaBtxAlexa488 yielded the expected three discrete photobleaching steps. We also performed positive control experiments with a nAChR containing enhanced green fluorescent protein in the gamma-subunit M3-M4 loop, which confirmed our nAChR beta19'Lys(BODIPYFL) measurements. Thus, we report on the cell-based single-molecule detection of nAChR beta19'Lys(BODIPYFL).

Energy and structure of the M2 helix in acetylcholine receptor-channel gating

We studied single-channel currents from neuromuscular acetylcholine receptor-channels with mutations in the pore-lining, M2 helix of the epsilon-subunit. Three parameters were quantified: 1), the diliganded gating equilibrium constant (E(2)), which reflects the energy difference between C(losed) and O(pen) conformations; 2), the correlation between the opening rate constant and E(2) on a log-log scale (Phi), which illuminates the energy character of the residue (C- versus O-like) within the C<-->O isomerization process; and 3), the open-channel current amplitude (i(0)), which reports whether a mutation alters the energetics of ion permeation. The largest E(2) changes were observed in the cytoplasmic half of epsilonM2 (5', 9', 12', 13', and 16'), with smaller changes apparent for residues > or =17'. Phi was approximately 0.54 for most epsilonM2 residues, but was approximately 0.32 at the positions that had largest E(2) changes. An arginine substitution reduced i(0) significantly at six positions, with the magnitude of the reduction increasing, 16'-->2'. The measurements suggest that the 9', 12', and 13' residues experience large and late free-energy changes in the channel-opening process. We speculate that in the gating isomerization the pore-facing residues >6' and <16' experience multiple energy perturbations associated with changes in protein structure and, perhaps, hydration.

Long-range coupling in an allosteric receptor revealed by mutant cycle analysis

The functional coupling of residues that are far apart in space is the quintessential property of allosteric proteins. For example, in Cys-loop receptors, the gating of an intrinsic ion channel is allosterically regulated by the binding of small molecule neurotransmitters 50-60 A from the channel gate. Some residues near the binding site must have as their primary function the communication of the binding event to the gating region. These gating pathway residues are essential to function, but their identification and characterization can be challenging. This work introduces a simple strategy, derived from mutant cycle analysis, for identifying gating pathway residues using macroscopic measurements alone. In the exemplar Cys-loop receptor, the nicotinic acetylcholine receptor, a well-characterized reporter mutation (betaL9'S) known to impact gating, was combined with mutations of target residues in the ligand-binding domain hypothesized or previously found to be functionally significant. A mutant cycle analysis of the macroscopic EC(50) measurements can then provide insights into the role of the target residue. This new method, elucidating long-range functional coupling in allosteric receptors, can be applied to several reporter mutations in a wide variety of receptors to identify previously characterized and novel mutations that impact the gating pathway. We support our interpretation of macroscopic data with single-channel studies. Elucidating long-range functional coupling in allosteric receptors should be broadly applicable to determining functional roles of residues in allosteric receptors.

Functional equivalence of the nicotinic acetylcholine receptor transmitter binding sites in the open state

The subunits of the muscle-type nicotinic acetylcholine receptor (AChR) are not uniformly oriented in the resting closed conformation: the two alpha subunits are rotated relative to its non-alpha subunits. In contrast, all the subunits overlay well with one another when agonist is bound to the AChR, suggesting that they are uniformly oriented in the open receptor. This gating-dependent increase in orientational uniformity due to rotation of the alpha subunits might affect the relative affinities of the two transmitter binding sites, making the two affinities dissimilar (functionally non-equivalent) in the initial ligand-bound closed state but similar (functionally equivalent) in the open state. To test this hypothesis, we measured single-channel activity of the alphaG153S gain-of-function mutant receptor evoked by choline, and estimated the resting closed-state and open-state affinities of the two transmitter binding sites. Both model-independent analyses and maximum-likelihood estimation of microscopic rate constants indicate that channel opening makes the binding sites' affinities more similar to each other. These results support the hypothesis that open-state affinities to the transmitter binding sites are primarily determined by the alpha subunits.

Probing protein packing surrounding the residues in and flanking the nicotinic acetylcholine receptor M2M3 loop

Nicotinic acetylcholine receptors (nAChR) are cation-selective, ligand-gated ion channels of the cysteine (Cys)-loop gene superfamily. The recent crystal structure of a bacterial homolog from Erwinia chrysanthemi (ELIC) agrees with previous structures of the N-terminal domain of AChBP (acetylcholine-binding protein) and of the electron-microscopy-derived Torpedo nAChR structure. However, the ELIC transmembrane domain is significantly more tightly packed than the corresponding region of the Torpedo nAChR. We investigated the tightness of protein packing surrounding the extracellular end of the M2 transmembrane segment and around the loop connecting the M2 and M3 segments using the substituted cysteine accessibility method. The M2 20' to 27' residues were highly water accessible and the variation in reaction rates were consistent with this region being alpha-helical. At all positions tested, the presence of ACh changed methanethiosulfonate ethylammonium (MTSEA) modification rates by <10-fold. In the presence of ACh, reaction rates for residues in the last extracellular alpha-helical turn of M2 and in the M2M3 loop increased, whereas rates in the penultimate alpha-helical turn of M2 decreased. Only three of eight M2M3 loop residues were accessible to MTSEA in both the presence and absence of ACh. We infer that the protein packing around the M2M3 loop is tight, consistent with its location at the interdomain interface where it is involved in the transduction of ligand binding in the extracellular domain to gating in the transmembrane domain. Our data indicate that the Torpedo nAChR transmembrane domain structure is a better model than the ELIC structure for eukaryotic Cys-loop receptors.

Mammalian nicotinic acetylcholine receptors: from structure to function

The classical studies of nicotine by Langley at the turn of the 20th century introduced the concept of a "receptive substance," from which the idea of a "receptor" came to light. Subsequent studies aided by the Torpedo electric organ, a rich source of muscle-type nicotinic receptors (nAChRs), and the discovery of alpha-bungarotoxin, a snake toxin that binds pseudo-irreversibly to the muscle nAChR, resulted in the muscle nAChR being the best characterized ligand-gated ion channel hitherto. With the advancement of functional and genetic studies in the late 1980s, the existence of nAChRs in the mammalian brain was confirmed and the realization that the numerous nAChR subtypes contribute to the psychoactive properties of nicotine and other drugs of abuse and to the neuropathology of various diseases, including Alzheimer's, Parkinson's, and schizophrenia, has since emerged. This review provides a comprehensive overview of these findings and the more recent revelations of the impact that the rich diversity in function and expression of this receptor family has on neuronal and nonneuronal cells throughout the body. Despite these numerous developments, our understanding of the contributions of specific neuronal nAChR subtypes to the many facets of physiology throughout the body remains in its infancy.

Gating at the mouth of the acetylcholine receptor channel: energetic consequences of mutations in the alphaM2-cap

Gating of nicotinic acetylcholine receptors from a C(losed) to an O(pen) conformation is the initial event in the postsynaptic signaling cascade at the vertebrate nerve-muscle junction. Studies of receptor structure and function show that many residues in this large, five-subunit membrane protein contribute to the energy difference between C and O. Of special interest are amino acids located at the two transmitter binding sites and in the narrow region of the channel, where C↔O gating motions generate a low↔high change in the affinity for agonists and in the ionic conductance, respectively. We have measured the energy changes and relative timing of gating movements for residues that lie between these two locations, in the C-terminus of the pore-lining M2 helix of the α subunit (‘αM2-cap’). This region contains a binding site for non-competitive inhibitors and a charged ring that influences the conductance of the open pore. αM2-cap mutations have large effects on gating but much smaller effects on agonist binding, channel conductance, channel block and desensitization. Three αM2-cap residues (αI260, αP265 and αS268) appear to move at the outset of channel-opening, about at the same time as those at the transmitter binding site. The results suggest that the αM2-cap changes its secondary structure to link gating motions in the extracellular domain with those in the channel that regulate ionic conductance.

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.

Acetylcholine receptor gating: movement in the alpha-subunit extracellular domain

Acetylcholine receptor channel gating is a brownian conformational cascade in which nanometer-sized domains (“Φ blocks”) move in staggering sequence to link an affinity change at the transmitter binding sites with a conductance change in the pore. In the α-subunit, the first Φ-block to move during channel opening is comprised of residues near the transmitter binding site and the second is comprised of residues near the base of the extracellular domain. We used the rate constants estimated from single-channel currents to infer the gating dynamics of Y127 and K145, in the inner and outer sheet of the β-core of the α-subunit. Y127 is at the boundary between the first and second Φ blocks, at a subunit interface. αY127 mutations cause large changes in the gating equilibrium constant and with a characteristic Φ-value (Φ = 0.77) that places this residue in the second Φ-block. We also examined the effect on gating of mutations in neighboring residues δI43 (Φ = 0.86), ɛN39 (complex kinetics), αI49 (no effect) and in residues that are homologous to αY127 on the ɛ, β, and δ subunits (no effect). The extent to which αY127 gating motions are coupled to its neighbors was estimated by measuring the kinetic and equilibrium constants of constructs having mutations in αY127 (in both α subunits) plus residues αD97 or δI43. The magnitude of the coupling between αD97 and αY127 depended on the αY127 side chain and was small for both H (0.53 kcal/mol) and C (−0.37 kcal/mol) substitutions. The coupling across the single α–δ subunit boundary was larger (0.84 kcal/mol). The Φ-value for K145 (0.96) indicates that its gating motion is correlated temporally with the motions of residues in the first Φ-block and is not synchronous with those of αY127. This suggests that the inner and outer sheets of the α-subunit β-core do not rotate as a rigid body.

Acetylcholine receptor gating at extracellular transmembrane domain interface: the "pre-M1" linker

Charged residues in the beta10-M1 linker region ("pre-M1") are important in the expression and function of neuromuscular acetylcholine receptors (AChRs). The perturbation of a salt bridge between pre-M1 residue R209 and loop 2 residue E45 has been proposed as being a principle event in the AChR gating conformational "wave." We examined the effects of mutations to all five residues in pre-M1 (positions M207-P211) plus E45 in loop 2 in the mouse alpha(1)-subunit. M207, Q208, and P211 mutants caused small (approximately threefold) changes in the gating equilibrium constant (K(eq)), but the changes for R209, L210, and E45 were larger. Of 19 different side chain substitutions at R209 on the wild-type background, only Q, K, and H generated functional channels, with the largest change in K(eq) (67-fold) from R209Q. Various R209 mutants were functional on different E45 backgrounds: H, Q, and K (E45A), H, A, N, and Q (E45R), and K, A, and N (E45L). Phi values for R209 (on the E45A background), L210, and E45 were 0.74, 0.35, and 0.80, respectively. Phi values for R209 on the wt and three other backgrounds could not be estimated because of scatter. The average coupling energy between 209/45 side chains (six different pairs) was only -0.33 kcal/mol (for both alpha subunits, combined). Pre-M1 residues are important for expression of functional channels and participate in gating, but the relatively modest changes in closed- vs. open-state energy caused mutations, the weak coupling energy between these residues and the functional activity of several unmatched-charge pairs are not consistent with the perturbation of a salt bridge between R209 and E45 playing the principle role in gating.

Acetylcholine receptor gating at extracellular transmembrane domain interface: the cys-loop and M2-M3 linker

Acetylcholine receptor channel gating is a propagated conformational cascade that links changes in structure and function at the transmitter binding sites in the extracellular domain (ECD) with those at a “gate” in the transmembrane domain (TMD). We used Φ-value analysis to probe the relative timing of the gating motions of α-subunit residues located near the ECD–TMD interface. Mutation of four of the seven amino acids in the M2–M3 linker (which connects the pore-lining M2 helix with the M3 helix), including three of the four residues in the core of the linker, changed the diliganded gating equilibrium constant (K_eq) by up to 10,000-fold (P272 > I274 > A270 > G275). The average Φ-value for the whole linker was ∼0.64. One interpretation of this result is that the gating motions of the M2–M3 linker are approximately synchronous with those of much of M2 (∼0.64), but occur after those of the transmitter binding site region (∼0.93) and loops 2 and 7 (∼0.77). We also examined mutants of six cys-loop residues (V132, T133, H134, F135, P136, and F137). Mutation of V132, H134, and F135 changed K_eq by 2800-, 10-, and 18-fold, respectively, and with an average Φ-value of 0.74, similar to those of other cys-loop residues. Even though V132 and I274 are close, the energetic coupling between I and V mutants of these positions was small (≤0.51 kcal mol⁻¹). The M2–M3 linker appears to be the key moving part that couples gating motions at the base of the ECD with those in TMD. These interactions are distributed along an ∼16-Å border and involve about a dozen residues.

Principles underlying energetic coupling along an allosteric communication trajectory of a voltage-activated K+ channel

The information flow between distal elements of a protein may rely on allosteric communication trajectories lying along the protein's tertiary or quaternary structure. To unravel the underlying features of energy parsing along allosteric pathways in voltage-gated K(+) channels, high-order thermodynamic coupling analysis was performed. We report that such allosteric trajectories are functionally conserved and delineated by well defined boundaries. Moreover, allosteric trajectories assume a hierarchical organization whereby increasingly stronger layers of cooperative residue interactions act to ensure efficient and cooperative long-range coupling between distal channel regions. Such long-range communication is brought about by a coupling of local and global conformational changes, suggesting that the allosteric trajectory also corresponds to a pathway of physical deformation. Supported by theoretical analyses and analogy to studies analyzing the contribution of long-range residue coupling to protein stability, we propose that such experimentally derived trajectory features are a general property of allosterically regulated proteins.

Conformational dynamics of the alphaM3 transmembrane helix during acetylcholine receptor channel gating

Muscle acetylcholine receptors are synaptic ion channels that "gate" between closed- and open-channel conformations. We used Phi-value analysis to probe the transition state of the diliganded gating reaction with regard to residues in the M3, membrane-spanning helix of the muscle acetylcholine receptor alpha-subunit. Phi (a fraction between 1 and 0) parameterizes the extent to which a mutation changes the opening versus the closing rate constant and, for a linear reaction mechanism, the higher the Phi-value, the "earlier" the gating motion. In the upper half of alphaM3 the gating motions of all five tested residues were temporally correlated (Phi approximately 0.30) and serve to link structural changes occurring at the middle of the M2, pore-lining helix with those occurring at the interface of the extracellular and transmembrane domains. alphaM3 belongs to a complex and diverse set of synchronously moving parts that change structure relatively late in the channel-opening process. The propagation of the gating Brownian conformational cascade has a complex spatial distribution in the transmembrane domain.

Closed-state cross-linking of adjacent beta1 subunits in alpha1beta1 GABAa receptors via introduced 6' cysteines

The pore structural changes associated with Cys-loop receptor gating are currently the subject of considerable interest. Several functional approaches have shown that surface exposure of pore-lining side chains does not change significantly during activation. However, a disulfide trapping study on alpha1(T6'C)beta1(T6'C) gamma-aminobutyric acid type A (GABA(A)) receptors (GABA(A)Rs), which showed that adjacent beta subunits cross-link in the open state only, concluded that channel gating is mediated by beta subunits contra-rotating through a summed angle of approximately 120 degrees. Such a large rotation is not easily reconciled with other evidence. The present study initially sought to investigate an observation that appeared inconsistent with the rotation model: namely that alpha1(T6'C)beta1(T6'C) GABA(A)Rs expressed in HEK293 cells form 6' cysteine-mediated disulfide bonds in the closed state. On the basis of electrophysiological and Western blotting experiments, we conclude that adjacent beta(T6'C) subunits dimerise efficiently in the closed state via cross-links between their respective 6' cysteines and that this locks the channels closed. This questions the beta subunit contra-rotation model of activation and raises the question of how the closed state cross-links form. We propose that beta subunit 6' cysteines move into sufficiently close proximity for disulfide formation via relatively large amplitude random thermal motions that appear to be a unique feature of beta subunits. Because dimerized channels are locked closed, we conclude either that the spontaneous beta subunit movements or asymmetries in the movements of adjacent beta subunits during activation are essential for pore opening. Our results identify a novel feature of GABA(A)R gating that may be important for understanding its activation mechanism.

A stepwise mechanism for acetylcholine receptor channel gating

Muscle contraction is triggered by the opening of acetylcholine receptors at the vertebrate nerve-muscle synapse. The M2 helix of this allosteric membrane protein lines the channel, and contains a 'gate' that regulates the flow of ions through the pore. We used single-molecule kinetic analysis to probe the transition state of the gating conformational change and estimate the relative timing of M2 motions in the alpha-subunit of the murine acetylcholine receptor. This analysis produces a 'Phi-value' for a given residue that reflects its open-like versus closed-like character at the transition state. Here we show that most of the residues throughout the length of M2 have a Phi-value of approximately 0.64 but that some near the middle have lower Phi-values of 0.52 or 0.31, suggesting that alphaM2 moves in three discrete steps. The core of the channel serves both as a gate that regulates ion flow and as a hub that directs the propagation of the gating isomerization through the membrane domain of the acetylcholine receptor.

Modulation of proton-induced current fluctuations in the human nicotinic acetylcholine receptor channel

The nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel that switches upon activation from a closed state to a full conducting state. We found that the mutation delta S268K, located at 12' position of the second transmembrane domain of the delta subunit of the human nAChR generates a long-lived intermediate conducting state, from which openings to a wild-type like conductance level occur on a submillisecond time scale. Aiming to understand the interplay between structural changes near the 12' position and channel gating, we investigated the influence of various parameters: different ligands (acetylcholine, choline and epibatidine), ligand concentrations, transmembrane voltages and both fetal and adult nAChRs. Since sojourns in the high conductance state are not fully resolved in time, spectral noise analysis was used as a complement to dwell time analysis to determine the gating rate constants. Open channel current fluctuations are described by a two-state Markov model. The characteristic time of the process is markedly influenced by the ligand and the receptor type, whereas the frequency of openings to the high conductance state increases with membrane hyperpolarization. Conductance changes are discussed with regard to reversible transfer reaction of single protons at the lysine 12' side chain.

Influence of agonists and antagonists on the segmental motion of residues near the agonist binding pocket of the acetylcholine-binding protein

Using the Lymnaea acetylcholine-binding protein as a surrogate of the extracellular domain of the nicotinic receptor, we combined site-directed labeling with fluorescence spectroscopy to assess possible linkages between ligand binding and conformational dynamics. Specifically, 2-[(5-fluoresceinyl)aminocarbonyl]ethyl methanethiosulfonate was conjugated to a free cysteine on loop C and to five substituted cysteines at strategic locations in the subunit sequence, and the backbone flexibility around each site of conjugation was measured with time-resolved fluorescence anisotropy. The sites examined were in loop C (Cys-188 using a C187S mutant), in the beta9 strand (T177C), in the beta10 strand (D194C), in the beta8-beta9 loop (N158C and Y164C), and in the beta7 strand (K139C). Conjugated fluorophores at these locations show distinctive anisotropy decay patterns indicating different degrees of segmental fluctuations near the agonist binding pocket. Ligand occupation and decay of anisotropy were assessed for one agonist (epibatidine) and two antagonists (alpha-bungarotoxin and d-tubocurarine). The Y164C and Cys-188 conjugates were also investigated with additional agonists (nicotine and carbamylcholine), partial agonists (lobeline and 4-hydroxy,2-methoxy-benzylidene anabaseine), and an antagonist (methyllycaconitine). With the exception of the T177C conjugate, both agonists and antagonists perturbed the backbone flexibility of each site; however, agonist-selective changes were only observed at Y164C in loop F where the agonists and partial agonists increased the range and/or rate of the fast anisotropy decay processes. The results reveal that agonists and antagonists produced distinctive changes in the flexibility of a portion of loop F.

Role of pairwise interactions between M1 and M2 domains of the nicotinic receptor in channel gating

The adult form of the nicotinic acetylcholine receptor (AChR) consists of five subunits (alpha(2)betaepsilondelta), each having four transmembrane domains (M1-M4). The atomic model of the nicotinic acetylcholine receptor shows that the pore-lining M2 domains make no extensive contacts with the rest of the transmembrane domains. However, there are several sites where close appositions between segments occur. It has been suggested that the pair alphaM1-F15' and alphaM2-L11' is one of the potential interactions between segments. To determine experimentally if these residues are interacting and to explore if this interhelical interaction is essential for channel gating, we combined mutagenesis with single-channel kinetic analysis. Mutations in alphaM1-F15' lead to profound changes in the opening rate and slighter changes in the closing rate. Channel gating is impaired as the volume of the residue increases. Rate-equilibrium linear free-energy relationship analysis reveals an approximately 70% open-state-like environment for alphaM1-F15' at the transition state of the gating reaction, suggesting that it moves early during the gating process. Replacing the residue at alphaM1-15' by that at alphaM2-11' and vice versa profoundly alters gating, but the combination of the two mutations restores gating to near normal, indicating that alphaM1-F15' and alphaM2-L11' are interchangeable. Double-mutant cycle analysis shows that these residues are energetically coupled. Thus, the interaction between M1 and M2 plays a key role in channel gating.