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.     

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.     

The role of loop 5 in acetylcholine receptor channel gating

Nicotinic acetylcholine receptor channel (AChR) gating is an organized sequence of molecular motions that couples a change in the affinity for ligands at the two transmitter binding sites with a change in the ionic conductance of the pore. Loop 5 (L5) is a nine-residue segment (mouse alpha-subunit 92-100) that links the beta4 and beta5 strands of the extracellular domain and that (in the alpha-subunit) contains binding segment A. Based on the structure of the acetylcholine binding protein, we speculate that in AChRs L5 projects from the transmitter binding site toward the membrane along a subunit interface. We used single-channel kinetics to quantify the effects of mutations to alphaD97 and other L5 residues with respect to agonist binding (to both open and closed AChRs), channel gating (for both unliganded and fully-liganded AChRs), and desensitization. Most alphaD97 mutations increase gating (up to 168-fold) but have little or no effect on ligand binding or desensitization. Rate-equilibrium free energy relationship analysis indicates that alphaD97 moves early in the gating reaction, in synchrony with the movement of the transmitter binding site (Phi = 0.93, which implies an open-like character at the transition state). alphaD97 mutations in the two alpha-subunits have unequal energetic consequences for gating, but their contributions are independent. We conclude that the key, underlying functional consequence of alphaD97 perturbations is to increase the unliganded gating equilibrium constant. L5 emerges as an important and early link in the AChR gating reaction which, in the absence of agonist, serves to increase the relative stability of the closed conformation of the protein.     

Photoactivation and dissociation of agonist molecules at the nicotinic acetylcholine receptor in voltage-clamped rat myoballs.

The photochemical properties of the azobenzene derivative, Bis-Q, were exploited to carry out an agonist concentration jump followed by a molecular rearrangement of bound agonist molecules at acetylcholine (ACh) receptor channels of voltage-clamped rat myoballs. Myoballs were bathed in solutions containing low concentrations of cis-Bis-Q, the inactive isomer. Whole-cell current relaxations were studied following a light flash that produced a concentration jump of agonist, trans-Bis-Q, followed by a second flash that produced net trans----cis photoisomerizations of Bis-Q molecules. The concentration-jump relaxation provided a measure of the mean burst duration for ACh receptor channels occupied by trans-Bis-Q (7.7 ms, 22 degrees C). The second current relaxation was a more rapid conductance decrease (phase 1, tau = 0.8 ms). Phase 1 may represent either the burst duration for receptors initially occupied by a single cis- and a single trans-Bis-Q molecule or that for unliganded receptors. Single-channel current recordings from excised outside-out membrane patches showed that single channels open following an agonist concentration jump comparable to that used in the whole-cell experiments; when many such records were averaged, a synthetic macroscopic relaxation was produced. Individual open channels closed faster following a flash that promoted trans----cis photoisomerizations of the bound ligand, thus confirming the whole-cell observations of phase 1.     

The intrinsic energy of the gating isomerization of a neuromuscular acetylcholine receptor channel

Nicotinic acetylcholine receptor (AChR) channels at neuromuscular synapses rarely open in the absence of agonists, but many different mutations increase the unliganded gating equilibrium constant (E0) to generate AChRs that are active constitutively. We measured E0 for two different sets of mutant combinations and by extrapolation estimated E0 for wild-type AChRs. The estimates were 7.6 and 7.8×10(-7) in adult-type mouse AChRs (-100 mV at 23°C). The values are in excellent agreement with one obtained previously by using a completely different method (6.5×10(-7), from monoliganded gating). E0 decreases with depolarization to the same extent as does the diliganded gating equilibrium constant, e-fold with ～60 mV. We estimate that at -100 mV the intrinsic energy of the unliganded gating isomerization is +8.4 kcal/mol (35 kJ/mol), and that in the absence of a membrane potential, the intrinsic chemical energy of this global conformational change is +9.4 kcal/mol (39 kJ/mol). Na+ and K+ in the extracellular solution have no measureable effect on E0, which suggests that unliganded gating occurs with only water occupying the transmitter binding sites. The results are discussed with regard to the energy changes in receptor activation and the competitive antagonism of ions in agonist binding.     

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.     

pH dependence of the acetylcholine receptor channel: a species variation

