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

The effect of unliganded gating on agonist response in nicotinic receptors

Understanding the agonist concentration-response curve (CRC) is the cornerstone in pharmacology. While CRC parameters, agonist potency (EC50) and efficacy (maximum response, Imax) are well-studied, the role of unliganded gating (minimum response, Imin) on CRC is often overlooked. This study explores the effect of unliganded gating on agonist response in muscle-type acetylcholine (ACh) receptors, focusing on the underexplored role of Imin in modulating EC50 and Imax. Three Gain-of-Function (GOF) mutations that increase, and two Loss-of-Function (LOF) mutations that decrease the unliganded gating equilibrium constant (L0) were studied using automated patch-clamp electrophysiology. GOF mutations enhanced agonist potency, whereas LOF mutations reduced it. The calculated CRC aligned well with empirical results, indicating that agonist CRC can be estimated from knowledge of L0. Reduction in agonist efficacy due to LOF mutations was calculated and subsequently validated using single-channel patch-clamp electrophysiology, a factor often obscured in normalized CRC. The study also evaluated the combined impact of mutations (L0) on CRC, confirming the predictive model. Further, no significant energetic coupling between distant residues (>15 Å) was found, indicating that the mutations' effects are localized and do not alter overall agonist affinity. These findings substantiate the role of unliganded gating in modulating agonist responses and establishes a predictive model for estimating CRC parameters from known changes in L0. The study highlights the importance of intrinsic activity in receptor theory.

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.

Agonist efficiency from concentration-response curves: Structural implications and applications

Agonists are evaluated by a concentration-response curve (CRC), with a midpoint (EC₅₀) that indicates potency, a high-concentration asymptote that indicates efficacy, and a low-concentration asymptote that indicates constitutive activity. A third agonist attribute, efficiency (η), is the fraction of binding energy that is applied to the conformational change that activates the receptor. We show that η can be calculated from EC₅₀ and the asymptotes of a CRC derived from either single-channel or whole-cell responses. For 20 agonists of skeletal muscle nicotinic receptors, the distribution of η-values is bimodal with population means at 51% (including acetylcholine, nornicotine, and dimethylphenylpiperazinium) and 40% (including epibatidine, varenicline, and cytisine). The value of η is related inversely to the size of the agonist's headgroup, with high- versus low-efficiency ligands having an average volume of 70 vs. 102 Å³. Most binding site mutations have only a small effect on acetylcholine efficiency, except for αY190A (35%), αW149A (60%), and those at αG153 (42%). If η is known, the EC₅₀ and high-concentration asymptote can be calculated from each other. Hence, an entire CRC can be estimated from the response to a single agonist concentration, and efficacy can be estimated from EC₅₀ of a CRC that has been normalized to 1. Given η, the level of constitutive activity can be estimated from a single CRC.

A mechanism for acetylcholine receptor gating based on structure, coupling, phi, and flip

Gupta et al. use single-channel electrophysiology to investigate the gating mechanism of acetylcholine receptor ion channels. They propose that channel opening starts at the M2–M3 linker and ligand-binding sites and proceeds through four brief intermediate conformations before ending with the collapse of a gate bubble.

Nicotinic acetylcholine receptors are allosteric proteins that generate membrane currents by isomerizing (“gating”) between resting and active conformations under the influence of neurotransmitters. Here, to explore the mechanisms that link the transmitter-binding sites (TBSs) with the distant gate, we use mutant cycle analyses to measure coupling between residue pairs, phi value analyses to sequence domain rearrangements, and current simulations to reproduce a microsecond shut component (“flip”) apparent in single-channel recordings. Significant interactions between amino acids separated by >15 Å are rare; an exception is between the αM2–M3 linkers and the TBSs that are ∼30 Å apart. Linker residues also make significant, local interactions within and between subunits. Phi value analyses indicate that without agonists, the linker is the first region in the protein to reach the gating transition state. Together, the phi pattern and flip component suggest that a complete, resting↔active allosteric transition involves passage through four brief intermediate states, with brief shut events arising from sojourns in all or a subset. We derive energy landscapes for gating with and without agonists, and propose a structure-based model in which resting→active starts with spontaneous rearrangements of the M2–M3 linkers and TBSs. These conformational changes stabilize a twisted extracellular domain to promote transmembrane helix tilting, gate dilation, and the formation of a “bubble” that collapses to initiate ion conduction. The energy landscapes suggest that twisting is the most energetically unfavorable step in the resting→active conformational change and that the rate-limiting step in the reverse process is bubble formation.

