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

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

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.

Agonist efficiency links binding and gating in a nicotinic receptor

Receptors signal by switching between resting (C) and active (O) shapes ('gating') under the influence of agonists. The receptor's maximum response depends on the difference in agonist binding energy, O minus C. In nicotinic receptors, efficiency (η) represents the fraction of agonist binding energy applied to a local rearrangement (an induced fit) that initiates gating. In this receptor, free energy changes in gating and binding can be interchanged by the conversion factor η. Efficiencies estimated from concentration-response curves (23 agonists, 53 mutations) sort into five discrete classes (%): 0.56 (17), 0.51(32), 0.45(13), 0.41(26), and 0.31(12), implying that there are 5 C versus O binding site structural pairs. Within each class efficacy and affinity are corelated linearly, but multiple classes hide this relationship. η unites agonist binding with receptor gating and calibrates one link in a chain of coupled domain rearrangements that comprises the allosteric transition of the protein.

Linear Free Energy Relationships for Enzymatic Reactions: Fresh Insight from a Venerable Probe

Linear free energy relationships (LFERs) for substituent effects on reactions that proceed through similar transition states provide insight into transition state structures. A classical approach to the analysis of LFERs showed that differences in the slopes of Brønsted correlations for addition of substituted alkyl alcohols to ring-substituted 1-phenylethyl carbocations and to the β-galactopyranosyl carbocation intermediate of reactions catalyzed by β-galactosidase provide evidence that the enzyme catalyst modifies the curvature of the energy surface at the saddle point for the transition state for nucleophile addition. We have worked to generalize the use of LFERs in the determination of enzyme mechanisms. The defining property of enzyme catalysts is their specificity for binding the transition state with a much higher affinity than the substrate. Triosephosphate isomerase (TIM), orotidine 5′-monophosphate decarboxylase (OMPDC), and glycerol 3-phosphate dehydrogenase (GPDH) show effective catalysis of reactions of phosphorylated substrates and strong phosphite dianion activation of reactions of phosphodianion truncated substrates, with rate constants k_cat/K_m (M^–1 s^–1) and k_cat/K_dK_HPi (M^–2 s^–1), respectively. Good linear logarithmic correlations, with a slope of 1.1, between these kinetic parameters determined for reactions catalyzed by five or more variant forms of each catalyst are observed, where the protein substitutions are mainly at side chains which function to stabilize the cage complex between the enzyme and substrate. This shows that the enzyme-catalyzed reactions of a whole substrate and substrate pieces proceed through transition states of similar structures. It provides support for the proposal that the dianion binding energy of whole phosphodianion substrates and of phosphite dianion is used to drive the conversion of these protein catalysts from flexible and entropically rich ground states to stiff and catalytically active Michaelis complexes that show the same activity toward catalysis of the reactions of whole and phosphodianion truncated substrates. There is a good linear correlation, with a slope of 0.73, between values of the dissociation constants log K_i for release of the transition state analog phosphoglycolate (PGA) trianion and log k_cat/K_m for isomerization of GAP for wild-type and variants of TIM. This correlation shows that the substituted amino acid side chains act to stabilize the complex between TIM and the PGA trianion and that ca. 70% of this stabilization is observed at the transition state for substrate deprotonation. The correlation provides evidence that these side chains function to enhance the basicity of the E165 side chain of TIM, which deprotonates the bound carbon acid substrate. There is a good linear correlation, with a slope of 0.74, between the values of ΔG^‡ and ΔG° determined by electron valence bond (EVB) calculations to model deprotonation of dihydroxyacetone phosphate (DHAP) in water and when bound to wild-type and variant forms of TIM to form the enediolate reaction intermediate. This correlation provides evidence that the stabilizing interactions of the transition state for TIM-catalyzed deprotonation of DHAP are optimized by placement of amino acid side chains in positions that provide for the maximum stabilization of the charged reaction intermediate, relative to the neutral substrate.

Chasing the open-state structure of pentameric ligand-gated ion channels

Members of the pentameric ligand-gated ion channel family have been crystallized in different conformations, including one in which the transmembrane pore is surprisingly wide. Gonzalez-Gutierrez et al. show that the open-channel conformation of animal members is more similar to the models with narrow pores.

Remarkable advances have been made toward the structural characterization of ion channels in the last two decades. However, the unambiguous assignment of well-defined functional states to the obtained structural models has proved challenging. In the case of the superfamily of nicotinic-receptor channels (also referred to as pentameric ligand-gated ion channels [pLGICs]), for example, two different types of model of the open-channel conformation have been proposed on the basis of structures solved to resolutions better than 4.0 Å. At the level of the transmembrane pore, the open-state models of the proton-gated pLGIC from Gloeobacter violaceus (GLIC) and the invertebrate glutamate-gated Cl^– channel (GluCl) are very similar to each other, but that of the glycine receptor (GlyR) is considerably wider. Indeed, the mean distances between the axis of ion permeation and the Cα atoms at the narrowest constriction of the pore (position −2′) differ by ∼2 Å in these two classes of model, a large difference when it comes to understanding the physicochemical bases of ion conduction and charge selectivity. Here, we take advantage of the extreme open-channel stabilizing effect of mutations at pore-facing position 9′. We find that the I9′A mutation slows down entry into desensitization of GLIC to the extent that macroscopic currents decay only slightly by the end of pH 4.5 solution applications to the extracellular side for several minutes. We crystallize (at pH 4.5) two variants of GLIC carrying this mutation and solve their structures to resolutions of 3.12 Å and 3.36 Å. Furthermore, we perform all-atom molecular dynamics simulations of ion permeation and picrotoxinin block, using the different open-channel structural models. On the basis of these results, we favor the notion that the open-channel structure of pLGICs from animals is much closer to that of the narrow models (of GLIC and GluCl) than it is to that of the GlyR.