The effects of pH changes on the miniature endplate current (mepc) and on endplate current fluctuations (acetylcholine [ACh] noise) were examined at the neuromuscular junction in vitro in two species of frogs. In Rana pipiens the relationship between the decay time constant of the mepc (tau') and pH had a symmetrical bell shape; the value of tau' being largest at pH 7 and decreasing at more acid or more alkaline pH. In acid pH the mepc amplitude (A) decreased relative to its value at pH 7, and in alkaline pH A increased. In Rana ridibunda a narrower and asymmetric bell-shaped dependence of tau' on pH, having a maximum of pH 5.5, was found. The mepc amplitude was again reduced in acid pH but had a peak at pH 5.5. Also, its value at pH 9 was larger than at pH 7. These results were obtained with a number of different buffers and were not found to be sensitive to the nature of the buffer chosen. By performing ACh-noise analysis we found that in Rana pipiens at acid pH (5.5-5.0), the single channel conductance (gamma) and the single channel open time (tau) were significantly reduced relative to their value at pH 7. However, in Rana ridibunda at acid pH (5.4) gamma was unchanged and tau was markedly increased relative to their values at pH 7. The results can be explained quantitatively by electrostatic interaction between two fixed and titratable ionic groups and a mobile charge in the receptor molecule. The model fits the data for groups having pKs approximately 4.8 and approximately 9.8 for Rana pipiens and approximately 4.6 and approximately 6.3 for Rana ridibunda. The groups can be tentatively identified as amino acid residues; glutamic or aspartic and lysine or tyrosine for Rana pipiens; glutamic or aspartic and histidine for Rana ridibunda. The difference in the fitted values of the other model parameters for these two species can be attributed to differences in the spatial configuration of the charged groups.     

An intersubunit hydrogen bond in the nicotinic acetylcholine receptor that contributes to channel gating

The muscle nicotinic acetylcholine receptor is a large, allosteric, ligand-gated ion channel with the subunit composition alpha2betagammadelta. Although much is now known about the structure of the binding site, relatively little is understood about how the binding event is communicated to the channel gate, causing the pore to open. Here we identify a key hydrogen bond near the binding site that is involved in the gating pathway. Using mutant cycle analysis with the novel unnatural residue alpha-hydroxyserine, we find that the backbone N-H of alphaSer-191 in loop C makes a hydrogen bond to an anionic side chain of the complementary subunit upon agonist binding. However, the anionic partner is not the glutamate predicted by the crystal structures of the homologous acetylcholine-binding protein. Instead, the hydrogen-bonding partner is the extensively researched aspartate gammaAsp-174/deltaAsp-180, which had originally been identified as a key binding residue for cationic agonists.     

An energy-efficient gating mechanism in the acetylcholine receptor channel suggested by molecular and Brownian dynamics

Acetylcholine receptors mediate electrical signaling between nerve and muscle by opening and closing a transmembrane ion conductive pore. Molecular and Brownian dynamics simulations are used to shed light on the location and mechanism of the channel gate. Four separate 5 ns molecular dynamics simulations are carried out on the imaged structure of the channel, a hypothetical open structure with a slightly wider pore and a mutant structure in which a central ring of hydrophobic residues is replaced by polar groups. Water is found to partially evacuate the pore during molecular simulations of the imaged structure, whereas ions face a large energy barrier and do not conduct through the channel in Brownian dynamics simulations. The pore appears to be in a closed configuration despite containing an unobstructed pathway across the membrane as a series of hydrophobic residues in the center of the channel provide an unfavorable home to water and ions. When the channel is widened slightly, water floods into the channel and ions conduct at a rate comparable to the currents measured experimentally in open channels. The pore remains permeable to ions provided the extracellular end of the pore-lining helix is restrained near the putative open configuration to mimic the presence of the ligand binding domain. Replacing some of the hydrophobic residues with polar ones decreases the barrier for ion permeation but does not result in significant currents. The channel is posited to utilize an energy efficient gating mechanism in which only minor conformational changes of the hydrophobic region of the pore are required to create macroscopic changes in conductance.     

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.     

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 alpha and non-alpha subunits, or homo-pentameric, composed of alpha7 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 alpha1 and epsilon or alpha7 subunits. The replacement of M3 in alpha1 by epsilonM3 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 epsilon chimeric subunit. The duration of the open state decreases with the increase in the number of alpha1M3 segments, indicating additive contributions of M3 of all subunits to channel closing. Each alpha1M3 segment decreases the energy barrier of the closing process by approximately 0.8 kcal/mol. Partial chimeric subunits show that small stretches of the M3 segment contribute additively to the open duration. The replacement of alpha1 sequence by alpha7 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.     

Structure of the transition state of gating in the acetylcholine receptor channel pore: a phi-value analysis

The gating mechanism of the acetylcholine receptor channel (AChR) was investigated by using rate equilibrium linear free energy relationships (LFERs) to probe the transition state between the closed and open conformations. The properties of the transition state of gating in the second transmembrane segment (M2) of the delta subunit, one of the five homologous pore-lining segments, was measured on a residue-by-residue basis. Series of point mutations were engineered at individual positions of this domain, and the corresponding constructs were characterized electrophysiologically, at the single-channel level. Fully liganded AChR opening and closing rate constants were estimated, and Phi-values (which are a measure of the extent of the conformational change realized at the transition state) were calculated for each reaction series as the slope of the Br?nsted relationship (log rate constant versus log equilibrium constant). Our results indicate that, at the transition state of gating, the extracellular half of deltaM2 partly resembles the open state (Phi-values between 0.24 and 0.38) while the intracellular half completely resembles the closed state (Phi-values between -0.18 and 0.03), with a break point near the middle of the M2 segment. This suggests that during gating the two halves of deltaM2 move asynchronously, with the rearrangement of the extracellular portion preceding (following) that of the intracellular part of deltaM2 during opening (closing). This particular sequence of molecular events indicates that the gating conformational change, which starts at the extracellular acetylcholine-binding sites (when opening), does not propagate exclusively along the primary sequence of the protein. In addition, our data are consistent with the deltaM2 segment bending or swiveling around its central residues during gating. We also elaborate on unsettled aspects of the analysis such as the accuracy of two-point LFERs, the physical interpretation of fractional Phi-values, and the existence of single versus parallel transition states for the gating reaction.     

Naturally occurring mutations at the acetylcholine receptor binding site independently alter ACh binding and channel gating

By defining functional defects in a congenital myasthenic syndrome (CMS), we show that two mutant residues, located in a binding site region of the acetylcholine receptor (AChR) epsilon subunit, exert opposite effects on ACh binding and suppress channel gating. Single channel kinetic analysis reveals that the first mutation, epsilon N182Y, increases ACh affinity for receptors in the resting closed state, which promotes sequential occupancy of the binding sites and discloses rate constants for ACh occupancy of the nonmutant alphadelta site. Studies of the analogous mutation in the delta subunit, deltaN187Y, disclose rate constants for ACh occupancy of the nonmutant alpha epsilon site. The second CMS mutation, epsilon D175N, reduces ACh affinity for receptors in the resting closed state; occupancy of the mutant site still promotes gating because a large difference in affinity is maintained between closed and open states. epsilon D175N impairs overall gating, however, through an effect independent of ACh occupancy. When mapped on a structural model of the AChR binding site, epsilon N182Y localizes to the interface with the alpha subunit, and epsilon D175 to the entrance of the ACh binding cavity. Both epsilon N182Y and epsilon D175 show state specificity in affecting closed relative to desensitized state affinities, suggesting that the protein chain harboring epsilon N182 and epsilon D175 rearranges in the course of receptor desensitization. The overall results show that key residues at the ACh binding site differentially stabilize the agonist bound to closed, open and desensitized states, and provide a set point for gating of the channel.     

Acetylcholine-induced current in perfused rat myoballs

Spherical "myoballs" were grown under tissue culture conditions from striated muscle of neonatal rat thighs. The myoballs were examined electrophysiologically with a suction pipette which was used to pass current and perfuse internally. A microelectrode was used to record membrane potential. Experiments were performed with approximately symmetrical (intracellular and extracellular) sodium aspartate solutions. The resting potential, acetylcholine (ACh) reversal potential, and sodium channel reversal potential were all approximately 0 mV. ACh-induced currents were examined by use of both voltage jumps and voltage ramps in the presence of iontophoretically applied agonist. The voltage-jump relaxations had a single exponential time-course. The time constant, tau, was exponentially related to membrane potential, increasing e-fold for 81 mV hyperpolarization. The equilibrium current- voltage relationship was also approximately exponential, from -120 to +81 mV, increasing e-fold for 104 mV hyperpolarization. The data are consistent with a first-order gating process in which the channel opening rate constant is slightly voltage dependent. The instantaneous current-voltage relationship was sublinear in the hyperpolarizing direction. Several models are discussed which can account for the nonlinearity. Evidence is presented that the "selectivity filter" for the ACh channel is located near the intracellular membrane surface.     

Plasticity of acetylcholine receptor gating motions via rate-energy relationships

Like other protein conformational changes, ion channel gating requires the protein to achieve a high-energy transition-state structure. It is not known whether ion channel gating takes place on a broad energy landscape on which many alternative transition state structures are accessible, or on a narrow energy landscape where only a few transition-state structures are possible. To address this question, we measured how rate-equilibrium free energy relationships (REFERs) for di-liganded and unliganded acetylcholine receptor gating vary as a function of the gating equilibrium constant. A large slope for the REFER plot indicates an openlike transition state, whereas a small slope indicates a closedlike transition state. Due to this relationship between REFERs and transition-state structure, the sensitivity of the REFER slope to mutation-induced energetic perturbations allows estimation of the breadth of the energy landscape underlying gating. The relatively large sensitivity of di-liganded REFER slopes to energetic perturbations suggests that the motions underlying di-liganded gating take place on a broad, shallow energy landscape where many alternative transition-state structures are accessible.     

Statistical methods for model discrimination. Applications to gating kinetics and permeation of the acetylcholine receptor channel.

Methods are described for discrimination of models of the gating kinetics and permeation of single ionic channels. Both maximum likelihood and regression procedures are discussed. In simple situations, where models are nested, standard hypothesis tests can be used. More commonly, however, non-nested models are of interest, and several procedures are described for model discrimination in these cases, including Monte Carlo methods, which allow the comparison of models at significance levels of choice. As an illustration, the methods are applied to single-channel data from acetylcholine receptor channels.     

FKBP12 binding modulates ryanodine receptor channel gating

The ryanodine receptor (RyR1)/calcium release channel on the sarcoplasmic reticulum of skeletal muscle is comprised of four 565,000-dalton RyR1s, each of which binds one FK506 binding protein (FKBP12). RyR1 is required for excitation-contraction coupling in skeletal muscle. FKBP12, a cis-trans peptidyl-prolyl isomerase, is required for the normal gating of the RyR1 channel. In the absence of FKBP12, RyR1 channels exhibit increased gating frequency, suggesting that FKBP12 "stabilizes" the channel in the open and closed states. We now show that substitution of a Gly, Glu, or Ile for Val2461 in RyR1 prevents FKBP12 binding to RyR1, resulting in channels with increased gating frequency. In the case of the V2461I mutant RyR1, normal channel function can be restored by adding FKBP12.6, an isoform of FKBP12. These data identify Val2461 as a critical residue required for FKBP12 binding to RyR1 and demonstrate the functional role for FKBP12 in the RyR1 channel complex.     

Nonequilibrium gating and voltage dependence of the ClC-0 Cl- channel

The gating of ClC-0, the voltage-dependent Cl- channel from Torpedo electric organ, is strongly influenced by Cl- ions in the external solution. Raising external Cl- over the range 1-600 mM favors the fast- gating open state and disfavors the slow-gating inactivated state. Analysis of purified single ClC-0 channels reconstituted into planar lipid bilayers was used to identify the role of Cl- ions in the channel's fast voltage-dependent gating process. External, but not internal, Cl- had a major effect on the channel's opening rate constant. The closing rate was more sensitive to internal Cl- than to external Cl-. Both opening and closing rates varied with voltage. A model was derived that postulates (a) that in the channel's closed state, Cl- is accessible to a site located at the outer end of the conduction pore, where it binds in a voltage-independent fashion, (b) that this closed conformation can open, whether liganded by Cl- or not, in a weakly voltage-dependent fashion, (c) that the Cl(-)-liganded closed channel undergoes a conformational change to a different closed state, such that concomitant with this change, Cl- ion moves inward, conferring voltage-dependence to this step, and (d) that this new Cl(-)- liganded closed state opens with a very high rate. According to this picture, Cl- movement within the pre-open channel is the major source of voltage dependence, and charge movement intrinsic to the channel protein contributes very little to voltage-dependent gating of ClC-0. Moreover, since the Cl- activation site is probably located in the ion conduction pathway, the fast gating of ClC-0 is necessarily coupled to ion conduction, a nonequilibrium process.