Mutations Causing Slow-Channel Myasthenia Reveal That a Valine Ring in the Channel Pore of Muscle AChR is Optimized for Stabilizing Channel Gating

We identify two novel mutations in acetylcholine receptor (AChR) causing a slow-channel congenital myasthenia syndrome (CMS) in three unrelated patients (Pts). Pt 1 harbors a heterozygous βV266A mutation (p.Val289Ala) in the second transmembrane domain (M2) of the AChR β subunit (CHRNB1). Pts 2 and 3 carry the same mutation at an equivalent site in the ε subunit (CHRNE), εV265A (p.Val285Ala). The mutant residues are conserved across all AChR subunits of all species and are components of a valine ring in the channel pore, which is positioned four residues above the leucine ring. Both βV266A and εV265A reduce the amino acid size and lengthen the channel opening bursts by fourfold by enhancing gating efficiency by approximately 30-fold. Substitution of alanine for valine at the corresponding position in the δ and α subunit prolongs the burst duration four- and eightfold, respectively. Replacing valine at ε codon 265 either by a still smaller glycine or by a larger leucine also lengthens the burst duration. Our analysis reveals that each valine in the valine ring contributes to channel kinetics equally, and the valine ring has been optimized in the course of evolution to govern channel gating.

Function of the M1 π-helix in endplate receptor activation and desensitization

KEY POINTS

A conserved proline in M1 causes a kink between α and π helical segments. The kink is under greater tension in the resting versus active conformation. The kink and the agonist do not interact directly. The π-helix separates the gating functions of the extracellular and transmembrane domains. Mutations of the conserved proline and propofol increase desensitization.

ABSTRACT

Nicotinic acetylcholine receptors (AChRs) switch on/off to generate transient membrane currents (C↔O; closed-open 'gating') and enter/recover from long-lived, refractory states (O↔D; 'desensitization'). The M1 transmembrane helix of the muscle endplate AChR is linked to a β-strand of the extracellular domain that extends to a neurotransmitter binding site. We used electrophysiology to measure the effects of mutations of amino acids that are located at a proline kink in M1 that separates π and α helices, in both α (N217, V218 and P221) and non-α subunits. In related receptors, the kink is straighter and more stable in O vs. C structures (gating is 'spring-loaded'). None of the AChR kink mutations had a measureable effect on agonist affinity but many influenced the allosteric gating constant substantially. Side chains in the M1 α-helix experience extraordinarily large energy differences between C and O structures, probably because of a ∼2 Å displacement and tilt of M2 relative to M1. There is a discrete break in the character of the gating transition state between αN217 and αV218, indicating that the π-helix is a border between extracellular- and transmembrane-domain function. Mutations of the conserved M1 proline, and the anaesthetic propofol, increase a rate constant for desensitization. The results suggest that straightening of the M1 proline kink triggers AChR desensitization.

Heterologous expression and nonsense suppression provide insights into agonist behavior at α6β2 nicotinic acetylcholine receptors

The α6-containing subtypes of the nicotinic acetylcholine receptor (nAChR) are localized to presynaptic terminals of the dopaminergic pathways of the central nervous system. Selective ligands for these nAChRs are potentially useful in both Parkinson's disease and addiction. For these and other goals, it is important to distinguish the binding behavior of agonists at the α6-β2 binding site versus other subtypes. To study this problem, we apply nonsense suppression-based non-canonical amino acid mutagenesis. We report a combination of four mutations in α6β2 that yield high-level heterologous expression in Xenopus oocytes. By varying mRNA injection ratios, two populations were observed with unique characteristics, likely due to differing stoichiometries. Responses to nine known nAChR agonists were analyzed at the receptor, and their corresponding EC50 values and efficacies are reported. The system is compatible with nonsense suppression, allowing structure-function studies between Trp149 - a conserved residue on loop B found to make a cation-π interaction at several nAChR subtypes - and several agonists. These studies reveal that acetylcholine forms a strong cation-π interaction with the conserved tryptophan, while nicotine and TC299423 do not, suggesting a unique pharmacology for the α6β2 nAChR.

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.