Pathways and Barriers for Ion Translocation through the 5-HT3A Receptor Channel

Pentameric ligand gated ion channels (pLGICs) are ionotropic receptors that mediate fast intercellular communications at synaptic level and include either cation selective (e.g., nAChR and 5-HT₃) or anion selective (e.g., GlyR, GABA_A and GluCl) membrane channels. Among others, 5-HT₃ is one of the most studied members, since its first cloning back in 1991, and a large number of studies have successfully pinpointed protein residues critical for its activation and channel gating. In addition, 5-HT₃ is also the target of a few pharmacological treatments due to the demonstrated benefits of its modulation in clinical trials. Nonetheless, a detailed molecular analysis of important protein features, such as the origin of its ion selectivity and the rather low conductance as compared to other channel homologues, has been unfeasible until the recent crystallization of the mouse 5-HT₃A receptor. Here, we present extended molecular dynamics simulations and free energy calculations of the whole 5-HT₃A protein with the aim of better understanding its ion transport properties, such as the pathways for ion permeation into the receptor body and the complex nature of the selectivity filter. Our investigation unravels previously unpredicted structural features of the 5-HT₃A receptor, such as the existence of alternative intersubunit pathways for ion translocation at the interface between the extracellular and the transmembrane domains, in addition to the one along the channel main axis. Moreover, our study offers a molecular interpretation of the role played by an arginine triplet located in the intracellular domain on determining the characteristic low conductance of the 5-HT₃A receptor, as evidenced in previous experiments. In view of these results, possible implications on other members of the superfamily are suggested.

The atypical cation-conduction and gating properties of ELIC underscore the marked functional versatility of the pentameric ligand-gated ion-channel fold

The superfamily of pentameric ligand-gated ion channels (pLGICs) is unique among ionotropic receptors in that the same overall structure has evolved to generate multiple members with different combinations of agonist specificities and permeant-ion charge selectivities. However, aside from these differences, pLGICs have been typically regarded as having several invariant functional properties. These include pore blockade by extracellular quaternary-ammonium cations in the micromolar-to-millimolar concentration range (in the case of the cation-selective members), and a gain-of-function phenotype, which manifests as a slower deactivation time course, as a result of mutations that reduce the hydrophobicity of the transmembrane pore lining. Here, we tested this notion on three distantly related cation-selective members of the pLGIC superfamily: the mouse muscle nicotinic acetylcholine receptor (nAChR), and the bacterial GLIC and ELIC channels. Remarkably, we found that, whereas low millimolar concentrations of TMA(+) and TEA(+) block the nAChR and GLIC, neither of these two quaternary-ammonium cations blocks ELIC at such concentrations; instead, both carry measurable inward currents when present as the only cations on the extracellular side. Also, we found that, whereas lidocaine binding speeds up the current-decay time courses of the nAChR and GLIC in the presence of saturating concentrations of agonists, the binding of lidocaine to ELIC slows this time course down. Furthermore, whereas mutations that reduce the hydrophobicity of the side chains at position 9' of the M2 α-helices greatly slowed the deactivation time course of the nAChR and GLIC, these mutations had little effect--or even sped up deactivation--when engineered in ELIC. Our data indicate that caution should be exercised when generalizing results obtained with ELIC to the rest of the pLGICs, but more intriguingly, they hint at the possibility that ELIC is a representative of a novel branch of the superfamily with markedly divergent pore properties despite a well-conserved three-dimensional architecture.

Functional anatomy of an allosteric protein

Synaptic receptors are allosteric proteins that switch on and off to regulate cell signalling. Here, we use single-channel electrophysiology to measure and map energy changes in the gating conformational change of a nicotinic acetylcholine receptor. Two separated regions in the α-subunits--the transmitter-binding sites and αM2-αM3 linkers in the membrane domain--have the highest ϕ-values (change conformation the earliest), followed by the extracellular domain, most of the membrane domain and the gate. Large gating-energy changes occur at the transmitter-binding sites, α-subunit interfaces, the αM1 helix and the gate. We hypothesize that rearrangements of the linkers trigger the global allosteric transition, and that the hydrophobic gate unlocks in three steps. The mostly local character of side-chain energy changes and the similarly high ϕ-values of separated domains, both with and without ligands, suggest that gating is not strictly a mechanical process initiated by the affinity change for the agonist.

Hydrophobic interactions between the voltage sensor and pore mediate inactivation in Kv11.1 channels

