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.     

Functional asymmetry of transmembrane segments in nicotinic acetylcholine receptors

Nicotinic acetylcholine receptors are heteropentameric ion channels that open upon activation to a single conducting state. The second transmembrane segments of each subunit were identified as channel-forming elements, but their respective contribution in the gating process remains unclear. Moreover, the detailed impact of variations of the membrane potential, such as occurring during an action potential, on the transmembrane domains, is unknown. Residues at the 12′ position, close to the center of each second transmembrane segment, play a key role in channel gating. We examined their functional symmetry by substituting a lysine to that position of each subunit and measuring the electrical activity of single channels. For 12′ lysines in the α, γ and δ subunits rapid transitions between an intermediate and large conductance appeared, which are interpreted as single lysine protonation events. From the kinetics of these transitions we calculated the pK _a values of respective lysines and showed that they vary differently with membrane hyperpolarization. Respective mutations in β or ε subunits gave receptors with openings of either intermediate or large conductance, suggesting extreme pK _a values in two open state conformations. The results demonstrate that these parts of the highly homologous transmembrane domains, as probed by the 12′ lysines, sense unequal microenvironments and are differently affected by physiologically relevant voltage changes. Moreover, observation of various gating events for mutants of α subunits suggests that the open channel pore exists in multiple conformations, which in turn supports the notion of functional asymmetry of the channel.     

Intramembrane Aromatic Interactions Influence the Lipid Sensitivities of Pentameric Ligand-gated Ion Channels

Although the Torpedo nicotinic acetylcholine receptor (nAChR) reconstituted into phosphatidylcholine (PC) membranes lacking cholesterol and anionic lipids adopts a conformation where agonist binding is uncoupled from channel gating, the underlying mechanism remains to be defined. Here, we examine the mechanism behind lipid-dependent uncoupling by comparing the propensities of two prokaryotic homologs, Gloebacter and Erwinia ligand-gated ion channel (GLIC and ELIC, respectively), to adopt a similar uncoupled conformation. Membrane-reconstituted GLIC and ELIC both exhibit folded structures in the minimal PC membranes that stabilize an uncoupled nAChR. GLIC, with a large number of aromatic interactions at the interface between the outermost transmembrane α-helix, M4, and the adjacent transmembrane α-helices, M1 and M3, retains the ability to flux cations in this uncoupling PC membrane environment. In contrast, ELIC, with a level of aromatic interactions intermediate between that of the nAChR and GLIC, does not undergo agonist-induced channel gating, although it does not exhibit the expected biophysical characteristics of the uncoupled state. Engineering new aromatic interactions at the M4-M1/M3 interface to promote effective M4 interactions with M1/M3, however, increases the stability of the transmembrane domain to restore channel function. Our data provide direct evidence that M4 interactions with M1/M3 are modulated during lipid sensing. Aromatic residues strengthen M4 interactions with M1/M3 to reduce the sensitivities of pentameric ligand-gated ion channels to their surrounding membrane environment.     

Kinetics of unliganded acetylcholine receptor channel gating.

Open- and closed-state lifetimes of unliganded acetylcholine receptor channel activity were analyzed by the method of likelihood maximazation. For both open times and closed times, the best-fitting density is most often a sum of two exponentials. These multiple open states cannot depend on the number of receptor binding sites occupied since they are observed in the absence of ligand. The rate of spontaneous opening and the faster decay constant of closing increased as the membrane was hyperpolarized. The voltage dependence of the rate of spontaneous opening is stronger than that for curare-liganded channels. Evidence that the acetylcholine receptor channel can open spontaneously in the absence of ligand has been presented previously (Sanchez et al, 1983; Brehm et al, 1984; Jackson, 1984). To add to this evidence, alpha-bungarotoxin was added to the patch electrode, causing the frequency of openings to decay with time. The rate constant determined from this decay is similar to rate constants reported for the binding of iodinated alpha-bungarotoxin to the acetylcholine receptor. The frequency of unliganded channel opening has been estimated as 2 X 10(-3) s-1 per receptor. A comparison of carbamylcholine-liganded and spontaneous gating transition rates suggests that ligand binding increases the rate of opening by a factor of 1.4 X 10(7). Carbamylcholine binding increases the mean open time by a factor of 5. Thus, a cholinergic agonist activates the acetylcholine receptor by destabilizing the closed state. The liganded and unliganded channel gating rates were used to analyze the energetics of ligand activation of the acetylcholine receptor channel, and to relate the open channel dissociation constant to the closed channel dissociation constant.     

Role of acetylcholine receptor domains in ion selectivity

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