Molecular interactions between mecamylamine enantiomers and the transmembrane domain of the human α4β2 nicotinic receptor

To characterize the binding sites of mecamylamine enantiomers on the transmembrane domain (TMD) of human (h) (α4)3(β2)2 and (α4)2(β2)3 nicotinic acetylcholine receptors (AChRs), we used nuclear magnetic resonance (NMR), molecular docking, and radioligand binding approaches. The interactions of (S)-(+)- and (R)-(-)-mecamylamine with several residues, determined by high-resolution NMR, within the hα4β2-TMD indicate different modes of binding at several luminal (L) and nonluminal (NL) sites. In general, the residues sensitive to each mecamylamine enantiomer are similar at both receptor stoichiometries. However, some differences were observed. The molecular docking experiments were crucial for delineating the location and orientation of each enantiomer in its binding site. In the (α4)2(β2)3-TMD, (S)-(+)-mecamylamine interacts with the L1 (i.e., between positions -3' and -5') and L2 (i.e., between positions 16' and 20') sites, whereas the β2-intersubunit (i.e., cytoplasmic end of two β2-TMDs) and α4/β2-intersubunit (i.e., cytoplasmic end of α4-TM1 and β2-TM3) sites are shared by both enantiomers. In the (α4)3(β2)2-TMD, both enantiomers bind with different orientations to the L1' (closer to ring 2') and α4-intrasubunit (i.e., at the cytoplasmic ends of α4-TM1 and α4-TM2) sites, but only (R)-(-)-mecamylamine interacts with the L2' (i.e., closer to ring 20') and α4-TM3-intrasubunit sites. Our findings are important because they provide, for the first time, a structural understanding of the allosteric modulation elicited by mecamylamine enantiomers at each hα4β2 stoichiometry. This advancement could be beneficial for the development of novel therapies for the treatment of several neurological disorders.

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

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

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

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

Action of nicotine and analogs on acetylcholine receptors having mutations of transmitter-binding site residue αG153

A primary target for nicotine is the acetylcholine receptor channel (AChR). Some of the ability of nicotine to activate differentially AChR subtypes has been traced to a transmitter-binding site amino acid that is glycine in lower affinity and lysine in higher affinity AChRs. We studied the effects of mutations of this residue (αG153) in neuromuscular AChRs activated by nicotine and eight other agonists including nornicotine and anabasine. All of the mutations increased the unliganded gating equilibrium constant. The affinity of the resting receptor (K_d) and the net binding energy from the agonist for gating (ΔG_B) were estimated by cross-concentration fitting of single-channel currents. In all but one of the agonist/mutant combinations there was a moderate decrease in K_d and essentially no change in ΔG_B. The exceptional case was nicotine plus lysine, which showed a large, >8,000-fold decrease in K_d but no change in ΔG_B. The extraordinary specificity of this combination leads us to speculate that AChRs with a lysine at position αG153 may be exposed to a nicotine-like compound in vivo.

Loop C and the mechanism of acetylcholine receptor-channel gating

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

Structural basis for ion permeation mechanism in pentameric ligand-gated ion channels

To understand the molecular mechanism of ion permeation in pentameric ligand-gated ion channels (pLGIC), we solved the structure of an open form of GLIC, a prokaryotic pLGIC, at 2.4 Å. Anomalous diffraction data were used to place bound anions and cations. This reveals ordered water molecules at the level of two rings of hydroxylated residues (named Ser6' and Thr2') that contribute to the ion selectivity filter. Two water pentagons are observed, a self-stabilized ice-like water pentagon and a second wider water pentagon, with one sodium ion between them. Single-channel electrophysiology shows that the side-chain hydroxyl of Ser6' is crucial for ion translocation. Simulations and electrostatics calculations complemented the description of hydration in the pore and suggest that the water pentagons observed in the crystal are important for the ion to cross hydrophobic constriction barriers. Simulations that pull a cation through the pore reveal that residue Ser6' actively contributes to ion translocation by reorienting its side chain when the ion is going through the pore. Generalization of these findings to the pLGIC family is proposed.

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.