Kv11.1 channels are critical for the maintenance of a normal heart rhythm. The flow of potassium ions through these channels is controlled by two voltage-regulated gates, termed "activation" and "inactivation," located at opposite ends of the pore. Crucially in Kv11.1 channels, inactivation gating occurs much more rapidly, and over a distinct range of voltages, compared with activation gating. Although it is clear that the fourth transmembrane segments (S4), within each subunit of the tetrameric channel, are important for controlling the opening and closing of the activation gate, their role during inactivation gating is much less clear. Here, we use rate equilibrium free energy relationship (REFER) analysis to probe the contribution of the S4 "voltage-sensor" helix during inactivation of Kv11.1 channels. Contrary to the important role that charged residues play during activation gating, it is the hydrophobic residues (Leu529, Leu530, Leu532, and Val535) that are the key molecular determinants of inactivation gating. Within the context of an interconnected multi-domain model of Kv11.1 inactivation gating, our REFER analysis indicates that the S4 helix and the S4-S5 linker undergo a conformational rearrangement shortly after that of the S5 helix and S5P linker, but before the S6 helix. Combining REFER analysis with double mutant cycle analysis, we provide evidence for a hydrophobic interaction between residues on the S4 and S5 helices. Based on a Kv11.1 channel homology model, we propose that this hydrophobic interaction forms the basis of an intersubunit coupling between the voltage sensor and pore domain that is an important mediator of inactivation gating.

Gating of the proton-gated ion channel from Gloeobacter violaceus at pH 4 as revealed by X-ray crystallography

Cryoelectron microscopy and X-ray crystallography have recently been used to generate structural models that likely represent the unliganded closed-channel conformation and the fully liganded open-channel conformation of different members of the nicotinic-receptor superfamily. To characterize the structure of the closed-channel conformation in its liganded state, we identified a number of positions in the loop between transmembrane segments 2 (M2) and 3 (M3) of a proton-gated ortholog from the bacterium Gloeobacter violaceus (GLIC) where mutations to alanine reduce the liganded-gating equilibrium constant, and solved the crystal structures of two such mutants (T25'A and Y27'A) at pH ~4.0. At the level of backbone atoms, the liganded closed-channel model presented here differs from the liganded open-channel structure of GLIC in the pre-M1 linker, the M3-M4 loop, and much more prominently, in the extracellular half of the pore lining, where the more pronounced tilt of the closed-channel M2 α-helices toward the pore's long axis narrows the permeation pathway. On the other hand, no differences between the liganded closed-channel and open-channel models could be detected at the level of the extracellular domain, where conformational changes are expected to underlie the low-to-high proton-affinity switch that drives gating of proton-bound channels. Thus, the liganded closed-channel model is nearly indistinguishable from the recently described "locally closed" structure. However, because cross-linking strategies (which could have stabilized unstable conformations) and mutations involving ionizable side chains (which could have affected proton-gated channel activation) were purposely avoided, we favor the notion that this structure represents one of the end states of liganded gating rather than an unstable intermediate.

Identification of a possible secondary picrotoxin-binding site on the GABA(A) receptor

The type A GABA receptors (GABARs) are ligand-gated ion channels (LGICs) found in the brain and are the major inhibitory neurotransmitter receptors. Upon binding of an agonist, the GABAR opens and increases the intraneuronal concentration of chloride ions, thus hyperpolarizing the cell and inhibiting the transmission of the nerve action potential. GABARs also contain many other modulatory binding pockets that differ from the agonist-binding site. The composition of the GABAR subunits can alter the properties of these modulatory sites. Picrotoxin is a noncompetitive antagonist for LGICs, and by inhibiting GABAR, picrotoxin can cause overstimulation and induce convulsions. We use addition of picrotoxin to probe the characteristics and possible mechanism of an additional modulatory pocket located at the interface between the ligand-binding domain and the transmembrane domain of the GABAR. Picrotoxin is widely regarded as a pore-blocking agent that acts at the cytoplasmic end of the channel. However, there are also data to suggest that there may be an additional, secondary binding site for picrotoxin. Through homology modeling, molecular docking, and molecular dynamics simulations, we show that binding of picrotoxin to this interface pocket correlates with these data, and negative modulation occurs at the pocket via a kinking of the pore-lining helices into a more closed orientation.

Pore helices play a dynamic role as integrators of domain motion during Kv11.1 channel inactivation gating

Proteins that form ion-selective pores in the membrane of cells are integral to many rapid signaling processes, including regulating the rhythm of the heartbeat. In potassium channels, the selectivity filter is critical for both endowing an exquisite selectivity for potassium ions, as well as for controlling the flow of ions through the pore. Subtle rearrangements in the complex hydrogen-bond network that link the selectivity filter to the surrounding pore helices differentiate conducting (open) from nonconducting (inactivated) conformations of the channel. Recent studies suggest that beyond the selectivity filter, inactivation involves widespread rearrangements of the channel protein. Here, we use rate equilibrium free energy relationship analysis to probe the structural changes that occur during selectivity filter gating in Kv11.1 channels, at near atomic resolution. We show that the pore helix plays a crucial dynamic role as a bidirectional interface during selectivity filter gating. We also define the molecular bases of the energetic coupling between the pore helix and outer helix of the pore domain that occurs early in the transition from open to inactivated states, as well as the coupling between the pore helix and inner helix late in the transition. Our data demonstrate that the pore helices are more than just static structural elements supporting the integrity of the selectivity filter; instead they play a crucial dynamic role during selectivity filter gating.

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.