Targeted molecular dynamics study of C-loop closure and channel gating in nicotinic receptors

The initial coupling between ligand binding and channel gating in the human α7 nicotinic acetylcholine receptor (nAChR) has been investigated with targeted molecular dynamics (TMD) simulation. During the simulation, eight residues at the tip of the C-loop in two alternating subunits were forced to move toward a ligand-bound conformation as captured in the crystallographic structure of acetylcholine binding protein (AChBP) in complex with carbamoylcholine. Comparison of apo- and ligand-bound AChBP structures shows only minor rearrangements distal from the ligand-binding site. In contrast, comparison of apo and TMD simulation structures of the nAChR reveals significant changes toward the bottom of the ligand-binding domain. These structural rearrangements are subsequently translated to the pore domain, leading to a partly open channel within 4 ns of TMD simulation. Furthermore, we confirmed that two highly conserved residue pairs, one located near the ligand-binding pocket (Lys145 and Tyr188), and the other located toward the bottom of the ligand-binding domain (Arg206 and Glu45), are likely to play important roles in coupling agonist binding to channel gating. Overall, our simulations suggest that gating movements of the α7 receptor may involve relatively small structural changes within the ligand-binding domain, implying that the gating transition is energy-efficient and can be easily modulated by agonist binding/unbinding.     

A unified view of cystic fibrosis transmembrane conductance regulator (CFTR) gating: combining the allosterism of a ligand-gated channel with the enzymatic activity of an ATP-binding cassette (ABC) transporter

The cystic fibrosis transmembrane conductance regulator (CFTR) is a unique ion channel in that its gating is coupled to an intrinsic enzymatic activity (ATP hydrolysis). This enzymatic activity derives from the evolutionary origin of CFTR as an ATP-binding cassette transporter. CFTR gating is distinct from that of a typical ligand-gated channel because its ligand (ATP) is usually consumed during the gating cycle. However, recent findings indicate that CFTR gating exhibits allosteric properties that are common to conventional ligand-gated channels (e.g. unliganded openings and constitutive mutations). Here, we provide a unified view of CFTR gating that combines the allosterism of a ligand-gated channel with its unique enzymatic activity.     

A hydrophobic gate in an ion channel: the closed state of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is the prototypic member of the 'Cys-loop' superfamily of ligand-gated ion channels which mediate synaptic neurotransmission, and whose other members include receptors for glycine, gamma-aminobutyric acid and serotonin. Cryo-electron microscopy has yielded a three-dimensional structure of the nAChR in its closed state. However, the exact nature and location of the channel gate remains uncertain. Although the transmembrane pore is constricted close to its center, it is not completely occluded. Rather, the pore has a central hydrophobic zone of radius about 3 A. Model calculations suggest that such a constriction may form a hydrophobic gate, preventing movement of ions through a channel. We present a detailed and quantitative simulation study of the hydrophobic gating model of the nicotinic receptor, in order to fully evaluate this hypothesis. We demonstrate that the hydrophobic constriction of the nAChR pore indeed forms a closed gate. Potential of mean force (PMF) calculations reveal that the constriction presents a barrier of height about 10 kT to the permeation of sodium ions, placing an upper bound on the closed channel conductance of 0.3 pS. Thus, a 3 A radius hydrophobic pore can form a functional barrier to the permeation of a 1 A radius Na+ ion. Using a united-atom force field for the protein instead of an all-atom one retains the qualitative features but results in differing conductances, showing that the PMF is sensitive to the detailed molecular interactions.     

Subunit-selective role of the M3 transmembrane domain of the nicotinic acetylcholine receptor in channel gating

The nicotinic acetylcholine receptor (AChR) can be either hetero-pentameric, composed of α and non-α subunits, or homo-pentameric, composed of α7 subunits. To explore the subunit-selective contributions of transmembrane domains to channel gating we analyzed single-channel activity of chimeric muscle AChRs. We exchanged M3 between α1 and ? or α7 subunits. The replacement of M3 in α1 by ?M3 significantly alters activation properties. Channel activity appears as bursts of openings whose durations are 20-fold longer than those of wild-type AChRs. In contrast, 7-fold briefer openings are observed in AChRs containing the reverse ? chimeric subunit. The duration of the open state decreases with the increase in the number of α1M3 segments, indicating additive contributions of M3 of all subunits to channel closing. Each α1M3 segment decreases the energy barrier of the closing process by ∼ 0.8 kcal/mol. Partial chimeric subunits show that small stretches of the M3 segment contribute additively to the open duration. The replacement of α1 sequence by α7 in M3 leads to 3-fold briefer openings whereas in M1 it leads to 10-fold prolonged openings, revealing that the subunit-selective role is unique to each transmembrane segment.     