An integrated catch-and-hold mechanism activates nicotinic acetylcholine receptors

In neuromuscular acetylcholine (ACh) receptor channels (AChRs), agonist molecules bind with a low affinity (LA) to two sites that can switch to high affinity (HA) and increase the probability of channel opening. We measured (by using single-channel kinetic analysis) the rate and equilibrium constants for LA binding and channel gating for several different agonists of adult-type mouse AChRs. Almost all of the variation in the equilibrium constants for LA binding was from differences in the association rate constants. These were consistently below the limit set by diffusion and were substantially different even though the agonists had similar sizes and the same charge. This suggests that binding to resting receptors is not by diffusion alone and, hence, that each binding site can undergo two conformational changes ("catch" and "hold") that connect three different structures (apo-, LA-bound, and HA-bound). Analyses of ACh-binding protein structures suggest that this binding site, too, may adopt three discrete structures having different degrees of loop C displacement ("capping"). For the agonists we tested, the logarithms of the equilibrium constants for LA binding and LA↔HA gating were correlated. Although agonist binding and channel gating have long been considered to be separate processes in the activation of ligand-gated ion channels, this correlation implies that the catch-and-hold conformational changes are energetically linked and together comprise an integrated process having a common structural basis. We propose that loop C capping mainly reflects agonist binding, with its two stages corresponding to the formation of the LA and HA complexes. The catch-and-hold reaction coordinate is discussed in terms of preopening states and thermodynamic cycles of activation.

Drying transition in the hydrophobic gate of the GLIC channel blocks ion conduction

The theoretical prediction of water drying transitions near nonpolar surfaces has stimulated an intensive search for biological processes exploiting this extreme form of hydrophobicity. Here we quantitatively demonstrate that drying of a hydrophobic constriction is the major determinant of ion conductance in the GLIC pentameric ion channel. Molecular-dynamics simulations show that in the closed state, the channel conductance is ∼12 orders-of-magnitude lower than in the open state. This large drop in conductance is remarkable because even in the functionally closed conformation the pore constriction remains wide enough for the passage of sodium ions, aided by a continuous bridge of ∼12 water molecules. However, we find that the free energy cost of hydrating the hydrophobic gate is large, accounting almost entirely for the energetic barrier blocking ion passage. The free energies of transferring a sodium ion into a prehydrated gate in functionally closed and open states differ by only 1.2 kcal/mol, compared to an 11 kcal/mol difference in the costs of hydrating the hydrophobic gate. Conversely, ion desolvation effects play only minor roles in GLIC ion channel gating. Our simulations help rationalize experiments probing the gating kinetics of the nicotinic acetylcholine receptor in response to mutations of pore-lining residues. The molecular character and phase behavior of water should thus be included in quantitative descriptions of ion channel gating.

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

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

Linear rate-equilibrium relations arising from ion channel-bilayer energetic coupling

Linear rate-equilibrium (RE) relations, also known as linear free energy relations, are widely observed in chemical reactions, including protein folding, enzymatic catalysis, and channel gating. Despite the widespread occurrence of linear RE relations, the principles underlying the linear relation between changes in activation and equilibrium energy in macromolecular reactions remain enigmatic. When examining amphiphile regulation of gramicidin channel gating in lipid bilayers, we noted that the gating process could be described by a linear RE relation with a simple geometric interpretation. This description is possible because the gating process provides a well-understood reaction, in which structural changes in a bilayer-embedded model protein can be studied at the single-molecule level. It is thus possible to obtain quantitative information about the energetics of the reaction transition state and its position on a spatial coordinate. It turns out that the linear RE relation for the gramicidin monomer-dimer reaction can be understood, and the quantitative relation between changes in activation energy and equilibrium energy can be interpreted, by considering the effects of amphiphiles on the changes in bilayer elastic energy associated with channel gating. We are not aware that a similar simple mechanistic explanation of a linear RE relation has been provided for a chemical reaction in a macromolecule. RE relations generally should be useful for examining how amphiphile-induced changes in bilayer properties modulate membrane protein folding and function, and for distinguishing between direct (e.g., due to binding) and indirect (bilayer-mediated) effects.

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.