Conformational transitions underlying pore opening and desensitization in membrane-embedded Gloeobacter violaceus ligand-gated ion channel (GLIC)

Direct structural insight into the mechanisms underlying activation and desensitization remain unavailable for the pentameric ligand-gated channel family. Here, we report the structural rearrangements underlying gating transitions in membrane-embedded GLIC, a prokaryotic homologue, using site-directed spin labeling and electron paramagnetic resonance (EPR) spectroscopy. We particularly probed the conformation of pore-lining second transmembrane segment (M2) under conditions that favor the closed and the ligand-bound desensitized states. The spin label mobility, intersubunit spin-spin proximity, and the solvent-accessibility parameters in the two states clearly delineate the underlying protein motions within M2. Our results show that during activation the extracellular hydrophobic region undergoes major changes involving an outward translational movement, away from the pore axis, leading to an increase in the pore diameter, whereas the lower end of M2 remains relatively immobile. Most notably, during desensitization, the intervening polar residues in the middle of M2 move closer to form a solvent-occluded barrier and thereby reveal the location of a distinct desensitization gate. In comparison with the crystal structure of GLIC, the structural dynamics of the channel in a membrane environment suggest a more loosely packed conformation with water-accessible intrasubunit vestibules penetrating from the extracellular end all the way to the middle of M2 in the closed state. These regions have been implicated to play a major role in alcohol and drug modulation. Overall, these findings represent a key step toward understanding the fundamentals of gating mechanisms in this class of channels.

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 (E₀) to generate AChRs that are active constitutively. We measured E₀ for two different sets of mutant combinations and by extrapolation estimated E₀ for wild-type AChRs. The estimates were 7.6 and 7.8 × 10⁻⁷ 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⁻⁷, from monoliganded gating). E₀ 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 E₀, 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.

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.

Mutations that stabilize the open state of the Erwinia chrisanthemi ligand-gated ion channel fail to change the conformation of the pore domain in crystals

The determination of structural models of the various stable states of an ion channel is a key step toward the characterization of its conformational dynamics. In the case of nicotinic-type receptors, different structures have been solved but, thus far, these different models have been obtained from different members of the superfamily. In the case of the bacterial member ELIC, a cysteamine-gated channel from Erwinia chrisanthemi, a structural model of the protein in the absence of activating ligand (and thus, conceivably corresponding to the closed state of this channel) has been previously generated. In this article, electrophysiological characterization of ELIC mutants allowed us to identify pore mutations that slow down the time course of desensitization to the extent that the channel seems not to desensitize at all for the duration of the agonist applications (>20 min). Thus, it seems reasonable to conclude that the probability of ELIC occupying the closed state is much lower for the ligand-bound mutants than for the unliganded wild-type channel. To gain insight into the conformation adopted by ELIC under these conditions, we solved the crystal structures of two of these mutants in the presence of a concentration of cysteamine that elicits an intracluster open probability of >0.9. Curiously, the obtained structural models turned out to be nearly indistinguishable from the model of the wild-type channel in the absence of bound agonist. Overall, our findings bring to light the limited power of functional studies in intact membranes when it comes to inferring the functional state of a channel in a crystal, at least in the case of the nicotinic-receptor superfamily.

Effects of lipid-analog detergent solubilization on the functionality and lipidic cubic phase mobility of the Torpedo californica nicotinic acetylcholine receptor

Over the past three decades, the Torpedo californica nicotinic acetylcholine receptor (nAChR) has been one of the most extensively studied membrane protein systems. However, the effects of detergent solubilization on nAChR stability and function are poorly understood. The use of lipid-analog detergents for nAChR solubilization has been shown to preserve receptor stability and functionality. The present study used lipid-analog detergents from phospholipid-analog and cholesterol-analog detergent families for solubilization and affinity purification of the receptor and probed nAChR ion channel function using planar lipid bilayers (PLBs) and stability using analytical size exclusion chromatography (A-SEC) in the detergent-solubilized state. We also examined receptor mobility on the lipidic cubic phase (LCP) by measuring the nAChR mobile fraction and diffusion coefficient through fluorescence recovery after photobleaching (FRAP) experiments using lipid-analog and non-lipid-analog detergents. Our results show that it is possible to isolate stable and functional nAChRs using lipid-analog detergents, with characteristic ion channel currents in PLBs and minimal aggregation as observed in A-SEC. Furthermore, fractional mobility and diffusion coefficient values observed in FRAP experiments were similar to the values observed for these parameters in the recently LCP-crystallized β(2)-adrenergic receptor. The overall results show that phospholipid-analog detergents with 16 carbon acyl-chains support nAChR stability, functionality and LCP mobility.

Charge substitution for a deep-pore residue reveals structural dynamics during BK channel gating