Characterization of the channel properties of a neuronal acetylcholine receptor reconstituted into planar lipid bilayers

An alpha-toxin-binding membrane protein, isolated from the head and thoracic ganglia of the locus (Locusta migratoria), was reconstituted into planar lipid bilayers. Cholinergic agonists such as acetylcholine, carbamylcholine, and suberyldicholine induced fluctuations of single channels, which suggests that the protein represents a functional cholinergic receptor channel. The antagonist d-tubocurarine blocked the activation of the channels, whereas hexamethonium had only a weak effect; similar properties have been described for nicotinic insect receptors in situ. The channel was selectively permeable to monovalent cations but was impermeable to anions. The conductance of the channel (75 pS in 100 mM NaCl) was independent of the type of agonist used to activate the receptor. Kinetic analysis of the channel gating revealed that, at high agonist concentrations (50 microM carbamylcholine), more than one closed state exists and that multiple gating events, bursting as well as fast flickering, appeared. At very high agonist concentrations (500 microM carbamylcholine), desensitization was observed. Channel kinetics were dependent on the transmembrane potential. Comparing the conductance, the kinetics, and the pharmacology of nicotinic acetylcholine receptor from insect ganglia and fish electroplax reconstituted into bilayers revealed obvious similarities but also significant differences.     

Computational study of the transmembrane domain of the acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel protein whose transmembrane domain (TM-domain) is believed to be responsible for channel gating via a hydrophobic effect. In this work, we perform molecular dynamics and Brownian dynamics simulations to investigate the effect of transmembrane potential on the conformation and water occupancy of TM-domain, and the resulting ion permeation events. The results show that the behavior of the hydrophobic gate is voltage-dependent. Large hyperpolarized membrane potential can change the conformation of TM-domain and water occupancy in this region, which may enable ion conduction. An electrostatic gating mechanism is also proposed from our simulations, which seems to play a role in addition to the well-known hydrophobic effect.     

Mechanics of channel gating of the nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is a key molecule involved in the propagation of signals in the central nervous system and peripheral synapses. Although numerous computational and experimental studies have been performed on this receptor, the structural dynamics of the receptor underlying the gating mechanism is still unclear. To address the mechanical fundamentals of nAChR gating, both conventional molecular dynamics (CMD) and steered rotation molecular dynamics (SRMD) simulations have been conducted on the cryo-electron microscopy (cryo-EM) structure of nAChR embedded in a dipalmitoylphosphatidylcholine (DPPC) bilayer and water molecules. A 30-ns CMD simulation revealed a collective motion amongst C-loops, M1, and M2 helices. The inward movement of C-loops accompanying the shrinking of acetylcholine (ACh) binding pockets induced an inward and upward motion of the outer β-sheet composed of β9 and β10 strands, which in turn causes M1 and M2 to undergo anticlockwise motions around the pore axis. Rotational motion of the entire receptor around the pore axis and twisting motions among extracellular (EC), transmembrane (TM), and intracellular MA domains were also detected by the CMD simulation. Moreover, M2 helices undergo a local twisting motion synthesized by their bending vibration and rotation. The hinge of either twisting motion or bending vibration is located at the middle of M2, possibly the gate of the receptor. A complementary twisting-to-open motion throughout the receptor was detected by a normal mode analysis (NMA). To mimic the pulsive action of ACh binding, nonequilibrium MD simulations were performed by using the SRMD method developed in one of our laboratories. The result confirmed all the motions derived from the CMD simulation and NMA. In addition, the SRMD simulation indicated that the channel may undergo an open-close (O ↔ C) motion. The present MD simulations explore the structural dynamics of the receptor under its gating process and provide a new insight into the gating mechanism of nAChR at the atomic level.     

Acetylcholine Receptor Gating is Influenced by the Polarity of Amino Acids at Position 9′ in the M2 Domain

Ligand-gated ion channels contain a conserved leucine at position 9′ (L9′) in the M2 transmembrane domain. We used multiple substitutions at this position in the γ subunit of the mouse acetylcholine receptor (AChR) (γL9′) to examine the role of residue polarity at this position in the gating process at both the macroscopic and single-channel levels. The midpoint of the macroscopic dose-response relationship (EC₅₀) and the channel closing rate constant, α, decreased as the polarity of the residue at that position increased, suggesting a stabilization of the open state of the channel. Both parameters showed similar dependencies on the polarity of the substituted residue. These data support the notion that during AChR gating, the amino acid at the 9′ position moves into a polar environment, and that interactions between this residue and the polar environment determine the stability of the open state. Since this residue is conserved in all other members of the ligand-gated ion channel family, we suggest that a similar mechanism applies to the other members of the family. Received: 17 September 1999/Revised: 15 December 1999     