Temperature dependence of acetylcholine receptor channels activated by different agonists

The temperature dependence of agonist binding and channel gating were measured for wild-type adult neuromuscular acetylcholine receptors activated by acetylcholine, carbamylcholine, or choline. With acetylcholine, temperature changed the gating rate constants (Q(10) ≈ 3.2) but had almost no effect on the equilibrium constant. The enthalpy change associated with gating was agonist-dependent, but for all three ligands it was approximately equal to the corresponding free-energy change. The equilibrium dissociation constant of the resting conformation (K(d)), the slope of the rate-equilibrium free-energy relationship (Φ), and the acetylcholine association and dissociation rate constants were approximately temperature-independent. In the mutant αG153S, the choline association and dissociation rate constants were temperature-dependent (Q(10) ≈ 7.4) but K(d) was not. By combining two independent mutations, we were able to compensate for the catalytic effect of temperature on the decay time constant of a synaptic current. At mouse body temperature, the channel-opening and -closing rate constants are ∼400 and 16 ms(-1). We hypothesize that the agonist dependence of the gating enthalpy change is associated with differences in ligand binding, specifically to the open-channel conformation of the protein.

Decrypting the sequence of structural events during the gating transition of pentameric ligand-gated ion channels based on an interpolated elastic network model

Despite many experimental and computational studies of the gating transition of pentameric ligand-gated ion channels (pLGICs), the structural basis of how ligand binding couples to channel gating remains unknown. By using a newly developed interpolated elastic network model (iENM), we have attempted to compute a likely transition pathway from the closed- to the open-channel conformation of pLGICs as captured by the crystal structures of two prokaryotic pLGICs. The iENM pathway predicts a sequence of structural events that begins at the ligand-binding loops and is followed by the displacements of two key loops (loop 2 and loop 7) at the interface between the extracellular and transmembrane domain, the tilting/bending of the pore-lining M2 helix, and subsequent movements of M4, M3 and M1 helices in the transmembrane domain. The predicted order of structural events is in broad agreement with the Φ-value analysis of α subunit of nicotinic acetylcholine receptor mutants, which supports a conserved core mechanism for ligand-gated channel opening in pLGICs. Further perturbation analysis has supported the critical role of certain intra-subunit and inter-subunit interactions in dictating the above sequence of events.

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

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

Pore opening and closing of a pentameric ligand-gated ion channel

Nerve signaling in humans and chemical sensing in bacteria both rely on the controlled opening and closing of the ion-conducting pore in pentameric ligand-gated ion channels. With the help of a multiscale simulation approach that combines mixed elastic network model calculations with molecular dynamics simulations, we study the opening and closing of the pore in Gloeobacter violaceus channel GLIC at atomic resolution. In our simulations of the GLIC transmembrane domain, we first verify that the two endpoints of the transition are open and closed to sodium ion conduction, respectively. We then show that a two-stage tilting of the pore-lining helices induces cooperative drying and iris-like closing of the channel pore. From the free energy profile of the gating transition and from unrestrained simulations, we conclude that the pore of the isolated GLIC transmembrane domain closes spontaneously. The mechanical work of opening the pore is performed primarily on the M2-M3 loop. Strong interactions of this short and conserved loop with the extracellular domain are therefore crucial to couple ligand binding to channel opening.

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

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

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

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

Relationships between structural dynamics and functional kinetics in oligomeric membrane receptors