The pore-lining amino acids of ion channel proteins reside on the interface between a polar (the pore) and a nonpolar environment (the rest of the protein). The structural dynamics of this region, which physically controls ionic flow, are essential components of channel gating. Using large-conductance, Ca²⁺-dependent K⁺ (BK) channels, we devised a systematic charge–substitution method to probe conformational changes in the pore region during channel gating. We identified a deep-pore residue (314 in hSlo1) as a marker of structural dynamics. We manipulated the charge states of this residue by substituting amino acids with different valence and pKa, and by adjusting intracellular pH. We found that the charged states of the 314 residues stabilized an open state of the BK channel. With models based on known structures of related channels, we postulate a dynamic rearrangement of the deep-pore region during BK channel opening/closing, which involves a change of the degree of pore exposure for 314.

Mapping heat exchange in an allosteric protein

Nicotinic acetylcholine receptors (AChRs) are synaptic ion channels that spontaneously isomerize (i.e., gate) between resting and active conformations. We used single-molecule electrophysiology to measure the temperature dependencies of mouse neuromuscular AChR gating rate and equilibrium constants. From these we estimated free energy, enthalpy, and entropy changes caused by mutations of amino acids located between the transmitter binding sites and the middle of the membrane domain. The range of equilibrium enthalpy change (13.4 kcal/mol) was larger than for free energy change (5.5 kcal/mol at 25°C). For two residues, the slope of the rate-equilibrium free energy relationship (Φ) was approximately constant with temperature. Mutant cycle analysis showed that both free energies and enthalpies are additive for energetically independent mutations. We hypothesize that changes in energy associated with changes in structure mainly occur close to the site of the mutation, and, hence, that it is possible to make a residue-by-residue map of heat exchange in the AChR gating isomerization. The structural correlates of enthalpy changes are discussed for 12 different mutations in the protein.

Mapping the sequence of conformational changes underlying selectivity filter gating in the K(v)11.1 potassium channel

The potassium channel selectivity filter both discriminates between K(+) and sodium ions and contributes to gating of ion flow. Static structures of conducting (open) and nonconducting (inactivated) conformations of this filter are known; however, the sequence of protein rearrangements that connect these two states is not. We show that closure of the selectivity filter gate in the human K(v)11.1 K(+) channel (also known as hERG, for ether-a-go-go-related gene), a key regulator of the rhythm of the heartbeat, is initiated by K(+) exit, followed in sequence by conformational rearrangements of the pore domain outer helix, extracellular turret region, voltage sensor domain, intracellular domains and pore domain inner helix. In contrast to the simple wave-like sequence of events proposed for opening of ligand-gated ion channels, a complex spatial and temporal sequence of widespread domain motions connect the open and inactivated states of the K(v)11.1 K(+) channel.

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.

Probing the transition state of the allosteric pathway of the Shaker Kv channel pore by linear free-energy relations

Long-range coupling between distant functional elements of proteins may rely on allosteric communication trajectories lying along the protein structure, as described in the case of the Shaker voltage-activated potassium (Kv) channel model allosteric system. Communication between the distant Kv channel activation and slow inactivation pore gates was suggested to be mediated by a network of local pairwise and higher-order interactions among the functionally unique residues that constitute the allosteric trajectory. The mechanism by which conformational changes propagate along the Kv channel allosteric trajectory to achieve pore opening, however, remains unclear. Such conformational changes may propagate in either a concerted or a sequential manner during the reaction coordinate of channel opening. Residue-level structural information on the transition state of channel gating is required to discriminate between these possibilities. Here, we combine patch-clamp electrophysiology recordings of Kv channel gating and analysis using linear free-energy relations, focusing on a select set of residues spanning the allosteric trajectory of the Kv channel pore. We show that all allosteric trajectory residues tested exhibit an open-like conformation in the transition state of channel opening, implying that coupling interactions occur along the trajectory break in a concerted manner upon moving from the closed to the open state. Energetic coupling between the Kv channel gates thus occurs in a concerted fashion in both the spatial and the temporal dimensions, strengthening the notion that such trajectories correspond to pathways of mechanical deformation along which conformational changes propagate.

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.

The gating isomerization of neuromuscular acetylcholine receptors

Acetylcholine receptor-channels are allosteric proteins that isomerize ('gate') between conformations that have a low vs. high affinity for the transmitter and conductance for ions. In order to comprehend the mechanism by which the affinity and conductance changes are linked it is of value to know the magnitude, timing and distribution of energy flowing through the system. Knowing both the di- and unliganded gating equilibrium constants (E(2) and E(0)) is a foundation for understanding the AChR gating mechanism and for engineering both the ligand and the protein to operate in predictable ways. In adult mouse neuromuscular receptors activated by acetylcholine, E(2) = 28 and E(0) approximately 6.5 x 10(7). At each (equivalent) transmitter binding site acetylcholine provides approximately 5.2 kcal mol(1) to motivate the isomerization. The partial agonist choline provides approximately 3.3 kcal mol(1). The relative time of a residue's gating energy change is revealed by the slope of its rate-equilibrium constant relationship. A map of this parameter suggests that energy propagates as a conformational cascade between the transmitter binding sites and the gate region. Although gating energy changes are widespread throughout the protein, some residues are particularly sensitive to perturbations. Several specific proposals for the structural events that comprise the gating conformational cascade are discussed.

Application of rate-equilibrium free energy relationship analysis to nonequilibrium ion channel gating mechanisms