An anionic ryanoid, 10-O-succinoylryanodol,provides insights into the mechanisms governing the interaction of ryanoids and the subsequent altered function of ryanodine-receptor channels

We have investigated the interactions of a novel anionic ryanoid, 10-O-succinoylryanodol, with individual mammalian cardiac muscle ryanodine receptor channels under voltage clamp conditions. As is the case for all ryanoids so far examined, the interaction of 10-O-succinoylryanodol with an individual RyR channel produces profound alterations in both channel gating and rates of ion translocation. In the continued presence of the ryanoid the channel fluctuates between periods of normal and modified gating, indicating a reversible interaction of the ligand with its receptor. Unlike the majority of ryanoids, we observe a range of different fractional conductance states of RyR in the presence of 10-O-succinoylryanodol. We demonstrate that 10-O-succinoylryanodol is a very flexible molecule and propose that each fractional conductance state arises from the interaction of a different conformer of the ryanoid molecule with the RyR channel. The probability of channel modification by 10-O-succinoylryanodol is dependent on the transmembrane holding potential. Comparison of the voltage dependence of channel modification by this novel anionic ryanoid with previous data obtained with cationic and neutral ryanoids reveals that the major influence of transmembrane potential on the probability of RyR channel modification by ryanoids results from an alteration in receptor affinity. These investigations also demonstrate that the charge of the ryanoid has a major influence on the rate of association of the ligand with its receptor indicating that ionic interactions are likely to be involved in this reaction.     

Fourier transform coupled tryptophan scanning mutagenesis identifies a bending point on the lipid-exposed δM3 transmembrane domain of the Torpedo californica nicotinic acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is a member of a family of ligand-gated ion channels that mediate diverse physiological functions, including fast synaptic transmission along the peripheral and central nervous systems. Several studies have made significant advances toward determining the structure and dynamics of the lipid-exposed domains of the nAChR. However, a high-resolution atomic structure of the nAChR still remains elusive. In this study, we extended the Fourier transform coupled tryptophan scanning mutagenesis (FT-TrpScanM) approach to gain insight into the secondary structure of the dM3 transmembrane domain of the Torpedo californica nAChR, to monitor conformational changes experienced by this domain during channel gating, and to identify which lipid-exposed positions are linked to the regulation of ion channel kinetics. The perturbations produced by periodic tryptophan substitutions along the dM3 transmembrane domain were characterized by two-electrode voltage clamp and ¹²⁵I-labeled a-bungarotoxin binding assays. The periodicity profiles and Fourier transform spectra of this domain revealed similar helical structures for the closed- and open-channel states. However, changes in the oscillation patterns observed between positions Val-299 and Val-304 during transition between the closed- and open-channel states can be explained by the structural effects caused by the presence of a bending point introduced by a Thr-Gly motif at positions 300-301. The changes in periodicity and localization of residues between the closed-and open-channel states could indicate a structural transition between helix types in this segment of the domain. Overall, the data further demonstrate a functional link between the lipid–exposed transmembrane domain and the nAChR gating machinery.     

Modeling Neuronal Nicotinic and GABA Receptors: Important Interface Salt-Links and Protein Dynamics

Protein motions in the Cys-loop ligand-gated ion receptors that govern the gating mechanism are still not well understood. The details as to how motions in the ligand-binding domain are translated to the transmembrane domain and how subunit rotations are linked to bring about the cooperative movements involved in gating are under investigation. Homology models of the α4β2 nicotinic acetylcholine (nACh) and β2α1γ2 GABA receptors were constructed based on the torpedo neuromuscular-like nicotinic receptor structure. The template constructed for the full electron microscopy structure must be considered more reliable for structure-function studies due to the preservation of the E45-R209 salt-link. Many other salt-links are seen to transiently form, including switching off of the E45-R209 link, within a network of potential salt-links at the binding domain to the transmembrane domain interface region. Several potentially important intersubunit salt-links form in both the nAChR and GABAR structures during the simulation and appear conserved across many subunit combinations, such as the salt-link between α4.E262 and β2.K255 in nAChR (β2.E262 and α1.K263 in GABAR), at the top of the pore-lining M2 helices, and the intersubunit link of R210 on the M1-linker to E168 on the β8-sheet of the adjacent subunit in the GABA receptor (E175-K46 being the structurally equivalent link in the nAChR, with reversed polarity). A network of other salt-links may be vital for transmitting the cooperative gating motions between subunits that become biased upon ligand binding. The changes seen in the simulations suggest that this network of salt-links helps to set limits and specific states for the conformational changes involved in gating of the receptor. We hope that these hypotheses will be tested experimentally in the near future.     