Recent efforts to broaden understanding of the molecular mechanisms of membrane receptors in signal transduction make use of rate-equilibrium free-energy relationships (REFERs), previously applied to chemical reactions, enzyme kinetics, and protein folding. For oligomeric membrane receptors, we distinguish between a), the Leffler parameter α_L, to characterize the global transition state for the interconversion between conformations; and b), the Fersht parameter, ϕ_F, to assign the degree of progression of individual residue positions at the transition state. For both α_L and ϕ_F, insights are achieved by using harmonic energy profiles to reflect the dynamic nature of proteins, as illustrated with single-channel results reported for normal and mutant nicotinic receptors. We also describe new applications of α_L based on published results. For large-conductance calcium-activated potassium channels, data are satisfactorily fit with an α_L value of 0.65, in accord with REFERs. In contrast, results reported for the flip conformational state of glycine and nicotinic receptors are in disaccord with REFERs, since they yield α_L values outside the usual limits of 0–1. Concerning published ϕ_F values underlying the conformational wave hypothesis for nicotinic receptors, we note that interpretations may be complicated by variations in the width of harmonic energy profiles.

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.

One-microsecond molecular dynamics simulation of channel gating in a nicotinic receptor homologue

Recently discovered bacterial homologues of eukaryotic pentameric ligand-gated ion channels, such as the Gloeobacter violaceus receptor (GLIC), are increasingly used as structural and functional models of signal transduction in the nervous system. Here we present a one-microsecond-long molecular dynamics simulation of the GLIC channel pH stimulated gating mechanism. The crystal structure of GLIC obtained at acidic pH in an open-channel form is equilibrated in a membrane environment and then instantly set to neutral pH. The simulation shows a channel closure that rapidly takes place at the level of the hydrophobic furrow and a progressively increasing quaternary twist. Two major events are captured during the simulation. They are initiated by local but large fluctuations in the pore, taking place at the top of the M2 helix, followed by a global tertiary relaxation. The two-step transition of the first subunit starts within the first 50 ns of the simulation and is followed at 450 ns by its immediate neighbor in the pentamer, which proceeds with a similar scenario. This observation suggests a possible two-step domino-like tertiary mechanism that takes place between adjacent subunits. In addition, the dynamical properties of GLIC described here offer an interpretation of the paradoxical properties of a permeable A13'F mutant whose crystal structure determined at 3.15 A shows a pore too narrow to conduct ions.

Gating transition of pentameric ligand-gated ion channels

Pentameric ligand-gated ion channels are an important family of membrane proteins and play key roles in physiological processes, including signal transduction at chemical synapses. Here, we study the conformational changes associated with the opening and closing of the channel pore. Based on recent crystal structures of two prokaryotic members of the family in open and closed states, respectively, mixed elastic network models are constructed for the transmembrane domain. To explore the conformational changes in the gating transition, a coarse-grained transition path is computed that smoothly connects the closed and open conformations of the channel. We find that the conformational transition involves no major rotations of the transmembrane helices, and is instead characterized by a concerted tilting of helices M2 and M3. In addition, helix M2 changes its bending state, which results in an early closure of the pore during the open-to-closed transition.

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.

Structural basis of activation of cys-loop receptors: the extracellular-transmembrane interface as a coupling region

Cys-loop receptors mediate rapid transmission throughout the nervous system by converting a chemical signal into an electric one. They are pentameric proteins with an extracellular domain that carries the transmitter binding sites and a transmembrane region that forms the ion pore. Their essential function is to couple the binding of the agonist at the extracellular domain to the opening of the ion pore. How the structural changes elicited by agonist binding are propagated through a distance of 50 A to the gate is therefore central for the understanding of the receptor function. A step forward toward the identification of the structures involved in gating has been given by the recently elucidated high-resolution structures of Cys-loop receptors and related proteins. The extracellular-transmembrane interface has attracted attention because it is a structural transition zone where beta-sheets from the extracellular domain merge with alpha-helices from the transmembrane domain. Within this zone, several regions form a network that relays structural changes from the binding site toward the pore, and therefore, this interface controls the beginning and duration of a synaptic response. In this review, the most recent findings on residues and pairwise interactions underlying channel gating are discussed, the main focus being on the extracellular-transmembrane interface.