Rate-equilibrium free energy relationship (REFER) analysis provides information on transition-state structures and has been applied to reveal the temporal sequence in which the different regions of an ion channel protein move during a closed-open conformational transition. To date, the theory used to interpret REFER relationships has been developed only for equilibrium mechanisms. Gating of most ion channels is an equilibrium process, but recently several ion channels have been identified to have retained nonequilibrium traits in their gating cycles, inherited from transporter-like ancestors. So far it has not been examined to what extent REFER analysis is applicable to such systems. By deriving the REFER relationships for a simple nonequilibrium mechanism, this paper addresses whether an equilibrium mechanism can be distinguished from a nonequilibrium one by the characteristics of their REFER plots, and whether information on the transition-state structures can be obtained from REFER plots for gating mechanisms that are known to be nonequilibrium cycles. The results show that REFER plots do not carry information on the equilibrium nature of the underlying gating mechanism. Both equilibrium and nonequilibrium mechanisms can result in linear or nonlinear REFER plots, and complementarity of REFER slopes for opening and closing transitions is a trivial feature true for any mechanism. Additionally, REFER analysis provides limited information about the transition-state structures for gating schemes that are known to be nonequilibrium cycles.

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.

Recent advances in understanding the structure of nicotinic acetylcholine receptors

Nicotinic acetylcholine receptors (nAChRs), members of the Cys-loop ligand-gated ion channels (LGICs) superfamily, are involved in signal transduction upon binding of the neurotransmitter acetylcholine or exogenous ligands, such as nicotine. nAChRs are pentameric assemblies of homologous subunits surrounding a central pore that gates cation flux, and are expressed at the neuromuscular junction and in the nervous system and several nonneuronal cell types. The 17 known nAChR subunits assemble into a variety of pharmacologically distinct receptor subtypes. nAChRs are implicated in a range of physiological functions and pathophysiological conditions related to muscle contraction, learning and memory, reward, motor control, arousal, and analgesia, and therefore present an important target for drug research. Such studies would be greatly facilitated by knowledge of the high-resolution structure of the nAChR. Although this information is far from complete, important progress has been made mainly based on electron microscopy studies of Torpedo nAChR and the high-resolution X-ray crystal structures of the homologous molluscan acetylcholine-binding proteins, the extracellular domain of the mouse nAChR alpha1 subunit, and two prokaryotic pentameric LGICs. Here, we review some of the latest advances in our understanding of nAChR structure and gating.

Mechanistic contributions of residues in the M1 transmembrane domain of the nicotinic receptor to channel gating