The extracellular linker of muscle acetylcholine receptor channels is a gating control element

We describe the functional consequences of mutations in the linker between the second and third transmembrane segments (M2-M3L) of muscle acetylcholine receptors at the single-channel level. Hydrophobic mutations (Ile, Cys, and Phe) placed near the middle of the linker of the alpha subunit (alphaS269) prolong apparent openings elicited by low concentrations of acetylcholine (ACh), whereas hydrophilic mutations (Asp, Lys, and Gln) are without effect. Because the gating kinetics of the alphaS269I receptor (a congenital myasthenic syndrome mutant) in the presence of ACh are too fast, choline was used as the agonist. This revealed an approximately 92-fold increased gating equilibrium constant, which is consistent with an approximately 10-fold decreased EC(50) in the presence of ACh. With choline, this mutation accelerates channel opening approximately 28-fold, slows channel closing approximately 3-fold, but does not affect agonist binding to the closed state. These ratios suggest that, with ACh, alphaS269I acetylcholine receptors open at a rate of approximately 1.4 x 10(6) s(-1) and close at a rate of approximately 760 s(-1). These gating rate constants, together with the measured duration of apparent openings at low ACh concentrations, further suggest that ACh dissociates from the diliganded open receptor at a rate of approximately 140 s(-1). Ile mutations at positions flanking alphaS269 impair, rather than enhance, channel gating. Inserting or deleting one residue from this linker in the alpha subunit increased and decreased, respectively, the apparent open time approximately twofold. Contrary to the alphaS269I mutation, Ile mutations at equivalent positions of the beta, straightepsilon, and delta subunits do not affect apparent open-channel lifetimes. However, in beta and straightepsilon, shifting the mutation one residue to the NH(2)-terminal end enhances channel gating. The overall results indicate that this linker is a control element whose hydrophobicity determines channel gating in a position- and subunit-dependent manner. Characterization of the transition state of the gating reaction suggests that during channel opening the M2-M3L of the alpha subunit moves before the corresponding linkers of the beta and straightepsilon subunits.     

Channel opening motion of alpha7 nicotinic acetylcholine receptor as suggested by normal mode analysis

The gating motion of the human nicotinic acetylcholine receptor (nAChR) alpha7 was investigated with normal mode analysis (NMA) of two homology models. The first model, referred to as model I, was built from both the Lymnaea stagnalis acetylcholine binding protein (AChBP) and the transmembrane (TM) domain of the Torpedo marmorata nAChR. The second model, referred to as model C, was based solely on the recent electron microscopy structure of the T. marmorata nAChR. Despite structural differences, both models exhibit nearly identical patterns of flexibility and correlated motions. In addition, both models show a similar global twisting motion that may represent channel gating. The similar results obtained for the two models indicate that NMA is most sensitive to the contact topology of the structure rather than its finer detail. The major difference between the low-frequency motions sampled for the two models is that a symmetrical pore-breathing motion, favoring channel opening, is present as the second most dominant motion in model I, whilst largely absent from model C. The absence of this mode in model C can be attributed to its less symmetrical architecture. Finally, as a further goal of the present study, an approximate open channel model, consistent with many experimental findings, has been produced.     

Single-Channel Current Through Nicotinic Receptor Produced by Closure of Binding Site C-Loop

We investigated the initial coupling of agonist binding to channel gating of the nicotinic acetylcholine receptor using targeted molecular-dynamics (TMD) simulation. After TMD simulation to accelerate closure of the C-loops at the agonist binding sites, the region of the pore that passes through the cell membrane expands. To determine whether the structural changes in the pore result in ion conduction, we used a coarse-grained ion conduction simulator, Biology Boltzmann transport Monte Carlo, and applied it to two structural frames taken before and after TMD simulation. The structural model before TMD simulation represents the channel in the proposed “resting” state, whereas the model after TMD simulation represents the channel in the proposed “active” state. Under external voltage biases, the channel in the “active” state was permeable to cations. Our simulated ion conductance approaches that obtained experimentally and recapitulates several functional properties characteristic of the nicotinic acetylcholine receptor. Thus, closure of the C-loop triggers a structural change in the channel sufficient to account for the open channel current. This approach of applying Biology Boltzmann transport Monte Carlo simulation can be used to further investigate the binding to gating transduction mechanism and the structural bases for ion selection and translocation.