The nicotinic receptor (AChR) is a pentamer of homologous subunits with an alpha(2)betaepsilondelta composition in adult muscle. Each subunit contains four transmembrane domains (M1-M4). Position 15' of the M1 domain is phenylalanine in alpha subunits while it is isoleucine in non-alpha subunits. Given this peculiar conservation pattern, we studied its contribution to muscle AChR activation by combining mutagenesis with single-channel kinetic analysis. AChRs containing the mutant alpha subunit (alphaF15'I) as well as those containing the reverse mutations in the non-alpha subunits (betaI15'F, deltaI15'F, and epsilonI15'F) show prolonged lifetimes of the diliganded open channel resulting from a slower closing rate with respect to wild-type AChRs. The kinetic changes are not equivalent among subunits, the beta subunit, being the one that produces the most significant stabilization of the open state. Kinetic analysis of betaI15'F of AChR channels activated by the low-efficacious agonist choline revealed a 10-fold decrease in the closing rate, a 2.5-fold increase in the opening rate, a 28-fold increase in the gating equilibrium constant in the diliganded receptor, and a significant increase opening in the absence of agonist. Mutations at betaI15' showed that the structural bases of its contribution to gating is complex. Rate-equilibrium linear free-energy relationships suggest an approximately 70% closed-state-like environment for the beta15' position at the transition state of gating. The overall results identify position 15' as a subunit-selective determinant of channel gating and add new experimental evidence that gives support to the involvement of the M1 domain in the operation of the channel gating apparatus.

Flip and flop: a molecular determinant for AMPA receptor channel opening

Alternative splicing in the extracellular ligand binding domain of the AMPA receptors generates two variants, i.e., flip and flop. The flop variant of the GluR2 AMPA receptor is known to desensitize faster than the flip counterpart, whereas the GluR1 flip and flop variants exhibit the same rate of desensitization. However, whether the alternative splicing affects the channel opening kinetic properties of these receptors is unknown. Using a laser-pulse photolysis technique, we have characterized the channel opening kinetic mechanism for the flip and flop channels of GluR1 and GluR2, respectively. We find that the flop variant of GluR2 opens the channel, following the binding of glutamate, with the same rate as the flip channel, but closes its channel more rapidly. The difference in the kinetic properties between the two receptor isoforms can be described by a model we proposed previously in which the channel closing rate, a measure of the stability of the open channel state, controls an apparent tendency of the channel to desensitize, most likely, through the closed channel state. Specifically, the flop sequence of GluR2 promotes the channel to close more rapidly and consequently to desensitize with a faster rate than the flip sequence. For GluR1, the alternative splicing does not seem to affect the channel opening kinetics, since the flip and flop variants of GluR1 have the same channel opening rate, and the same channel closing rate. As expected and indeed observed, the flop variant desensitizes with the same rate as the flip variant does. On the basis of these results, we hypothesize that the flip/flop sequence cassette of AMPA receptors, in a sequence-dependent manner, regulates the rate of the channel closing process, in the microsecond time domain, through which it further regulates the channel desensitization in the millisecond time region.

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.

In silico models for the human alpha4beta2 nicotinic acetylcholine receptor

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is one of the most widely expressed nAChR subtypes in the brain. Its subunits have high sequence identity (54 and 46% for alpha4 and beta2, respectively) with alpha and beta subunits in Torpedo nAChR. Using the known structure of the Torpedo nAChR as a template, the closed-channel structure of the alpha4beta2 nAChR was constructed through homology modeling. Normal-mode analysis was performed on this closed structure and the resulting lowest frequency mode was applied to it for a "twist-to-open" motion, which increased the minimum pore radius from 2.7 to 3.4 A and generated an open-channel model. Nicotine could bind to the predicted agonist binding sites in the open-channel model but not in the closed one. Both models were subsequently equilibrated in a ternary lipid mixture via extensive molecular dynamics (MD) simulations. Over the course of 11 ns MD simulations, the open channel remained open with filled water, but the closed channel showed a much lower water density at its hydrophobic gate comprised of residues alpha4-V259 and alpha4-L263 and their homologous residues in the beta2 subunits. Brownian dynamics simulations of Na+ permeation through the open channel demonstrated a current-voltage relationship that was consistent with experimental data on the conducting state of alpha4beta2 nAChR. Besides establishment of the well-equilibrated closed- and open-channel alpha4beta2 structural models, the MD simulations on these models provided valuable insights into critical factors that potentially modulate channel gating. Rotation and tilting of TM2 helices led to changes in orientations of pore-lining residue side chains. Without concerted movement, the reorientation of one or two hydrophobic side chains could be enough for channel opening. The closed- and open-channel structures exhibited distinct patterns of electrostatic interactions at the interface of extracellular and transmembrane domains that might regulate the signal propagation of agonist binding to channel opening. A potential prominent role of the beta2 subunit in channel gating was also elucidated in the study.

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.

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.

The protein folding problem

The "protein folding problem" consists of three closely related puzzles: (a) What is the folding code? (b) What is the folding mechanism? (c) Can we predict the native structure of a protein from its amino acid sequence? Once regarded as a grand challenge, protein folding has seen great progress in recent years. Now, foldable proteins and nonbiological polymers are being designed routinely and moving toward successful applications. The structures of small proteins are now often well predicted by computer methods. And, there is now a testable explanation for how a protein can fold so quickly: A protein solves its large global optimization problem as a series of smaller local optimization problems, growing and assembling the native structure from peptide fragments, local structures first.

Conformational variability of the glycine receptor M2 domain in response to activation by different agonists

Models describing the structural changes mediating Cys loop receptor activation generally give little attention to the possibility that different agonists may promote activation via distinct M2 pore-lining domain structural rearrangements. We investigated this question by comparing the effects of different ligands on the conformation of the external portion of the homomeric alpha1 glycine receptor M2 domain. Conformational flexibility was assessed by tethering a rhodamine fluorophore to cysteines introduced at the 19' or 22' positions and monitoring fluorescence and current changes during channel activation. During glycine activation, fluorescence of the label attached to R19'C increased by approximately 20%, and the emission peak shifted to lower wavelengths, consistent with a more hydrophobic fluorophore environment. In contrast, ivermectin activated the receptors without producing a fluorescence change. Although taurine and beta-alanine were weak partial agonists at the alpha1R19'C glycine receptor, they induced large fluorescence changes. Propofol, which drastically enhanced these currents, did not induce a glycine-like blue shift in the spectral emission peak. The inhibitors strychnine and picrotoxin elicited fluorescence and current changes as expected for a competitive antagonist and an open channel blocker, respectively. Glycine and taurine (or beta-alanine) also produced an increase and a decrease, respectively, in the fluorescence of a label attached to the nearby L22'C residue. Thus, results from two separate labeled residues support the conclusion that the glycine receptor M2 domain responds with distinct conformational changes to activation by different agonists.

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

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

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.

An intersubunit trigger of channel gating in the muscle nicotinic receptor

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

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.

Channel opening by anesthetics and GABA induces similar changes in the GABAA receptor M2 segment

For many general anesthetics, their molecular basis of action involves interactions with GABA(A) receptors. Anesthetics produce concentration-dependent effects on GABA(A) receptors. Low concentrations potentiate submaximal GABA-induced currents. Higher concentrations directly activate the receptors. Functional effects of anesthetics have been characterized, but little is known about the conformational changes they induce. We probed anesthetic-induced conformational changes in the M2 membrane-spanning, channel-lining segment using disulfide trapping between engineered cysteines. Previously, we showed that oxidation by copper phenanthroline in the presence of GABA of the M2 6' cysteine mutants, alpha(1)T261Cbeta(1)T256C and alpha(1)beta(1)T256C resulted in formation of an intersubunit disulfide bond between the adjacent beta-subunits that significantly increased the channels' spontaneous open probability. Oxidation in GABA's absence had no effect. We examined the effect on alpha(1)T261Cbeta(1)T256C and on alpha(1)beta(1)T256C of oxidation by copper phenanthroline in the presence of potentiating and directly activating concentrations of the general anesthetics propofol, pentobarbital, and isoflurane. Oxidation in the presence of potentiating concentration of anesthetics had little effect. Oxidation in the presence of directly activating anesthetic concentrations significantly increased the channels' spontaneous open probability. We infer that activation by anesthetics and GABA induces a similar conformational change at the M2 segment 6' position that is related to channel opening.

Tryptophan-scanning mutagenesis in the alphaM3 transmembrane domain of the muscle-type acetylcholine receptor. A spring model revealed

Membrane proteins constitute a large fraction of all proteins, yet very little is known about their structure and conformational transitions. A fundamental question that remains obscure is how protein domains that are in direct contact with the membrane lipids move during the conformational change of the membrane protein. Important structural and functional information of several lipid-exposed transmembrane domains of the acetylcholine receptor (AChR) and other ion channel membrane proteins have been provided by the tryptophan-scanning mutagenesis. Here, we use the tryptophan-scanning mutagenesis to monitor the conformational change of the alphaM3 domain of the muscle-type AChR. The perturbation produced by the systematic tryptophan substitution along the alphaM3 domain were characterized through two-electrode voltage clamp and 125I-labeled alpha-bungarotoxin binding. The periodicity profiles of the changes in AChR expression (closed state) and ACh EC50 (open-channel state) disclose two different helical structures; a thinner-elongated helix for the closed state and a thicker-shrunken helix for the open-channel state. The existence of two different helical structures suggest that the conformational transition of the alphaM3 domain between both states resembles a spring motion and reveals that the lipid-AChR interface plays a key role in the propagation of the conformational wave evoked by agonist binding. In addition, the present study also provides evidence about functional and structural differences between the alphaM3 domains of the Torpedo and muscle-type receptors AChR.

Agonist-activated ion channels

This paper looks at ion channels as an example of the pharmacologist's stock in trade, the action of an agonist on a receptor to produce a response. Looked at in this way, ion channels have been helpful because they are still the only system which is simple enough for quantitative investigation of transduction mechanisms. A short history is given of attempts to elucidate what happens between the time when agonist first binds, and the time when the channel opens.

Spontaneous mobility of GABAA receptor M2 extracellular half relative to noncompetitive antagonist action

The gamma-aminobutyric acid type A receptor beta(3) homopentamer is spontaneously open and highly sensitive to many noncompetitive antagonists(NCAs) and Zn(2+). Our earlier study of the M2 cytoplasmic half (-1' to 10') established a model in which NCAs bind at pore-lining residues Ala(2)', Thr(6)', and Leu(9)'. To further define transmembrane 2 (M2) structure relative to NCA action, we extended the Cys scanning to the extra cellular half of the beta(3) homopentamer (11' to 20'). Spontaneous disulfides formed with T13'C, L18'C, and E20'C from M2/M2 cross-linking and with I14'C (weak), H17'C, and R19'Con bridging M2/M3 intersubunits, based on single (M2 Cys only) and dual (M2 Cys plus M3 C289S) mutations. Induced disulfides also formed with T16'C, but there were few or none with M11'C, T12'C, and N15'C. These findings show conformational flexibility/mobility in the M2 extracellular half 17' to 20' region interpreted as a deformed beta-like conformation in the open channel. The NCA radioligands used were [(3)H]1-(4-ethynylphenyl)-4-n-propyl-2,6,7-trioxabicyclo[2.2.2]octane ([(3)H]EBOB) and [(3)H]3,3-bis-trifluoromethylbicyclo[2.2.1]heptane-2,2-dicarbonitrile with essentially the same results. NCA binding was disrupted by individual Cys substitutions at 13',14',16',17', and 19'. The inactivity of T13'C/T13'S may have been due to disturbance of the channel gate; I14'S and T16'S showed much better binding activity than their Cys counterparts, and the low activities of H17'C and R19'C were reversed by dithiothreitol. Zn(2+) potency for inhibition of [(3)H]EBOB binding was lowered 346-fold by the mutation H17'A. We propose that NCAs enter their binding site both directly, through the channel pore, and indirectly, through the water cavity of adjacent subunits.

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.

Recent advances in Cys-loop receptor structure and function

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

Do proteins with similar folds have similar transition state structures? A diffuse transition state of the 169 residue apoflavodoxin

Apoflavodoxin from Anabaena PCC 7119 is a 169 residue globular protein of known structure and energetics. Here, we present a comprehensive Phi-value analysis to characterize the structure of its transition state. A total of 34 non-disruptive mutations are made throughout the structure and a range of Phi-values from zero to one are observed. In addition, a small set of eight aliphatic small-to-large mutations have been introduced in the hydrophobic core of the protein and they have been analyzed to investigate the feasibility of stabilizing the unfolding transition state by creating new non-native interactions. We find that the transition state of apoflavodoxin (so far the largest protein subjected to Phi-analysis) is diffuse and that it can be stabilized by unspecific hydrophobic interactions that can speed up the folding reaction. The data gathered on the apoflavodoxin transition state are compared with results from experimental studies in other proteins to revisit the relationship between the native state topology and transition state structure.